1
|
Liu L, Liu Y, Xiong Y, Hu XJ. Cox regression of clustered event times with covariates missing not at random. Scand Stat Theory Appl 2019. [DOI: 10.1111/sjos.12409] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Li Liu
- School of Mathematics and StatisticsWuhan University Wuhan China
| | - Yanyan Liu
- School of Mathematics and StatisticsWuhan University Wuhan China
| | - Yi Xiong
- Department of Statistics and Actuarial ScienceSimon Fraser University Burnaby Canada
| | - X. Joan Hu
- Department of Statistics and Actuarial ScienceSimon Fraser University Burnaby Canada
| |
Collapse
|
2
|
Park K, Qiu P. Comparing crossing hazard rate functions by joint modelling of survival and longitudinal data. J STAT COMPUT SIM 2019. [DOI: 10.1080/00949655.2019.1668392] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
Affiliation(s)
- Kayoung Park
- Department of Mathematics and Statistics, Old Dominion University, Norfolk, VA, USA
| | - Peihua Qiu
- Department of Biostatistics, University of Florida, Gainesville, FL, USA
| |
Collapse
|
3
|
Affiliation(s)
- Lu You
- Department of Biostatistics, University of Florida, Gainesville, FL
| | - Peihua Qiu
- Department of Biostatistics, University of Florida, Gainesville, FL
| |
Collapse
|
4
|
Ahn S, Lim J, Paik MC, Sacco RL, Elkind MS. Cox model with interval-censored covariate in cohort studies. Biom J 2018; 60:797-814. [PMID: 29775990 DOI: 10.1002/bimj.201700090] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/30/2017] [Revised: 12/19/2017] [Accepted: 02/27/2018] [Indexed: 11/07/2022]
Abstract
In cohort studies the outcome is often time to a particular event, and subjects are followed at regular intervals. Periodic visits may also monitor a secondary irreversible event influencing the event of primary interest, and a significant proportion of subjects develop the secondary event over the period of follow-up. The status of the secondary event serves as a time-varying covariate, but is recorded only at the times of the scheduled visits, generating incomplete time-varying covariates. While information on a typical time-varying covariate is missing for entire follow-up period except the visiting times, the status of the secondary event are unavailable only between visits where the status has changed, thus interval-censored. One may view interval-censored covariate of the secondary event status as missing time-varying covariates, yet missingness is partial since partial information is provided throughout the follow-up period. Current practice of using the latest observed status produces biased estimators, and the existing missing covariate techniques cannot accommodate the special feature of missingness due to interval censoring. To handle interval-censored covariates in the Cox proportional hazards model, we propose an available-data estimator, a doubly robust-type estimator as well as the maximum likelihood estimator via EM algorithm and present their asymptotic properties. We also present practical approaches that are valid. We demonstrate the proposed methods using our motivating example from the Northern Manhattan Study.
Collapse
Affiliation(s)
- Soohyun Ahn
- Department of Mathematics, Ajou University, Suwon, Korea
| | - Johan Lim
- Department of Statistics, Seoul National University, Seoul, Korea
| | | | - Ralph L Sacco
- Department of Neurology, Miller School of Medicine, University of Miami, Miami, FL, USA
| | | |
Collapse
|
5
|
Portier F, El Ghouch A, Van Keilegom I. Efficiency and bootstrap in the promotion time cure model. BERNOULLI 2017. [DOI: 10.3150/16-bej852] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
|
6
|
Su YR, Wang JL. SEMIPARAMETRIC EFFICIENT ESTIMATION FOR SHARED-FRAILTY MODELS WITH DOUBLY-CENSORED CLUSTERED DATA. Ann Stat 2016. [PMID: 29527068 DOI: 10.1214/15-aos1406] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
Abstract
In this paper, we investigate frailty models for clustered survival data that are subject to both left- and right-censoring, termed "doubly-censored data". This model extends current survival literature by broadening the application of frailty models from right-censoring to a more complicated situation with additional left censoring. Our approach is motivated by a recent Hepatitis B study where the sample consists of families. We adopt a likelihood approach that aims at the nonparametric maximum likelihood estimators (NPMLE). A new algorithm is proposed, which not only works well for clustered data but also improve over existing algorithm for independent and doubly-censored data, a special case when the frailty variable is a constant equal to one. This special case is well known to be a computational challenge due to the left censoring feature of the data. The new algorithm not only resolves this challenge but also accommodate the additional frailty variable effectively. Asymptotic properties of the NPMLE are established along with semi-parametric efficiency of the NPMLE for the finite-dimensional parameters. The consistency of Bootstrap estimators for the standard errors of the NPMLE is also discussed. We conducted some simulations to illustrate the numerical performance and robustness of the proposed algorithm, which is also applied to the Hepatitis B data.
Collapse
Affiliation(s)
- Yu-Ru Su
- Biostatistics and Biomathematics, Public Health Science Division, Fred Hutchinson Cancer Research Center, Seattle, 98103, U.S.A
| | - Jane-Ling Wang
- Department of Statistics, University of California, Davis, California, 95616, U.S.A
| |
Collapse
|
7
|
Cao H, Churpek MM, Zeng D, Fine JP. Analysis of the Proportional Hazards Model with Sparse Longitudinal Covariates. J Am Stat Assoc 2015; 110:1187-1196. [PMID: 26576066 PMCID: PMC4643320 DOI: 10.1080/01621459.2014.957289] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
Abstract
Regression analysis of censored failure observations via the proportional hazards model permits time-varying covariates which are observed at death times. In practice, such longitudinal covariates are typically sparse and only measured at infrequent and irregularly spaced follow-up times. Full likelihood analyses of joint models for longitudinal and survival data impose stringent modelling assumptions which are difficult to verify in practice and which are complicated both inferentially and computationally. In this article, a simple kernel weighted score function is proposed with minimal assumptions. Two scenarios are considered: half kernel estimation in which observation ceases at the time of the event and full kernel estimation for data where observation may continue after the event, as with recurrent events data. It is established that these estimators are consistent and asymptotically normal. However, they converge at rates which are slower than the parametric rates which may be achieved with fully observed covariates, with the full kernel method achieving an optimal convergence rate which is superior to that of the half kernel method. Simulation results demonstrate that the large sample approximations are adequate for practical use and may yield improved performance relative to last value carried forward approach and joint modelling method. The analysis of the data from a cardiac arrest study demonstrates the utility of the proposed methods.
Collapse
Affiliation(s)
- Hongyuan Cao
- Department of Statistics, University of Missouri-Columbia, Columbia, MO, 65201
| | - Mathew M Churpek
- Department of Health Studies and Department of Medicine, University of Chicago, Chicago, IL, 60637
| | - Donglin Zeng
- Department of Biostatistics, University of North Carolina at Chapel Hill, Chapel Hill, NC, 27514
| | - Jason P Fine
- Department of Biostatistics, University of North Carolina at Chapel Hill, Chapel Hill, NC, 27514
| |
Collapse
|
8
|
Park KY, Qiu P. Model selection and diagnostics for joint modeling of survival and longitudinal data with crossing hazard rate functions. Stat Med 2014; 33:4532-46. [PMID: 25043230 DOI: 10.1002/sim.6259] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/22/2013] [Revised: 05/09/2014] [Accepted: 06/10/2014] [Indexed: 11/11/2022]
Abstract
Comparison of two hazard rate functions is important for evaluating treatment effect in studies concerning times to some important events. In practice, it may happen that the two hazard rate functions cross each other at one or more unknown time points, representing temporal changes of the treatment effect. Also, besides survival data, there could be longitudinal data available regarding some time-dependent covariates. When jointly modeling the survival and longitudinal data in such cases, model selection and model diagnostics are especially important to provide reliable statistical analysis of the data, which are lacking in the literature. In this paper, we discuss several criteria for assessing model fit that have been used for model selection and apply them to the joint modeling of survival and longitudinal data for comparing two crossing hazard rate functions. We also propose hypothesis testing and graphical methods for model diagnostics of the proposed joint modeling approach. Our proposed methods are illustrated by a simulation study and by a real-data example concerning two early breast cancer treatments.
Collapse
Affiliation(s)
- Ka Young Park
- Department of Biostatistics, University of Florida, Gainesville, FL 32610, U.S.A
| | | |
Collapse
|
9
|
|
10
|
Su YR, Wang JL. MODELING LEFT-TRUNCATED AND RIGHT-CENSORED SURVIVAL DATA WITH LONGITUDINAL COVARIATES. Ann Stat 2012; 40:1465-1488. [PMID: 29479122 DOI: 10.1214/12-aos996] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
Abstract
There is a surge in medical follow-up studies that include longitudinal covariates in the modeling of survival data. So far, the focus has been largely on right censored survival data. We consider survival data that are subject to both left truncation and right censoring. Left truncation is well known to produce biased sample. The sampling bias issue has been resolved in the literature for the case which involves baseline or time-varying covariates that are observable. The problem remains open however for the important case where longitudinal covariates are present in survival models. A joint likelihood approach has been shown in the literature to provide an effective way to overcome those difficulties for right censored data, but this approach faces substantial additional challenges in the presence of left truncation. Here we thus propose an alternative likelihood to overcome these difficulties and show that the regression coefficient in the survival component can be estimated unbiasedly and efficiently. Issues about the bias for the longitudinal component are discussed. The new approach is illustrated numerically through simulations and data from a multi-center AIDS cohort study.
Collapse
Affiliation(s)
- Yu-Ru Su
- Graduate Group in Biostatistics, University of California, Davis, California, 95616, U.S.A
| | - Jane-Ling Wang
- Department of Statistics, University of California, Davis, California, 95616, U.S.A
| |
Collapse
|
11
|
Hernandez-Quintero A, Dupuy JF, Escarela G. Analysis of a semiparametric mixture model for competing risks. ANN I STAT MATH 2009. [DOI: 10.1007/s10463-009-0229-1] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
|
12
|
Estimation of the mean quality-adjusted survival using a multistate model for the sojourn times. J Stat Plan Inference 2008. [DOI: 10.1016/j.jspi.2007.10.010] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
|