AUTHOR=Kang Guolian , Mirzaei Sedigheh S. , Zhang Hui , Zhu Liang , Rai Shesh N. , Srivastava Deo Kumar 

TITLE=Robust Behrens–Fisher Statistic Based on Trimmed Means and Its Usefulness in Analyzing High-Throughput Data

JOURNAL=Frontiers in Systems Biology

VOLUME=2

YEAR=2022

URL=https://www.frontiersin.org/journals/systems-biology/articles/10.3389/fsysb.2022.877601

DOI=10.3389/fsysb.2022.877601

ISSN=2674-0702

ABSTRACT=<p>In the context of high-throughput data, the differences in continuous markers between two groups are usually assessed by ordering the p-values obtained from the two-sample pooled <italic>t</italic>-test or Wilcoxon–Mann–Whitney test and choosing a stringent cutoff such as 10<sup>–8</sup> to control the family-wise error rate <inline-formula id="inf1"><mml:math id="m1" xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mrow><mml:mo>(</mml:mo><mml:mi>F</mml:mi><mml:mi>W</mml:mi><mml:mi>E</mml:mi><mml:mi>R</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:math></inline-formula> or false discovery rate <inline-formula id="inf2"><mml:math id="m2" xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mi>F</mml:mi><mml:mi>D</mml:mi><mml:mi>R</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>.</mml:mo></mml:mrow></mml:math></inline-formula> All markers with p-values below the cutoff are declared to be significantly associated with the phenotype. This inherently assumes that the test procedure provides valid type I error estimates in extreme tails of the null distribution. The aforementioned tests assume homoscedasticity in the two groups, and the <italic>t</italic>-test further assumes underlying distributions to be normally distributed. Cao et al. (Biometrika, 2013, 100, 495–502) have shown that in the context of multiple hypotheses testing the approach based on <inline-formula id="inf3"><mml:math id="m3" xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mrow><mml:mi>F</mml:mi><mml:mi>D</mml:mi><mml:mi>R</mml:mi></mml:mrow></mml:math></inline-formula> may not be valid under non-normality and/or heteroscedasticity. Therefore, having a test statistic that is robust to these violations is needed. In this study, we propose a robust analog of Behrens–Fisher statistic based on trimmed means, conduct an extensive simulation study to compare its performance with other competing approaches, and demonstrate its usefulness by applying it to DNA methylation data used by Teschendorff et al. (Genome Res., 2010, 20, 440–446). An R program to implement the proposed method is provided in the Supplementary Material.</p>