Generalized gamma frailty model

The quasi-xgamma frailty model with survival analysis under heterogeneity problem, validation testing, and risk analysis for emergency care data

Frailty models are important for survival data because they allow for the possibility of unobserved heterogeneity problem. The problem of heterogeneity can be existed due to a variety of factors, such as genetic predisposition, environmental factors, or lifestyle choices. Frailty models can help to identify these factors and to better understand their impact on survival. In this study, we suggest a novel quasi xgamma frailty (QXg-F) model for the survival analysis. In this work, the test of Rao-Robson and Nikulin is employed to test the validity and suitability of the probabilistic model, we examine the distribution's properties and evaluate its performance in comparison with many relevant cox-frailty models. To show how well the QXg-F model captures heterogeneity and enhances model fit, we use simulation studies and real data applications, including a fresh dataset gathered from an emergency hospital in Algeria. According to our research, the QXg-F model is a viable replacement for the current frailty modeling distributions and has the potential to improve the precision of survival analyses in a number of different sectors, including emergency care. Moreover, testing the ability and the importance of the new QXg-F model in insurance is investigated using simulations via different methods and application to insurance data.

A two-way flexible generalized gamma transformation cure rate model

We propose a two-way flexible cure rate model. The first flexibility is provided by considering a family of Box-Cox transformation cure models that include the commonly used cure models as special cases. The second flexibility is provided by proposing the wider class of generalized gamma distributions to model the associated lifetime. The advantage of this two-way flexibility is that it allows us to carry out tests of hypotheses to select an adequate cure model (within the family of Box-Cox transformation cure models) and a suitable lifetime distribution (within the wider class of generalized gamma distributions) that jointly provides the best fit to a given data. First, we study the maximum likelihood estimation of the generalized gamma Box-Cox transformation (GGBCT) model parameters. Then, we use the flexibility of our proposed model to carry out power studies to demonstrate the power of likelihood ratio test in rejecting mis-specified models. Furthermore, we study the bias and efficiency of the estimators of the cure rates under model mis-specification. Our findings strongly suggest the importance of selecting a correct lifetime distribution and a correct cure rate model, which can be achieved through the proposed two-way flexible model. Finally, we illustrate the applicability of our proposed model using a data from a breast cancer study and show that our model provides a better fit than the existing semiparametric Box-Cox transformation cure model with piecewise exponential approximation to the lifetime distribution.

Left-censored recurrent event analysis in epidemiological studies: a proposal for when the number of previous episodes is unknown

BACKGROUND

When dealing with recurrent events in observational studies it is common to include subjects who became at risk before follow-up. This phenomenon is known as left censoring, and simply ignoring these prior episodes can lead to biased and inefficient estimates. We aimed to propose a statistical method that performs well in this setting.

METHODS

Our proposal was based on the use of models with specific baseline hazards. In this, the number of prior episodes were imputed when unknown and stratified according to whether the subject had been at risk of presenting the event before t = 0. A frailty term was also used. Two formulations were used for this "Specific Hazard Frailty Model Imputed" based on the "counting process" and "gap time." Performance was then examined in different scenarios through a comprehensive simulation study.

RESULTS

The proposed method performed well even when the percentage of subjects at risk before follow-up was very high. Biases were often below 10% and coverages were around 95%, being somewhat conservative. The gap time approach performed better with constant baseline hazards, whereas the counting process performed better with non-constant baseline hazards.

CONCLUSIONS

The use of common baseline methods is not advised when knowledge of prior episodes experienced by a participant is lacking. The approach in this study performed acceptably in most scenarios in which it was evaluated and should be considered an alternative in this context. It has been made freely available to interested researchers as R package miRecSurv.

Penalized variable selection for cause-specific hazard frailty models with clustered competing-risks data

Competing risks data usually arise when an occurrence of an event precludes other types of events from being observed. Such data are often encountered in a clustered clinical study such as a multi-center clinical trial. For the clustered competing-risks data which are correlated within a cluster, competing-risks models allowing for frailty terms have been recently studied. To the best of our knowledge, however, there is no literature on variable selection methods for cause-specific hazard frailty models. In this article, we propose a variable selection procedure for fixed effects in cause-specific competing risks frailty models using a penalized h-likelihood (HL). Here, we study three penalty functions, LASSO, SCAD, and HL. Simulation studies demonstrate that the proposed procedure using the HL penalty works well, providing a higher probability of choosing the true model than LASSO and SCAD methods without losing prediction accuracy. The proposed method is illustrated by using two kinds of clustered competing-risks cancer data sets.

Weighted Lindley frailty model: estimation and application to lung cancer data

In this paper, we propose a novel frailty model for modeling unobserved heterogeneity present in survival data. Our model is derived by using a weighted Lindley distribution as the frailty distribution. The respective frailty distribution has a simple Laplace transform function which is useful to obtain marginal survival and hazard functions. We assume hazard functions of the Weibull and Gompertz distributions as the baseline hazard functions. A classical inference procedure based on the maximum likelihood method is presented. Extensive simulation studies are further performed to verify the behavior of maximum likelihood estimators under different proportions of right-censoring and to assess the performance of the likelihood ratio test to detect unobserved heterogeneity in different sample sizes. Finally, to demonstrate the applicability of the proposed model, we use it to analyze a medical dataset from a population-based study of incident cases of lung cancer diagnosed in the state of São Paulo, Brazil.

Nonproportional hazards model with a frailty term for modeling subgroups with evidence of long-term survivors: Application to a lung cancer dataset

With advancements in medical treatments for cancer, an increase in the life expectancy of patients undergoing new treatments is expected. Consequently, the field of statistics has evolved to present increasingly flexible models to explain such results better. In this paper, we present a lung cancer dataset with some covariates that exhibit nonproportional hazards (NPHs). Besides, the presence of long-term survivors is observed in subgroups. The proposed modeling is based on the generalized time-dependent logistic model with each subgroup's effect time and a random term effect (frailty). In practice, essential covariates are not observed for several reasons. In this context, frailty models are useful in modeling to quantify the amount of unobservable heterogeneity. The frailty distribution adopted was the weighted Lindley distribution, which has several interesting properties, such as the Laplace transform function on closed form, flexibility in the probability density function, among others. The proposed model allows for NPHs and long-term survivors in subgroups. Parameter estimation was performed using the maximum likelihood method, and Monte Carlo simulation studies were conducted to evaluate the estimators' performance. We exemplify this model's use by applying data of patients diagnosed with lung cancer in the state of São Paulo, Brazil.

A systematic review of the quality of conduct and reporting of survival analyses of tuberculosis outcomes in Africa

BACKGROUND

Survival analyses methods (SAMs) are central to analysing time-to-event outcomes. Appropriate application and reporting of such methods are important to ensure correct interpretation of the data. In this study, we systematically review the application and reporting of SAMs in studies of tuberculosis (TB) patients in Africa. It is the first review to assess the application and reporting of SAMs in this context.

METHODS

Systematic review of studies involving TB patients from Africa published between January 2010 and April 2020 in English language. Studies were eligible if they reported use of SAMs. Application and reporting of SAMs were evaluated based on seven author-defined criteria.

RESULTS

Seventy-six studies were included with patient numbers ranging from 56 to 182,890. Forty-three (57%) studies involved a statistician/epidemiologist. The number of published papers per year applying SAMs increased from two in 2010 to 18 in 2019 (P = 0.004). Sample size estimation was not reported by 67 (88%) studies. A total of 22 (29%) studies did not report summary follow-up time. The survival function was commonly presented using Kaplan-Meier survival curves (n = 51, (67%) studies) and group comparisons were performed using log-rank tests (n = 44, (58%) studies). Sixty seven (91%), 3 (4.1%) and 4 (5.4%) studies reported Cox proportional hazard, competing risk and parametric survival regression models, respectively. A total of 37 (49%) studies had hierarchical clustering, of which 28 (76%) did not adjust for the clustering in the analysis. Reporting was adequate among 4.0, 1.3 and 6.6% studies for sample size estimation, plotting of survival curves and test of survival regression underlying assumptions, respectively. Forty-five (59%), 52 (68%) and 73 (96%) studies adequately reported comparison of survival curves, follow-up time and measures of effect, respectively.

CONCLUSION

The quality of reporting survival analyses remains inadequate despite its increasing application. Because similar reporting deficiencies may be common in other diseases in low- and middle-income countries, reporting guidelines, additional training, and more capacity building are needed along with more vigilance by reviewers and journal editors.

A general frailty model to accommodate individual heterogeneity in the acquisition of multiple infections: An application to bivariate current status data

The analysis of multivariate time-to-event (TTE) data can become complicated due to the presence of clustering, leading to dependence between multiple event times. For a long time, (conditional) frailty models and (marginal) copula models have been used to analyze clustered TTE data. In this article, we propose a general frailty model employing a copula function between the frailty terms to construct flexible (bivariate) frailty distributions with the application to current status data. The model has the advantage to impose a less restrictive correlation structure among latent frailty variables as compared to traditional frailty models. Specifically, our model uses a copula function to join the marginal distributions of the frailty vector. In this article, we considered different copula functions, and we relied on marginal gamma distributions due to their mathematical convenience. Based on a simulation study, our novel model outperformed the commonly used additive correlated gamma frailty model, especially in the case of a negative association between the frailties. At the end of the article, the new methodology is illustrated on real-life data applications entailing bivariate serological survey data.

Semiparametric Regression Estimation for Recurrent Event Data with Errors in Covariates under Informative Censoring

Recurrent event data arise frequently in many longitudinal follow-up studies. Hence, evaluating covariate effects on the rates of occurrence of such events is commonly of interest. Examples include repeated hospitalizations, recurrent infections of HIV, and tumor recurrences. In this article, we consider semiparametric regression methods for the occurrence rate function of recurrent events when the covariates may be measured with errors. In contrast to the existing works, in our case the conventional assumption of independent censoring is violated since the recurrent event process is interrupted by some correlated events, which is called informative drop-out. Further, some covariates may be measured with errors. To accommodate for both informative censoring and measurement error, the occurrence of recurrent events is modelled through an unspecified frailty distribution and accompanied with a classical measurement error model. We propose two corrected approaches based on different ideas, and we show that they are numerically identical when estimating the regression parameters. The asymptotic properties of the proposed estimators are established, and the finite sample performance is examined via simulations. The proposed methods are applied to the Nutritional Prevention of Cancer trial for assessing the effect of the plasma selenium treatment on the recurrence of squamous cell carcinoma.

Likelihood inference for COM-Poisson cure rate model with interval-censored data and Weibull lifetimes

In this paper, we consider a competing cause scenario and assume the number of competing causes to follow a Conway–Maxwell Poisson distribution which can capture both over and under dispersion that is usually encountered in discrete data. Assuming the population of interest having a component cure and the form of the data to be interval censored, as opposed to the usually considered right-censored data, the main contribution is in developing the steps of the expectation maximization algorithm for the determination of the maximum likelihood estimates of the model parameters of the flexible Conway–Maxwell Poisson cure rate model with Weibull lifetimes. An extensive Monte Carlo simulation study is carried out to demonstrate the performance of the proposed estimation method. Model discrimination within the Conway–Maxwell Poisson distribution is addressed using the likelihood ratio test and information-based criteria to select a suitable competing cause distribution that provides the best fit to the data. A simulation study is also carried out to demonstrate the loss in efficiency when selecting an improper competing cause distribution which justifies the use of a flexible family of distributions for the number of competing causes. Finally, the proposed methodology and the flexibility of the Conway–Maxwell Poisson distribution are illustrated with two known data sets from the literature: smoking cessation data and breast cosmesis data.

Bayesian cure rate models induced by frailty in survival analysis

Frailty models provide a convenient way of modeling unobserved dependence and heterogeneity in survival data which, if not accounted for duly, would result incorrect inference. Gamma frailty models are commonly used for this purpose, but alternative continuous distributions are possible as well. However, with cure rate being present in survival data, these continuous distributions may not be appropriate since individuals with long-term survival times encompass zero frailty. So, we propose here a flexible probability distribution induced by a discrete frailty, and then present some special discrete probability distributions. We specifically focus on a special hyper-Poisson distribution and then develop the corresponding Bayesian simulation, influence diagnostics and an application to real dataset by means of intensive Markov chain Monte Carlo algorithm. These illustrate the usefulness of the proposed model as well as the inferential results developed here.

Power and sample size calculation for paired recurrent events data based on robust nonparametric tests

The purpose of this paper is to develop a formula for calculating the required sample size for paired recurrent events data. The developed formula is based on robust non-parametric tests for comparing the marginal mean function of events between paired samples. This calculation can accommodate the associations among a sequence of paired recurrent event times with a specification of correlated gamma frailty variables for a proportional intensity model. We evaluate the performance of the proposed method with comprehensive simulations including the impacts of paired correlations, homogeneous or nonhomogeneous processes, marginal hazard rates, censoring rate, accrual and follow-up times, as well as the sensitivity analysis for the assumption of the frailty distribution. The use of the formula is also demonstrated using a premature infant study from the neonatal intensive care unit of a tertiary center in southern Taiwan. Copyright © 2017 John Wiley & Sons, Ltd.

A frailty model for recurrent events during alternating restraint and non-restraint time periods

We consider recurrent events of the same type that occur during alternating restraint and non-restraint time periods. This research is motivated by a study on juvenile recidivism, where the probationers were followed for re-offenses during alternating placement periods and free-time periods. During the placement periods, the probationers were under a restricted environment with direct supervision of the probation officers. During the free-time periods, the probationers were released to home and not under direct supervision. Although re-offenses can occur during both types of time periods, the intensities of the re-offenses are very different. Thus, these two types of time periods should be modeled differently. The same data structure also arises in many biomedical settings, as exemplified by tumor metastases during chemotherapy and chemo-free periods. In this paper, we propose a joint modeling framework that explicitly accounts for the different types of time periods, as well as the within-subject dependence during the same type and between different types of time periods. The estimation procedure is implemented in SAS and is easily accessible to practical investigators. We evaluate the proposed method through simulation studies under several realistic scenarios and demonstrate the feasibility of the proposed method by applying it to the juvenile recidivism dataset. Copyright © 2016 John Wiley & Sons, Ltd.

Analyzing recurrent events when the history of previous episodes is unknown or not taken into account: proceed with caution

OBJECTIVE

Researchers in public health are often interested in examining the effect of several exposures on the incidence of a recurrent event. The aim of the present study is to assess how well the common-baseline hazard models perform to estimate the effect of multiple exposures on the hazard of presenting an episode of a recurrent event, in presence of event dependence and when the history of prior-episodes is unknown or is not taken into account.

METHODS

Through a comprehensive simulation study, using specific-baseline hazard models as the reference, we evaluate the performance of common-baseline hazard models by means of several criteria: bias, mean squared error, coverage, confidence intervals mean length and compliance with the assumption of proportional hazards.

RESULTS

Results indicate that the bias worsen as event dependence increases, leading to a considerable overestimation of the exposure effect; coverage levels and compliance with the proportional hazards assumption are low or extremely low, worsening with increasing event dependence, effects to be estimated, and sample sizes.

CONCLUSIONS

Common-baseline hazard models cannot be recommended when we analyse recurrent events in the presence of event dependence. It is important to have access to the history of prior-episodes per subject, it can permit to obtain better estimations of the effects of the exposures.

Expectation maximization-based likelihood inference for flexible cure rate models with Weibull lifetimes

Recently, a flexible cure rate survival model has been developed by assuming the number of competing causes of the event of interest to follow the Conway–Maxwell–Poisson distribution. This model includes some of the well-known cure rate models discussed in the literature as special cases. Data obtained from cancer clinical trials are often right censored and expectation maximization algorithm can be used in this case to efficiently estimate the model parameters based on right censored data. In this paper, we consider the competing cause scenario and assuming the time-to-event to follow the Weibull distribution, we derive the necessary steps of the expectation maximization algorithm for estimating the parameters of different cure rate survival models. The standard errors of the maximum likelihood estimates are obtained by inverting the observed information matrix. The method of inference developed here is examined by means of an extensive Monte Carlo simulation study. Finally, we illustrate the proposed methodology with a real data on cancer recurrence.

Safety Profile of Nifurtimox and Treatment Interruption for Chronic Chagas Disease in Colombian Adults

Nifurtimox (NFX) is one of the approved drugs used to treat Chagas disease. Safety profile studies and models on risk factors for treatment interruption in adults are scarce in Latin America. This study evaluated retrospectively the medical records of adult Chagas disease patients treated with NFX between 2007 and 2012 in Bogotá, Colombia. An accelerated failure time model was used, and associations were expressed as time ratio (TR). In total, 76 adult patients with NFX were included: 60 (79.0%) completed 60 days of treatment, 61 (80.3%) presented adverse drug reactions (ADRs), and 16 (21.0%) required treatment interruption. The predominant symptoms were epigastric pain (23.7%), nauseas (18.4%), sleep disturbances (18.4%), loss of appetite (17.1%), and temporary loss of memory (15.2%). ADRs were classified as mild (64.5%), moderate (30.4%), and severe (5.1%). Time of treatment was significantly longer when presenting ≤ 3 ADRs (TR: 1.78; 95% CI: 1.04–3.03), presence of non-severe ADRs (TR: 6.52; 95% CI: 3.24–13.1), doses of NFX ≤ 8 mg/kg/day (TR: 1.78; 95% CI: 0.90–3.49), and age < 48 years (TR: 1.57; 95% CI: 0.90–2.74). Treatment with NFX in adults caused a high frequency of ADRs, but most of the cases were mild and did not require treatment interruption. Severity and number of ADRs were the main predictors for treatment interruption.

The Cox shared frailty model with log-skew-normal frailties

Clustered survival data arise when groups of failure times share a common ingredient; typically, they refer to the same individual or individuals with a common factor. When the association between failure times within the same cluster is of interest, statistical methods called frailty models have been used. The frailty is an unobserved random component which affects the risk level, changing from cluster to cluster but shared by all observations within the same cluster. Various probability distributions have been proposed for the frailty term, with special emphasis on the gamma and log-normal distribution. Since adequate modelling of the frailty distribution is essential to properly investigate the dependence structure, we introduce a new class of frailty models with a flexible distribution form. Specifically, we adopt the skew-normal distribution for the log-transformed frailty, leading to an extension of the log-normal model. After presenting the methodology connected to this choice, we illustrate it with a case study of multiple myeloma patients with autologous stem cells transplantation.