marginal.LR: Calculate the integrated likelihood ratio.

In lunge111/intSEQ: Differential Expression Analysis of RNA-seq Data with Integrated Likelihood Method

Description

A function to calculate the integrated likelihood ratio statistics.

Usage

marginal.LR(thetas, disp.fit, data,  thet, fit.thet, lambda1, 
lambda2, meanmeth, is.smooth = TRUE, m0, ma,const,w1,w2)

Arguments

thetas	Selected dispersions that located in the neighborhood of objective gene.
disp.fit	Selected smoothed dispersions that located in the neighborhood of objective gene.
data	The count data of objective gene
thet	The estimated dispersion of objective gene.
fit.thet	The smoothed dispersion of objective gene.
lambda1	The first parameter of prior variance. See details of intSEQ.
lambda2	The second parameter of prior variance. See details of intSEQ.
meanmeth	Use the estimated dispersion or the local mean as the mean of normal prior.
is.smooth	logical. Whether to use the smoothed dispersion or the estimator itself as the mean of normal prior.
m0	A numeric vector indicating the mean of objective gene under null model.
ma	A numeric vector indicating the mean of objective gene under full model.
const	A multiplying constant to make the product of likelihoods non-zero.
w1	See Details in intSEQ.
w2	See Details in intSEQ.

Details

The marginal.LR returns the integrated likelihood ratio statistic for a objective gene. See intSEQ for detail.

Value

A scalar indicating the integrated likelihood ratio statistic

Author(s)

Yilun Zhang, David Rocke

References

our paper

See Also

intSEQ

lunge111/intSEQ documentation built on May 20, 2019, 9:38 a.m.