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.

The negative binomial-beta Weibull regression model to predict the cure of prostate cancer

In this article, for the first time, we propose the negative binomial-beta Weibull (BW) regression model for studying the recurrence of prostate cancer and to predict the cure fraction for patients with clinically localized prostate cancer treated by open radical prostatectomy. The cure model considers that a fraction of the survivors are cured of the disease. The survival function for the population of patients can be modeled by a cure parametric model using the BW distribution. We derive an explicit expansion for the moments of the recurrence time distribution for the uncured individuals. The proposed distribution can be used to model survival data when the hazard rate function is increasing, decreasing, unimodal and bathtub shaped. Another advantage is that the proposed model includes as special sub-models some of the well-known cure rate models discussed in the literature. We derive the appropriate matrices for assessing local influence on the parameter estimates under different perturbation schemes. We analyze a real data set for localized prostate cancer patients after open radical prostatectomy.

A log-linear regression model for the beta-Birnbaum-Saunders distribution with censored data

The beta-Birnbaum-Saunders (Cordeiro and Lemonte, 2011) and Birnbaum-Saunders (Birnbaum and Saunders, 1969a) distributions have been used quite effectively to model failure times for materials subject to fatigue and lifetime data. We define the log-beta-Birnbaum-Saunders distribution by the logarithm of the beta-Birnbaum-Saunders distribution. Explicit expressions for its generating function and moments are derived. We propose a new log-beta-Birnbaum-Saunders regression model that can be applied to censored data and be used more effectively in survival analysis. We obtain the maximum likelihood estimates of the model parameters for censored data and investigate influence diagnostics. The new location-scale regression model is modified for the possibility that long-term survivors may be presented in the data. Its usefulness is illustrated by means of two real data sets. (C) 2011 Elsevier B.V. All rights reserved.

Competing regression models for longitudinal data

The choice of an appropriate family of linear models for the analysis of longitudinal data is often a matter of concern for practitioners. To attenuate such difficulties, we discuss some issues that emerge when analyzing this type of data via a practical example involving pretestposttest longitudinal data. In particular, we consider log-normal linear mixed models (LNLMM), generalized linear mixed models (GLMM), and models based on generalized estimating equations (GEE). We show how some special features of the data, like a nonconstant coefficient of variation, may be handled in the three approaches and evaluate their performance with respect to the magnitude of standard errors of interpretable and comparable parameters. We also show how different diagnostic tools may be employed to identify outliers and comment on available software. We conclude by noting that the results are similar, but that GEE-based models may be preferable when the goal is to compare the marginal expected responses.

Mosquito (Diptera: Culicidae) assemblages associated with Nidularium and Vriesea bromeliads in Serra do Mar, Atlantic Forest, Brazil

Background: The most substantial and best preserved area of Atlantic Forest is within the biogeographical sub-region of Serra do Mar. The topographic complexity of the region creates a diverse array of microclimates, which can affect species distribution and diversity inside the forest. Given that Atlantic Forest includes highly heterogeneous environments, a diverse and medically important Culicidae assemblage, and possible species co-occurrence, we evaluated mosquito assemblages from bromeliad phytotelmata in Serra do Mar (southeastern Brazil). Methods: Larvae and pupae were collected monthly from Nidularium and Vriesea bromeliads between July 2008 and June 2009. Collection sites were divided into landscape categories (lowland, hillslope and hilltop) based on elevation and slope. Correlations between bromeliad mosquito assemblage and environmental variables were assessed using multivariate redundancy analysis. Differences in species diversity between bromeliads within each category of elevation were explored using the Renyi diversity index. Univariate binary logistic regression analyses were used to assess species co-occurrence. Results: A total of 2,024 mosquitoes belonging to 22 species were collected. Landscape categories (pseudo-F value = 1.89, p = 0.04), bromeliad water volume (pseudo-F = 2.99, p = 0.03) and bromeliad fullness (Pseudo-F = 4.47, p < 0.01) influenced mosquito assemblage structure. Renyi diversity index show that lowland possesses the highest diversity indices. The presence of An. homunculus was associated with Cx. ocellatus and the presence of An. cruzii was associated with Cx. neglectus, Cx. inimitabilis fuscatus and Cx. worontzowi. Anopheles cruzii and An. homunculus were taken from the same bromeliad, however, the co-occurrence between those two species was not statistically significant. Conclusions: One of the main findings of our study was that differences in species among mosquito assemblages were influenced by landscape characteristics. The bromeliad factor that influenced mosquito abundance and assemblage structure was fullness. The findings of the current study raise important questions about the role of An. homunculus in the transmission of Plasmodium in Serra do Mar, southeastern Atlantic Forest.

Accuracy of positron emission tomography/computed tomography and clinical assessment in the detection of complete rectal tumor regression after neoadjuvant chemoradiation Long-Term Results of a Prospective Trial (National Clinical Trial 00254683)

BACKGROUND: Neoadjuvant chemoradiation (CRT) therapy may result in significant tumor regression in patients with rectal cancer. Patients who develop complete tumor regression have been managed by treatment strategies that are alternatives to standard total mesorectal excision. Therefore, assessment of tumor response with positron emission tomography/computed tomography (PET/CT) after neoadjuvant treatment may offer relevant information for the selection of patients to receive alternative treatment strategies. METHODS: Patients with clinical T2 (cT2) through cT4NxM0 rectal adenocarcinoma were included prospectively. Neoadjuvant therapy consisted of 54 grays of radiation and 5-fluorouracil-based chemotherapy. Baseline PET/CT studies were obtained before CRT followed by PET/CT studies at 6 weeks and 12 weeks after the completion of CRT. Clinical assessment was performed at 12 weeks after CRT completion. PET/CT results were compared with clinical and pathologic data. RESULTS: In total, 99 patients were included in the study. Twenty-three patients were complete responders (16 had a complete clinical response, and 7 had a complete pathologic response). The PET/CT response evaluation at 12 weeks indicated that 18 patients had a complete response, and 81 patients had an incomplete response. There were 5 false-negative and 10 false-positive PET/CT results. PET/CT for the detection of residual cancer had 93% sensitivity, 53% specificity, a 73% negative predictive value, an 87% positive predictive value, and 85% accuracy. Clinical assessment alone resulted in an accuracy of 91%. PET/CT information may have detected misdiagnoses made by clinical assessment alone, improving overall accuracy to 96%. CONCLUSIONS: Assessment of tumor response at 12 weeks after CRT completion with PET/CT imaging may provide a useful additional tool with good overall accuracy for the selection of patients who may avoid unnecessary radical resection after achieving a complete clinical response. Cancer 2012;35013511. (C) 2011 American Cancer Society.

Binary and ternary inclusion complexes of dapsone in cyclodextrins and polymers: preparation, characterization and evaluation

Dapsone (DAP) is a synthetic sulfone drug with bacteriostatic activity, mainly against Mycobacterium leprae. In this study we have investigated the interactions of DAP with cyclodextrins, 2-hydroxypropyl-beta-cyclodextrin (HP beta CD) and beta-cyclodextrin (beta CD), in the presence and absence of water-soluble polymers, in order to improve its solubility and bioavailability. Solid systems DAP/HP beta CD and DAP/beta CD, in the presence or absence of polyvinylpyrrolidone (PVP K30) or hydroxypropyl methylcellulose (HPMC), were prepared. The binary and ternary systems were evaluated and characterized by SEM, DSC, XRD and NMR analysis as well as phase solubility assays, in order to investigate the interactions between DAP and the excipients in aqueous solution. This study revealed that inclusion complexes of DAP and cyclodextrins (HP beta CD and beta CD) can be produced in order to improve DAP solubility and bioavailability in the presence or absence of polymers (PVP K30 and HPMC). The more stable inclusion complex was obtained with HP beta CD, and consequently HP beta CD was more efficient in improving DAP solubility than beta CD, and the addition of polymers had no influence on DAP solubility or on the stability of the DAP/CDs complexes.

A general class of zero-or-one inflated beta regression models

This paper proposes a general class of regression models for continuous proportions when the data contain zeros or ones. The proposed class of models assumes that the response variable has a mixed continuous-discrete distribution with probability mass at zero or one. The beta distribution is used to describe the continuous component of the model, since its density has a wide range of different shapes depending on the values of the two parameters that index the distribution. We use a suitable parameterization of the beta law in terms of its mean and a precision parameter. The parameters of the mixture distribution are modeled as functions of regression parameters. We provide inference, diagnostic, and model selection tools for this class of models. A practical application that employs real data is presented. (C) 2011 Elsevier B.V. All rights reserved.

Bayesian analysis of random regression models using B-splines to model test-day milk yield of Holstein cattle in Brazil

The objective of this paper is to model variations in test-day milk yields of first lactations of Holstein cows by RR using B-spline functions and Bayesian inference in order to fit adequate and parsimonious models for the estimation of genetic parameters. They used 152,145 test day milk yield records from 7317 first lactations of Holstein cows. The model established in this study was additive, permanent environmental and residual random effects. In addition, contemporary group and linear and quadratic effects of the age of cow at calving were included as fixed effects. Authors modeled the average lactation curve of the population with a fourth-order orthogonal Legendre polynomial. They concluded that a cubic B-spline with seven random regression coefficients for both the additive genetic and permanent environment effects was to be the best according to residual mean square and residual variance estimates. Moreover they urged a lower order model (quadratic B-spline with seven random regression coefficients for both random effects) could be adopted because it yielded practically the same genetic parameter estimates with parsimony. (C) 2012 Elsevier B.V. All rights reserved.

Suicide rates and trends in Sao Paulo, Brazil, according to gender, age and demographic aspects: a joinpoint regression analysis

Objective: To evaluate suicide rates and trends in Sao Paulo by sex, age-strata, and methods. Methods: Data was collected from State registry from 1996 to 2009. Population was estimated using the National Census. We utilized joinpoint regression analysis to explore temporal trends. We also evaluated marital status, ethnicity, birthplace and methods for suicide. Results: In the period analyzed, 6,002 suicides were accrued with a rate of 4.6 per 100,000 (7.5 in men and 2.0 in women); the male-to-female ratio was around 3.7. Trends for men presented a significant decline of 5.3% per year from 1996 to 2002, and a significant increase of 2.5% from 2002 onwards. Women did not present significant changes. For men, the elderly (> 65 years) had a significant reduction of 2.3% per year, while younger men (25-44 years) presented a significant increase of 8.6% from 2004 onwards. Women did not present significant trend changes according to age. Leading suicide methods were hanging and poisoning for men and women, respectively. Other analyses showed an increased suicide risk ratio for singles and foreigners. Conclusions: Specific epidemiological trends for suicide in the city of Sao Paulo that warrant further investigation were identified. High-risk groups - such as immigrants - could benefit from targeted strategies of suicide prevention.

Random regression analyses using B-spline functions to model growth of Nellore cattle

The objective of this study was to estimate (co)variance components using random regression on B-spline functions to weight records obtained from birth to adulthood. A total of 82 064 weight records of 8145 females obtained from the data bank of the Nellore Breeding Program (PMGRN/Nellore Brazil) which started in 1987, were used. The models included direct additive and maternal genetic effects and animal and maternal permanent environmental effects as random. Contemporary group and dam age at calving (linear and quadratic effect) were included as fixed effects, and orthogonal Legendre polynomials of age (cubic regression) were considered as random covariate. The random effects were modeled using B-spline functions considering linear, quadratic and cubic polynomials for each individual segment. Residual variances were grouped in five age classes. Direct additive genetic and animal permanent environmental effects were modeled using up to seven knots (six segments). A single segment with two knots at the end points of the curve was used for the estimation of maternal genetic and maternal permanent environmental effects. A total of 15 models were studied, with the number of parameters ranging from 17 to 81. The models that used B-splines were compared with multi-trait analyses with nine weight traits and to a random regression model that used orthogonal Legendre polynomials. A model fitting quadratic B-splines, with four knots or three segments for direct additive genetic effect and animal permanent environmental effect and two knots for maternal additive genetic effect and maternal permanent environmental effect, was the most appropriate and parsimonious model to describe the covariance structure of the data. Selection for higher weight, such as at young ages, should be performed taking into account an increase in mature cow weight. Particularly, this is important in most of Nellore beef cattle production systems, where the cow herd is maintained on range conditions. There is limited modification of the growth curve of Nellore cattle with respect to the aim of selecting them for rapid growth at young ages while maintaining constant adult weight.

Prediction of enzymatic infarct size in ST-segment elevation myocardial infarction

Objectives Predictors of adverse outcomes following myocardial infarction (MI) are well established; however, little is known about what predicts enzymatically estimated infarct size in patients with acute ST-elevation MI. The Complement And Reduction of INfarct size after Angioplasty or Lytics trials of pexelizumab used creatine kinase (CK)-MB area under the curve to determine infarct size in patients treated with primary percutaneous coronary intervention (PCI) or fibrinolysis. Methods Prediction of infarct size was carried out by measuring CK-MB area under the curve in patients with ST-segment elevation MI treated with reperfusion therapy from January 2000 to April 2002. Infarct size was calculated in 1622 patients (PCI=817; fibrinolysis=805). Logistic regression was used to examine the relationship between baseline demographics, total ST-segment elevation, index angiographic findings (PCI group), and binary outcome of CK-MB area under the curve greater than 3000 ng/ml. Results Large infarcts occurred in 63% (515) of the PCI group and 69% (554) of the fibrinolysis group. Independent predictors of large infarcts differed depending on mode of reperfusion. In PCI, male sex, no prior coronary revascularization and diabetes, decreased systolic blood pressure, sum of ST-segment elevation, total (angiographic) occlusion, and nonright coronary artery culprit artery were independent predictors of larger infarcts (C index=0.73). In fibrinolysis, younger age, decreased heart rate, white race, no history of arrhythmia, increased time to fibrinolytic therapy in patients treated up to 2 h after symptom onset, and sum of ST-segment elevation were independently associated with a larger infarct size (C index=0.68). Conclusion Clinical and patient data can be used to predict larger infarcts on the basis of CK-MB quantification. These models may be helpful in designing future trials and in guiding the use of novel pharmacotherapies aimed at limiting infarct size in clinical practice. Coron Artery Dis 23:118-125 (C) 2012 Wolters Kluwer Health | Lippincott Williams & Wilkins.

Incorporating label dependency into the binary relevance framework for multi-label classification

In multi-label classification, examples can be associated with multiple labels simultaneously. The task of learning from multi-label data can be addressed by methods that transform the multi-label classification problem into several single-label classification problems. The binary relevance approach is one of these methods, where the multi-label learning task is decomposed into several independent binary classification problems, one for each label in the set of labels, and the final labels for each example are determined by aggregating the predictions from all binary classifiers. However, this approach fails to consider any dependency among the labels. Aiming to accurately predict label combinations, in this paper we propose a simple approach that enables the binary classifiers to discover existing label dependency by themselves. An experimental study using decision trees, a kernel method as well as Naive Bayes as base-learning techniques shows the potential of the proposed approach to improve the multi-label classification performance.

Predicting smear negative pulmonary tuberculosis with classification trees and logistic regression: a cross-sectional study

Abstract Background Smear negative pulmonary tuberculosis (SNPT) accounts for 30% of pulmonary tuberculosis cases reported yearly in Brazil. This study aimed to develop a prediction model for SNPT for outpatients in areas with scarce resources. Methods The study enrolled 551 patients with clinical-radiological suspicion of SNPT, in Rio de Janeiro, Brazil. The original data was divided into two equivalent samples for generation and validation of the prediction models. Symptoms, physical signs and chest X-rays were used for constructing logistic regression and classification and regression tree models. From the logistic regression, we generated a clinical and radiological prediction score. The area under the receiver operator characteristic curve, sensitivity, and specificity were used to evaluate the model's performance in both generation and validation samples. Results It was possible to generate predictive models for SNPT with sensitivity ranging from 64% to 71% and specificity ranging from 58% to 76%. Conclusion The results suggest that those models might be useful as screening tools for estimating the risk of SNPT, optimizing the utilization of more expensive tests, and avoiding costs of unnecessary anti-tuberculosis treatment. Those models might be cost-effective tools in a health care network with hierarchical distribution of scarce resources.

Automatic lumen segmentation in IVOCT images using binary morphological reconstruction

Abstract Background Atherosclerosis causes millions of deaths, annually yielding billions in expenses round the world. Intravascular Optical Coherence Tomography (IVOCT) is a medical imaging modality, which displays high resolution images of coronary cross-section. Nonetheless, quantitative information can only be obtained with segmentation; consequently, more adequate diagnostics, therapies and interventions can be provided. Since it is a relatively new modality, many different segmentation methods, available in the literature for other modalities, could be successfully applied to IVOCT images, improving accuracies and uses. Method An automatic lumen segmentation approach, based on Wavelet Transform and Mathematical Morphology, is presented. The methodology is divided into three main parts. First, the preprocessing stage attenuates and enhances undesirable and important information, respectively. Second, in the feature extraction block, wavelet is associated with an adapted version of Otsu threshold; hence, tissue information is discriminated and binarized. Finally, binary morphological reconstruction improves the binary information and constructs the binary lumen object. Results The evaluation was carried out by segmenting 290 challenging images from human and pig coronaries, and rabbit iliac arteries; the outcomes were compared with the gold standards made by experts. The resultant accuracy was obtained: True Positive (%) = 99.29 ± 2.96, False Positive (%) = 3.69 ± 2.88, False Negative (%) = 0.71 ± 2.96, Max False Positive Distance (mm) = 0.1 ± 0.07, Max False Negative Distance (mm) = 0.06 ± 0.1. Conclusions In conclusion, by segmenting a number of IVOCT images with various features, the proposed technique showed to be robust and more accurate than published studies; in addition, the method is completely automatic, providing a new tool for IVOCT segmentation.