tests/testthat/test_QR.R

set.seed(1)
dds <- makeExampleDESeqDataSet(n=100,betaSD=1)
dds <- DESeq(dds)
ddsNoQR <- nbinomWaldTest(dds,useQR=FALSE)
res <- results(dds)
resNoQR <- results(ddsNoQR)
expect_equal(res$log2FoldChange, resNoQR$log2FoldChange, tolerance=1e-6)
nlhuong/ZeroInflatedDESeq2 documentation built on May 23, 2019, 9:06 p.m.