tests/testthat/test_LRT_prior.R
In nlhuong/ZeroInflatedDESeq2: Differential gene expression analysis based on the negative binomial distribution
dds <- makeExampleDESeqDataSet(n=100)
colData(dds)$condition <- factor(rep(1:3,each=4))
colData(dds)$group <- factor(rep(1:2,6))
design(dds) <- ~ group + condition
dds <- DESeq(dds,test="LRT",reduced=~ group,betaPrior=TRUE)
expect_true(any(attr(dds,"betaPriorVar") < 1e6))
res <- results(dds)
expect_true(grepl("LRT",mcols(res)$description[colnames(res) == "stat"]))
design(dds) <- ~ group * condition
dds <- DESeq(dds,test="LRT",reduced=~ group,betaPrior=TRUE)
res <- results(dds)
nlhuong/ZeroInflatedDESeq2 documentation built on May 23, 2019, 9:06 p.m.
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
dds <- makeExampleDESeqDataSet(n=100)
colData(dds)$condition <- factor(rep(1:3,each=4))
colData(dds)$group <- factor(rep(1:2,6))
design(dds) <- ~ group + condition
dds <- DESeq(dds,test="LRT",reduced=~ group,betaPrior=TRUE)
expect_true(any(attr(dds,"betaPriorVar") < 1e6))
res <- results(dds)
expect_true(grepl("LRT",mcols(res)$description[colnames(res) == "stat"]))
design(dds) <- ~ group * condition
dds <- DESeq(dds,test="LRT",reduced=~ group,betaPrior=TRUE)
res <- results(dds)
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
Embedding an R snippet on your website
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.