data-raw/03_dds_small.R

In roryk/bcbioRnaseq: RNA-Seq Utilities

library(pryr)
library(devtools)
library(DESeq2)
load_all()

dds_small <- as(bcb_small, "DESeqDataSet")
design(dds_small) <- ~ treatment
dds_small <- DESeq(dds_small)
validObject(dds_small)

stopifnot(design(dds_small) == ~ treatment)
stopifnot(identical(
    resultsNames(dds_small),
    c("Intercept", "treatment_folic_acid_vs_control")
))

rld_small <- rlog(dds_small)
validObject(rld_small)

resultsNames(dds_small)
res_small <- results(
    dds_small,
    contrast = c(
        factor = "treatment",
        numerator = "folic_acid",
        denominator = "control"
    )
)
# Shrink log2 fold changes
res_small <- lfcShrink(
    dds = dds_small,
    coef = 2L,
    res = res_small
)
validObject(res_small)

use_data(
    dds_small, rld_small, res_small,
    overwrite = TRUE, compress = "xz"
)