Lifetime Data Analysis最新文献

Multi-state models are used to describe how individuals transition through different states over time. The distribution of the time spent in different states, referred to as 'length of stay', is often of interest. Methods for estimating expected length of stay in a given state are well established. The focus of this paper is on the distribution of the time spent in different states conditional on the complete pathway taken through the states, which we call 'conditional length of stay'. This work is motivated by questions about length of stay in hospital wards and intensive care units among patients hospitalised due to Covid-19. Conditional length of stay estimates are useful as a way of summarising individuals' transitions through the multi-state model, and also as inputs to mathematical models used in planning hospital capacity requirements. We describe non-parametric methods for estimating conditional length of stay distributions in a multi-state model in the presence of censoring, including conditional expected length of stay (CELOS). Methods are described for an illness-death model and then for the more complex motivating example. The methods are assessed using a simulation study and shown to give unbiased estimates of CELOS, whereas naive estimates of CELOS based on empirical averages are biased in the presence of censoring. The methods are applied to estimate conditional length of stay distributions for individuals hospitalised due to Covid-19 in the UK, using data on 42,980 individuals hospitalised from March to July 2020 from the COVID19 Clinical Information Network.

Over the last decades, the challenges in survival models have been changing considerably and full probabilistic modeling is crucial in many medical applications. Motivated from a new biological interpretation of cancer metastasis, we introduce a general method for obtaining more flexible cure rate models. The proposal model extended the promotion time cure rate model. Furthermore, it includes several well-known models as special cases and defines many new special models. We derive several properties of the hazard function for the proposed model and establish mathematical relationships with the promotion time cure rate model. We consider a frequentist approach to perform inferences, and the maximum likelihood method is employed to estimate the model parameters. Simulation studies are conducted to evaluate its performance with a discussion of the obtained results. A real dataset from population-based study of incident cases of melanoma diagnosed in the state of São Paulo, Brazil, is discussed in detail.

For clinical trial design and analysis, there has been extensive work related to using joint models for longitudinal and time-to-event data without a cure fraction (i.e., when all patients are at risk for the event of interest), but comparatively little treatment has been given to design and analysis of clinical trials using joint models that incorporate a cure fraction. In this paper, we develop a Bayesian clinical trial design methodology focused on evaluating the treatment's effect on a time-to-event endpoint using a promotion time cure rate model, where the longitudinal process is incorporated into the hazard model for the promotion times. A piecewise linear hazard model for the period after assessment of the longitudinal measure ends is proposed as an alternative to extrapolating the longitudinal trajectory. This may be advantageous in scenarios where the period of time from the end of longitudinal measurements until the end of observation is substantial. Inference for the time-to-event endpoint is based on a novel estimand which combines the treatment's effect on the probability of cure and its effect on the promotion time distribution, mediated by the longitudinal outcome. We propose an approach for sample size determination such that the design has a high power and a well-controlled type I error rate with both operating characteristics defined from a Bayesian perspective. We demonstrate the methodology by designing a breast cancer clinical trial with a primary time-to-event endpoint where longitudinal outcomes are measured periodically during follow up.

We propose an inferential procedure for additive hazards regression with high-dimensional survival data, where the covariates are prone to measurement errors. We develop a double bias correction method by first correcting the bias arising from measurement errors in covariates through an estimating function for the regression parameter. By adopting the convex relaxation technique, a regularized estimator for the regression parameter is obtained by elaborately designing a feasible loss based on the estimating function, which is solved via linear programming. Using the Neyman orthogonality, we propose an asymptotically unbiased estimator which further corrects the bias caused by the convex relaxation and regularization. We derive the convergence rate of the proposed estimator and establish the asymptotic normality for the low-dimensional parameter estimator and the linear combination thereof, accompanied with a consistent estimator for the variance. Numerical experiments are carried out on both simulated and real datasets to demonstrate the promising performance of the proposed double bias correction method.

The proportional hazards (PH) model is, arguably, the most popular model for the analysis of lifetime data arising from epidemiological studies, among many others. In such applications, analysts may be faced with censored outcomes and/or studies which institute enrollment criterion leading to left truncation. Censored outcomes arise when the event of interest is not observed but rather is known relevant to an observation time(s). Left truncated data occur in studies that exclude participants who have experienced the event prior to being enrolled in the study. If not accounted for, both of these features can lead to inaccurate inferences about the population under study. Thus, to overcome this challenge, herein we propose a novel unified PH model that can be used to accommodate both of these features. In particular, our approach can seamlessly analyze exactly observed failure times along with interval-censored observations, while aptly accounting for left truncation. To facilitate model fitting, an expectation-maximization algorithm is developed through the introduction of carefully structured latent random variables. To provide modeling flexibility, a monotone spline representation is used to approximate the cumulative baseline hazard function. The performance of our methodology is evaluated through a simulation study and is further illustrated through the analysis of two motivating data sets; one that involves child mortality in Nigeria and the other prostate cancer.

Quantile residual lifetime (QRL) is of significant interest in many clinical studies as an easily interpretable quantity compared to other summary measures of survival distributions. In cancer or other chronic diseases, treatments are often compared based on the distributions or quantiles of the residual lifetime. Thus a common problem of interest is to test the equality of the QRL between two populations. In this paper, we propose two classes of tests to compare two QRLs; one class is based on the difference between two estimated QRLs, and the other is based on the estimating function of the QRL, where the estimated QRL from one sample is plugged into the QRL-estimating-function of the other sample. We outline the asymptotic properties of these test statistics. Simulation studies demonstrate that the proposed tests produced Type I errors closer to the nominal level and are superior to some existing tests based on both Type I error and power. Our proposed test statistics are also computationally less intensive and more straightforward compared to tests based on the confidence intervals. We applied the proposed methods to a randomized multicenter phase III trial for breast cancer patients.

The incubation period is a key characteristic of an infectious disease. In the outbreak of a novel infectious disease, accurate evaluation of the incubation period distribution is critical for designing effective prevention and control measures . Estimation of the incubation period distribution based on limited information from retrospective inspection of infected cases is highly challenging due to censoring and truncation. In this paper, we consider a semiparametric regression model for the incubation period and propose a sieve maximum likelihood approach for estimation based on the symptom onset time, travel history, and basic demographics of reported cases. The approach properly accounts for the pandemic growth and selection bias in data collection. We also develop an efficient computation method and establish the asymptotic properties of the proposed estimators. We demonstrate the feasibility and advantages of the proposed methods through extensive simulation studies and provide an application to a dataset on the outbreak of COVID-19.

We develop an efficient algorithm to compute the likelihood of the phase-type ageing model. The proposed algorithm uses the uniformisation method to stabilise the numerical calculation. It also utilises a vectorised formula to only calculate the necessary elements of the probability distribution. Our algorithm, with an error's upper bound, could be adjusted easily to tackle the likelihood calculation of the Coxian models. Furthermore, we compare the speed and the accuracy of the proposed algorithm with those of the traditional method using the matrix exponential. Our algorithm is faster and more accurate than the traditional method in calculating the likelihood. Based on our experiments, we recommend using 20 sets of randomly-generated initial values for the optimisation to get a reliable estimate for which the evaluated likelihood is close to the maximum likelihood.

Recurrent event and failure time data arise frequently in many clinical and observational studies. In this article, we propose a joint modeling of generalized scale-change models for the recurrent event process and the failure time, and allow the two processes to be correlated through a shared frailty. The proposed joint model is flexible in that it requires neither the Poisson assumption for the recurrent event process nor a parametric assumption on the frailty distribution. Estimating equation approaches are developed for parameter estimation, and the asymptotic properties of the resulting estimators are established. Simulation studies are conducted to evaluate the finite sample performances of the proposed method. An application to a medical cost study of chronic heart failure patients is provided.

An important complexity in censored data is that only partial information on the variables of interest is observed. In recent years, a large family of asymmetric distributions and maximum likelihood estimation for the parameters in that family has been studied, in the complete data case. In this paper, we exploit the appealing family of quantile-based asymmetric distributions to obtain flexible distributions for modelling right censored survival data. The flexible distributions can be generated using a variety of symmetric distributions and monotonic link functions. The interesting feature of this family is that the location parameter coincides with an index-parameter quantile of the distribution. This family is also suitable to characterize different shapes of the hazard function (constant, increasing, decreasing, bathtub and upside-down bathtub or unimodal shapes). Statistical inference is done for the whole family of distributions. The parameter estimation is carried out by optimizing a non-differentiable likelihood function. The asymptotic properties of the estimators are established. The finite-sample performance of the proposed method and the impact of censorship are investigated via simulations. Finally, the methodology is illustrated on two real data examples (times to weaning in breast-fed data and German Breast Cancer data).