SensitivityPlot: Generate sensitivity plots

In TE: Insertion/Deletion Dynamics for Transposable Elements

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

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

TE documentation built on May 1, 2019, 10:13 p.m.