Bias in the estimation of false discovery rate in microarray studies.

@article{citeulike:1122674, abstract = {MOTIVATION: The false discovery rate (FDR) provides a key statistical assessment for microarray studies. Its value depends on the proportion pi(0) of non-differentially expressed (non-DE) genes. In most microarray studies, many genes have small effects not easily separable from non-DE genes. As a result, current methods often overestimate pi(0) and FDR, leading to unnecessary loss of power in the overall analysis. METHODS: For the common two-sample comparison we derive a natural mixture model of the test statistic and an explicit bias formula in the standard estimation of pi(0). We suggest an improved estimation of pi(0) based on the mixture model and describe a practical likelihood-based procedure for this purpose. RESULTS: The analysis shows that a large bias occurs when pi(0) is far from 1 and when the non-centrality parameters of the distribution of the test statistic are near zero. The theoretical result also explains substantial discrepancies between non-parametric and model-based estimates of pi(0). Simulation studies indicate mixture-model estimates are less biased than standard estimates. The method is applied to breast cancer and lymphoma data examples. AVAILABILITY: An R-package OCplus containing functions to compute pi(0) based on the mixture model, the resulting FDR and other operating characteristics of microarray data, is freely available at http://www.meb.ki.se/\~{}yudpaw CONTACT: yudi.pawitan@meb.ki.se and alexander.ploner@meb.ki.se.}, address = {Department of Medical Epidemiology and Biostatistics, Karolinska Institutet, Stockholm, Sweden. yudi.pawitan@meb.ki.se}, author = {Pawitan, Y. and Murthy, K. R. and Michiels, S. and Ploner, A. }, citeulike-article-id = {1122674}, issn = {1367-4803}, journal = {Bioinformatics}, keywords = {false\_discovery\_rate, microarray}, month = {October}, number = {20}, pages = {3865--3872}, posted-at = {2007-10-15 12:12:48}, priority = {5}, title = {Bias in the estimation of false discovery rate in microarray studies.}, url = {http://view.ncbi.nlm.nih.gov/pubmed/16105901}, volume = {21}, year = {2005} }

TY - JOUR ID - citeulike:1122674 L3 - citeulike-article-id:1122674 TI - Bias in the estimation of false discovery rate in microarray studies. JF - Bioinformatics VL - 21 IS - 20 SP - 3865 EP - 3872 AD - Department of Medical Epidemiology and Biostatistics, Karolinska Institutet, Stockholm, Sweden. yudi.pawitan@meb.ki.se SN - 1367-4803 N2 - MOTIVATION: The false discovery rate (FDR) provides a key statistical assessment for microarray studies. Its value depends on the proportion pi(0) of non-differentially expressed (non-DE) genes. In most microarray studies, many genes have small effects not easily separable from non-DE genes. As a result, current methods often overestimate pi(0) and FDR, leading to unnecessary loss of power in the overall analysis. METHODS: For the common two-sample comparison we derive a natural mixture model of the test statistic and an explicit bias formula in the standard estimation of pi(0). We suggest an improved estimation of pi(0) based on the mixture model and describe a practical likelihood-based procedure for this purpose. RESULTS: The analysis shows that a large bias occurs when pi(0) is far from 1 and when the non-centrality parameters of the distribution of the test statistic are near zero. The theoretical result also explains substantial discrepancies between non-parametric and model-based estimates of pi(0). Simulation studies indicate mixture-model estimates are less biased than standard estimates. The method is applied to breast cancer and lymphoma data examples. AVAILABILITY: An R-package OCplus containing functions to compute pi(0) based on the mixture model, the resulting FDR and other operating characteristics of microarray data, is freely available at http://www.meb.ki.se/~yudpaw CONTACT: yudi.pawitan@meb.ki.se and alexander.ploner@meb.ki.se. KW - false_discovery_rate KW - microarray AU - Pawitan, Y AU - Murthy, KR AU - Michiels, S AU - Ploner, A PY - 2005/10/15/ UR - http://view.ncbi.nlm.nih.gov/pubmed/16105901 ER -