plot.segratioMCMC {polySegratioMM}R Documentation

MCMC plots for segregation ratio mixture models

Description

Standard MCMC trace and density plots for specified mixure model parameters and posterior probability distributions for specified markers

Usage

plot.segratioMCMC(x, ..., row.index = c(1:10), var.index = c(1:6),
marker.index = c(1:8))

Arguments

x object of class segratioMCMC
... extra options for printing
row.index which rows to print (Default: first 10)
var.index which mixture model variable to summarise (Default: all)
marker.index which markers to summarise (Default: 1:8)

Value

None.

Author(s)

Peter Baker p.baker1@uq.edu.au

See Also

dosagesMCMC readJags

Examples

## simulate small autooctaploid data set
a1 <- sim.autoMarkers(8,c(0.7,0.2,0.1),n.markers=100,n.individuals=50)
##print(a1)
sr <-  segregationRatios(a1$markers)
x <- setModel(3,8)

## Not run: 
## fit simple model in one hit and summarise

x.run <- runSegratioMM(sr, x, burn.in=200, sample=500)
plot(x.run$mcmc.mixture)
## End(Not run)

[Package polySegratioMM version 0.5-2 Index]