TE (version 0.3-0)

SensitivityPlot: Generate sensitivity plots

Description

Create sensitivity plots of a few families to investigate different death rate scenarios

Usage

SensitivityPlot(resList, col, xMax, markHalfPeak = FALSE,
  famLegend = TRUE, rLegend = names(resList), ...)

Arguments

resList

A list of families returned by EstDynamics

col

A vector of colors

xMax

The maximum of the x-axis

markHalfPeak

Whether to mark the time points with half-intensity

famLegend

Whether to create legend for families

rLegend

Text for the legend for families

...

Passed into matplot

Examples

# NOT RUN {
data(AetLTR)
copia3 <- subset(AetLTR, GroupID == 3 & !is.na(Chr))
copia9 <- subset(AetLTR, GroupID == 9 & !is.na(Chr))
res3 <- EstDynamics(copia3$Mismatch, copia3$UngapedLen)
res9 <- EstDynamics(copia9$Mismatch, copia9$UngapedLen)
SensitivityPlot(list(`Copia 3`=res3, `Copia 9`=res9))
# }