MasterBayes (version 2.58)

modeG: Posterior Mode of Genotypes

Description

Finds the mode of the posterior marginal distribution of genotypes

Usage

modeG(postG, threshold=0)

Arguments

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

Value

G

list of genotype objects

id

id vector

References

Hadfield J.D. et al, Molecular Ecology

See Also

MCMCped, genotype

Examples

# NOT RUN {
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]])
# }