optimPlot {DoubleExpSeq} | R Documentation |
The MLE solution is signified on a plot as the minimum of the negative log likelihood of the generalized beta prime distribution, parameterized in terms of the weight parameter.
optimPlot(y, m, groups, contrast=c(1,2), use.all.groups=TRUE,...)
y |
numeric matrix of inclusion counts. |
m |
numeric matrix of total counts: inclusion + exclusion. |
groups |
vector or factor giving the experimental group/condition for each sample/library. |
contrast |
numeric vector of length 2 specifying which levels of the "groups" factor should be compared. This is only relevant if "use.all.groups" is FALSE. |
use.all.groups |
logical. If TRUE, all data in "y" is used to estimate dispersions. If FALSE, only the 2 groups given in "contrasts" are used to estimate dispersions. Only makes a difference if "y" contains more than 2 groups. Default is TRUE. |
... |
further arguments passes to plot() |
The MLE estimate of the (transformed) weight parameter in the WEB-Seq method is checked to be a true global minimum of the negative log likelihood of the generalized beta prime distribution. The weight parameter is transformed from an infinite range to the (0,1) range before optimization and this is the range on which the estimate is checked.
A plot to the current device
Sean Ruddy
data(exon)
# If all groups were used to estimate dispersions
optimPlot(counts, offsets, groups)
# If only the 2 groups being compared were used to estimate dispersions
optimPlot(counts, offsets, groups, contrast=c(1,3), use.all.groups=FALSE)