Plot estimates for insertion and deletion rates

Share:

Description

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

Usage

1
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.

Author(s)

Utkarsh J. Dang and G. Brian Golding

udang@mcmaster.ca

See Also

See also plot.indelmiss.

Examples

1
2
3
4
5
6
7
8
#indel <- indelrates(datasource = "simulation", seed = 1, taxa = 5,
#             mu = 1, nu = 5, phyl = 5000,  nmiss = 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")