Estimação clássica e Bayesiana em modelos de sobrevida com fração de cura

In Survival Analysis, long duration models allow for the estimation of the healing fraction, which represents a portion of the population immune to the event of interest. Here we address classical and Bayesian estimation based on mixture models and promotion time models, using different distributions (exponential, Weibull and Pareto) to model failure time. The database used to illustrate the implementations is described in Kersey et al. (1987) and it consists of a group of leukemia patients who underwent a certain type of transplant. The specific implementations used were numeric optimization by BFGS as implemented in R (base::optim), Laplace approximation (own implementation) and Gibbs sampling as implemented in Winbugs. We describe the main features of the models used, the estimation methods and the computational aspects. We also discuss how different prior information can affect the Bayesian estimates

Analysis of ensemble models in the medium term hydropower scheduling

The medium term hydropower scheduling (MTHS) problem involves an attempt to determine, for each time stage of the planning period, the amount of generation at each hydro plant which will maximize the expected future benefits throughout the planning period, while respecting plant operational constraints. Besides, it is important to emphasize that this decision-making has been done based mainly on inflow earliness knowledge. To perform the forecast of a determinate basin, it is possible to use some intelligent computational approaches. In this paper one considers the Dynamic Programming (DP) with the inflows given by their average values, thus turning the problem into a deterministic one which the solution can be obtained by deterministic DP (DDP). The performance of the DDP technique in the MTHS problem was assessed by simulation using the ensemble prediction models. Features and sensitivities of these models are discussed. © 2012 IEEE.

Serum Metalloproteinases 2 and 9 as Predictors of Gait Status, Pressure Ulcer and Mortality after Hip Fracture

Introduction: The aim of this study is to evaluate the serum activity of metalloproteinases (MMPs) -2 and -9 as predictors of pressure ulcer (PU), gait status and mortality 6 months after hip fracture. Methods: Eighty-seven patients over the age of 65 admitted to the orthopedic unit from January to December 2010 with hip fracture were prospectively evaluated. Upon admission, patient demographic information, including age, gender and concomitant diseases, was recorded. Blood samples were taken for analysis of MMP -2 and -9 activity by gel zymography and for biochemical examination within the first 72 hours of the patient's admission, after clinical stabilization. The fracture pattern (neck, trochanteric or subtrochanteric), time from admission to surgery, surgery duration and length of hospital stay were also recorded. Results: Two patients were excluded due to the presence of pathological fractures (related to cancer), and three patients were excluded due to the presence of PU before admission. Eighty-two patients, with a mean age of 80.4 ± 7.3 years, were included in the analysis. Among these patients, 75.6% were female, 59.8% had PU, and 13.4% died 6 months after hip fracture. All patients underwent hip fracture repair. In a univariate analysis, there were no differences in serum MMP activity between hip fracture patients with or without PU. In addition, the multiple logistic regression analysis models, which were adjusted by age, gender, length of hospital stay and C-reactive protein, showed that the pro-MMP-9 complexed with neutrophil gelatinase-associated lipocalin form (130 kDa) was associated with gait status recovery 6 months after hip fracture. Conclusions: In conclusion, serum pro-MMP-9 is a predictor of gait status recovery 6 months after hip fracture. © 2013 Gumieiro et al.

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, 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)

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.

The impact of psychiatric diagnosis on treatment adherence and duration among victimized children and adolescents in Sao Paulo, Brazil

OBJECTIVE: Despite the high prevalence of substance abuse and mood disorders among victimized children and adolescents, few studies have investigated the association of these disorders with treatment adherence, represented by numbers of visits per month and treatment duration. We aimed to investigate the effects of substance abuse and mood disorders on treatment adherence and duration in a special programfor victimized children in Sao Paulo, Brazil. METHODS: A total of 351 participants were evaluated for psychiatric disorders and classified into one of five groups: mood disorders alone; substance abuse disorders alone; mood and substance abuse disorders; other psychiatric disorders; no psychiatric disorders. The associations between diagnostic classification and adherence to treatment and the duration of program participation were tested with logistic regression and survival analysis, respectively. RESULTS: Children with mood disorders alone had the highest rate of adherence (79.5%); those with substance abuse disorders alone had the lowest (40%); and those with both disorders had an intermediate rate of adherence (50%). Those with other psychiatric disorders and no psychiatric disorders also had high rates of adherence (75.6% and 72.9%, respectively). Living with family significantly increased adherence for children with substance abuse disorders but decreased adherence for those with no psychiatric disorders. The diagnostic correlates of duration of participation were similar to those for adherence. CONCLUSIONS: Mood and substance abuse disorders were strong predictive factors for treatment adherence and duration, albeit in opposite directions. Living with family seems to have a positive effect on treatment adherence for patients with substance abuse disorders. More effective treatment is needed for victimized substance-abusing youth.

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.

The impact of psychiatric diagnosis on treatment adherence and duration among victimized children and adolescents in São Paulo, Brazil

OBJECTIVE: Despite the high prevalence of substance abuse and mood disorders among victimized children and adolescents, few studies have investigated the association of these disorders with treatment adherence, represented by numbers of visits per month and treatment duration. We aimed to investigate the effects of substance abuse and mood disorders on treatment adherence and duration in a special program for victimized children in São Paulo, Brazil. METHODS: A total of 351 participants were evaluated for psychiatric disorders and classified into one of five groups: mood disorders alone; substance abuse disorders alone; mood and substance abuse disorders; other psychiatric disorders; no psychiatric disorders. The associations between diagnostic classification and adherence to treatment and the duration of program participation were tested with logistic regression and survival analysis, respectively. RESULTS: Children with mood disorders alone had the highest rate of adherence (79.5%); those with substance abuse disorders alone had the lowest (40%); and those with both disorders had an intermediate rate of adherence (50%). Those with other psychiatric disorders and no psychiatric disorders also had high rates of adherence (75.6% and 72.9%, respectively). Living with family significantly increased adherence for children with substance abuse disorders but decreased adherence for those with no psychiatric disorders. The diagnostic correlates of duration of participation were similar to those for adherence. CONCLUSIONS: Mood and substance abuse disorders were strong predictive factors for treatment adherence and duration, albeit in opposite directions. Living with family seems to have a positive effect on treatment adherence for patients with substance abuse disorders. More effective treatment is needed for victimized substance-abusing youth

The Integrated European Project "GEHA - GEnetics of Healthy Aging": recruitment, health status assessment and survival of the Italian 90+ sibpairs

The present study is part of the EU Integrated Project “GEHA – Genetics of Healthy Aging” (Franceschi C et al., Ann N Y Acad Sci. 1100: 21-45, 2007), whose aim is to identify genes involved in healthy aging and longevity, which allow individuals to survive to advanced age in good cognitive and physical function and in the absence of major age-related diseases. Aims The major aims of this thesis were the following: 1. to outline the recruitment procedure of 90+ Italian siblings performed by the recruiting units of the University of Bologna (UNIBO) and Rome (ISS). The procedures related to the following items necessary to perform the study were described and commented: identification of the eligible area for recruitment, demographic aspects related to the need of getting census lists of 90+siblings, mail and phone contact with 90+ subjects and their families, bioethics aspects of the whole procedure, standardization of the recruitment methodology and set-up of a detailed flow chart to be followed by the European recruitment centres (obtainment of the informed consent form, anonimization of data by using a special code, how to perform the interview, how to collect the blood, how to enter data in the GEHA Phenotypic Data Base hosted at Odense). 2. to provide an overview of the phenotypic characteristics of 90+ Italian siblings recruited by the recruiting units of the University of Bologna (UNIBO) and Rome (ISS). The following items were addressed: socio-demographic characteristics, health status, cognitive assessment, physical conditions (handgrip strength test, chair-stand test, physical ability including ADL, vision and hearing ability, movement ability and doing light housework), life-style information (smoking and drinking habits) and subjective well-being (attitude towards life). Moreover, haematological parameters collected in the 90+ sibpairs as optional parameters by the Bologna and Rome recruiting units were used for a more comprehensive evaluation of the results obtained using the above mentioned phenotypic characteristics reported in the GEHA questionnaire. 3. to assess 90+ Italian siblings as far as their health/functional status is concerned on the basis of three classification methods proposed in previous studies on centenarians, which are based on: • actual functional capabilities (ADL, SMMSE, visual and hearing abilities) (Gondo et al., J Gerontol. 61A (3): 305-310, 2006); • actual functional capabilities and morbidity (ADL, ability to walk, SMMSE, presence of cancer, ictus, renal failure, anaemia, and liver diseases) (Franceschi et al., Aging Clin Exp Res, 12:77-84, 2000); • retrospectively collected data about past history of morbidity and age of disease onset (hypertension, heart disease, diabetes, stroke, cancer, osteopororis, neurological diseases, chronic obstructive pulmonary disease and ocular diseases) (Evert et al., J Gerontol A Biol Sci Med Sci. 58A (3): 232-237, 2003). Firstly these available models to define the health status of long-living subjects were applied to the sample and, since the classifications by Gondo and Franceschi are both based on the present functional status, they were compared in order to better recognize the healthy aging phenotype and to identify the best group of 90+ subjects out of the entire studied population. 4. to investigate the concordance of health and functional status among 90+ siblings in order to divide sibpairs in three categories: the best (both sibs are in good shape), the worst (both sibs are in bad shape) and an intermediate group (one sib is in good shape and the other is in bad shape). Moreover, the evaluation wanted to discover which variables are concordant among siblings; thus, concordant variables could be considered as familiar variables (determined by the environment or by genetics). 5. to perform a survival analysis by using mortality data at 1st January 2009 from the follow-up as the main outcome and selected functional and clinical parameters as explanatory variables. Methods A total of 765 90+ Italian subjects recruited by UNIBO (549 90+ siblings, belonging to 258 families) and ISS (216 90+ siblings, belonging to 106 families) recruiting units are included in the analysis. Each subject was interviewed according to a standardized questionnaire, comprising extensively utilized questions that have been validated in previous European studies on elderly subjects and covering demographic information, life style, living conditions, cognitive status (SMMSE), mood, health status and anthropometric measurements. Moreover, subjects were asked to perform some physical tests (Hand Grip Strength test and Chair Standing test) and a sample of about 24 mL of blood was collected and then processed according to a common protocol for the preparation and storage of DNA aliquots. Results From the analysis the main findings are the following: - a standardized protocol to assess cognitive status, physical performances and health status of European nonagenarian subjects was set up, in respect to ethical requirements, and it is available as a reference for other studies in this field; - GEHA families are enriched in long-living members and extreme survival, and represent an appropriate model for the identification of genes involved in healthy aging and longevity; - two simplified sets of criteria to classify 90+ sibling according to their health status were proposed, as operational tools for distinguishing healthy from non healthy subjects; - cognitive and functional parameters have a major role in categorizing 90+ siblings for the health status; - parameters such as education and good physical abilities (500 metres walking ability, going up and down the stairs ability, high scores at hand grip and chair stand tests) are associated with a good health status (defined as “cognitive unimpairment and absence of disability”); - male nonagenarians show a more homogeneous phenotype than females, and, though far fewer in number, tend to be healthier than females; - in males the good health status is not protective for survival, confirming the male-female health survival paradox; - survival after age 90 was dependent mainly on intact cognitive status and absence of functional disabilities; - haemoglobin and creatinine levels are both associated with longevity; - the most concordant items among 90+ siblings are related to the functional status, indicating that they contain a familiar component. It is still to be investigated at what level this familiar component is determined by genetics or by environment or by the interaction between genetics, environment and chance (and at what level). Conclusions In conclusion, we could state that this study, in accordance with the main objectives of the whole GEHA project, represents one of the first attempt to identify the biological and non biological determinants of successful/unsuccessful aging and longevity. Here, the analysis was performed on 90+ siblings recruited in Northern and Central Italy and it could be used as a reference for others studies in this field on Italian population. Moreover, it contributed to the definition of “successful” and “unsuccessful” aging and categorising a very large cohort of our most elderly subjects into “successful” and “unsuccessful” groups provided an unrivalled opportunity to detect some of the basic genetic/molecular mechanisms which underpin good health as opposed to chronic disability. Discoveries in the topic of the biological determinants of healthy aging represent a real possibility to identify new markers to be utilized for the identification of subgroups of old European citizens having a higher risk to develop age-related diseases and disabilities and to direct major preventive medicine strategies for the new epidemic of chronic disease in the 21st century.

Evaluating Prediction Rules for t-Year Survivors With Censored Regression Models

Suppose that we are interested in establishing simple, but reliable rules for predicting future t-year survivors via censored regression models. In this article, we present inference procedures for evaluating such binary classification rules based on various prediction precision measures quantified by the overall misclassification rate, sensitivity and specificity, and positive and negative predictive values. Specifically, under various working models we derive consistent estimators for the above measures via substitution and cross validation estimation procedures. Furthermore, we provide large sample approximations to the distributions of these nonsmooth estimators without assuming that the working model is correctly specified. Confidence intervals, for example, for the difference of the precision measures between two competing rules can then be constructed. All the proposals are illustrated with two real examples and their finite sample properties are evaluated via a simulation study.

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.