On a new piecewise regression model with cure rate: Diagnostics and application to medical data

In this article, we discuss an extension of the classical negative binomial cure rate model with piecewise exponential distribution of the time to event for concurrent causes, which enables the modeling of monotonic and non-monotonic hazard functions (ie, the shape of the hazard function is not assumed as in traditional parametric models). This approach produces a flexible cure rate model, depending on the choice of time partition. We discuss local influence on this negative binomial power piecewise exponential model. We report on Monte Carlo simulation studies and application of the model to real melanoma and leukemia datasets.

Dynamic Bayesian adjustment of anticipatory covariates in retrospective data: application to the effect of education on divorce risk

We address a problem in inference from retrospective studies where the value of a variable is measured at the date of the survey but is used as covariate to events that have occurred long before the survey. This causes problem because the value of the current-date (anticipatory) covariate does not follow the temporal order of events. We propose a dynamic Bayesian approach for modelling jointly the anticipatory covariate and the event of interest, and allowing the effects of the anticipatory covariate to vary over time. The issues are illustrated with data on the effects of education attained by the survey-time on divorce risks among Swedish men. The overall results show that failure to adjust for the anticipatory nature of education leads to elevated relative risks of divorce across educational levels. The results are partially in accordance with previous findings based on analyses of the same data set. More importantly, our findings provide new insights in that the bias due to anticipatory covariates varies over marriage duration.

Flexible Bayesian survival modeling with semiparametric time-dependent and shape-restricted covariate effects

Presently, there are few options with available software to perform a fully Bayesian analysis of time-to-event data wherein the hazard is estimated semi- or non-parametrically. One option is the piecewise exponential model, which requires an often unrealistic assumption that the hazard is piecewise constant over time. The primary aim of this paper is to construct a tractable semiparametric alternative to the piecewise exponential model that assumes the hazard is continuous, and to provide modifiable, user-friendly software that allows the use of these methods in a variety of settings. To accomplish this aim, we use a novel model formulation for the log-hazard based on a low-rank thin plate linear spline that readily facilitates adjustment for covariates with time-dependent and proportional hazards effects, possibly subject to shape restrictions. We investigate the performance of our model choices via simulation. We then analyze colorectal cancer data from a clinical trial comparing the effectiveness of two novel treatment regimes relative to the standard of care for overall survival. We estimate a time-dependent hazard ratio for each novel regime relative to the standard of care while adjusting for the effect of aspartate transaminase, a biomarker of liver function, that is subject to a non-decreasing shape restriction.

Application of the Bayesian dynamic survival model in medicine

The Bayesian dynamic survival model (BDSM), a time-varying coefficient survival model from the Bayesian prospective, was proposed in early 1990s but has not been widely used or discussed. In this paper, we describe the model structure of the BDSM and introduce two estimation approaches for BDSMs: the Markov Chain Monte Carlo (MCMC) approach and the linear Bayesian (LB) method. The MCMC approach estimates model parameters through sampling and is computationally intensive. With the newly developed geoadditive survival models and software BayesX, the BDSM is available for general applications. The LB approach is easier in terms of computations but it requires the prespecification of some unknown smoothing parameters. In a simulation study, we use the LB approach to show the effects of smoothing parameters on the performance of the BDSM and propose an ad hoc method for identifying appropriate values for those parameters. We also demonstrate the performance of the MCMC approach compared with the LB approach and a penalized partial likelihood method available in software R packages. A gastric cancer trial is utilized to illustrate the application of the BDSM.

Cox Models With Smooth Functional Effect of Covariates Measured With Error

We propose, develop, and implement a fully Bayesian inferential approach for the Cox model when the log hazard function contains unknown smooth functions of the variables measured with error. Our approach is to model nonparametrically both the log-baseline hazard and the smooth components of the log-hazard functions using low-rank penalized splines. Careful implementation of the Bayesian inferential machinery is shown to produce remarkably better results than the naive approach. Our methodology was motivated by and applied to the study of progression time to chronic kidney disease as a function of baseline kidney function and applied to the Atherosclerosis Risk in Communities study, a large epidemiological cohort study. This article has supplementary material online.

Time-varying effects of prognostic factors associated with disease-free survival in breast cancer

Early detection and effective treatments have dramatically improved breast cancer survivorship, yet the risk of relapse persists even 15 years after the initial diagnosis. It is important to identify prognostic factors for late breast cancer events. The authors investigated time-varying effects of tumor characteristics on breast-cancer-free survival using data on 3,088 breast cancer survivors from 4 US states who participated in a randomized dietary intervention trial in 1995-2006, with maximum follow-up through 15 years (median, 9 years). A piecewise constant penalized spline approach incorporating time-varying coefficients was adopted, allowing for deviations from the proportional hazards assumption. This method is more flexible than standard approaches, provides direct estimates of hazard ratios across time intervals, and is computationally tractable. Having a stage II or III tumor was associated with a 3-fold higher hazard of breast cancer than having a stage I tumor during the first 2.5 years after diagnosis; this hazard ratio decreased to 2.1 after 7.7 years, but higher tumor stage remained a significant risk factor. Similar diminishing effects were found for poorly differentiated tumors. Interestingly, having a positive estrogen receptor status was protective up to 4 years after diagnosis but detrimental after 7.7 years (hazard ratio = 1.5). These results emphasize the importance of careful statistical modeling allowing for possibly time-dependent effects in long-term survivorship studies.

Regression splines in the time-dependent coefficient rates model for recurrent event data

Many epidemiologic studies involve the occurrence of recurrent events and much attention has been given for the development of modeling techniques that take into account the dependence structure of multiple event data. This paper presents a time-dependent coefficient rates model that incorporates regression splines in its estimation procedure. Such methods would be appropriate in situations where the effect of an exposure or covariates changes over time in recurrent event data settings. The finite sample properties of the estimators are studied via simulation. Using data from a randomized community trial that was designed to evaluate the effect of vitamin A supplementation on recurrent diarrheal episodes in small children, we model the functional form of the treatment effect on the time to the occurrence of diarrhea. The results describe how this effect varies over time. In summary, we observed a major impact of the vitamin A supplementation on diarrhea after 2 months of the dosage, with the effect diminishing after the third dosage. The proposed method can be viewed as a flexible alternative to the marginal rates model with constant effect in situations where the effect of interest may vary over time.

A semiparametric Bayesian proportional hazards model for interval censored data with frailty effects

BACKGROUND

Multivariate analysis of interval censored event data based on classical likelihood methods is notoriously cumbersome. Likelihood inference for models which additionally include random effects are not available at all. Developed algorithms bear problems for practical users like: matrix inversion, slow convergence, no assessment of statistical uncertainty.

METHODS

MCMC procedures combined with imputation are used to implement hierarchical models for interval censored data within a Bayesian framework.

RESULTS

Two examples from clinical practice demonstrate the handling of clustered interval censored event times as well as multilayer random effects for inter-institutional quality assessment. The software developed is called survBayes and is freely available at CRAN.

CONCLUSION

The proposed software supports the solution of complex analyses in many fields of clinical epidemiology as well as health services research.

Estimating the grid of time-points for the piecewise exponential model

One of the greatest challenges related to the use of piecewise exponential models (PEMs) is to find an adequate grid of time-points needed in its construction. In general, the number of intervals in such a grid and the position of their endpoints are ad-hoc choices. We extend previous works by introducing a full Bayesian approach for the piecewise exponential model in which the grid of time-points (and, consequently, the endpoints and the number of intervals) is random. We estimate the failure rates using the proposed procedure and compare the results with the non-parametric piecewise exponential estimates. Estimates for the survival function using the most probable partition are compared with the Kaplan-Meier estimators (KMEs). A sensitivity analysis for the proposed model is provided considering different prior specifications for the failure rates and for the grid. We also evaluate the effect of different percentage of censoring observations in the estimates. An application to a real data set is also provided. We notice that the posteriors are strongly influenced by prior specifications, mainly for the failure rates parameters. Thus, the priors must be fairly built, say, really disclosing the expert prior opinion.

The separation of timescales in Bayesian survival modeling of the time-varying effect of a time-dependent exposure

In this paper, we apply flexible Bayesian survival analysis methods to investigate the risk of lymphoma associated with kidney transplantation among patients with end-stage renal disease. Of key interest is the potentially time-varying effect of a time-dependent exposure: transplant status. Bayesian modeling of the baseline hazard and the effect of transplant requires consideration of 2 timescales: time since study start and time since transplantation, respectively. Previous related work has not dealt with the separation of multiple timescales. Using a hierarchical model for the hazard function, both timescales are incorporated via conditionally independent stochastic processes; smoothing of each process is specified via intrinsic conditional Gaussian autoregressions. Features of the corresponding posterior distribution are evaluated from draws obtained via a Metropolis-Hastings-Green algorithm.

Approaches in modelling long-term survival: an application to breast cancer

Several modelling techniques have been proposed for non-proportional hazards. In this work we consider different models which can be classified into three wide categories: models with time-varying effects of the covariates; frailty models and cure rate models. We present those different extensions of the proportional hazards model on an application of 2433 breast cancer patients with a long follow-up. We comment on the differences and similarities among the models and evaluate their performance using survival and hazard plots, Brier scores and pseudo-observations.

Bayesian semiparametric dynamic frailty models for multiple event time data

Many biomedical studies collect data on times of occurrence for a health event that can occur repeatedly, such as infection, hospitalization, recurrence of disease, or tumor onset. To analyze such data, it is necessary to account for within-subject dependency in the multiple event times. Motivated by data from studies of palpable tumors, this article proposes a dynamic frailty model and Bayesian semiparametric approach to inference. The widely used shared frailty proportional hazards model is generalized to allow subject-specific frailties to change dynamically with age while also accommodating nonproportional hazards. Parametric assumptions on the frailty distribution are avoided by using Dirichlet process priors for a shared frailty and for multiplicative innovations on this frailty. By centering the semiparametric model on a conditionally conjugate dynamic gamma model, we facilitate posterior computation and lack-of-fit assessments of the parametric model. Our proposed method is demonstrated using data from a cancer chemoprevention study.

Bayesian dynamic models for survival data with a cure fraction

In this paper, we propose a new class of semi-parametric cure rate models. Specifically, we construct dynamic models for piecewise hazard functions over a finite partition of the time axis. Allowing the size of partition and the levels of baseline hazard to be random, our proposed models provide a great flexibility in controlling the degree of parametricity in the right tail of the survival distribution and the amount of correlations among the log-baseline hazard levels. Several properties of the proposed models are derived, and propriety of the implied posteriors with improper noninformative priors for regression coefficients based on the proposed models is established for the fixed partition of the time axis. In addition, an efficient reversible jump computational algorithm is developed for carrying out posterior computation. A real data set from a melanoma clinical trial is analyzed in detail to further demonstrate the proposed methodology.

Dynamic survival models with spatial frailty

In many survival studies, covariates effects are time-varying and there is presence of spatial effects. Dynamic models can be used to cope with the variations of the effects and spatial components are introduced to handle spatial variation. This paper proposes a methodology to simultaneously introduce these components into the model. A number of specifications for the spatial components are considered. Estimation is performed via a Bayesian approach through Markov chain Monte Carlo methods. Models are compared to assess relevance of their components. Analysis of a real data set is performed, showing the relevance of both time-varying covariate effects and spatial components. Extensions to the methodology are proposed along with concluding remarks.

Reduced-rank hazard regression for modelling non-proportional hazards

The Cox proportional hazards model is the most common method to analyse survival data. However, the proportional hazards assumption might not hold. The natural extension of the Cox model is to introduce time-varying effects of the covariates. For some covariates such as (surgical)treatment non-proportionality could be expected beforehand. For some other covariates the non-proportionality only becomes apparent if the follow-up is long enough. It is often observed that all covariates show similar decaying effects over time. Such behaviour could be explained by the popular (gamma-) frailty model. However, the (marginal) effects of covariates in frailty models are not easy to interpret. In this paper we propose the reduced-rank model for time-varying effects of covariates. Starting point is a Cox model with p covariates and time-varying effects modelled by q time functions (constant included), leading to a pxq structure matrix that contains the regression coefficients for all covariate by time function interactions. By reducing the rank of this structure matrix a whole range of models is introduced, from the very flexible full-rank model (identical to a Cox model with time-varying effects) to the very rigid rank one model that mimics the structure of a gamma-frailty model, but is easier to interpret. We illustrate these models with an application to ovarian cancer patients.

A class of parametric dynamic survival models

A class of parametric dynamic survival models are explored in which only limited parametric assumptions are made, whilst avoiding the assumption of proportional hazards. Both the log-baseline hazard and covariate effects are modelled by piecewise constant and correlated processes. The method of estimation is to use Markov chain Monte Carlo simulations: Gibbs sampling with a Metropolis-Hastings step. In addition to standard right censored data sets, extensions to accommodate interval censoring and random effects are included. The model is applied to two well known and illustrative data sets, and the dynamic variability of covariate effects investigated.

Meta-analysis of polytetrafluoroethylene bypass grafts to infrapopliteal arteries

CONTEXT

Reports of polytetrafluoroethylene (PTFE) bypass grafting to the infrapopliteal arteries have often used survival analysis of acceptable quality to describe a wide range of long-term results. In theory, these results may be combined if variability between series and time intervals is considered.

OBJECTIVE

Meta-analysis was performed to gain insight into long-term graft patency and foot preservation after PTFE bypass grafting to infrapopliteal arteries.

DATA SOURCE

Studies published from 1982 through 2001 were identified from the PubMed database and pertinent original articles.

STUDY SELECTION

Three investigators selected 43 studies that used survival analysis, reported 2-year patency rates, and included at least 15 bypass procedures. Data extraction and transformation: Based on standard life-tables or survivor curves, an interval success rate was calculated for each month in each series. The monthly success rates were combined across series, enabling construction of pooled survivor curves.

DATA SYNTHESIS

Random-effects meta-analysis yielded 5-year pooled estimates (SE) of 30.5% (7.6%) for primary graft patency, 39.7% (5.5%) for secondary graft patency, and 55.7% (5.0%) for foot preservation. During the entire follow-up, pooled estimates were slightly higher for series of PTFE grafts with adjunctive procedures compared with series of PTFE grafts only. Sensitivity analysis: A simulation using only unfavorable assumptions showed a decrease of less than 5% at 5 years for all outcomes, and smaller differences at subgroup meta-analysis. Funnel plots suggested that publication bias was unlikely.

CONCLUSION

This meta-analysis indicated moderate success for PTFE bypass grafts to infrapopliteal arteries, but the role of adjunctive procedures at the distal anastomosis remains uncertain.

Bayesian modeling of time-varying and waning exposure effects

In epidemiologic studies, there is often interest in assessing the association between exposure history and disease incidence. For many diseases, incidence may depend not only on cumulative exposure, but also on the ages at which exposure occurred. This article proposes a flexible Bayesian approach for modeling age-varying and waning exposure effects. The Cox model is generalized to allow the hazard of disease to depend on an integral, across the exposed ages, of a piecewise polynomial function of age, multiplied by an exponential decay term. Linearity properties of the model facilitate posterior computation via a Gibbs sampler, which generalizes previous algorithms for Cox regression with time-dependent covariates. The approach is illustrated by an application to the study of protective effects of breastfeeding on incidence of childhood asthma.

Bayesian methods for a three-state model for rodent carcinogenicity studies

The objective of a chronic rodent bioassay is to assess the impact of a chemical compound on the development of tumors. However, most tumor types are not observable prior to necropsy, making direct estimation of the tumor incidence rate problematic. In such cases, estimation can proceed only if the study incorporates multiple interim sacrifices or we make use of simplified parametric or nonparametric models. In addition, it is widely accepted that other factors, such as weight, can be related to both dose level and tumor onset, confounding the association of interest. However, there is not typically enough information in the current study to assess such effects. The addition of historical data can help alleviate this problem. In this article, we propose a novel Bayesian semiparametric model for the analysis of data from rodent carcinogenicity studies. We develop informative prior distributions for covariate effects through the use of historical control data and outline a Gibbs sampling scheme. We implement the model by analyzing data from a National Toxicology Program chronic rodent bioassay.

Bayesian modeling of incidence and progression of disease from cross-sectional data

In the absence of longitudinal data, the current presence and severity of disease can be measured for a sample of individuals to investigate factors related to disease incidence and progression. In this article, Bayesian discrete-time stochastic models are developed for inference from cross-sectional data consisting of the age at first diagnosis, the current presence of disease, and one or more surrogates of disease severity. Semiparametric models are used for the age-specific hazards of onset and diagnosis, and a normal underlying variable approach is proposed for modeling of changes with latency time in disease severity. The model accommodates multiple surrogates of disease severity having different measurement scales and heterogeneity among individuals in disease progression. A Markov chain Monte Carlo algorithm is described for posterior computation, and the methods are applied to data from a study of uterine leiomyoma.

Time-varying exposure and the impact of stressful life events on onset of affective disorder

Stressful life events are now established as risk factors for the onset of affective disorder but few studies have investigated time-varying exposure effects. Discrete (grouped) time survival methods provide a flexible framework for evaluating multiple time-dependent covariates and time-varying covariate effects. Here, we use these methods to investigate the time-varying influence of life events on the onset of affective disorder. Various straightforward time-varying exposure models are compared, involving one or more (stepped) time-dependent covariates and time-dependent covariates constructed or estimated according to exponential decay. These models are applied to data from two quite different studies. The first, a small scale interviewer-based longitudinal study (n = 180) concerned with affective disorder onset following loss (or threat of loss) event experiences. The second, a questionnaire assessment as part of an ongoing population study (n = 3353), provides a history of marital loss events and of depressive disorder onset. From the first study the initial impact of loss events was found to decay with a half-life of 5 weeks. Psychological coping strategy was found to modify vulnerability to the adverse effects of these events. The second study revealed that while men had a lower immediate risk of disorder onset following loss event experience their risk period was greater than for women. Time-varying exposure effects were well described by the appropriate use of simple time-dependent covariates.