Testing equality of means in partially paired data with incompleteness in single response

Abstract

In testing differentially expressed genes between tumor and healthy tissues, data are usually collected in paired form. However, incomplete paired data often occur. While extensive statistical researches exist for paired data with incompleteness in both arms, hardly any recent work can be found on paired data with incompleteness in single arm. This paper aims to fill this gap by proposing some new methods, namely, P-value pooling methods and a nonparametric combination test. Simulation studies are conducted to investigate the performance of the proposed methods in terms of type I error and power at small to moderate sample sizes. A real data set from The Cancer Genome Atlas (TCGA) breast cancer study is analyzed using the proposed methods.

Keywords

Differentially expressed genes RNA sequencing meta-analysis nonparametric test

Get full access to this article

View all access options for this article.

References

Wang

Gerstein

Snyder

. RNA-Seq: a revolutionary tool for transcriptomics. Nat Rev Genet 2009; 10: 57–63.

Mortazavi

Williams

McCue

et al.

Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Meth 2008; 5: 621–628.

Davey

Hohenlohe

Etter

et al.

Genome-wide genetic marker discovery and genotyping using next-generation sequencing. Nat Rev Genet 2011; 12: 499–510.

Wilson

Ramos

Villaseor

et al.

Localization of human BRCA1 and its loss in high-grade, non-inherited breast carcinomas. Nat Genet 1999; 21: 236–240.

Kuriakose

Chen

et al.

Selection and validation of differentially expressed genes in head and neck cancer. Cell Mol Life Sci 2004; 61: 1372–1383.

Lapointe

Higgins

et al.

Gene expression profiling identifies clinically relevant subtypes of prostate cancer. Proc Natl Acad Sci U S A 2004; 101: 811–816.

Feng

Hawes

Stern

et al.

DNA methylation in tumor and matched normal tissues from non-small cell lung cancer patients. Cancer Epidemiol Biomarker Prevent 2008; 17: 645–654.

Finak

Bertos

Pepin

et al.

Stromal gene expression predicts clinical outcome in breast cancer. Nat Med 2008; 14: 518–527.

Kim

Lee

et al.

Statistical methods of translating microarray data into clinically relevant diagnostic information in colorectal cancer. Bioinformatics 2005; 21: 517–528.

10.

Lin

Stivers

. On difference of means with incomplete data. Biometrika 1974; 61: 325–334.

11.

Ekbohm

. On comparing means in the paired case with incomplete data on both responses. Biometrika 1976; 63: 299–304.

12.

Bhoj

. Testing equality of means of correlated variates with missing observations on both responses. Biometrika 1978; 65: 225–228.

13.

Kuan

Huang

. A simple and robust method for partially matched samples using the pvalues pooling approach. Stat Med 2013; 32: 3247–3259.

14.

Samawi

Vogel

. On some nonparametric tests for partially observed correlated data: Proposing new tests. J Stat Theory Appl 2015; 14: 131–155.

15.

Rempala

Looney

. Asymptotic properties of a two sample randomized test for partially dependent data. J Stat Plan Inference 2006; 136: 68–89.

16.

Lim

Liang

et al.

Permutation test for incomplete paired data with application to cDNA microarray data. Computat Stat Data Analysis 2012; 56: 510–521.

17.

Amro

Pauly

. Permuting incomplete paired data: a novel exact and asymptotic correct randomization test. J Stat Comput Simul 2017; 87: 1148–1159.

18.

Guo

Yuan

. A comparative review of methods for comparing means using partially paired data. Stat Meth Med Res 2017; 26: 1323–1340.

19.

Huang

Stern

Zhao

. Transcriptional profiles from paired normal samples offer complementary information on cancer patient survival evidence from TCGA pan-cancer data. Scientific Reports 2016; 6: 20567–20567.

20.

The Cancer Genome Atlas. About TCGA Data, https://tcga-data.nci.nih.gov/docs/publications/tcga/about.html (accessed 10 June 2017)

21.

Ritchie

Phipson

et al.

Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucl Acids Res 2015; 43: e47–e47.

22.

Love

Huber

Anders

. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol 2014; 15: 550–550.

23.

Mehta

Gurland

. Testing equality of means in the presence of correlation. Biometrika 1969; 56: 119–126.

24.

Lin

. Procedures for testing the difference of means with incomplete data. J Am Stat Assoc 1973; 68: 699–703.

25.

Morrison

. A test for equality of means of correlated variates with missing data on one response. Biometrika 1973; 60: 101–105.

26.

Little

RJA

. Inference about means from incomplete multivariate data. Biometrika 1976; 63: 593–604.

27.

Ekbohm

. Comparing means in the paired case with missing data on one response. Biometrika 1976; 63: 169–172.

28.

Mehta

Gurland

. A test for equality of means in the presence of correlation and missing values. Biometrika 1973; 60: 211–213.

29.

Lin

Stivers

. Testing for equality of means with incomplete data on one variable: a Monte Carlo study. J Am Stat Assoc 1975; 70: 190–193.

30.

DeSantis

Bryan

et al.

Breast cancer statistics, 2013. CA: Cancer J Clinicians 2014; 64: 52–62.

31.

Ford

Easton

Stratton

et al.

Genetic heterogeneity and penetrance analysis of the BRCA1 and BRCA2 genes in breast cancer families. Am J Human Genet 1998; 62: 676–689.

32.

Rahman

Seal

Thompson

et al.

PALB2, which encodes a BRCA2-interacting protein, is a breast cancer susceptibility gene. Nat Genet 2007; 39: 165–167.

33.

Harari

Yarden

. Molecular mechanisms underlying ErbB2/HER2 action in breast cancer. Oncogene 2000; 19: 6102–6114.

34.

Munoz

Henderson

Haber

et al.

Role of the MRP1/ABCC1 multidrug transporter protein in cancer. IUBMB Life 2007; 59: 752–757.

35.

De Jong

Nolte

Te Meerman

et al.

Genes other than BRCA1 and BRCA2 involved in breast cancer susceptibility. J Med Genet 2002; 39: 225–242.

36.

Rosenthal

. Combining results of independent studies. Psychol Bull 1978; 85: 185–193.

37.

Bailey

Gribskov

. Combining evidence using p-values: application to sequence homology searches. Bioinformatics 1998; 14: 48–54.

38.

Hartung

. A note on combining dependent tests of significance. Biometric J 1999; 41: 849–855.

39.

Brown

. A method for combining non-independent, one-sided tests of significance. Biometrics 1975; 31: 987–992.

40.

Kost

McDermott

. Combining dependent p-values. Stat Probabil Lett 2002; 60: 183–190.

41.

Hou

. A simple approximation for the distribution of the weighted combination of non-independent or independent probabilities. Stat Probab Lett 2005; 73: 179–187.

42.

Zaykin

. Optimally weighted Z-test is a powerful method for combining probabilities in meta-analysis. J Evolution Biol 2011; 24: 1836–1841.