indelmiss (version 1.0.10)

plotrates: Plot estimates for insertion and deletion rates

Description

Plotting command for use on an object of class "indelmiss".

Usage

plotrates(x, model, ci = TRUE, cil = 95, ...)

Value

Plot estimates for insertion and deletion rates. Called by plot.indelmiss.

Arguments

x

An object of class "indelmiss".

model

One of "M1", "M2", "M3", or "M4".

ci

TRUE plots confidence intervals around the estimates.

cil

Confidence interval level.

...

Any further commands to plot.

Author

Utkarsh J. Dang and G. Brian Golding

See Also

See also plot.indelmiss.

Examples

indel <- indelrates(datasource = "simulation", seed = 1, taxa = 5,
             mu = 1, nu = 5, phyl = 5000,  pmiss = 0, toi = 1, 
             bgtype="ancestornodes", bg = c(7, 9), 
             zerocorrection = TRUE, 
             modelnames = c("M1", "M2", "M3", "M4"),
             optmethod = "nlminb")
print(indel)
plotrates(indel, model="M4")