# NOT RUN {
# NOTE: Requires the edgeR package
# Add normalized counts and log2CPM as additional "assay" items in the DGEobj
dgeObj <- readRDS(system.file("exampleObj.RDS", package = "DGEobj"))
dgeList <- edgeR::calcNormFactors(edgeR::DGEList(dgeObj$counts), method="TMM")
log2cpm <- edgeR::cpm(dgeList, log = TRUE)
dgeObj <- addItems(dgeObj,
itemList = list(newDgelist = dgeList, Log2CPM = log2cpm),
itemTypes = list("assay", "assay"),
parents = list("counts", "newDgelist")
)
inventory(dgeObj)
# }
# NOT RUN {
# }
Run the code above in your browser using DataLab