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.

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.

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.

Exploratory Approach to the use of Techniques of Regression Analysis and Pooled Analysis in Studies of Consumer Behavior

This paper is part of a large study to assess the adequacy of the use of multivariate statistical techniques in theses and dissertations of some higher education institutions in the area of marketing with theme of consumer behavior from 1997 to 2006. The regression and conjoint analysis are focused on in this paper, two techniques with great potential of use in marketing studies. The objective of this study was to analyze whether the employement of these techniques suits the needs of the research problem presented in as well as to evaluate the level of success in meeting their premisses. Overall, the results suggest the need for more involvement of researchers in the verification of all the theoretical precepts of application of the techniques classified in the category of investigation of dependence among variables.

Multicollinearity and financial constraint in investment decisions: a Bayesian generalized ridge regression

This paper addresses the investment decisions considering the presence of financial constraints of 373 large Brazilian firms from 1997 to 2004, using panel data. A Bayesian econometric model was used considering ridge regression for multicollinearity problems among the variables in the model. Prior distributions are assumed for the parameters, classifying the model into random or fixed effects. We used a Bayesian approach to estimate the parameters, considering normal and Student t distributions for the error and assumed that the initial values for the lagged dependent variable are not fixed, but generated by a random process. The recursive predictive density criterion was used for model comparisons. Twenty models were tested and the results indicated that multicollinearity does influence the value of the estimated parameters. Controlling for capital intensity, financial constraints are found to be more important for capital-intensive firms, probably due to their lower profitability indexes, higher fixed costs and higher degree of property diversification.

Tests for the hysteresis hypothesis in Brazilian industrialized exports: A threshold cointegration analysis

This paper examines the hysteresis hypothesis in the Brazilian industrialized exports using a time series analysis. This hypothesis finds an empirical representation into the nonlinear adjustments of the exported quantity to relative price changes. Thus, the threshold cointegration analysis proposed by Balke and Fomby [Balke, N.S. and Fomby, T.B. Threshold Cointegration. International Economic Review, 1997; 38; 627-645.] was used for estimating models with asymmetric adjustment of the error correction term. Amongst sixteen industrial sectors selected, there was evidence of nonlinearities in the residuals of long-run relationships of supply or demand for exports in nine of them. These nonlinearities represent asymmetric and/or discontinuous responses of exports to different representative measures of real exchange rates, in addition to other components of long-run demand or supply equations. (C) 2007 Elsevier B.V. All rights reserved.

Hematological parameters and anaerobic threshold in Brazilian soccer players throughout a training program

We assessed the responses of hematological parameters and their relationship to the anaerobic threshold of Brazilian soccer players during a training program. Twelve athletes were evaluated at the beginning (week 0, T1), in the middle (week 6, T2), and at the end (week 12, T3) of the soccer training program. On the first day at 7:30 AM, before collecting the blood sample at rest for the determination of the hematological parameters, the athletes were conducted to the anthropometric evaluation. On the second day at 8:30 AM, the athletes had their anaerobic threshold measured. Analysis of variance with Newman-Keuls`post hoc was used for statistical comparisons between the parameters measured during the soccer training program. Correlations between the parameters analyzed were determined using the Pearson`s correlation coefficient. Erythrocytes concentration, hemoglobin, and hematocrit were significantly increased from T1 to T2. The specific soccer training program led to a rise in erythrocytes, hemoglobin, and hematocrit from T1 to T2. We assumed that these results occurred due to the plasma volume reduction and may be explained by the soccer training program characteristics. Furthermore, we did not observe any correlation between the anaerobic threshold and the hematological parameters.

Physical function interfering with pain and symptoms in fibromyalgia patients

Objectives. The aim of this study was to assess the relationship between variables of physical assessment - muscular strength, flexibility and dynamic balance - with pain, pain threshold, and fibromyalgia symptoms (FM). Methods. Our sample consists of 55 women, with age ranging from 30 to 55 years (mean of 46.5, (standard deviation, SD=6.6)), mean body mass index (BMI) of 28.7(3.8) and diagnosed for FM according to the American College of Rheumatology criteria. Pain intensity was measured using a visual analogue scale (VAS) and pain threshold (PT) using Fisher`s dolorimeter. FM symptoms were assessed by the Fibromyalgia Impact Questionnaire (FIQ); flexibility by the third finger to floor test (3FF); the muscular strength index (MSI) by the maximum volunteer isometric contraction at flexion and extension of right knee and elbow using a force transducer, dynamic balance by the time to get up and go (TUG) test and the functional reach test (FRT). Data were analysed using Pearson`s correlation, as well as simple and multivariate regression tests, with significance level of 5%. Results. PT and FIQ were weakly but significantly correlated with the TUG, MSI and 3FF as well as VAS with the TUG and MSI (p<0.05). VAS, PT and FIQ was not correlated with FRT. Simple regression suggests that, alone, TUG, FR, MSI and 3FF are low predictors of VAS, PT and FIQ. For the VAS, the best predictive model includes TUG and MSI, explaining 12.6% of pain. variability. For TP and total symptoms, as obtained by the FIQ, most predictive model includes 3FF and MSI, which respectively respond by 30% and 21% of the variability. Conclusion. Muscular strength, flexibility and balance are associated with pain, pain threshold, and symptoms in FM patients.

Risk Factors for Visual Field Progression in Treated Glaucoma

Objective: To determine intraocular pressure (IOP)-dependent and IOP-independent variables associated with visual field (VF) progression in treated glaucoma. Design: Retrospective cohort of the Glaucoma Progression Study. Methods: Consecutive, treated glaucoma patients with repeatable VF loss who had 8 or more VF examinations of either eye, using the Swedish Interactive Threshold Algorithm (24-2 SITA-Standard, Humphrey Field Analyzer II; Carl Zeiss Meditec, Inc, Dublin, California), during the period between January 1999 and September 2009 were included. Visual field progression was evaluated using automated pointwise linear regression. Evaluated data included age, sex, race, central corneal thickness, baseline VF mean deviation, mean follow-up IOP, peak IOP, IOP fluctuation, a detected disc hemorrhage, and presence of beta-zone parapapillary atrophy. Results: We selected 587 eyes of 587 patients (mean [SD] age, 64.9 [13.0] years). The mean (SD) number of VFs was 11.1 (3.0), spanning a mean (SD) of 6.4 (1.7) years. In the univariable model, older age (odds ratio [OR], 1.19 per decade; P = .01), baseline diagnosis of exfoliation syndrome (OR, 1.79; P = .01), decreased central corneal thickness (OR, 1.38 per 40 mu m thinner; P < .01), a detected disc hemorrhage (OR, 2.31; P < .01), presence of beta-zone parapapillary atrophy (OR, 2.17; P < .01), and all IOP parameters (mean follow-up, peak, and fluctuation; P < .01) were associated with increased risk of VF progression. In the multivariable model, peak IOP (OR, 1.13; P < .01), thinner central corneal thickness (OR, 1.45 per 40 mu m thinner; P < .01), a detected disc hemorrhage (OR, 2.59; P < .01), and presence of beta-zone parapapillary atrophy (OR, 2.38; P < .01) were associated with VF progression. Conclusions: IOP-dependent and IOP-independent risk factors affect disease progression in treated glaucoma. Peak IOP is a better predictor of progression than is IOP mean or fluctuation.

POISSON DISTRIBUTION TO ANALYZE NEAR-THRESHOLD MOTOR EVOKED POTENTIALS

Motor unit action potentials (MUAPs) evoked by repetitive, low-intensity transcranial magnetic stimulation can be modeled as a Poisson process. A mathematical consequence of such a model is that the ratio of the variance to the mean of the amplitudes of motor evoked potentials (MEPs) should provide an estimate of the mean size of the individual MUAPs that summate to generate each MEP. We found that this is, in fact, the case. Our finding thus supports the use of the Poisson distribution to model MEP generation and indicates that this model enables characterization of the motor unit population that contributes to near-threshold MEPs. Muscle Nerve 42: 825-828, 2010

Regression of glomerular injury by losartan in experimental diabetic nephropathy

Many features of chronic kidney disease may be reversed, but it is unclear whether advanced lesions, such as adhesions of sclerotic glomerular tufts to Bowman`s capsule (synechiae), can resolve during treatment. We previously showed, using a renal ablation model, that the renoprotective effect of the AT-1 receptor blocker, losartan, is dose-dependent. Here we determined if moderate and advanced glomerular lesions, associated with streptozotocin-induced diabetes, regress with conventional or high-dose losartan treatment. Using daily insulin injection for 10 months, we maintained diabetic adult male Munich-Wistar rats in a state of moderate hyperglycemia. Following this period, some rats continued to receive insulin with or without conventional or high-dose losartan for an additional 2 months. Diabetic rats pretreated with insulin for 10 months and age-matched non-diabetic rats served as controls. Mesangial expansion was found in the control diabetic rats and was exacerbated in those rats maintained on only insulin for an additional 2 months. Conventional and high-dose losartan treatments reduced this mesangial expansion and the severity of synechiae lesions below that found prior to treatment; however, the frequency of the latter was unchanged. There was no dose-response effect of losartan. Our results show that regression of mesangial expansion and contraction of sclerotic lesions is feasible in the treatment of diabetes, but complete resolution of advanced glomerulosclerosis may be hard to achieve.

Effect of low-dose combined oral contraceptive on aerobic capacity and anaerobic threshold level in active and sedentary young women

Background: The purpose of this study was to evaluate the effect of long-term use of oral contraceptives (DC) containing 0.20 mg of ethinylestradiol (EE) combined with 0.15 mg of gestodene (GEST) on the peak aerobic capacity and at the anaerobic threshold (AT) level in active and sedentary young women. Study Design: Eighty-eight women (23 +/- 2.1 years old) were divided into four groups active-OC (G1), active-NOC (G2), sedentary-OC (G3) and sedentary-NOC (G4) and were submitted to a continuous ergospirometric incremental test on a cycloergometer with 20 to 25 W min(-1) increments. Data were analyzed by two-way ANOVA with Tukey post hoc test. Level of significance was set at 5%. Results: The OC use effect for the variables relative and absolute oxygen uptake VO(2) mL kg(-1) min(-1); VO(2), L min(-1), respectively), carbon dioxide output (VCO(2), L min(-1)), ventilation (VE, L min(-1)), heart rate (HR, bpm), respiratory exchange ratio (RER) and power output (W) data, as well as the interaction between OC use and exercise effect on the peak of test and at the AT level did not differ significantly between the active groups (G1 and G2) and the sedentary groups (G3 and G4). As to the exercise effect, for all variables studied, it was noted that the active groups presented higher values for the variables VO(2), VCO(2), VE and power output (p<.05) than the sedentary groups. The RER and HR were similar (p>.05) at the peak and at the AT level between G1 vs. G3 and G2 vs. G4. Conclusions: Long-term use of OC containing EE 0.20 mg plus GEST 0.15 mg does not affect aerobic capacity at the peak and at the AT level of exercise tests. (C) 2010 Elsevier Inc. All rights reserved.