Development and validation of a risk model for prediction of hazardous alcohol consumption in general practice attendees : the PredictAL study

Background: Little is known about the risk of progression to hazardous alcohol use in people currently drinking at safe limits. We aimed to develop a prediction model (predictAL) for the development of hazardous drinking in safe drinkers. Methods: A prospective cohort study of adult general practice attendees in six European countries and Chile followed up over 6 months. We recruited 10,045 attendees between April 2003 to February 2005. 6193 European and 2462 Chilean attendees recorded AUDIT scores below 8 in men and 5 in women at recruitment and were used in modelling risk. 38 risk factors were measured to construct a risk model for the development of hazardous drinking using stepwise logistic regression. The model was corrected for over fitting and tested in an external population. The main outcome was hazardous drinking defined by an AUDIT score >= 8 in men and >= 5 in women. Results: 69.0% of attendees were recruited, of whom 89.5% participated again after six months. The risk factors in the final predictAL model were sex, age, country, baseline AUDIT score, panic syndrome and lifetime alcohol problem. The predictAL model's average c-index across all six European countries was 0.839 (95% CI 0.805, 0.873). The Hedge's g effect size for the difference in log odds of predicted probability between safe drinkers in Europe who subsequently developed hazardous alcohol use and those who did not was 1.38 (95% CI 1.25, 1.51). External validation of the algorithm in Chilean safe drinkers resulted in a c-index of 0.781 (95% CI 0.717, 0.846) and Hedge's g of 0.68 (95% CI 0.57, 0.78). Conclusions: The predictAL risk model for development of hazardous consumption in safe drinkers compares favourably with risk algorithms for disorders in other medical settings and can be a useful first step in prevention of alcohol misuse.

Modeling in-Hospital Patient Survival During the First 28 Days After Intensive Care Unit Admission: a Prognostic Model for Clinical Trials in General Critically Ill Patients

OBJECTIVE: The objective of the study was to develop a model for estimating patient 28-day in-hospital mortality using 2 different statistical approaches. DESIGN: The study was designed to develop an outcome prediction model for 28-day in-hospital mortality using (a) logistic regression with random effects and (b) a multilevel Cox proportional hazards model. SETTING: The study involved 305 intensive care units (ICUs) from the basic Simplified Acute Physiology Score (SAPS) 3 cohort. PATIENTS AND PARTICIPANTS: Patients (n = 17138) were from the SAPS 3 database with follow-up data pertaining to the first 28 days in hospital after ICU admission. INTERVENTIONS: None. MEASUREMENTS AND RESULTS: The database was divided randomly into 5 roughly equal-sized parts (at the ICU level). It was thus possible to run the model-building procedure 5 times, each time taking four fifths of the sample as a development set and the remaining fifth as the validation set. At 28 days after ICU admission, 19.98% of the patients were still in the hospital. Because of the different sampling space and outcome variables, both models presented a better fit in this sample than did the SAPS 3 admission score calibrated to vital status at hospital discharge, both on the general population and in major subgroups. CONCLUSIONS: Both statistical methods can be used to model the 28-day in-hospital mortality better than the SAPS 3 admission model. However, because the logistic regression approach is specifically designed to forecast 28-day mortality, and given the high uncertainty associated with the assumption of the proportionality of risks in the Cox model, the logistic regression approach proved to be superior.

DISTINCT CELLULAR MIGRATION INDUCED BY Leishmania infantum chagasi AND SALIVA FROM Lutzomyia longipalpis IN A HEMORRHAGIC POOL MODEL

Recruitment of a specific cell population after Leishmania infection can influence the outcome of the disease. Cellular migration in response to Leishmania or vector saliva has been reported in air pouch model, however, cellular migration induced by Leishmania associated with host's blood and vector saliva in this model has not been described. Herein we investigated cellular migration into air pouch of hamster after stimulation with combination of L. chagasi and host's blood and Lutzomyia longipalpis saliva. Migration induced by saliva was 3-fold more than those induced by L. chagasi alone. Additionally, L. chagasi associated with blood and saliva induced significantly even more leukocytes into air pouch than Leishmania alone. L. chagasi recruited a diverse cell population; however, most of these cells seem to have not migrated to the inflammatory exudate, remaining in the pouch lining tissue. These results indicate that L. chagasi can reduce leukocyte accumulation to the initial site of infection, and when associated with vector saliva in the presence of blood components, increase the influx of more neutrophils than macrophages, suggesting that the parasite has developed a strategy to minimize the initial inflammatory response, allowing an unlimited progression within the host. This work reinforces the importance of studies on the salivary components of sand fly vectors of leishmaniasis in the transmission process and the establishment of the infection.

Acute Kidney Injury after Pediatric Cardiac Surgery: Risk Factors and Outcomes. Proposal for a Predictive Model

Objectives: To characterize the epidemiology and risk factors for acute kidney injury (AKI) after pediatric cardiac surgery in our center, to determine its association with poor short-term outcomes, and to develop a logistic regression model that will predict the risk of AKI for the study population. Methods: This single-center, retrospective study included consecutive pediatric patients with congenital heart disease who underwent cardiac surgery between January 2010 and December 2012. Exclusion criteria were a history of renal disease, dialysis or renal transplantation. Results: Of the 325 patients included, median age three years (1 day---18 years), AKI occurred in 40 (12.3%) on the first postoperative day. Overall mortality was 13 (4%), nine of whom were in the AKI group. AKI was significantly associated with length of intensive care unit stay, length of mechanical ventilation and in-hospital death (p<0.01). Patients’ age and postoperative serum creatinine, blood urea nitrogen and lactate levels were included in the logistic regression model as predictor variables. The model accurately predicted AKI in this population, with a maximum combined sensitivity of 82.1% and specificity of 75.4%. Conclusions: AKI is common and is associated with poor short-term outcomes in this setting. Younger age and higher postoperative serum creatinine, blood urea nitrogen and lactate levels were powerful predictors of renal injury in this population. The proposed model could be a useful tool for risk stratification of these patients.

Validation of Two US Risk Scores for Percutaneous Coronary Intervention in a Single-Center Portuguese Population of Patients with Acute Coronary Syndrome

INTRODUCTION: New scores have been developed and validated in the US for in-hospital mortality risk stratification in patients undergoing coronary angioplasty: the National Cardiovascular Data Registry (NCDR) risk score and the Mayo Clinic Risk Score (MCRS). We sought to validate these scores in a European population with acute coronary syndrome (ACS) and to compare their predictive accuracy with that of the GRACE risk score. METHODS: In a single-center ACS registry of patients undergoing coronary angioplasty, we used the area under the receiver operating characteristic curve (AUC), a graphical representation of observed vs. expected mortality, and net reclassification improvement (NRI)/integrated discrimination improvement (IDI) analysis to compare the scores. RESULTS: A total of 2148 consecutive patients were included, mean age 63 years (SD 13), 74% male and 71% with ST-segment elevation ACS. In-hospital mortality was 4.5%. The GRACE score showed the best AUC (0.94, 95% CI 0.91-0.96) compared with NCDR (0.87, 95% CI 0.83-0.91, p=0.0003) and MCRS (0.85, 95% CI 0.81-0.90, p=0.0003). In model calibration analysis, GRACE showed the best predictive power. With GRACE, patients were more often correctly classified than with MCRS (NRI 78.7, 95% CI 59.6-97.7; IDI 0.136, 95% CI 0.073-0.199) or NCDR (NRI 79.2, 95% CI 60.2-98.2; IDI 0.148, 95% CI 0.087-0.209). CONCLUSION: The NCDR and Mayo Clinic risk scores are useful for risk stratification of in-hospital mortality in a European population of patients with ACS undergoing coronary angioplasty. However, the GRACE score is still to be preferred.

A SARIMA forecasting model to predict the number of cases of dengue in Campinas, State of São Paulo, Brazil

INTRODUCTION: Forecasting dengue cases in a population by using time-series models can provide useful information that can be used to facilitate the planning of public health interventions. The objective of this article was to develop a forecasting model for dengue incidence in Campinas, southeast Brazil, considering the Box-Jenkins modeling approach. METHODS: The forecasting model for dengue incidence was performed with R software using the seasonal autoregressive integrated moving average (SARIMA) model. We fitted a model based on the reported monthly incidence of dengue from 1998 to 2008, and we validated the model using the data collected between January and December of 2009. RESULTS: SARIMA (2,1,2) (1,1,1)12 was the model with the best fit for data. This model indicated that the number of dengue cases in a given month can be estimated by the number of dengue cases occurring one, two and twelve months prior. The predicted values for 2009 are relatively close to the observed values. CONCLUSIONS: The results of this article indicate that SARIMA models are useful tools for monitoring dengue incidence. We also observe that the SARIMA model is capable of representing with relative precision the number of cases in a next year.

Harvest control rule vs optimal harvesting of an age-structured population: The case of the Ibero-Atlantic Sardine Fishery

Research Masters

Population Dynamics Modeling of Arapaima gigas

Pirarucu (Arapaima gigas) has been of the most important natural fishing resources of the Amazon region. Due to its economic importance, and the necessity to preserve the species hand, field research concerning the habits and behavior of the pirarucu has been increasing for the last 20 years. The aim of this paper is to present a mathematical model for the pirarucu population dynamics considering the species peculiarities, particularly the male parental care over the offspring. The solution of the dynamical systems indicates three possible equilibrium points for the population. The first corresponds to extinction; the third corresponds to a stable population close to the environmental carrying capacity. The second corresponds to an unstable equilibrium located between extinction and full use of the carrying capacity. It is shown that lack of males’ parental care closes the gap between the point corresponding to the unstable equilibrium and the point of stable non-trivial equilibrium. If guarding failure reaches a critical point the two points coincide and the population tends irreversibly to extinction. If some event tends to destabilize the population equilibrium, as for instance inadequate parental care, the model responds in such a way as to restore the trajectory towards the stable equilibrium point avoiding the route to extinction. The parameters introduced to solve the system of equations are partially derived from limited but reliable field data collected at the Mamirauá Sustainable Development Reserve (MSDR) in the Brazilian Amazonian Region.

Risk factors for atherosclerosis in an elderly out patient population in the City of São Paulo

OBJECTIVE: To analyze in out clinic elderly patients of both sexes for the prevalence of risk factors for atherosclerosis and study their association with the complications of atherosclerosis. METHODS: Five hundred and sixteen outpatients, 152 men and 364 women, 60 years or older, were studied. The prevalences of hypertension, dyslipidemia, diabetes mellitus, cigarette smoking and obesity were determined in both sexes and compared using the chi-square test. The association between these factors and the presence of atherosclerotic complications was analyzed by logistic regression. RESULTS: The comparative analysis of the factors in both sexes showed that hypertension, total cholesterol > or = 240mg/dL, LDL-cholesterol > or = 160mg/dL, and body mass index >27.5 were more frequent among women, but HDL-cholesterol <35mg/dL and cigarette smoking were more frequent among men, and no difference occurred between sexes in relation to the frequency of triglycerides > or = 250mg/dL and diabetes mellitus. After adjustment of the variables in the regression model, we observed that in the total of elderly patients, risk factors for complications of atherosclerosis were: triglycerides > or = 250mg/dL, hypertension, and male sex. Among men, the risk factors were: LDL-cholesterol > or = 160mg/dL, diabetes mellitus, HDL-cholesterol <35mg/dL and hypertension. Among women, the risk factors were: tryglicerides > or = 250mg/dL and hypertension. CONCLUSION: The results showed that, in the elderly, the risk factors for atherosclerosis persist, but with different behaviors between men and women. The study suggests that the relative importance of the risk factors can change with the aging process.

Estimation of statistical cure from cancer using population-based data

Dissertação de mestrado em Estatística

Some predictors of cardiovascular mortality among the elderly population of Botucatu (SP)

OBJECTIVE: To detect factors associated with cardiovascular mortality in the elderly of Botucatu. METHODS: We evaluated 29 variables of interest in a cohort of patients aged ³60 using data from a survey conducted between 1983/84. The elderly cohort was analyzed in 1992 to detect the occurrence of cardiovascular deaths. Survival analysis was performed using the Kaplan-Meier method, the log-rank test, and Cox regression analysis. Three models were adapted for each group of variables, and a final model was chosen from those variables selected from each group. RESULTS: We identified predictor for cardiovascular death according to age for elderly males not supporting the family, not possessing a vehicle, and previous cardiovascular disease. In elderly females, the predictor variables were previous cardiovascular disease and diabetes mellitus. CONCLUSION: Socioeconomic indicators (family heading and vehicle ownerrship) may be added to well stabilished medical factors (diabete mellitus and hypertension to select target groups for programs intended to reduce deaths due to cardiovascular diseases in elderly people.

Validation study of an automated wrist monitor, omron model HEM-608, compared with the standard methods for blood pressure measurement

OBJECTIVE - The aim of our study was to assess the profile of a wrist monitor, the Omron Model HEM-608, compared with the indirect method for blood pressure measurement. METHODS - Our study population consisted of 100 subjects, 29 being normotensive and 71 being hypertensive. Participants had their blood pressure checked 8 times with alternate techniques, 4 by the indirect method and 4 with the Omron wrist monitor. The validation criteria used to test this device were based on the internationally recognized protocols. RESULTS - Our data showed that the Omron HEM-608 reached a classification B for systolic and A for diastolic blood pressure, according to the one protocol. The mean differences between blood pressure values obtained with each of the methods were -2.3 +7.9mmHg for systolic and 0.97+5.5mmHg for diastolic blood pressure. Therefore, we considered this type of device approved according to the criteria selected. CONCLUSION - Our study leads us to conclude that this wrist monitor is not only easy to use, but also produces results very similar to those obtained by the standard indirect method.

An Equation for the Prediction of Oxygen Consumption in a Brazilian Population

Background: The equations predicting maximal oxygen uptake (VO2max or peak) presently in use in cardiopulmonary exercise testing (CPET) softwares in Brazil have not been adequately validated. These equations are very important for the diagnostic capacity of this method. Objective: Build and validate a Brazilian Equation (BE) for prediction of VO2peak in comparison to the equation cited by Jones (JE) and the Wasserman algorithm (WA). Methods: Treadmill evaluation was performed on 3119 individuals with CPET (breath by breath). The construction group (CG) of the equation consisted of 2495 healthy participants. The other 624 individuals were allocated to the external validation group (EVG). At the BE (derived from a multivariate regression model), age, gender, body mass index (BMI) and physical activity level were considered. The same equation was also tested in the EVG. Dispersion graphs and Bland-Altman analyses were built. Results: In the CG, the mean age was 42.6 years, 51.5% were male, the average BMI was 27.2, and the physical activity distribution level was: 51.3% sedentary, 44.4% active and 4.3% athletes. An optimal correlation between the BE and the CPET measured VO2peak was observed (0.807). On the other hand, difference came up between the average VO2peak expected by the JE and WA and the CPET measured VO2peak, as well as the one gotten from the BE (p = 0.001). Conclusion: BE presents VO2peak values close to those directly measured by CPET, while Jones and Wasserman differ significantly from the real VO2peak.

Experimental design for mixed models with application to population pharmacokinetic studies

experimental design, mixed model, random coefficient regression model, population pharmacokinetics, approximate design

Modeling and simulation of population balances for particulate processes

This work focuses on the modeling and numerical approximations of population balance equations (PBEs) for the simulation of different phenomena occurring in process engineering. The population balance equation (PBE) is considered to be a statement of continuity. It tracks the change in particle size distribution as particles are born, die, grow or leave a given control volume. In the population balance models the one independent variable represents the time, the other(s) are property coordinate(s), e.g., the particle volume (size) in the present case. They typically describe the temporal evolution of the number density functions and have been used to model various processes such as granulation, crystallization, polymerization, emulsion and cell dynamics. The semi-discrete high resolution schemes are proposed for solving PBEs modeling one and two-dimensional batch crystallization models. The schemes are discrete in property coordinates but continuous in time. The resulting ordinary differential equations can be solved by any standard ODE solver. To improve the numerical accuracy of the schemes a moving mesh technique is introduced in both one and two-dimensional cases ...