Displaying 2 results from an estimated 2 matches for "codom".
Did you mean:
odom
2010 May 31
0
miss.loc function in MCMC Geneland: can't make it work
...the command I enter:
library(Geneland)
coord<-read.table("path.txt")
geno<-read.table("path.txt")
nullal<-read.table("path.txt")
dim(nullal)
dim(coord)
dim(geno)
plot(coord,xlab="east",ylab="north",asp=1)
MCMC(coordinates=coord,geno.dip.codom=geno,miss.loc=nullal,delta.coord=0,varnpop=TRUE,npopmax=40,spatial=TRUE,freq.model="Uncorrelated",nit=100000,thinning=10,path.mcmc="path_whichever
it is/")
Olivier BAGGIANO
PhD Candidate
Australian Rivers Institute
Griffith School of Environment
Griffith University, Nathan Q...
2011 Apr 04
0
Multithreading of Geneland
Hi all,
I would like to multithread 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 18...