plotrates {indelmiss}

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, ...)

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.

Value

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

Author(s)

Utkarsh J. Dang and G. Brian Golding

[e-mail]

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