Full_Name: Gilles Guillot
Version: 2.6.0
OS: windows XP professional
Submission from: (NULL) (129.240.88.50)
This sequence of command does not work in the R gui 2.6.0:
library(Geneland)
set.seed(1)
data <- simdata(nindiv=200,
coord.lim=c(0,1,0,1) ,
number.nuclei=5 ,
allele.numbers=rep(10,20),
IBD=FALSE,
npop=2,
give.tess.grid=FALSE)
geno <- data$genotypes
coord <- t(data$coord.indiv)
path.mcmc <- paste(tempdir(),"/",sep="")
set.seed(1)
mcmcFmodel(coordinates=coord,
genotypes=geno,
path.mcmc=path.mcmc,
rate.max=10,
delta.coord=0,
npopmin=1,
npopinit=5,
npopmax=5,
nb.nuclei.max=50,
nit=500,
thinning=1,
freq.model="Dirichlet",
varnpop=FALSE,
spatial=TRUE)
the call to mcmcFmodel freezes R
The same sequence of command works in the R command line of R 2.6.0
and also in the R GUI of 2.5.1
______________________________________________
R-devel@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel
