modeG {MasterBayes} | R Documentation |
Finds the mode of the posterior marginal distribution of genotypes
modeG(postG, threshold=0, ...)
postG |
posterior distribution of genotypes from an MCMCped model with argument write_postG=TRUE |
threshold |
threshold probability under which ML genotypes are replaced by NA |
... |
further arguments to be passed |
G |
list of genotype objects |
id |
id vector |
Jarrod Hadfield j.hadfield@ed.ac.uk
Hadfield J.D. et al, Molecular Ecology
MCMCped
, genotype
data(WarblerP) data(WarblerG) GdP<-GdataPed(WarblerG) var1<-expression(varPed(c("lat", "long"), gender="Male", relational="OFFSPRING")) # paternity is to be modelled as a function of distance # between offspring and male territories res1<-expression(varPed("offspring", restrict=0)) # indivdiuals from the offspring generation are excluded as parents res2<-expression(varPed("terr", gender="Female", relational="OFFSPRING", restrict="==")) # mothers not from the offspring territory are excluded PdP<-PdataPed(formula=list(var1,res1,res2), data=WarblerP, USsire=FALSE) tP<-tunePed(beta=30) model1<-MCMCped(PdP=PdP, GdP=GdP, tP=tP, nitt=3000, thin=2, burnin=1000, write_postG=TRUE) G<-modeG(model1$G)$G summary(G[[1]])