Gordon Smyth
August 16, 2005
Walter and Eliza Hall Institute of Medical Research
Aims
This laboratory explores some of the features of the limma package for assessing differential expression in microarray experiments. Examples are included of cDNA two-color microarrays and Affymetrix one-channel microarrays. Some pre-processing issues are also discussed for two-color arrays.
Lab exercises
Exercise | Platform | Design | Topics covered |
apoAI data | cDNA | Two group comparison with common reference | Introduction to linear models. Obtaining empirical Bayes statistics. Getting lists of differentially expressed genes. |
integrin beta7 data | cDNA | Direct comparisons with dye-swaps | Data entry for two color data. Highlighting control probes. Exploring different background correction methods. Allowing for genewise dye effects. |
Estrogen data | Affymetrix | 2x2 Factorial | More on linear models. Use of contrasts. Venn diagrams. Linking gene lists to annotation information on the internet. Gene set tests. |
Drosophila embryogenesis dataset | Affymetrix | Time course with series-level replication | Time course analysis using linear models and moderated F-statistics. |
Datasets used in the exercises
Please check whether you already have the "Drosophila Embryo" and "Estrogen" packages from the Required Software. (These data sets are stored within R packages.)
- ApoAI Knockout Data (1.2 MB)
- Integrin beta7 data (14.6 MB)
- Estrogen Data (23.6 MB)
- Drosophila Embryo Data (as an R package): Windows - Source (Mac or Linux) (3.5 MB)
Lab Handouts
The handouts in pdf can be accessed here: lab
handouts
Required R packages
Please install these packages before attempting to
repeat the lab exercises. Note: if you don't have write
permission to the system library directory of your R
installation, you can use the .libPaths() function with
something like .libPaths("C:/mylibdir")
before you run install.packages()
(or
equivalent) to install the packages in a customized
directory location.
A local cache of packages is here: Local Package Cache
Package |
Windows |
MacOS X |
Source |
limma_2.0.4 |
|||
statmod_1.2.0 |
|||
affy_1.6.7 |
|||
Biobase_1.5.12 |
|||
hgu95av2_1.8.4 |
|||
hgu95av2cdf_1.5.1 |
|||
xtable_1.2-5 |
Getting started
You should be running R 2.1.0 or 2.1.1 and limma 2.0.4. A good way to get started is to open up the Limma User's Guide:
library(limma)
If you're using Windows, just use the drop-down menu "Vignettes". Otherwise, type
limmaUsersGuide()
References
- Smyth, G. K., Thorne, N. P. and Wettenhall J. (2005) limma: Linear Models for Microarray Data User's Guide. http://bioinf.wehi.edu.au/limma (Included as part of the limma package.)
- Smyth, G. K. (2005). Limma: linear models for microarray data. In: Bioinformatics and Computational Biology Solutions using R and Bioconductor, R. Gentleman, V. Carey, S. Dudoit, R. Irizarry, W. Huber (eds.), Springer, New York.
- Smyth, G. K. (2004). Linear models and empirical Bayes methods for assessing differential expression in microarray experiments. Statistical Applications in Genetics and Molecular Biology 3, No. 1, Article 3. http://www.bepress.com/sagmb/vol3/iss1/art3/