| 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@sheffield.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=TRUE))
# 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]])