Multivariate frailty models for two types of recurrent events with a dependent terminal event: Application to breast cancer data

Individuals may experience more than one type of recurrent event and a terminal event during the life course of a disease. Follow‐up may be interrupted for several reasons, including the end of a study, or patients lost to follow‐up, which are noninformative censoring events. Death could also stop the follow‐up, hence, it is considered as a dependent terminal event. We propose a multivariate frailty model that jointly analyzes two types of recurrent events with a dependent terminal event. Two estimation methods are proposed: a semiparametrical approach using penalized likelihood estimation where baseline hazard functions are approximated by M‐splines, and another one with piecewise constant baseline hazard functions. Finally, we derived martingale residuals to check the goodness‐of‐fit. We illustrate our proposals with a real dataset on breast cancer. The main objective was to model the dependency between the two types of recurrent events (locoregional and metastatic) and the terminal event (death) after a breast cancer.     

Analysis of interval‐censored recurrent event processes subject to resolution

Interval‐censored recurrent event data arise when the event of interest is not readily observed but the cumulative event count can be recorded at periodic assessment times. In some settings, chronic disease processes may resolve, and individuals will cease to be at risk of events at the time of disease resolution. We develop an expectation‐maximization algorithm for fitting a dynamic mover‐stayer model to interval‐censored recurrent event data under a Markov model with a piecewise‐constant baseline rate function given a latent process. The model is motivated by settings in which the event times and the resolution time of the disease process are unobserved. The likelihood and algorithm are shown to yield estimators with small empirical bias in simulation studies. Data are analyzed on the cumulative number of damaged joints in patients with psoriatic arthritis where individuals experience disease remission.     

Joint regression analysis for survival data in the presence of two sets of semi‐competing risks

In many clinical trials, multiple time‐to‐event endpoints including the primary endpoint (e.g., time to death) and secondary endpoints (e.g., progression‐related endpoints) are commonly used to determine treatment efficacy. These endpoints are often biologically related. This work is motivated by a study of bone marrow transplant (BMT) for leukemia patients, who may experience the acute graft‐versus‐host disease (GVHD), relapse of leukemia, and death after an allogeneic BMT. The acute GVHD is associated with the relapse free survival, and both the acute GVHD and relapse of leukemia are intermediate nonterminal events subject to dependent censoring by the informative terminal event death, but not vice versa, giving rise to survival data that are subject to two sets of semi‐competing risks. It is important to assess the impacts of prognostic factors on these three time‐to‐event endpoints. We propose a novel statistical approach that jointly models such data via a pair of copulas to account for multiple dependence structures, while the marginal distribution of each endpoint is formulated by a Cox proportional hazards model. We develop an estimation procedure based on pseudo‐likelihood and carry out simulation studies to examine the performance of the proposed method in finite samples. The practical utility of the proposed method is further illustrated with data from the motivating example.     

Dynamic clinical prediction models for discrete time‐to‐event data with competing risks—A case study on the OUTCOMEREA database

The development of clinical prediction models requires the selection of suitable predictor variables. Techniques to perform objective Bayesian variable selection in the linear model are well developed and have been extended to the generalized linear model setting as well as to the Cox proportional hazards model. Here, we consider discrete time‐to‐event data with competing risks and propose methodology to develop a clinical prediction model for the daily risk of acquiring a ventilator‐associated pneumonia (VAP) attributed to P. aeruginosa (PA) in intensive care units. The competing events for a PA VAP are extubation, death, and VAP due to other bacteria. Baseline variables are potentially important to predict the outcome at the start of ventilation, but may lose some of their predictive power after a certain time. Therefore, we use a landmark approach for dynamic Bayesian variable selection where the set of relevant predictors depends on the time already spent at risk. We finally determine the direct impact of a variable on each competing event through cause‐specific variable selection.     

A penalized likelihood approach for a progressive three-state model with censored and truncated data: application to AIDS

We consider the estimation of the intensity and survival functions for a continuous time progressive three-state semi-Markov model with intermittently observed data. The estimator of the intensity function is defined nonparametrically as the maximum of a penalized likelihood. We thus obtain smooth estimates of the intensity and survival functions. This approach can accommodate complex observation schemes such as truncation and interval censoring. The method is illustrated with a study of hemophiliacs infected by HIV. The intensity functions and the cumulative distribution functions for the time to infection and for the time to AIDS are estimated. Covariates can easily be incorporated into the model.     

Zero‐modified Poisson model: Bayesian approach,influence diagnostics,and an application to a Brazilian leptospirosis notification data

In this paper, a Bayesian method for inference is developed for the zero‐modified Poisson (ZMP) regression model. This model is very flexible for analyzing count data without requiring any information about inflation or deflation of zeros in the sample. A general class of prior densities based on an information matrix is considered for the model parameters. A sensitivity study to detect influential cases that can change the results is performed based on the Kullback–Leibler divergence. Simulation studies are presented in order to illustrate the performance of the developed methodology. Two real datasets on leptospirosis notification in Bahia State (Brazil) are analyzed using the proposed methodology for the ZMP model.     

Parameter estimation of a plant uptake model for cyanide: application to hydroponic data

A plant uptake model is applied to describe free cyanide and ferrocyanide transport and fate in willow (Salix eriocephala var. Michaux) grown in hydroponics. The model is applied to experimental data to determine best-fit parameter values, their associated uncertainty, and their relative importance to field-scale phytoremediation applications. The fitted model results, using least-squares optimization of the observed log concentrations, indicate that free cyanide volatilization from leaf tissue and free cyanide cell wall adsorption were negligible. The free cyanide maximum uptake rate and assimilate (noncyanide 15N) first-order leaf loss rate were the only coefficients that significantly affected the model goodness of fit and were concurrently sensitive to data uncertainty in the parameter optimization. Saturation kinetics may be applicable for free cyanide uptake into plants, but not for ferrocyanide uptake, which may occur via preferential protein-mediated or inefficient transpiration stream uptake. Within the free cyanide system, the relative magnitudes of the saturation uptake parameters and the demonstration of an active role for plants in uptake relative to transpiration suggest the potential importance of preferential diffusion through the cell membranes as reported in the literature, rather than protein-mediated uptake. The fitted 13-parameter model matched the observed data well except for the predicted stem and leaf tissue assimilate concentrations, which were significantly underestimated, particularly in the free cyanide system. These low predicted values, combined with the slightly underestimated solution free cyanide removal, suggest that noncyanide 15N redistribution in phloem should be considered.     

A penalized likelihood approach for an illness-death model with interval-censored data: application to age-specific incidence of dementia

We consider the problem of estimating the intensity functions for a continuous time 'illness-death' model with intermittently observed data. In such a case, it may happen that a subject becomes diseased between two visits and dies without being observed. Consequently, there is an uncertainty about the precise number of transitions. Estimating the intensity of transition from health to illness by survival analysis (treating death as censoring) is biased downwards. Furthermore, the dates of transitions between states are not known exactly. We propose to estimate the intensity functions by maximizing a penalized likelihood. The method yields smooth estimates without parametric assumptions. This is illustrated using data from a large cohort study on cerebral ageing. The age-specific incidence of dementia is estimated using an illness-death approach and a survival approach.     

Analysis of interval‐grouped data in weed science: The binnednp Rcpp package

相似文献   

Examining methodology to identify patterns of consulting in primary care for different groups of patients before a diagnosis of cancer: An exemplar applied to oesophagogastric cancer

BackgroundCurrent methods for estimating the timeliness of cancer diagnosis are not robust because dates of key defining milestones, for example first presentation, are uncertain. This is exacerbated when patients have other conditions (multimorbidity), particularly those that share symptoms with cancer. Methods independent of this uncertainty are needed for accurate estimates of the timeliness of cancer diagnosis, and to understand how multimorbidity impacts the diagnostic process.MethodsParticipants were diagnosed with oesophagogastric cancer between 2010 and 2019. Controls were matched on year of birth, sex, general practice and multimorbidity burden calculated using the Cambridge Multimorbidity Score. Primary care data (Clinical Practice Research Datalink) was used to explore population-level consultation rates for up to two years before diagnosis across different multimorbidity burdens. Five approaches were compared on the timing of the consultation frequency increase, the inflection point for different multimorbidity burdens, different aggregated time-periods and sample sizes.ResultsWe included 15,410 participants, of which 13,328 (86.5 %) had a measurable multimorbidity burden. Our new maximum likelihood estimation method found evidence that the inflection point in consultation frequency varied with multimorbidity burden, from 154 days (95 %CI 131.8–176.2) before diagnosis for patients with no multimorbidity, to 126 days (108.5–143.5) for patients with the greatest multimorbidity burden. Inflection points identified using alternative methods were closer to diagnosis for up to three burden groups. Sample size reduction and changing the aggregation period resulted in inflection points closer to diagnosis, with the smallest change for the maximum likelihood method.DiscussionExisting methods to identify changes in consultation rates can introduce substantial bias which depends on sample size and aggregation period. The direct maximum likelihood method was less prone to this bias than other methods and offers a robust, population-level alternative for estimating the timeliness of cancer diagnosis.