search for: postprocesschain

...t;,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 posterior probability + + lpd <- rep(NA,nrun) + for(irun in 1:nrun) + { + path.lpd <- paste(path.mcmc,"log.posterior.density.txt",sep="") + lpd[irun] <- me...

Hello, anybody know the procedure to plot the geneland result in a map? thank you ********************************************************** Vincenzo Landi Post Doctorate student Animal genomic and breeding cell:0039/3395388713 Fax. 075-5857122 [[alternative HTML version deleted]]

...thread that script, to detect structure from multilocus genetic data : >library(Geneland) > >geno = read.table("cot966gen_test.txt") #the file is show after >MCMC(geno.dip.codom = geno, varnpop=T, npopmax=20, spatial = F, nit=100000, thinnin=100, path.mcmc="./") >PostProcessChain(path.mcmc="./", nxdom=100, nydom=100, burnin=200) I have an 8 cores computer, and since I have to compute that bit on a thousand of line like that : 209 209 217 217 180 180 154 154 181 181 192 192 -9 -9 211 211 -9 -9 160 160 -9 -9 -9 -9 -9 -9 254 254 140 140 181 181 -9 -9 -9 -9 211 211...