Survival Analysis of Patients with Heart Failure: Implications of Time-Varying Regression Effects in Modeling Mortality

Background: Several models have been designed to predict survival of patients with heart failure. These, while available and widely used for both stratifying and deciding upon different treatment options on the individual level, have several limitations. Specifically, some clinical variables that may influence prognosis may have an influence that change over time. Statistical models that include such characteristic may help in evaluating prognosis. The aim of the present study was to analyze and quantify the impact of modeling heart failure survival allowing for covariates with time-varying effects known to be independent predictors of overall mortality in this clinical setting. Methodology: Survival data from an inception cohort of five hundred patients diagnosed with heart failure functional class III and IV between 2002 and 2004 and followed-up to 2006 were analyzed by using the proportional hazards Cox model and variations of the Cox's model and also of the Aalen's additive model. Principal Findings: One-hundred and eighty eight (188) patients died during follow-up. For patients under study, age, serum sodium, hemoglobin, serum creatinine, and left ventricular ejection fraction were significantly associated with mortality. Evidence of time-varying effect was suggested for the last three. Both high hemoglobin and high LV ejection fraction were associated with a reduced risk of dying with a stronger initial effect. High creatinine, associated with an increased risk of dying, also presented an initial stronger effect. The impact of age and sodium were constant over time. Conclusions: The current study points to the importance of evaluating covariates with time-varying effects in heart failure models. The analysis performed suggests that variations of Cox and Aalen models constitute a valuable tool for identifying these variables. The implementation of covariates with time-varying effects into heart failure prognostication models may reduce bias and increase the specificity of such models.

MATHEMATICAL MODELS AND POLYHEDRAL STUDIES FOR INTEGRAL SHEET METAL DESIGN

We deal with the optimization of the production of branched sheet metal products. New forming techniques for sheet metal give rise to a wide variety of possible profiles and possible ways of production. In particular, we show how the problem of producing a given profile geometry can be modeled as a discrete optimization problem. We provide a theoretical analysis of the model in order to improve its solution time. In this context we give the complete convex hull description of some substructures of the underlying polyhedron. Moreover, we introduce a new class of facet-defining inequalities that represent connectivity constraints for the profile and show how these inequalities can be separated in polynomial time. Finally, we present numerical results for various test instances, both real-world and academic examples.

A Bayesian analysis of the Conway-Maxwell-Poisson cure rate model

The purpose of this paper is to develop a Bayesian analysis for the right-censored survival data when immune or cured individuals may be present in the population from which the data is taken. In our approach the number of competing causes of the event of interest follows the Conway-Maxwell-Poisson distribution which generalizes the Poisson distribution. Markov chain Monte Carlo (MCMC) methods are used to develop a Bayesian procedure for the proposed model. Also, some discussions on the model selection and an illustration with a real data set are considered.

Survival analysis and iniquities in older Brazilians: a six year follow up survey in São Paulo, Brazil

Maria Lucia Lebrão is the Coordinator of the SABE study. Jair LF Santos and Yeda AO Duarte receive support from National Council of Research (CNPq). The SABE study is supported by The São Paulo Research Foundation (FAPESP).

An in-depth analysis of the HIV-1/AIDS dynamics by comprehensive mathematical modeling

This work presents major results from a novel dynamic model intended to deterministically represent the complex relation between HIV-1 and the human immune system. The novel structure of the model extends previous work by representing different host anatomic compartments under a more in-depth cellular and molecular immunological phenomenology. Recently identified mechanisms related to HIV-1 infection as well as other well known relevant mechanisms typically ignored in mathematical models of HIV-1 pathogenesis and immunology, such as cell-cell transmission, are also addressed. (C) 2011 Elsevier Ltd. All rights reserved.

Survival, Causes of Death, and Prognostic Factors in Systemic Sclerosis: Analysis of 947 Brazilian Patients

Objective. To analyze survival, prognostic factors, and causes of death in a large cohort of patients with systemic sclerosis (SSc). Methods. From 1991 to 2010, 947 patients with SSc were treated at 2 referral university centers in Brazil. Causes of death were considered SSc-related and non-SSc-related. Multiple logistic regression analysis was used to identify prognostic factors. Survival at 5 and 10 years was estimated using the Kaplan-Meier method. Results. One hundred sixty-eight patients died during the followup. Among the 110 deaths considered related to SSc, there was predominance of lung (48.1%) and heart (24.5%) involvement. Most of the 58 deaths not related to SSc were caused by infection, cardiovascular or cerebrovascular disease, and cancer. Male sex, modified Rodnan skin score (mRSS) > 20, osteoarticular involvement, lung involvement, and renal crisis were the main prognostic factors associated to death. Overall survival rate was 90% for 5 years and 84% for 10 years. Patients presented worse prognosis if they had diffuse SSc (85% vs 92% at 5 yrs, respectively, and 77% vs 87% at 10 yrs, compared to limited SSc), male sex (77% vs 90% at 5 yrs and 64% vs 86% at 10 yrs, compared to female sex), and mRSS > 20 (83% vs 90% at 5 yrs and 66% vs 86% at 10 yrs, compared to mRSS <20). Conclusion. Survival was worse in male patients with diffuse SSc, and lung and heart involvement represented the main causes of death in this South American series of patients with SSc. (First Release Aug 15 2012; J Rheumatol 2012;39:1971-8; doi:10.3899/jrheum.111582)

Comparison of two experimental models of pulmonary hypertension

Objective: To compare two models of pulmonary hypertension (monocrotaline and monocrotaline+pneumonectomy) regarding hemodynamic severity, structure of pulmonary arteries, inflammatory markers (IL-1 and PDGF), and 45-day survival. Methods: We used 80 Sprague-Dawley rats in two study protocols: structural analysis; and survival analysis. The rats were divided into four groups: control; monocrotaline (M), pneumonectomy (P), and monocrotaline+pneumonectomy (M+P). In the structural analysis protocol, 40 rats (10/group) were catheterized for the determination of hemodynamic variables, followed by euthanasia for the removal of heart and lung tissue. The right ventricle (RV) was dissected from the interventricular septum (IS), and the ratio between RV weight and the weight of the left ventricle (LV) plus IS (RV/LV+IS) was taken as the index of RV hypertrophy. In lung tissues, we performed histological analyses, as well as using ELISA to determine IL-1 and PDGF levels. In the survival protocol, 40 animals (10/group) were followed for 45 days. Results: The M and M+P rats developed pulmonary hypertension, whereas the control and P rats did not. The RV/LV+IS ratio was significantly higher in M+P rats than in M rats, as well as being significantly higher in M and M+P rats than in control and P rats. There were no significant differences between the M and M+P rats regarding the area of the medial layer of the pulmonary arteries; IL-1 and PDGF levels; or survival. Conclusions: On the basis of our results, we cannot conclude that the monocrotaline+pneumonectomy model is superior to the monocrotaline model.

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.

Efficacy of brief alcohol screening intervention for college students (BASICS): a meta-analysis of randomized controlled trials

Background: Many studies reported that brief interventions are effective in reducing excessive drinking. This study aimed to assess the efficacy of a protocol of brief intervention for college students (BASICS), delivered face-to-face, to reduce risky alcohol consumption and negative consequences. Methods: A systematic review with meta-analysis was performed by searching for randomized controlled trials (RCTs) in Medline, PsycInfo, Web of Science and Cochrane Library databases. A quality assessment of RCTs was made by using a validated scale. Combined mean effect sizes, using meta-analysis random-effects models, were calculated. Results: 18 studies were included in the review. The sample sizes ranged from 54 to 1275 (median = 212). All studies presented a good evaluation of methodological quality and four were found to have excellent quality. After approximately 12 months of follow-up, students receiving BASICS showed a significant reduction in alcohol consumption (difference between means = -1.50 drinks per week, 95% CI: -3.24 to -0.29) and alcohol-related problems (difference between means = -0.87, 95% CI: -1.58 to -0.20) compared to controls. Conclusions: Overall, BASICS lowered both alcohol consumption and negative consequences in college students. Gender and peer factors seem to play an important role as moderators of behavior change in college drinking. Characteristics of BASICS procedure have been evaluated as more favorable and acceptable by students in comparison with others interventions or control conditions. Considerations for future researches were discussed.

Latent residual analysis in binary regression with skewed link

Model diagnostics is an integral part of model determination and an important part of the model diagnostics is residual analysis. We adapt and implement residuals considered in the literature for the probit, logistic and skew-probit links under binary regression. New latent residuals for the skew-probit link are proposed here. We have detected the presence of outliers using the residuals proposed here for different models in a simulated dataset and a real medical dataset.

The polysurvival model with long-term survivors

Long-term survival models have historically been considered for analyzing time-to-event data with long-term survivors fraction. However, situations in which a fraction (1 - p) of systems is subject to failure from independent competing causes of failure, while the remaining proportion p is cured or has not presented the event of interest during the time period of the study, have not been fully considered in the literature. In order to accommodate such situations, we present in this paper a new long-term survival model. Maximum likelihood estimation procedure is discussed as well as interval estimation and hypothesis tests. A real dataset illustrates the methodology.

The log-exponentiated generalized gamma regression model for censored data

For the first time, we introduce a generalized form of the exponentiated generalized gamma distribution [Cordeiro et al. The exponentiated generalized gamma distribution with application to lifetime data, J. Statist. Comput. Simul. 81 (2011), pp. 827-842.] that is the baseline for the log-exponentiated generalized gamma regression model. The new distribution can accommodate increasing, decreasing, bathtub- and unimodal-shaped hazard functions. A second advantage is that it includes classical distributions reported in the lifetime literature as special cases. We obtain explicit expressions for the moments of the baseline distribution of the new regression model. The proposed model can be applied to censored data since it includes as sub-models several widely known regression models. It therefore can be used more effectively in the analysis of survival data. We obtain maximum likelihood estimates for the model parameters by considering censored data. We show that our extended regression model is very useful by means of two applications to real data.

A 15-year study on the treatment of undernourished children at a nutrition rehabilitation centre (CREN), Brazil

Objective: To build a life table and determine the factors related to the time of treatment of undernourished children at a nutrition rehabilitation centre (CREN), Sao Paulo, Brazil. Design: Nutritional status was assessed from weight-for-age, height-for-age and BMI-for-age Z-scores, while neuropsychomotor development was classified according to the milestones of childhood development. Life tables, Kaplan-Meier survival curves and Cox multiple regression models were employed in data analysis. Setting: CREN (Centre of Nutritional Recovery and Education), Sao Paulo, Brazil. Subjects: Undernourished children (n 228) from the southern slums of Sao Paulo who had received treatment at CREN under a day-hospital regime between the years 1994 and 2009. Results: The Kaplan-Meier curves of survival analysis showed statistically significant differences in the periods of treatment at CREN between children presenting different degrees of neuropsychomotor development (log-rank = 6.621; P = 0.037). Estimates based on the multivariate Cox model revealed that children aged >= 24 months at the time of admission exhibited a lower probability of nutritional rehabilitation (hazard ratio (HR) = 0.49; P = 0.046) at the end of the period compared with infants aged up 12 months. Children presenting slow development were better rehabilitated in comparison with those exhibiting adequate evolution (HR = 4.48; P = 0.023). No significant effects of sex, degree of undernutrition or birth weight on the probability of nutritional rehabilitation were found. Conclusions: Age and neuropsychomotor developmental status at the time of admission to CREN are critical factors in determining the duration of treatment.

Influence diagnostics for elliptical semiparametric mixed models

In this paper we extend semiparametric mixed linear models with normal errors to elliptical errors in order to permit distributions with heavier and lighter tails than the normal ones. Penalized likelihood equations are applied to derive the maximum penalized likelihood estimates (MPLEs) which appear to be robust against outlying observations in the sense of the Mahalanobis distance. A reweighed iterative process based on the back-fitting method is proposed for the parameter estimation and the local influence curvatures are derived under some usual perturbation schemes to study the sensitivity of the MPLEs. Two motivating examples preliminarily analyzed under normal errors are reanalyzed considering some appropriate elliptical errors. The local influence approach is used to compare the sensitivity of the model estimates.

The tangential variation of a localized flux-type eigenvalue problem

In this work the differentiability of the principal eigenvalue lambda = lambda(1)(Gamma) to the localized Steklov problem -Delta u + qu = 0 in Omega, partial derivative u/partial derivative nu = lambda chi(Gamma)(x)u on partial derivative Omega, where Gamma subset of partial derivative Omega is a smooth subdomain of partial derivative Omega and chi(Gamma) is its characteristic function relative to partial derivative Omega, is shown. As a key point, the flux subdomain Gamma is regarded here as the variable with respect to which such differentiation is performed. An explicit formula for the derivative of lambda(1) (Gamma) with respect to Gamma is obtained. The lack of regularity up to the boundary of the first derivative of the principal eigenfunctions is a further intrinsic feature of the problem. Therefore, the whole analysis must be done in the weak sense of H(1)(Omega). The study is of interest in mathematical models in morphogenesis. (C) 2011 Elsevier Inc. All rights reserved.