PlotFamilies {TE}

Plot the age distributions or insertion rates for multiple families.

Description

Plot the age distributions or insertion rates for multiple families.

Usage

PlotFamilies(resList, type = c("insRt", "ageDist"), ...)

Arguments

resList

A list of TEfit/TEfit2 objects, which can be mixed

type

Whether to plot the insertion rates ('insRt') or the age distributions ('ageDist').

...

Passed into plotting functions.

Value

A list of line data (plotDat) and peak locations (peakDat).

Examples

data(AetLTR)
copia3 <- subset(AetLTR, GroupID == 3 & !is.na(Chr))
gypsy24 <- subset(AetLTR, GroupID == 24 & !is.na(Chr))
res3 <- EstDynamics(copia3$Mismatch, copia3$UngapedLen)
res24 <- EstDynamics2(gypsy24$Mismatch, gypsy24$UngapedLen)

# Plot insertion rates
PlotFamilies(list(`Copia 3`=res3, `Gypsy 24`=res24))

# Plot age distributions
PlotFamilies(list(`Copia 3`=res3, `Gypsy 24`=res24), type='ageDist')