Bioinformatics Advance Access published online on March 14, 2008
Bioinformatics, doi:10.1093/bioinformatics/btn100
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Unequal group variances in microarray data analyses
1 Department of Statistics, University of Örebro, Sweden; 2Department of Medical Epidemiology and Biostatistics, Karolinska Institutet, Stockholm, Sweden; 3Department of Biomedical Sciences and Biotechnology, University of Brescia, Italy
To whom correspondence should be addressed. Prof. Yudi Pawitan, E-mail: yudi.pawitan{at}ki.se
| Abstract |
|---|
Motivation: In searching for differentially expressed (DE) genes in microarray data, we often observe a fraction of the genes to have unequal variability between groups. This is not an issue in large samples, where a valid test exists that uses individual variances separately. The problem arises in the small-sample setting, where the approximately-valid Welch test lacks sensitivity, while the more sensitive moderated t test assumes equal variance.
Methods: We introduce a moderated Welch test (MWT) that allows nequal variance between groups. It is based on (i) weighting of pooled and unpooled standard errors and (ii) improved estimation of the gene-level variance that exploits the information from across the genes.
Results: When a nontrivial proportion of genes has unequal variability, FDR estimates based on the standard t and moderated t tests are often too optimistic, while the standard Welch test has low sensitivity. The MWT is shown to (i) perform better than the standard t, the standard Welch and the moderated t tests when the variances are unequal between groups, and (ii) perform similarly to the moderated t, and better than the standard t andWelch tests when the group variances are equal. These results mean that MWT is more reliable than other existing tests over wider range of data conditions.
Availability: R package to perform MWT is available at http://www.meb.ki.se/~yudpaw
Contact: yudi.pawitan{at}ki.se
Associate Editor: Prof. Martin Bishop
Received on January 21, 2008; revised on February 23, 2008; accepted on March 12, 2008