Residuals for log-Burr XII regression models in survival analysis

In this paper, we compare three residuals to assess departures from the error assumptions as well as to detect outlying observations in log-Burr XII regression models with censored observations. These residuals can also be used for the log-logistic regression model, which is a special case of the log-Burr XII regression model. For different parameter settings, sample sizes and censoring percentages, various simulation studies are performed and the empirical distribution of each residual is displayed and compared with the standard normal distribution. These studies suggest that the residual analysis usually performed in normal linear regression models can be straightforwardly extended to the modified martingale-type residual in log-Burr XII regression models with censored data.

A Brazilian registry of juvenile dermatomyositis: onset features and classification of 189 cases

Objective To describe onset features, classification and treatment of juvenile dermatomyositis (JDM) and juvenile polymyositis (JPM) from a multicentre registry. Methods Inclusion criteria were onset age lower than 18 years and a diagnosis of any idiopathic inflammatory myopathy (IIM) by attending physician. Bohan & Peter (1975) criteria categorisation was established by a scoring algorithm to define JDM and JPM based oil clinical protocol data. Results Of the 189 cases included, 178 were classified as JDM, 9 as JPM (19.8: 1) and 2 did not fit the criteria; 6.9% had features of chronic arthritis and connective tissue disease overlap. Diagnosis classification agreement occurred in 66.1%. Medial? onset age was 7 years, median follow-up duration was 3.6 years. Malignancy was described in 2 (1.1%) cases. Muscle weakness occurred in 95.8%; heliotrope rash 83.5%; Gottron plaques 83.1%; 92% had at least one abnormal muscle enzyme result. Muscle biopsy performed in 74.6% was abnormal in 91.5% and electromyogram performed in 39.2% resulted abnormal in 93.2%. Logistic regression analysis was done in 66 cases with all parameters assessed and only aldolase resulted significant, as independent variable for definite JDM (OR=5.4, 95%CI 1.2-24.4, p=0.03). Regarding treatment, 97.9% received steroids; 72% had in addition at least one: methotrexate (75.7%), hydroxychloroquine (64.7%), cyclosporine A (20.6%), IV immunoglobulin (20.6%), azathioprine (10.3%) or cyclophosphamide (9.6%). In this series 24.3% developed calcinosis and mortality rate was 4.2%. Conclusion Evaluation of predefined criteria set for a valid diagnosis indicated aldolase as the most important parameter associated with de, methotrexate combination, was the most indicated treatment.

A log-extended Weibull regression model

A bathtub-shaped failure rate function is very useful in survival analysis and reliability studies. The well-known lifetime distributions do not have this property. For the first time, we propose a location-scale regression model based on the logarithm of an extended Weibull distribution which has the ability to deal with bathtub-shaped failure rate functions. We use the method of maximum likelihood to estimate the model parameters and some inferential procedures are presented. We reanalyze a real data set under the new model and the log-modified Weibull regression model. We perform a model check based on martingale-type residuals and generated envelopes and the statistics AIC and BIC to select appropriate models. (C) 2009 Elsevier B.V. All rights reserved.

Different Approaches for Modeling Grouped Survival Data: A Mango Tree Study

Interval-censored survival data, in which the event of interest is not observed exactly but is only known to occur within some time interval, occur very frequently. In some situations, event times might be censored into different, possibly overlapping intervals of variable widths; however, in other situations, information is available for all units at the same observed visit time. In the latter cases, interval-censored data are termed grouped survival data. Here we present alternative approaches for analyzing interval-censored data. We illustrate these techniques using a survival data set involving mango tree lifetimes. This study is an example of grouped survival data.

Log-modified Weibull regression models with censored data: Sensitivity and residual analysis

This paper proposes a regression model considering the modified Weibull distribution. This distribution can be used to model bathtub-shaped failure rate functions. Assuming censored data, we consider maximum likelihood and Jackknife estimators for the parameters of the model. We derive the appropriate matrices for assessing local influence on the parameter estimates under different perturbation schemes and we also present some ways to perform global influence. Besides, for different parameter settings, sample sizes and censoring percentages, various simulations are performed and the empirical distribution of the modified deviance residual is displayed and compared with the standard normal distribution. These studies suggest that the residual analysis usually performed in normal linear regression models can be straightforwardly extended for a martingale-type residual in log-modified Weibull regression models with censored data. Finally, we analyze a real data set under log-modified Weibull regression models. A diagnostic analysis and a model checking based on the modified deviance residual are performed to select appropriate models. (c) 2008 Elsevier B.V. All rights reserved.

On estimation and influence diagnostics for zero-inflated negative binomial regression models

The zero-inflated negative binomial model is used to account for overdispersion detected in data that are initially analyzed under the zero-Inflated Poisson model A frequentist analysis a jackknife estimator and a non-parametric bootstrap for parameter estimation of zero-inflated negative binomial regression models are considered In addition an EM-type algorithm is developed for performing maximum likelihood estimation Then the appropriate matrices for assessing local influence on the parameter estimates under different perturbation schemes and some ways to perform global influence analysis are derived In order to study departures from the error assumption as well as the presence of outliers residual analysis based on the standardized Pearson residuals is discussed The relevance of the approach is illustrated with a real data set where It is shown that zero-inflated negative binomial regression models seems to fit the data better than the Poisson counterpart (C) 2010 Elsevier B V All rights reserved

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.

Does using stepwise variable selection to build sequential path analysis models make sense?

Causal inference methods - mainly path analysis and structural equation modeling - offer plant physiologists information about cause-and-effect relationships among plant traits. Recently, an unusual approach to causal inference through stepwise variable selection has been proposed and used in various works on plant physiology. The approach should not be considered correct from a biological point of view. Here, it is explained why stepwise variable selection should not be used for causal inference, and shown what strange conclusions can be drawn based upon the former analysis when one aims to interpret cause-and-effect relationships among plant traits.

Contribution of soil organic carbon to the ion exchange capacity of tropical soils

Highly weathered soils represent about 3 billion ha of the tropical region. Oxisols represent about 60% of the Brazilian territory (more than 5 million km 2), in areas of great agricultural importance. Soil organic carbon (SOC) can be responsible for more than 80% of the cation exchange capacity (CEC) of highly weathered soils, such as Oxisols and Ultisols. The objective of this study was to estimate the contribution of the SOC to the CEC of Brazilian soils from different orders. Surface samples (0.0 to 0.2 m) of 30 uncultivated soils (13 Oxisols, 6 Ultisols, 5 Alfisols, 3 Entisols, I Histosol, 1 Inceptisol. and I Molisol), under native forests and from reforestation sites from Sao Paulo State, Brazil, were collected in order to obtain a large variation of (electro)chemical, physical, and mineralogical soil attributes. Total content of SOC was quantified by titulometric and colorimetric methods. Effective cation exchange capacity (ECEC) was obtained by two methods: the indirect method-summation-estimated the ECECi from the sum of basic cations (Ca+ Mg+ K+ Na) and exchangeable Al; and the direct ECECd obtained by the compulsive exchange method, using unbuffered BaCl2 solution. The contribution of SOC to the soil CEC was estimated by the Bennema statistical method. The amount of SOC var ied from 6.6 g kg(-1) to 213.4 g kg(-1). while clay contents varied from 40 g kg(-1) to 716 g kg(-1). Soil organic carbon contents were strongly associated to the clay contents, suggesting that clay content was the primary variable in controling the variability of SOC contents in the samples. Cation exchange capacity varied from 7.0 mmol(c) kg(-1) to 137.8 mmol(c) kg(-1) and had a positive Correlation with SOC. The mean contribution (per grain) of the SOC (1.64 mmol(c)) for the soil CEC was more than 44 times higher than the contribution of the clay fraction (0.04 mmol(c),). A regression model that considered the SOC content as the only significant variable explained 60% of the variation in the soil total CEC. The importance of SOC was related to soil pedogenetic process, since its contribution to the soil CEC was more evident in Oxisols with predominance of Fe and Al (oxihydr)oxides in the mineral fraction or in Ultisols, that presented illuviated clay. The influence of SOC in the sign and in the magnitude of the net charge of soils reinforce the importance of agricultural management systems that preserve high levels of SOC, in order to improve their sustainability.

Parametric correlation functions to model the structure of permanent environmental (co)variances in milk yield random regression models

The objective of the present study was to estimate milk yield genetic parameters applying random regression models and parametric correlation functions combined with a variance function to model animal permanent environmental effects. A total of 152,145 test-day milk yields from 7,317 first lactations of Holstein cows belonging to herds located in the southeastern region of Brazil were analyzed. Test-day milk yields were divided into 44 weekly classes of days in milk. Contemporary groups were defined by herd-test-day comprising a total of 2,539 classes. The model included direct additive genetic, permanent environmental, and residual random effects. The following fixed effects were considered: contemporary group, age of cow at calving (linear and quadratic regressions), and the population average lactation curve modeled by fourth-order orthogonal Legendre polynomial. Additive genetic effects were modeled by random regression on orthogonal Legendre polynomials of days in milk, whereas permanent environmental effects were estimated using a stationary or nonstationary parametric correlation function combined with a variance function of different orders. The structure of residual variances was modeled using a step function containing 6 variance classes. The genetic parameter estimates obtained with the model using a stationary correlation function associated with a variance function to model permanent environmental effects were similar to those obtained with models employing orthogonal Legendre polynomials for the same effect. A model using a sixth-order polynomial for additive effects and a stationary parametric correlation function associated with a seventh-order variance function to model permanent environmental effects would be sufficient for data fitting.

Random regression models to estimate test-day milk yield genetic parameters Holstein cows in Southeastern Brazil

A total of 152,145 weekly test-day milk yield records from 7317 first lactations of Holstein cows distributed in 93 herds in southeastern Brazil were analyzed. Test-day milk yields were classified into 44 weekly classes of DIM. The contemporary groups were defined as herd-year-week of test-day. The model included direct additive genetic, permanent environmental and residual effects as random and fixed effects of contemporary group and age of cow at calving as covariable, linear and quadratic effects. Mean trends were modeled by a cubic regression on orthogonal polynomials of DIM. Additive genetic and permanent environmental random effects were estimated by random regression on orthogonal Legendre polynomials. Residual variances were modeled using third to seventh-order variance functions or a step function with 1, 6,13,17 and 44 variance classes. Results from Akaike`s and Schwarz`s Bayesian information criterion suggested that a model considering a 7th-order Legendre polynomial for additive effect, a 12th-order polynomial for permanent environment effect and a step function with 6 classes for residual variances, fitted best. However, a parsimonious model, with a 6th-order Legendre polynomial for additive effects and a 7th-order polynomial for permanent environmental effects, yielded very similar genetic parameter estimates. (C) 2008 Elsevier B.V. All rights reserved.

A Log-Linear Regression Model for the Beta-Weibull Distribution

We introduce the log-beta Weibull regression model based on the beta Weibull distribution (Famoye et al., 2005; Lee et al., 2007). We derive expansions for the moment generating function which do not depend on complicated functions. The new regression model represents a parametric family of models that includes as sub-models several widely known regression models that can be applied to censored survival data. We employ a frequentist analysis, a jackknife estimator, and a parametric bootstrap for the parameters of the proposed model. We derive the appropriate matrices for assessing local influences on the parameter estimates under different perturbation schemes and present some ways to assess global influences. Further, for different parameter settings, sample sizes, and censoring percentages, several simulations are performed. In addition, the empirical distribution of some modified residuals are displayed and compared with the standard normal distribution. These studies suggest that the residual analysis usually performed in normal linear regression models can be extended to a modified deviance residual in the proposed regression model applied to censored data. We define martingale and deviance residuals to evaluate the model assumptions. The extended regression model is very useful for the analysis of real data and could give more realistic fits than other special regression models.

Test procedure for variable rate fertilizer on coffee

The objective was to develop and test a procedure for applying variable rates of fertilizers and evaluate yield response in coffee (Coffea arabica L.) with regard to the application of phosphorus and potassium. The work was conducted during the 2004 season in a 6.4 ha field located in central Sao Paulo state. Two treatments were applied with alternating strips of fixed and variable rates during the whole season: one following the fertilizing procedures recommended locally, and the other based on a grid soil sampling. A prototype pneumatic fertilizer applicator was used, carrying two conveyor belts, one for each row. Harvesting was done with a commercial harvester equipped with a customized volumetric yield monitor, separating the two treatments. Data were analyzed based on geostatistics, correlations and regressions. The procedure showed to be feasible and effective. The area that received fertilizer applications at a variable rate showed a 34% yield increase compared to the area that received a fixed rate. The variable rate fertilizer resulted in a savings of 23% in phosphate fertilizer and a 13% increase in potassium fertilizer, when compared to fixed rate fertilizer. Yield in 2005, the year after the variable rate treatments, still presented residual effect from treatments carried out during the previous cycle.

Paclitaxel associated with cholesterol-rich nanoemulsions promotes atherosclerosis regression in the rabbit

A cholesterol-rich nanoemulsion (LDE) that resembles LDL binds to the LDL receptors and after injection into the blood stream may concentrate in cells with LDL receptor overexpression, as occurs in neoplasias and other proliferative processes. Thus, LDE can be used as vehicle to target drugs against those cells. The current study was designed to verify in rabbits whether LDE concentrates in the lesioned rabbit artery and whether a paclitaxel derivative, paclitaxel oleate, associated to LDE could reduce the atherosclerotic lesions. Sixteen male New Zealand rabbits were fed a 1% cholesterol diet for 60 days. Starting from day 30 under cholesterol feeding, eight animals were treated with four weekly intravenous injections of LDE-paclitaxel (4 mg/kg) and eight with four weekly intravenous saline solution injections for additional 30 days. On day 60, the animals were sacrificed for analysis. The uptake of LDE labeled with [C-14]-cholesteryl oleate by the aortic arch of cholesterol-fed rabbits was twice as much that observed in animals fed only regular chow. LDE-paclitaxel reduced the lesion areas of cholesterol-fed animals by 60% and intima-media ratio fourfold and inhibited the macrophage migration and the smooth muscle cell proliferation and invasion of the intima. LDE-paclitaxel treatment had no toxicity. In conclusion, LDE-paclitaxel produced pronounced atherosclerosis regression without toxicity and has shown remarkable potential in cardiovascular therapeutics. (c) 2008 Published by Elsevier Ireland Ltd.

Stability of clavulanic acid under variable pH, ionic strength and temperature conditions. A new kinetic approach

Clavulanic acid (CA) is a beta-lactam antibiotic that alone exhibits only weak antibacterial activity, but is a potent inhibitor of beta-lactamases enzymes. For this reason it is used as a therapeutic in conjunction with penicillins and cephalosporins. However, it is a well-known fact that it is unstable not only during its production phase, but also during downstream processing. Therefore, the main objective of this study was the evaluation of CA long-term stability under different conditions of pH and temperature, in the presence of variable levels of different salts, so as to suggest the best conditions to perform its simultaneous production and recovery by two-phase polymer/salt liquid-liquid extractive fermentation. To this purpose, the CA stability was investigated at different values of pH (4.0-8.0) and temperature (20-45 degrees C), and the best conditions were met at a pH 6.0-7.2 and 20 degrees C. Its stability was also investigated at 30 degrees C in the presence of NaCl, Na(2)SO(4), CaCl(2) and MgSO(4) at concentrations of 0.1 and 0.5 M in Mcllvaine buffer (pH 6.5). All salts led to increased CA instability with respect to the buffer alone, and this effect decreased in following sequence: Na(2)SO(4) > MgSO(4) > CaCl(2) > NaCl. Kinetic and thermodynamic parameters of CA degradation were calculated adopting a new model that took into consideration the equilibrium between the active and a reversibly inactivated form of CA after long-time degradation. (C) 2009 Elsevier B.V. All rights reserved.