Pre-process features

In the video exercise, you saw that the sample distributions for the doxorubicin study were extremely right-skewed. Thus, the first step you need to take is to pre-process the features: log-transform, normalize, and filter.

Este ejercicio forma parte del curso

Differential Expression Analysis with limma in R

Ver curso

Instrucciones del ejercicio

The ExpressionSet object eset_raw with the raw data has been loaded in your workspace. The limma package is loaded.

Log transform the measurements. Use plotDensities to visualize. Label the samples by their genotype.
Quantile normalize the measurements with normalizeBetweenArrays and re-visualize.
Use rowMeans to determine which genes have a mean expression level greater than 0.
Filter the genes (i.e. rows) with the logical vector keep and re-visualize.

Ejercicio interactivo práctico

Prueba este ejercicio y completa el código de muestra.

# Create a new ExpressionSet to store the processed data
eset <- eset_raw
# Log transform
exprs(eset) <- ___(exprs(eset))
___(eset,  group = pData(eset)[___], legend = "topright")
# Quantile normalize
exprs(eset) <- ___(exprs(eset))
___(eset,  group = pData(eset)[___], legend = "topright")
# Determine the genes with mean expression level greater than 0
keep <- ___(exprs(eset)) > ___
sum(keep)
# Filter the genes
eset <- eset[___]
___(eset, group = pData(eset)[___], legend = "topright")

Editar y ejecutar código