Probability fold change: a robust computational approach for identifying differentially expressed gene lists.
Identifying genes that are differentially expressed under different experimental conditions is a fundamental task in microarray studies. However, different ranking methods generate very different gene lists, and this could profoundly impact follow-up analyses and biological interpretation. Therefore, developing improved ranking methods are critical in microarray data analysis. We developed a new algorithm, the probabilistic fold change (PFC), which ranks genes based on a confidence interval estimate of fold change. We performed extensive testing using multiple benchmark data sources including the MicroArray Quality Control (MAQC) data sets. We corroborated our observations with MAQC data sets using qRT-PCR data sets and Latin square spike-in data sets. Along with PFC, we tested six other popular ranking algorithms including Mean Fold Change (FC), SAM, t-statistic (T), Bayesian-t (BAYT), Intensity-Conditional Fold Change (CFC), and Rank Product (RP). PFC achieved reproducibility and accuracy that are consistently among the best of the seven ranking algorithms while other ranking algorithms would show weakness in some cases. Contrary to common belief, our results demonstrated that statistical accuracy will not translate to biological reproducibility and therefore both quality aspects need to be evaluated.
Transcriptional Genomics Core, Cedars-Sinai Medical Center, David Geffen School of Medicine at UCLA, Los Angeles, CA 90048, USA. dengx@ucla.edu
Xutao
X
Deng
Jun
J
Xu
James
J
Hui
Charles
C
Wang
D000465:Algorithms; D000818:Animals; D001499:Bayes Theorem; D001699:Biometry; D016001:Confidence Intervals; D003627:Data Interpretation, Statistical; D030541:Databases, Genetic; D015870:Gene Expression; D006801:Humans; D008957:Models, Genetic; D015233:Models, Statistical; D020411:Oligonucleotide Array Sequence Analysis / Q000592:standards / Q000706:statistics & numerical data*; D011786:Quality Control; D015203:Reproducibility of Results; D020133:Reverse Transcriptase Polymerase Chain Reaction / Q000706:statistics & numerical data; D018401:Sample Size