Sample size calculation for microarray experiments with blocked one-way design
1 Department of Biostatistics and Bioinformatics, Duke University Medical Center, North Carolina 27710, USA
2 Center for Fibroid Biology and Therapy, Duke University Medical Center, North Carolina 27710, USA
BMC Bioinformatics 2009, 10:164 doi:10.1186/1471-2105-10-164Published: 28 May 2009
One of the main objectives of microarray analysis is to identify differentially expressed genes for different types of cells or treatments. Many statistical methods have been proposed to assess the treatment effects in microarray experiments.
In this paper, we consider discovery of the genes that are differentially expressed among K (> 2) treatments when each set of K arrays consists of a block. In this case, the array data among K treatments tend to be correlated because of block effect. We propose to use the blocked one-way ANOVA F-statistic to test if each gene is differentially expressed among K treatments. The marginal p-values are calculated using a permutation method accounting for the block effect, adjusting for the multiplicity of the testing procedure by controlling the false discovery rate (FDR). We propose a sample size calculation method for microarray experiments with a blocked one-way design. With FDR level and effect sizes of genes specified, our formula provides a sample size for a given number of true discoveries.
The calculated sample size is shown via simulations to provide an accurate number of true discoveries while controlling the FDR at the desired level.