runLimma | R Documentation |
This function runs DEG analysis with Limma at user defined FDR and LFC cutoff by providing CMAP02 normalized expression values and annotation of treatment v.s. control samples. For more context of the CMAP02 database, please consult the vignette of this package.
runLimma(df, comp_list, fdr = 0.05, foldchange = 1, verbose = TRUE)
df |
data.frame containing normalized intensity values of CMAP02 samples |
comp_list |
list of CEL file ids of treatment and control samples for
each compound treatment. The list for CMAP02 data is generated from
the |
fdr |
cutoff of false discovery rate (FDR) for defining DEGs |
foldchange |
cutoff of log2 fold change (LFC) for defining DEGs |
verbose |
TRUE or FALSE |
list containing DEGs and log2FC matrix
path <- system.file("extdata", "cmap_instances_02.txt", package="signatureSearchData")
cmap_inst <- read.delim(path, check.names=FALSE)
comp_list <- sampleList(cmap_inst, myby="CMP_CELL")
df <- as.data.frame(matrix(runif(70), ncol=7))
colnames(df) <- unlist(comp_list[1])
degList <- runLimma(df, comp_list[1])
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.