EstimateWEBDisp {DoubleExpSeq} | R Documentation |
Calculation of shrunken dispersion estimates via weighted likelihood where the weight parameter is estimated using an empirical bayes strategy.
EstimateWEBDisp(y, m, groups, neff = NULL, S = NULL)
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. |
neff |
numeric vector of length equal to the number of rows of "y" where each value is the effective sample size for the event. Default is NULL in which case the effective sample size is calculated within the function. |
S |
numeric vector of length equal to the number of rows of "y" where each value is the random variable for each event whose distribution across exons is gamma. Default is NULL in which case the vector is calculated internally. |
Shrunken dispersion estimates are obtained by maximizing the weighted sum of the likelihood for a given event and the sum of likelihoods for all events, the common likelihood. The weight given to the common likelihood is estimated via empirical bayes.
vector of length equal to the number of rows of "y" where each value is the estimate of dispersion.
Sean Ruddy
data(exon)
dispersions <- EstimateWEBDisp( counts, offsets, groups)