Displaying 2 results from an estimated 2 matches for "npopinit".
Did you mean:
noinit
2009 Jan 13
3
problem whit Geneland
...- 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)
But at the end R say to me that do had found the function "mcmcFmodel".
Any sugestion?
thanks...
2010 May 20
1
Geneland error on unix: Error in MCMC(........ :, unused argument(s) (ploidy = 2, genotypes = geno)
...coord<-read.table(XY.TABLE)
+
+ ## Loop for multiple runs
+
+ for(irun in 1:nrun)
+ {
+ ## define path to MCMC directory
+
+ path.mcmc <- paste(path.sp,irun,"/",sep="")
+ dir.create(path.mcmc)
+ MCMC(coordinates=coord, ploidy=2, genotypes=geno, varnpop=T,
npopmax=MAXPOP, npopinit=MINPOP, spatial=T, freq.model="Correlated",
nit=NITS, thinning=500, rate.max=nrow(geno), delta.coord=100,
path.mcmc=path.mcmc)
+ ## MCMC postprocessing
+ PostProcessChain(coordinates=coord,path.mcmc=path.mcmc,genotypes=geno,
nxdom=50,nydom=50,burnin=burnin)
+ }
+ ## Computing average p...