nbinomGLMTest {DESeq} | R Documentation |
For each gene, the function calculates a chi-square p value by simply calculating:
1 - pchisq(resReduced$deviance - resFull$deviance, attr(resReduced,
"df.residual") - attr(resFull, "df.residual"))
nbinomGLMTest(resFull, resReduced)
resFull, resReduced |
GLM fit data frames, as returned by |
a vector of p values
Simon Anders, anders@embl.de
cds <- makeExampleCountDataSet()[ 1:100, ] cds <- estimateSizeFactors( cds ) cds <- estimateDispersions( cds, method="pooled" ) fit1 <- fitNbinomGLMs( cds, count ~ condition ) fit0 <- fitNbinomGLMs( cds, count ~ 1 ) nbinomGLMTest( fit1, fit0 )