modeG {MasterBayes} | R Documentation

## 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
`...` | further arguments to be passed

### Value

`G` | list of `genotype` objects
`id` | id vector

### Author(s)

Jarrod Hadfield j.hadfield@sheffield.ac.uk

### References

Hadfield J.D. et al, Molecular Ecology

### See Also

`MCMCped`, `genotype`

### Examples

    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]])

[Package MasterBayes version 2.0 Index]