Parameter estimation in Cox models with missing failure indicators and the OPPERA study

In a prospective cohort study, examining all participants for incidence of the condition of interest may be prohibitively expensive. For example, the "gold standard" for diagnosing temporomandibular disorder (TMD) is a physical examination by a trained clinician. In large studies, examining all participants in this manner is infeasible. Instead, it is common to use questionnaires to screen for incidence of TMD and perform the "gold standard" examination only on participants who screen positively. Unfortunately, some participants may leave the study before receiving the "gold standard" examination. Within the framework of survival analysis, this results in missing failure indicators. Motivated by the Orofacial Pain: Prospective Evaluation and Risk Assessment (OPPERA) study, a large cohort study of TMD, we propose a method for parameter estimation in survival models with missing failure indicators. We estimate the probability of being an incident case for those lacking a "gold standard" examination using logistic regression. These estimated probabilities are used to generate multiple imputations of case status for each missing examination that are combined with observed data in appropriate regression models. The variance introduced by the procedure is estimated using multiple imputation. The method can be used to estimate both regression coefficients in Cox proportional hazard models as well as incidence rates using Poisson regression. We simulate data with missing failure indicators and show that our method performs as well as or better than competing methods. Finally, we apply the proposed method to data from the OPPERA study.

[1]  Roderick J. A. Little,et al.  Statistical Analysis with Missing Data: Little/Statistical Analysis with Missing Data , 2002 .

[2]  R. Ohrbach,et al.  Study protocol, sample characteristics, and loss to follow-up: the OPPERA prospective cohort study. , 2013, The journal of pain : official journal of the American Pain Society.

[3]  Eric Bair,et al.  Preclinical episodes of orofacial pain symptoms and their association with health care behaviors in the OPPERA prospective cohort study , 2013, PAIN®.

[4]  Eric Bair,et al.  Potential genetic risk factors for chronic TMD: genetic associations from the OPPERA case control study. , 2011, The journal of pain : official journal of the American Pain Society.

[5]  Ping Chen,et al.  Regression analysis of right-censored failure time data with missing censoring indicators , 2009 .

[6]  D. Rubin INFERENCE AND MISSING DATA , 1975 .

[7]  Sundarraman Subramanian,et al.  Efficient estimation of regression coefficients and baseline hazard under proportionality of conditional hazards , 2000 .

[8]  Robert Tibshirani,et al.  An Introduction to the Bootstrap , 1994 .

[9]  Amalia S Magaret,et al.  Incorporating validation subsets into discrete proportional hazards models for mismeasured outcomes , 2008, Statistics in medicine.

[10]  William Maixner,et al.  Orofacial pain prospective evaluation and risk assessment study--the OPPERA study. , 2011, The journal of pain : official journal of the American Pain Society.

[11]  Ralf Bender,et al.  Generating survival times to simulate Cox proportional hazards models , 2005, Statistics in medicine.

[12]  S. Dworkin,et al.  Research diagnostic criteria for temporomandibular disorders: review, criteria, examinations and specifications, critique. , 1992, Journal of craniomandibular disorders : facial & oral pain.

[13]  Eric Bair,et al.  Study methods, recruitment, sociodemographic findings, and demographic representativeness in the OPPERA study. , 2011, The journal of pain : official journal of the American Pain Society.

[14]  L. Magder,et al.  Logistic regression when the outcome is measured with uncertainty. , 1997, American journal of epidemiology.

[15]  Eric Bair,et al.  Potential autonomic risk factors for chronic TMD: descriptive data and empirically identified domains from the OPPERA case-control study. , 2011, The journal of pain : official journal of the American Pain Society.

[16]  R. Ohrbach,et al.  Clinical findings and pain symptoms as potential risk factors for chronic TMD: descriptive data and empirically identified domains from the OPPERA case-control study. , 2011, The journal of pain : official journal of the American Pain Society.

[17]  Sundarraman Subramanian,et al.  Product‐limit Estimators and Cox Regression with Missing Censoring Information , 1998 .

[18]  D. Rubin Multiple Imputation After 18+ Years , 1996 .

[19]  Suman Bhattacharya,et al.  An Audit Strategy for Progression‐Free Survival , 2011, Biometrics.

[20]  R. Ohrbach,et al.  Potential psychosocial risk factors for chronic TMD: descriptive data and empirically identified domains from the OPPERA case-control study. , 2011, The journal of pain : official journal of the American Pain Society.

[21]  Eric Bair,et al.  Pain sensitivity risk factors for chronic TMD: descriptive data and empirically identified domains from the OPPERA case control study. , 2011, The journal of pain : official journal of the American Pain Society.

[22]  Zhiliang Ying,et al.  Non- and semi-parametric analysis of failure time data with missing failure indicators , 2007 .

[23]  D. Rubin,et al.  Statistical Analysis with Missing Data. , 1989 .

[24]  Michael R. Kosorok,et al.  Analysis of Time-to-Event Data With Incomplete Event Adjudication , 2004 .