Outlier effects on robust joint modelling of longitudinal and survival data

Robust joint modelling is an emerging field of research. Through the advancements in electronic patient healthcare records, the popularly of joint modelling approaches has grown rapidly in recent years providing simultaneous analysis of longitudinal and survival data. This research advances previous work through the development of a novel robust joint modelling methodology for one of the most common types of standard joint models, that which links a linear mixed model with a Cox proportional hazards model. Through t-distributional assumptions, longitudinal outliers are accommodated with their detrimental impact being down weighed and thus providing more efficient and reliable estimates. The robust joint modelling technique and its major benefits are showcased through the analysis of Northern Irish end stage renal disease patients. With an ageing population and growing prevalence of chronic kidney disease within the United Kingdom, there is a pressing demand to investigate the detrimental relationship between the changing haemoglobin levels of haemodialysis patients and their survival. As outliers within the NI renal data were found to have significantly worse survival, identification of outlying individuals through robust joint modelling may aid nephrologists to improve patient's survival. A simulation study was also undertaken to explore the difference between robust and standard joint models in the presence of increasing proportions and extremity of longitudinal outliers. More efficient and reliable estimates were obtained by robust joint models with increasing contrast between the robust and standard joint models when a greater proportion of more extreme outliers are present. Through illustration of the gains in efficiency and reliability of parameters when outliers exist, the potential of robust joint modelling is evident. The research presented in this thesis highlights the benefits and stresses the need to utilise a more robust approach to joint modelling in the presence of longitudinal outliers.

Influence of Remote Monitoring on Long-Term Cardiovascular Outcomes after Cardioverter-Defibrillator Implantation

AIMS: Device-based remote monitoring (RM) has been linked to improved clinical outcomes at short to medium-term follow-up. Whether this benefit extends to long-term follow-up is unknown. We sought to assess the effect of device-based RM on long-term clinical outcomes in recipients of implantable cardioverter-defibrillators (ICD). METHODS: We performed a retrospective cohort study of consecutive patients who underwent ICD implantation for primary prevention. RM was initiated with patient consent according to availability of RM hardware at implantation. Patients with concomitant cardiac resynchronization therapy were excluded. Data on hospitalizations, mortality and cause of death were systematically assessed using a nationwide healthcare platform. A Cox proportional hazards model was employed to estimate the effect of RM on mortality and a composite endpoint of cardiovascular mortality and hospital admission due to heart failure (HF). RESULTS: 312 patients were included with a median follow-up of 37.7months (range 1 to 146). 121 patients (38.2%) were under RM since the first outpatient visit post-ICD and 191 were in conventional follow-up. No differences were found regarding age, left ventricular ejection fraction, heart failure etiology or NYHA class at implantation. Patients under RM had higher long-term survival (hazard ratio [HR] 0.50, CI 0.27-0.93, p=0.029) and lower incidence of the composite outcome (HR 0.47, CI 0.27-0.82, p=0.008). After multivariate survival analysis, overall survival was independently associated with younger age, higher LVEF, NYHA class lower than 3 and RM. CONCLUSION: RM was independently associated with increased long-term survival and a lower incidence of a composite endpoint of hospitalization for HF or cardiovascular mortality.

Impact on hip fracture mortality after the establishment of an orthogeriatric care program in a colombian hospital

El objetivo del estudio es evaluar la mortalidad a un año en pacientes con fractura de cadera, mayores de 65 años tratados en un programa establecido de orto-geriatría. 298 se trataron de acuerdo al protocolo de orto-geriatría, se calculo la mortalidad a un año, se establecieron los predictores de mortalidad orto-geriátrico. La sobrevida anual se incremento de 80% a 89% (p = .039) durante los cuatro años de seguimiento del programa y disminuyo el riesgo de mortalidad anual postoperatorio (Hazard Ratio = 0.54, p = .049). La enfermedad cardiaca y la edad maor a 85 años fueron predictores positivos para mortalidad.

A semi-Markov multistate model for estimation of the mean quality-adjusted survival for non-progressive processes

We discuss the estimation of the expected value of the quality-adjusted survival, based on multistate models. We generalize an earlier work, considering the sojourn times in health states are not identically distributed, for a given vector of covariates. Approaches based on semiparametric and parametric (exponential and Weibull distributions) methodologies are considered. A simulation study is conducted to evaluate the performance of the proposed estimator and the jackknife resampling method is used to estimate the variance of such estimator. An application to a real data set is also included.

Estimation of the mean quality-adjusted survival using a multistate model for the sojourn times

In clinical trials, it may be of interest taking into account physical and emotional well-being in addition to survival when comparing treatments. Quality-adjusted survival time has the advantage of incorporating information about both survival time and quality-of-life. In this paper, we discuss the estimation of the expected value of the quality-adjusted survival, based on multistate models for the sojourn times in health states. Semiparametric and parametric (with exponential distribution) approaches are considered. A simulation study is presented to evaluate the performance of the proposed estimator and the jackknife resampling method is used to compute bias and variance of the estimator. (C) 2007 Elsevier B.V. All rights reserved.

A general threshold stress hybrid hazard model for lifetime data

In this paper we propose a hybrid hazard regression model with threshold stress which includes the proportional hazards and the accelerated failure time models as particular cases. To express the behavior of lifetimes the generalized-gamma distribution is assumed and an inverse power law model with a threshold stress is considered. For parameter estimation we develop a sampling-based posterior inference procedure based on Markov Chain Monte Carlo techniques. We assume proper but vague priors for the parameters of interest. A simulation study investigates the frequentist properties of the proposed estimators obtained under the assumption of vague priors. Further, some discussions on model selection criteria are given. The methodology is illustrated on simulated and real lifetime data set.

Independent predictors and a prognostic model for surgical outcome in refractory frontal lobe epilepsy

Purpose: Refractory frontal lobe epilepsy (FLE) remains one of the most challenging surgically remediable epilepsy syndromes. Nevertheless, definition of independent predictors and predictive models of postsurgical seizure outcome remains poorly explored in FLE. Methods: We retrospectively analyzed data from 70 consecutive patients with refractory FLE submitted to surgical treatment at our center from July 1994 to December 2006. Univariate results were submitted to logistic regression models and Cox proportional hazards regression to identify isolated risk factors for poor surgical results and to construct predictive models for surgical outcome in FLE. Results: From 70 patients submitted to surgery, 45 patients (64%) had favorable outcome and 37 (47%) became seizure free. Isolated risk factors for poor surgical outcome are expressed in hazard ratio (H.R.) and were time of epilepsy (H.R.=4.2; 95% C.I.=.1.5-11.7; p=0.006), ictal EEG recruiting rhythm (H.R. = 2.9; 95% C.I. = 1.1-7.7; p=0.033); normal MRI (H.R. = 4.8; 95% C.I. = 1.4-16.6; p = 0.012), and MRI with lesion involving eloquent cortex (H.R. = 3.8; 95% C.I. = 1.2-12.0; p = 0.021). Based on these variables and using a logistic regression model we constructed a model that correctly predicted long-term surgical outcome in up to 80% of patients. Conclusion: Among independent risk factors for postsurgical seizure outcome, epilepsy duration is a potentially modifiable factor that could impact surgical outcome in FLE. Early diagnosis, presence of an MRI lesion not involving eloquent cortex, and ictal EEG without recruited rhythm independently predicted favorable outcome in this series. (C) 2011 Elsevier B.V. All rights reserved.

A Cox Model for Biostatistics of the Future

Professor Sir David R. Cox (DRC) is widely acknowledged as among the most important scientists of the second half of the twentieth century. He inherited the mantle of statistical science from Pearson and Fisher, advanced their ideas, and translated statistical theory into practice so as to forever change the application of statistics in many fields, but especially biology and medicine. The logistic and proportional hazards models he substantially developed, are arguably among the most influential biostatistical methods in current practice. This paper looks forward over the period from DRC's 80th to 90th birthdays, to speculate about the future of biostatistics, drawing lessons from DRC's contributions along the way. We consider "Cox's model" of biostatistics, an approach to statistical science that: formulates scientific questions or quantities in terms of parameters gamma in probability models f(y; gamma) that represent in a parsimonious fashion, the underlying scientific mechanisms (Cox, 1997); partition the parameters gamma = theta, eta into a subset of interest theta and other "nuisance parameters" eta necessary to complete the probability distribution (Cox and Hinkley, 1974); develops methods of inference about the scientific quantities that depend as little as possible upon the nuisance parameters (Barndorff-Nielsen and Cox, 1989); and thinks critically about the appropriate conditional distribution on which to base infrences. We briefly review exciting biomedical and public health challenges that are capable of driving statistical developments in the next decade. We discuss the statistical models and model-based inferences central to the CM approach, contrasting them with computationally-intensive strategies for prediction and inference advocated by Breiman and others (e.g. Breiman, 2001) and to more traditional design-based methods of inference (Fisher, 1935). We discuss the hierarchical (multi-level) model as an example of the future challanges and opportunities for model-based inference. We then consider the role of conditional inference, a second key element of the CM. Recent examples from genetics are used to illustrate these ideas. Finally, the paper examines causal inference and statistical computing, two other topics we believe will be central to biostatistics research and practice in the coming decade. Throughout the paper, we attempt to indicate how DRC's work and the "Cox Model" have set a standard of excellence to which all can aspire in the future.

Bayesian and survival model for tumor relapse status and disease-specific survival, with applications to breast cancer

Breast cancer is the most common non-skin cancer and the second leading cause of cancer-related death in women in the United States. Studies on ipsilateral breast tumor relapse (IBTR) status and disease-specific survival will help guide clinic treatment and predict patient prognosis.^ After breast conservation therapy, patients with breast cancer may experience breast tumor relapse. This relapse is classified into two distinct types: true local recurrence (TR) and new ipsilateral primary tumor (NP). However, the methods used to classify the relapse types are imperfect and are prone to misclassification. In addition, some observed survival data (e.g., time to relapse and time from relapse to death)are strongly correlated with relapse types. The first part of this dissertation presents a Bayesian approach to (1) modeling the potentially misclassified relapse status and the correlated survival information, (2) estimating the sensitivity and specificity of the diagnostic methods, and (3) quantify the covariate effects on event probabilities. A shared frailty was used to account for the within-subject correlation between survival times. The inference was conducted using a Bayesian framework via Markov Chain Monte Carlo simulation implemented in softwareWinBUGS. Simulation was used to validate the Bayesian method and assess its frequentist properties. The new model has two important innovations: (1) it utilizes the additional survival times correlated with the relapse status to improve the parameter estimation, and (2) it provides tools to address the correlation between the two diagnostic methods conditional to the true relapse types.^ Prediction of patients at highest risk for IBTR after local excision of ductal carcinoma in situ (DCIS) remains a clinical concern. The goals of the second part of this dissertation were to evaluate a published nomogram from Memorial Sloan-Kettering Cancer Center, to determine the risk of IBTR in patients with DCIS treated with local excision, and to determine whether there is a subset of patients at low risk of IBTR. Patients who had undergone local excision from 1990 through 2007 at MD Anderson Cancer Center with a final diagnosis of DCIS (n=794) were included in this part. Clinicopathologic factors and the performance of the Memorial Sloan-Kettering Cancer Center nomogram for prediction of IBTR were assessed for 734 patients with complete data. Nomogram for prediction of 5- and 10-year IBTR probabilities were found to demonstrate imperfect calibration and discrimination, with an area under the receiver operating characteristic curve of .63 and a concordance index of .63. In conclusion, predictive models for IBTR in DCIS patients treated with local excision are imperfect. Our current ability to accurately predict recurrence based on clinical parameters is limited.^ The American Joint Committee on Cancer (AJCC) staging of breast cancer is widely used to determine prognosis, yet survival within each AJCC stage shows wide variation and remains unpredictable. For the third part of this dissertation, biologic markers were hypothesized to be responsible for some of this variation, and the addition of biologic markers to current AJCC staging were examined for possibly provide improved prognostication. The initial cohort included patients treated with surgery as first intervention at MDACC from 1997 to 2006. Cox proportional hazards models were used to create prognostic scoring systems. AJCC pathologic staging parameters and biologic tumor markers were investigated to devise the scoring systems. Surveillance Epidemiology and End Results (SEER) data was used as the external cohort to validate the scoring systems. Binary indicators for pathologic stage (PS), estrogen receptor status (E), and tumor grade (G) were summed to create PS+EG scoring systems devised to predict 5-year patient outcomes. These scoring systems facilitated separation of the study population into more refined subgroups than the current AJCC staging system. The ability of the PS+EG score to stratify outcomes was confirmed in both internal and external validation cohorts. The current study proposes and validates a new staging system by incorporating tumor grade and ER status into current AJCC staging. We recommend that biologic markers be incorporating into revised versions of the AJCC staging system for patients receiving surgery as the first intervention.^ Chapter 1 focuses on developing a Bayesian method to solve misclassified relapse status and application to breast cancer data. Chapter 2 focuses on evaluation of a breast cancer nomogram for predicting risk of IBTR in patients with DCIS after local excision gives the statement of the problem in the clinical research. Chapter 3 focuses on validation of a novel staging system for disease-specific survival in patients with breast cancer treated with surgery as the first intervention. ^

A joint hazard and time scaling model to compare survival curves.

To provide a more general method for comparing survival experience, we propose a model that independently scales both hazard and time dimensions. To test the curve shape similarity of two time-dependent hazards, h1(t) and h2(t), we apply the proposed hazard relationship, h12(tKt)/ h1(t) = Kh, to h1. This relationship doubly scales h1 by the constant hazard and time scale factors, Kh and Kt, producing a transformed hazard, h12, with the same underlying curve shape as h1. We optimize the match of h12 to h2 by adjusting Kh and Kt. The corresponding survival relationship S12(tKt) = [S1(t)]KtKh transforms S1 into a new curve S12 of the same underlying shape that can be matched to the original S2. We apply this model to the curves for regional and local breast cancer contained in the National Cancer Institute's End Results Registry (1950-1973). Scaling the original regional curves, h1 and S1 with Kt = 1.769 and Kh = 0.263 produces transformed curves h12 and S12 that display congruence with the respective local curves, h2 and S2. This similarity of curve shapes suggests the application of the more complete curve shapes for regional disease as templates to predict the long-term survival pattern for local disease. By extension, this similarity raises the possibility of scaling early data for clinical trial curves according to templates of registry or previous trial curves, projecting long-term outcomes and reducing costs. The proposed model includes as special cases the widely used proportional hazards (Kt = 1) and accelerated life (KtKh = 1) models.

A diszkrét kiválasztási modell becslése Cox-regresszióval (Parameter estimation for multinomial discrete choice model using Cox-regression)

Considering the so-called "multinomial discrete choice" model the focus of this paper is on the estimation problem of the parameters. Especially, the basic question arises how to carry out the point and interval estimation of the parameters when the model is mixed i.e. includes both individual and choice-specific explanatory variables while a standard MDC computer program is not available for use. The basic idea behind the solution is the use of the Cox-proportional hazards method of survival analysis which is available in any standard statistical package and provided a data structure satisfying certain special requirements it yields the MDC solutions desired. The paper describes the features of the data set to be analysed.

Regression models for grouped survival data: Estimation and sensitivity analysis

In this study, regression models are evaluated for grouped survival data when the effect of censoring time is considered in the model and the regression structure is modeled through four link functions. The methodology for grouped survival data is based on life tables, and the times are grouped in k intervals so that ties are eliminated. Thus, the data modeling is performed by considering the discrete models of lifetime regression. The model parameters are estimated by using the maximum likelihood and jackknife methods. To detect influential observations in the proposed models, diagnostic measures based on case deletion, which are denominated global influence, and influence measures based on small perturbations in the data or in the model, referred to as local influence, are used. In addition to those measures, the local influence and the total influential estimate are also employed. Various simulation studies are performed and compared to the performance of the four link functions of the regression models for grouped survival data for different parameter settings, sample sizes and numbers of intervals. Finally, a data set is analyzed by using the proposed regression models. (C) 2010 Elsevier B.V. All rights reserved.

Long-term survival of women with breast cancer in New South Wales

Several long-term studies of breast cancer survival have shown continued excess mortality from breast cancer up to 20-40 years following treatment. The purpose of this report was to investigate temporal trends in long-term survival from breast cancer in all New South Wales (NSW) women. Breast cancer cases incident in 1972-1996 (54,228) were derived from the NSW Central Cancer Registry a population-based registry which began in 1972. All cases of breast cancer not known to be dead were matched against death records. The expected survival for NSW women was derived from published annual life tables. Relative survival analysis compared the survival of cancer cases with the age, sex and period matched mortality of the total population. Cases were considered alive at the end of 1996, except when known to be dead. Proportional hazards regression was employed to model survival on age, period and degree of spread at diagnosis. Survival at 5, 10, 15, 20 and 25 years of follow-up was 76 per cent, 65 per cent, 60 per cent, 57 per cent and 56 per cent. The annual hazard rate for excess mortality was 4.3 per cent in year 1, maximal at 6.5 per cent in year 3, declining to 4.7 per cent in year 5, 2.7 per cent in year 10, 1.4 per cent in year 15, 1.0 per cent for years 16-20, and 0.4 per cent for years 20-25 of follow-up. Relative survival was highest in 40-49 year-olds. Cases diagnosed most recently (1992-1996) had the highest survival, compared with cases diagnosed in previous periods. Five-year survival improved over time, especially from the late 1980s for women in the screening age group (50-69 years). Survival was highest for those with localised cancer at diagnosis: 88.4 per cent, 79.1 per cent, 74.6 per cent, 72.7 per cent and 72.8 per cent at 5, 10, 15, 20 and 25 years follow-up (excluding those aged greater than or equal to 70 years). There was no significant difference between the survival of the breast cancer cases and the general population at 20-25 years follow-up. Degree of spread was less predictive of survival 5-20 years after diagnosis, compared with 0-5 years after diagnosis, and was not significant at 20-25 years of follow-up. Relative survival from breast cancer in NSW women continues to decrease to 25 years after diagnosis, but there is little excess mortality after 15 years follow-up, especially for those with localised cancer at diagnosis, and the minimal excess mortality at 20-25 years of follow-up is not statistically significant. (C) 2002 Elsevier Science Ltd. All rights reserved.

Modelling the Distribution of Ischaemic Stroke-Specific Survival Time Using an EM-based Mixture Approach with Random Effects Adjustment

A two-component survival mixture model is proposed to analyse a set of ischaemic stroke-specific mortality data. The survival experience of stroke patients after index stroke may be described by a subpopulation of patients in the acute condition and another subpopulation of patients in the chronic phase. To adjust for the inherent correlation of observations due to random hospital effects, a mixture model of two survival functions with random effects is formulated. Assuming a Weibull hazard in both components, an EM algorithm is developed for the estimation of fixed effect parameters and variance components. A simulation study is conducted to assess the performance of the two-component survival mixture model estimators. Simulation results confirm the applicability of the proposed model in a small sample setting. Copyright (C) 2004 John Wiley Sons, Ltd.

Impact of number of vessels disease on outcome of patients with stable coronary artery disease: 5-year follow-up of the medical, angioplasty, and bypass surgery study (MASS)

Objective: To evaluate whether the number of vessels disease has an impact on clinical outcomes as well as on therapeutic results accordingly to medical, percutaneous, or surgery treatment in chronic coronary artery disease. Methods: We evaluated 825 individuals enrolled in MASS study, a randomized study to compare treatment options for single or multivessel coronary artery disease with preserved left ventricular function, prospectively followed during 5 years. The incidence of overall mortality and the composite end-point of death, myocardial infarction, and refractory angina were compared in three groups: single vessel disease (SVD n = 214), two-vessel disease (2VD n = 253) and three-vessel disease (3VD n = 358). The relationship between baseline variables and the composite end-point was assessed using a Cox proportional hazards survival model. Results: Most baseline characteristics were similar among groups, except age (younger in SVD and older in 3VD, p < 0.001), lower incidence of hypertension in SVD (p < 0.0001), and lower levels of total and LDL-cholesterol in 3VD (p = 0.004 and p = 0.005, respectively). There were no statistical differences in composite end-point in 5 years among groups independent of the kind of treatment; however, there was a higher mortality rate in 3VD (p < 0.001). When we stratified our analysis for each treatment option, bypass surgery was associated with a tower number of composite end-point in all groups (SVD p < 0.001, 2VD p = 0.002, 3VD p < 0.001). In multivariate analysis, we found higher mortality risk in 3VD comparing to SVD (p = 0.005, HR 3.14, 95%Cl 1.4-7.0). Conclusion: Three-vessel disease was associated with worse prognosis compared to single-or two-vessel disease in patients with stable coronary disease and preserved ventricular function at 5-year follow-up. In addition, event-free survival rates were higher after bypass surgery, independent of the number of vessels diseased in these subsets of patients. (c) 2008 European Association for Cardio-Thoracic Surgery. Published by Elsevier B.V. All rights reserved.