Displaying 1 result from an estimated 1 matches for "nullal".
Did you mean:
ullal
2010 May 31
0
miss.loc function in MCMC Geneland: can't make it work
...e MCMC is stopped. I only get this message when entering the 'null
allele' matrix and using the miss.loc function in the MCMC.
Any idea what this is about? Here is 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,thinnin...