Visualize batch effects
In the olfactory stem cell experiment, there were 7 treatments with 4 replicates each for a total of 28 samples. However, these 28 samples were processed in 4 separate batches. The effect of the treatments is of biological interest, but the effect of the batches is technical noise. Using dimension reduction, determine which of these two effects had a larger impact on the gene expression data.
Diese Übung ist Teil des Kurses
Differential Expression Analysis with limma in R
Anleitung zur Übung
The ExpressionSet object eset with the olfactory stem cell data has been loaded in your workspace.
Use
plotMDSto plot the principal components. Label the samples by the treatment they received, and setgene.selectionto"common".Re-visualize the principal components, labeling the samples by their batch.
Interaktive Übung
Versuche dich an dieser Übung, indem du diesen Beispielcode vervollständigst.
# Load package
library(limma)
# Plot principal components labeled by treatment
___(eset, labels = pData(eset)[, "___"], gene.selection = "___")
# Plot principal components labeled by batch
___(eset, labels = pData(eset)[, "___"], gene.selection = "___")