Lifetime Data Analysis最新文献

We propose a new goodness-of-fit procedure designed to verify the reciprocal property that characterizes the fatigue-life or Birnbaum-Saunders (BS) distribution. Under this property, scaling a random variable that takes positive values by its median results in the same distribution as its reciprocal, a feature frequently encountered in reliability and survival studies. Our procedure employs total time on test (TTT) curves to compare the behavior of the observed data and its reciprocal counterpart, capturing both local and global discrepancies through supremum- and area-based statistics. We establish the theoretical validity of these statistics under mild assumptions, showing that they deliver accurate inference for moderate to large samples. Simulation evidence indicates that our TTT-based procedures are sensitive to subtle departures from log-symmetry, particularly when the distribution underlying the data has heavier or lighter tails than the assumed one. Illustrative real data examples further reveal how overlooking deviations from the reciprocal property can distort reliability estimates and predictions of failure times, showing the practical importance of the new goodness-of-fit procedure. Overall, our findings strengthen the BS framework and provide robust tools for model validation and selection when log-symmetric modeling assumptions are in place.

Double truncation commonly occurs in astronomy, epidemiology and economics. Compared to one-sided truncation, double truncation, which combines both left and right truncation, is more challenging to handle and the methods for analyzing doubly truncated data are limited. For the situation, a common approach is to perform conditional analysis conditional on truncation times, which is simple but may not be efficient. Corresponding to this, we propose a pairwise pseudo-likelihood approach that aims to recover some information missed in the conditional methods and can yield more efficient estimation. The resulting estimator is shown to be consistent and asymptotically normal. An extensive simulation study indicates that the proposed procedure works well in practice and is indeed more efficient than the conditional approach. The proposed methodology applied to an AIDS study.

Here we revisit a goodness-of-fit testing problem for randomly right-censored data in the presence of cured subjects, i.e. the population consists of two parts: the cured or non-susceptible group, who will never experience the event of interest versus those who will undergo the event of interest when followed up sufficiently long. We consider the modifications of proposed characterization-based goodness-of-fit tests for the exponential distribution constructed via the inverse probability of censoring weighted U- or V-approach. We present their asymptotic properties and extend our discussion to encompass suitable generalizations applicable to a variety of tests formulated using the same methodology. A comparative power study of these proposed tests against a recent CvM-based competitor and the modifications of the most prominent competitors identified in prior studies that did not consider the presence of cured subjects, demonstrates good finite sample performance. Novel tests are illustrated on a real dataset related to leukemia relapse.

In longitudinal observational studies with time-to-event outcomes, a common objective in causal analysis is to estimate the causal survival curve under hypothetical intervention scenarios. The g-formula is a useful tool for this analysis. To enhance the traditional parametric g-formula, we developed an alternative g-formula estimator, which incorporates the Bayesian Additive Regression Trees into the modeling of the time-evolving generative components, aiming to mitigate the bias due to model misspecification. We focus on binary time-varying treatments and introduce a general class of g-formulas for discrete survival data that can incorporate longitudinal balancing scores. The minimum sufficient formulation of these longitudinal balancing scores is linked to the nature of treatment strategies, i.e., static or dynamic. For each type of treatment strategy, we provide posterior sampling algorithms. We conducted simulations to illustrate the empirical performance of the proposed method and demonstrate its practical utility using data from the Yale New Haven Health System's electronic health records.

Survival data is doubly truncated when only participants who experience an event during a random interval are included in the sample. Existing methods typically correct for double truncation bias in Cox regression through inverse probability weighting via the nonparametric maximum likelihood estimate (NPMLE) of the selection probabilities. This approach relies on two key assumptions, quasi-independent truncation and positivity of the sampling probabilities, yet there are no methods available to thoroughly assess these assumptions in the regression context. Furthermore, these estimators can be particularly sensitive to extreme event times. Finally, current double truncation methods rely on bootstrapping for variance estimation. Aside from the unnecessary computational burden, there are often identifiability issues with the NPMLE during bootstrap resampling. To address these limitations of current methods, we propose a class of robust Cox regression coefficient estimators with time-varying inverse probability weights and extend these estimators to conduct sensitivity analysis regarding possible non-positivity of the sampling probabilities. Also, we develop a nonparametric test and graphical diagnostic for verifying the quasi-independent truncation assumption. Finally, we provide closed-form standard errors for the NPMLE as well as for the proposed estimators. The proposed estimators are evaluated through extensive simulations and illustrated using an AIDS study.

The published version of the manuscript (D´iaz, Hoffman, Hejazi Lifetime Data Anal 30, 213-236, 2024) contained an error (We would like to thank Kara Rudolph for pointing out an issue that led to uncovering the error)) in the definition of the outcome that had cascading effects and created errors in the definition of multiple objects in the paper. We correct those errors here. For completeness, we reproduce the entire manuscript, underlining places where we made a correction.Longitudinal modified treatment policies (LMTP) have been recently developed as a novel method to define and estimate causal parameters that depend on the natural value of treatment. LMTPs represent an important advancement in causal inference for longitudinal studies as they allow the non-parametric definition and estimation of the joint effect of multiple categorical, ordinal, or continuous treatments measured at several time points. We extend the LMTP methodology to problems in which the outcome is a time-to-event variable subject to a competing event that precludes observation of the event of interest. We present identification results and non-parametric locally efficient estimators that use flexible data-adaptive regression techniques to alleviate model misspecification bias, while retaining important asymptotic properties such as $\sqrt{n}$ -consistency. We present an application to the estimation of the effect of the time-to-intubation on acute kidney injury amongst COVID- 19 hospitalized patients, where death by other causes is taken to be the competing event.

We propose a two-stage estimation procedure for a copula-based model with semi-competing risks data, where the non-terminal event is subject to dependent censoring by the terminal event, and both events are subject to independent censoring. With a copula-based model, the marginal survival functions of individual event times are specified by semiparametric transformation models, and the dependence between the bivariate event times is specified by a parametric copula function. For the estimation procedure, in the first stage, the parameters associated with the marginal of the terminal event are estimated using only the corresponding observed outcomes, and in the second stage, the marginal parameters for the non-terminal event time and the copula parameter are estimated together via maximizing a pseudo-likelihood function based on the joint distribution of the bivariate event times. We derived the asymptotic properties of the proposed estimator and provided an analytic variance estimator for inference. Through simulation studies, we showed that our approach leads to consistent estimates with less computational cost and more robustness than the one-stage procedure developed in Chen YH (Lifetime Data Anal 18:36-57, 2012), where all parameters were estimated simultaneously. In addition, our approach demonstrates more desirable finite-sample performances over another existing two-stage estimation method proposed in Zhu H et al., (Commu Statistics-Theory Methods 51(22):7830-7845, 2021) . An R package PMLE4SCR is developed to implement our proposed method.

Recurrent event data with a terminal event arise in follow-up studies. The current literature has primarily focused on the effect of covariates on the recurrent event process using marginal estimating equation approaches or joint modeling approaches via frailties. In this article, we propose a conditional model for recurrent event data with a terminal event, which provides an intuitive interpretation of the effect of the terminal event: at an early time, the rate of recurrent events is nearly independent of the terminal event, but the dependence gets stronger as time goes close to the terminal event time. A two-stage likelihood-based approach is proposed to estimate parameters of interest. Asymptotic properties of the estimators are established. The finite-sample behavior of the proposed method is examined through simulation studies. A real data of colorectal cancer is analyzed by the proposed method for illustration.

Recurrent events are common in medical practice or epidemiologic studies when each subject experiences a particular event repeatedly over time. In some long-term observations of recurrent events, a terminal event such as death may exist in recurrent event data. Meanwhile, some inspected subjects will withdraw from a study for some time for various reasons and then resume, which may happen more than once. The period between the subject leaving and returning to the study is called an intermittent gap. One naive method typically ignores gaps and treats the events as usual recurrent events, which could result in misleading estimation results. In this article, we consider a semiparametric proportional rates model for recurrent event data with intermittent gaps and a terminal event. An estimation procedure is developed for the model parameters, and the asymptotic properties of the resulting estimators are established. Simulation studies demonstrate that the proposed estimators perform satisfactorily compared to the naive method that ignores gaps. A diabetes study further shows the utility of the proposed method.

Based on the expectile loss function and the adaptive LASSO penalty, the paper proposes and studies the estimation methods for the accelerated failure time (AFT) model. In this approach, we need to estimate the survival function of the censoring variable by the Kaplan-Meier estimator. The AFT model parameters are first estimated by the expectile method and afterwards, when the number of explanatory variables can be large, by the adaptive LASSO expectile method which directly carries out the automatic selection of variables. We also obtain the convergence rate and asymptotic normality for the two estimators, while showing the sparsity property for the censored adaptive LASSO expectile estimator. A numerical study using Monte Carlo simulations confirms the theoretical results and demonstrates the competitive performance of the two proposed estimators. The usefulness of these estimators is illustrated by applying them to three survival data sets.