Inference on Survival Data with Covariate Measurement Error - An Imputation-based Approach

We propose a new method for fitting proportional hazards models with error-prone covariates. Regression coefficients are estimated by solving an estimating equation that is the average of the partial likelihood scores based on imputed true covariates. For the purpose of imputation, a linear spline model is assumed on the baseline hazard. We discuss consistency and asymptotic normality of the resulting estimators, and propose a stochastic approximation scheme to obtain the estimates. The algorithm is easy to implement, and reduces to the ordinary Cox partial likelihood approach when the measurement error has a degenerative distribution. Simulations indicate high efficiency and robustness. We consider the special case where error-prone replicates are available on the unobserved true covariates. As expected, increasing the number of replicate for the unobserved covariates increases efficiency and reduces bias. We illustrate the practical utility of the proposed method with an Eastern Cooperative Oncology Group clinical trial where a genetic marker, c-myc expression level, is subject to measurement error.

Age-specific and sex-specific weight gain norms to monitor antiretroviral therapy in children in low-income and middle-income countries.

BACKGROUND Viral load and CD4% are often not available in resource-limited settings for monitoring children's responses to antiretroviral therapy (ART). We aimed to construct normative curves for weight gain at 6, 12, 18, and 24 months following initiation of ART in children, and to assess the association between poor weight gain and subsequent responses to ART. DESIGN Analysis of data from HIV-infected children younger than 10 years old from African and Asian clinics participating in the International epidemiologic Databases to Evaluate AIDS. METHODS The generalized additive model for location, scale, and shape was used to construct normative percentile curves for weight gain at 6, 12, 18, and 24 months following ART initiation. Cox proportional models were used to assess the association between lower percentiles (< 50th) of weight gain distribution at the different time points and subsequent death, virological suppression, and virological failure. RESULTS Among 7173 children from five regions of the world, 45% were underweight at baseline. Weight gain below the 50th percentile at 6, 12, 18, and 24 months of ART was associated with increased risk of death, independent of baseline characteristics. Poor weight gain was not associated with increased hazards of virological suppression or virological failure. CONCLUSION Monitoring weight gain on ART using age-specific and sex-specific normative curves specifically developed for HIV-infected children on ART is a simple, rapid, sustainable tool that can aid in the identification of children who are at increased risk of death in the first year of ART.

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.

Factors determining disease duration in Alzheimer's disease:a postmortem study of 103 cases using the Kaplan-Meier estimator and Cox regression

Factors associated with duration of dementia in a consecutive series of 103 Alzheimer's disease (AD) cases were studied using the Kaplan-Meier estimator and Cox regression analysis (proportional hazard model). Mean disease duration was 7.1 years (range: 6 weeks-30 years, standard deviation = 5.18); 25% of cases died within four years, 50% within 6.9 years, and 75% within 10 years. Familial AD cases (FAD) had a longer duration than sporadic cases (SAD), especially cases linked to presenilin (PSEN) genes. No significant differences in duration were associated with age, sex, or apolipoprotein E (Apo E) genotype. Duration was reduced in cases with arterial hypertension. Cox regression analysis suggested longer duration was associated with an earlier disease onset and increased senile plaque (SP) and neurofibrillary tangle (NFT) pathology in the orbital gyrus (OrG), CA1 sector of the hippocampus, and nucleus basalis of Meynert (NBM). The data suggest shorter disease duration in SAD and in cases with hypertensive comorbidity. In addition, degree of neuropathology did not influence survival, but spread of SP/NFT pathology into the frontal lobe, hippocampus, and basal forebrain was associated with longer disease duration. © 2014 R. A. Armstrong.

Aplicação do escore MELD em pacientes submetidos a transplante de fígado: análise retrospectiva da sobrevida e dos fatores preditivos a curto e longo prazo

BACKGROUND: The model for end-stage liver disease (MELD) was developed to predict short-term mortality in patients with cirrhosis. There are few reports studying the correlation between MELD and long-term posttransplantation survival. AIM: To assess the value of pretransplant MELD in the prediction of posttransplant survival. METHODS: The adult patients (age >18 years) who underwent liver transplantation were examined in a retrospective longitudinal cohort of patients, through the prospective data base. We excluded acute liver failure, retransplantation and reduced or split-livers. The liver donors were evaluated according to: age, sex, weight, creatinine, bilirubin, sodium, aspartate aminotransferase, personal antecedents, brain death cause, steatosis, expanded criteria donor number and index donor risk. The recipients' data were: sex, age, weight, chronic hepatic disease, Child-Turcotte-Pugh points, pretransplant and initial MELD score, pretransplant creatinine clearance, sodium, cold and warm ischemia times, hospital length of stay, blood requirements, and alanine aminotransferase (ALT >1,000 UI/L = liver dysfunction). The Kaplan-Meier method with the log-rank test was used for the univariable analyses of posttransplant patient survival. For the multivariable analyses the Cox proportional hazard regression method with the stepwise procedure was used with stratifying sodium and MELD as variables. ROC curve was used to define area under the curve for MELD and Child-Turcotte-Pugh. RESULTS: A total of 232 patients with 10 years follow up were available. The MELD cutoff was 20 and Child-Turcotte-Pugh cutoff was 11.5. For MELD score > 20, the risk factors for death were: red cell requirements, liver dysfunction and donor's sodium. For the patients with hyponatremia the risk factors were: negative delta-MELD score, red cell requirements, liver dysfunction and donor's sodium. The regression univariated analyses came up with the following risk factors for death: score MELD > 25, blood requirements, recipient creatinine clearance pretransplant and age donor >50. After stepwise analyses, only red cell requirement was predictive. Patients with MELD score < 25 had a 68.86%, 50,44% and 41,50% chance for 1, 5 and 10-year survival and > 25 were 39.13%, 29.81% and 22.36% respectively. Patients without hyponatremia were 65.16%, 50.28% and 41,98% and with hyponatremia 44.44%, 34.28% and 28.57% respectively. Patients with IDR > 1.7 showed 53.7%, 27.71% and 13.85% and index donor risk <1.7 was 63.62%, 51.4% and 44.08%, respectively. Age donor > 50 years showed 38.4%, 26.21% and 13.1% and age donor <50 years showed 65.58%, 26.21% and 13.1%. Association with delta-MELD score did not show any significant difference. Expanded criteria donors were associated with primary non-function and severe liver dysfunction. Predictive factors for death were blood requirements, hyponatremia, liver dysfunction and donor's sodium. CONCLUSION: In conclusion MELD over 25, recipient's hyponatremia, blood requirements, donor's sodium were associated with poor survival.

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.

A Dynamic Model of Hand-and-Foot Syndrome in Patients Receiving Capecitabine

For the purpose of developing a longitudinal model to predict hand-and-foot syndrome (HFS) dynamics in patients receiving capecitabine, data from two large phase III studies were used. Of 595 patients in the capecitabine arms, 400 patients were randomly selected to build the model, and the other 195 were assigned for model validation. A score for risk of developing HFS was modeled using the proportional odds model, a sigmoidal maximum effect model driven by capecitabine accumulation as estimated through a kinetic-pharmacodynamic model and a Markov process. The lower the calculated creatinine clearance value at inclusion, the higher was the risk of HFS. Model validation was performed by visual and statistical predictive checks. The predictive dynamic model of HFS in patients receiving capecitabine allows the prediction of toxicity risk based on cumulative capecitabine dose and previous HFS grade. This dose-toxicity model will be useful in developing Bayesian individual treatment adaptations and may be of use in the clinic.

Myocardial Perfusion Imaging Is a Strong Predictor of Death in Women

OBJECTIVES We sought to assess the prognostic value and risk classification improvement using contemporary single-photon emission computed tomography myocardial perfusion imaging (SPECT-MPI) to predict all-cause mortality. BACKGROUND Myocardial perfusion is a strong estimator of prognosis. Evidence published to date has not established the added prognostic value of SPECT-MPI nor defined an approach to detect improve classification of risk in women from a developing nation. METHODS A total of 2,225 women referred for SPECT-MPI were followed by a mean period of 3.7 +/- 1.4 years. SPECT-MPI results were classified as abnormal on the presence of any perfusion defect. Abnormal scans were further classified as with mild/moderate reversible, severe reversible, partial reversible, or fixed perfusion defects. Risk estimates for incident mortality were categorized as <1%/year, 1% to 2%/year, and >2%/year using Cox proportional hazard models. Risk-adjusted models incorporated clinical risk factors, left ventricular ejection fraction (LVEF), and perfusion variables. RESULTS All-cause death occurred in 139 patients. SPECT-MPI significantly risk stratified the population; patients with abnormal scans had significantly higher death rates compared with patients with normal scans, 13.1% versus 4.0%, respectively (p < 0.001). Cox analysis demonstrated that after adjusting for clinical risk factors and LVEF, SPECT-MPI improved the model discrimination (integrated discrimination index = 0.009; p = 0.02), added significant incremental prognostic information (global chi-square increased from 87.7 to 127.1; p < 0.0001), and improved risk prediction (net reclassification improvement = 0.12; p = 0.005). CONCLUSIONS SPECT-MPI added significant incremental prognostic information to clinical and left ventricular functional variables while enhancing the ability to classify this Brazilian female population into low-and high-risk categories of all-cause mortality. (J Am Coll Cardiol Img 2011;4:880-8) (C) 2011 by the American College of Cardiology Foundation

Coronary Revascularization Induces a Shift From Cardiac Toward Noncardiac Mortality Without Improving Survival in Vascular Surgery Patient

OBJECTIVE: Although evidence has shown that ischemic heart disease (IHD) in vascular surgery patients has a negative impact on the prognosis after surgery, it is unclear whether directed treatment of IHD may influence cause-specific and overall mortality. The objective of this study was to determine the prognostic implication of coronary revascularization (CR) on overall and cause-specific mortality in vascular surgery patients. METHODS: Patients undergoing surgery for abdominal aortic aneurysm, carotid artery stenosis, or peripheral artery disease in a university hospital in The Netherlands between January 2003 and December 2011 were retrospectively included. Survival estimates were obtained by Kaplan-Meier and Cox regression analysis. RESULTS: A total of 1104 patients were included. Adjusted survival analyses showed that IHD significantly increased the risk of overall mortality (hazard ratio [HR], 1.50; 95% confidence interval, 1.21-1.87) and cardiovascular death (HR, 1.93; 95% confidence interval, 1.35-2.76). Compared with those without CR, patients previously undergoing CR had similar overall mortality (HR, 1.38 vs 1.62; P = .274) and cardiovascular mortality (HR, 1.83 vs 2.02; P = .656). Nonrevascularized IHD patients were more likely to die of IHD (6.9% vs 35.7%), whereas revascularized IHD patients more frequently died of cardiovascular causes unrelated to IHD (39.1% vs 64.3%; P = .018). CONCLUSIONS: This study confirms the significance of IHD for postoperative survival of vascular surgery patients. CR was associated with lower IHD-related death rates. However, it failed to provide an overall survival benefit because of an increased rate of cardiovascular mortality unrelated to IHD. Intensification of secondary prevention regimens may be required to prevent this shift toward non-IHD-related death and thereby improve life expectancy.

Adalimumab in acute sciatica reduces the long-term need for surgery: a 3-year follow-up of a randomised double-blind placebo-controlled trial.

INTRODUCTION: Two subcutaneous injections of adalimumab in severe acute sciatica significantly reduced the number of back operations in a short-term randomised controlled clinical trial. OBJECTIVE: To determine in a 3-year follow-up study whether the short-term benefit of adalimumab in sciatica is sustained over a longer period of time. METHODS: The primary outcome of this analysis was incident discectomy. Three years after randomisation, information on surgery could be retrieved in 56/61 patients (92%).A multivariate Cox proportional hazard models, adjusted for potential confounders, was used to determine factors predisposing to surgery. RESULTS: Twenty-three (41%) patients had back surgery within 3 years, 8/29 (28%) in the adalimumab group and 15/27 (56%) in the placebo group, p=0.04. Adalimumab injections reduced the need for back surgery by 61% (HR)=0.39 (95% CI 0.17 to 0.92). In a multivariate model, treatment with a tumour necrosis factor-α antagonist remained the strongest protective factor (HR=0.17, p=0.002). Other significant predictors of surgery were a good correlation between symptoms and MRI findings (HR=11.6, p=0.04), baseline intensity of leg pain (HR=1.3, p=0.06), intensity of back pain (HR=1.4, p=0.03) and duration of sickness leave (HR=1.01 per day, p=0.03). CONCLUSION: A short course of adalimumab in patients with severe acute sciatica significantly reduces the need for back surgery.

Improving appropriateness of antibiotic therapy: randomized trial of an intervention to foster reassessment of prescription after 3 days.

OBJECTIVES: Reassessment of ongoing antibiotic therapy is an important step towards appropriate use of antibiotics. This study was conducted to evaluate the impact of a short questionnaire designed to encourage reassessment of intravenous antibiotic therapy after 3 days. PATIENTS AND METHODS: Patients hospitalized on the surgical and medical wards of a university hospital and treated with an intravenous antibiotic for 3-4 days were randomly allocated to either an intervention or control group. The intervention consisted of mailing to the physician in charge of the patient a three-item questionnaire referring to possible adaptation of the antibiotic therapy. The primary outcome was the time elapsed from randomization until a first modification of the initial intravenous antibiotic therapy. It was compared within both groups using Cox proportional-hazard modelling. RESULTS: One hundred and twenty-six eligible patients were randomized in the intervention group and 125 in the control group. Time to modification of intravenous antibiotic therapy was 14% shorter in the intervention group (adjusted hazard ratio for modification 1.28, 95% CI 0.99-1.67, P = 0.06). It was significantly shorter in the intervention group compared with a similar group of 151 patients observed during a 2 month period preceding the study (adjusted hazard ratio 1.17, 95% CI 1.03-1.32, P = 0.02). CONCLUSION: The results suggest that a short questionnaire, easily adaptable to automatization, has the potential to foster reassessment of antibiotic therapy.

Increased expression of urokinase-type plasminogen activator mRNA determines adverse prognosis in ErbB2-positive primary breast cancer.

PURPOSE: To evaluate and validate mRNA expression markers capable of identifying patients with ErbB2-positive breast cancer associated with distant metastasis and reduced survival. PATIENTS AND METHODS: Expression of 60 genes involved in breast cancer biology was assessed by quantitative real-time PCR (qrt-PCR) in 317 primary breast cancer patients and correlated with clinical outcome data. Results were validated subsequently using two previously published and publicly available microarray data sets with different patient populations comprising 295 and 286 breast cancer samples, respectively. RESULTS: Of the 60 genes measured by qrt-PCR, urokinase-type plasminogen activator (uPA or PLAU) mRNA expression was the most significant marker associated with distant metastasis-free survival (MFS) by univariate Cox analysis in patients with ErbB2-positive tumors and an independent factor in multivariate analysis. Subsequent validation in two microarray data sets confirmed the prognostic value of uPA in ErbB2-positive tumors by both univariate and multivariate analysis. uPA mRNA expression was not significantly associated with MFS in ErbB2-negative tumors. Kaplan-Meier analysis showed in all three study populations that patients with ErbB2-positive/uPA-positive tumors exhibited significantly reduced MFS (hazard ratios [HR], 4.3; 95% CI, 1.6 to 11.8; HR, 2.7; 95% CI, 1.2 to 6.2; and, HR, 2.8; 95% CI, 1.1 to 7.1; all P < .02) as compared with the group with ErbB2-positive/uPA-negative tumors who exhibited similar outcome to those with ErbB2-negative tumors, irrespective of uPA status. CONCLUSION: After evaluation of 898 breast cancer patients, uPA mRNA expression emerged as a powerful prognostic indicator in ErbB2-positive tumors. These results were consistent among three independent study populations assayed by different techniques, including qrt-PCR and two microarray platforms.

Anti-TNF therapy in acute sciatica reduces the long-term need for surgery: A three-year follow-up of a randomized double blind placebo controlled trial

Introduction: Two subcutaneous injections of adalimumab in severeacute sciatica have demonstrated a significant benefit on the numberof back surgeries in a short-term randomized controlled clinical trial[1]. This 3-year follow-up study aimed to determine whether theshort-term benefit was sustained over a longer period of time.Methods: Information on surgery was retrieved in 56/61 patients(93%). We used a Cox proportional hazard models to determinefactors predisposing to surgery.Results: Twenty-three (41%) patients had back surgery within 3 years,8/29 (28%) in the adalimumab group and 15/ 27 (56%) in the placebogroup, p = 0.038. Adalimumab injections reduced the need for backsurgery by 61% (Hazard Ratio (HR): 0.39 (95% CI: 0.17-0.92). In amultivariate model, treatment with a TNF-α antagonist remained thestrongest protective factor (HR 0.17, p = 0.002). Other significantpredictors of surgery were a good correlation between symptomsand MRI findings (HR = 11.6, p = 0.04), baseline intensity of leg pain(HR = 1.3, p = 0.06), intensity of back pain (HR = 1.4, p = 0.03)and duration of sickness leave (HR = 1.01 per day, p = 0.03).Conclusion: A short course of adalimumab in patients with severeacute sciatica significantly reduces the need for back surgery.