Influence diagnostics in beta regression

We consider the issue of assessing influence of observations in the class of beta regression models, which is useful for modelling random variables that assume values in the standard unit interval and are affected by independent variables. We propose a Cook-like distance and also measures of local influence under different perturbation schemes. Applications using real data are presented. (c) 2008 Elsevier B.V.. All rights reserved.

Improved maximum-likelihood estimation in a regression model with general parametrization

We analyse the finite-sample behaviour of two second-order bias-corrected alternatives to the maximum-likelihood estimator of the parameters in a multivariate normal regression model with general parametrization proposed by Patriota and Lemonte [A. G. Patriota and A. J. Lemonte, Bias correction in a multivariate regression model with genereal parameterization, Stat. Prob. Lett. 79 (2009), pp. 1655-1662]. The two finite-sample corrections we consider are the conventional second-order bias-corrected estimator and the bootstrap bias correction. We present the numerical results comparing the performance of these estimators. Our results reveal that analytical bias correction outperforms numerical bias corrections obtained from bootstrapping schemes.

Covariance matrix formula for Birnbaum-Saunders regression models

The Birnbaum-Saunders regression model is commonly used in reliability studies. We derive a simple matrix formula for second-order covariances of maximum-likelihood estimators in this class of models. The formula is quite suitable for computer implementation, since it involves only simple operations on matrices and vectors. Some simulation results show that the second-order covariances can be quite pronounced in small to moderate sample sizes. We also present empirical applications.

Influence diagnostics in Birnbaum-Saunders nonlinear regression models

We consider the issue of assessing influence of observations in the class of Birnbaum-Saunders nonlinear regression models, which is useful in lifetime data analysis. Our results generalize those in Galea et al. [8] which are confined to Birnbaum-Saunders linear regression models. Some influence methods, such as the local influence, total local influence of an individual and generalized leverage are discussed. Additionally, the normal curvatures for studying local influence are derived under some perturbation schemes. We also give an application to a real fatigue data set.

Asymptotic skewness in Birnbaum-Saunders nonlinear regression models

The family of distributions proposed by Birnbaum and Saunders (1969) can be used to model lifetime data and it is widely applicable to model failure times of fatiguing materials. We give a simple matrix formula of order n(-1/2), where n is the sample size, for the skewness of the distributions of the maximum likelihood estimates of the parameters in Birnbaum-Saunders nonlinear regression models, recently introduced by Lemonte and Cordeiro (2009). The formula is quite suitable for computer implementation, since it involves only simple operations on matrices and vectors, in order to obtain closed-form skewness in a wide range of nonlinear regression models. Empirical and real applications are analyzed and discussed. (C) 2010 Elsevier B.V. All rights reserved.

Some restriction tests in a new class of regression models for proportions

The main purpose of this work is to study the behaviour of Skovgaard`s [Skovgaard, I.M., 2001. Likelihood asymptotics. Scandinavian journal of Statistics 28, 3-32] adjusted likelihood ratio statistic in testing simple hypothesis in a new class of regression models proposed here. The proposed class of regression models considers Dirichlet distributed observations, and the parameters that index the Dirichlet distributions are related to covariates and unknown regression coefficients. This class is useful for modelling data consisting of multivariate positive observations summing to one and generalizes the beta regression model described in Vasconcellos and Cribari-Neto [Vasconcellos, K.L.P., Cribari-Neto, F., 2005. Improved maximum likelihood estimation in a new class of beta regression models. Brazilian journal of Probability and Statistics 19,13-31]. We show that, for our model, Skovgaard`s adjusted likelihood ratio statistics have a simple compact form that can be easily implemented in standard statistical software. The adjusted statistic is approximately chi-squared distributed with a high degree of accuracy. Some numerical simulations show that the modified test is more reliable in finite samples than the usual likelihood ratio procedure. An empirical application is also presented and discussed. (C) 2009 Elsevier B.V. All rights reserved.

Birnbaum-Saunders nonlinear regression models

We introduce, for the first time, a new class of Birnbaum-Saunders nonlinear regression models potentially useful in lifetime data analysis. The class generalizes the regression model described by Rieck and Nedelman [Rieck, J.R., Nedelman, J.R., 1991. A log-linear model for the Birnbaum-Saunders distribution. Technometrics 33, 51-60]. We discuss maximum-likelihood estimation for the parameters of the model, and derive closed-form expressions for the second-order biases of these estimates. Our formulae are easily computed as ordinary linear regressions and are then used to define bias corrected maximum-likelihood estimates. Some simulation results show that the bias correction scheme yields nearly unbiased estimates without increasing the mean squared errors. Two empirical applications are analysed and discussed. Crown Copyright (C) 2009 Published by Elsevier B.V. All rights reserved.

Bias correction in a multivariate normal regression model with general parameterization

This paper derives the second-order biases Of maximum likelihood estimates from a multivariate normal model where the mean vector and the covariance matrix have parameters in common. We show that the second order bias can always be obtained by means of ordinary weighted least-squares regressions. We conduct simulation studies which indicate that the bias correction scheme yields nearly unbiased estimators. (C) 2009 Elsevier B.V. All rights reserved.

Aminoacetone, a putative endogenous source of methylglyoxal, causes oxidative stress and death to insulin-producing RINm5f cells

Aminoacetone (AA), triose phosphates, and acetone are putative endogenous sources of potentially cytotoxic and genotoxic methylglyoxal (MG), which has been reported to be augmented in the plasma of diabetic patients. In these patients, accumulation of MG derived from aminoacetone, a threonine and glycine catabolite, is inferred from the observed concomitant endothelial overexpression of circulating semicarbazide-sensitive amine oxidases. These copper-dependent enzymes catalyze the oxidation of primary amines, such as AA and methylamine, by molecular oxygen, to the corresponding aldehydes, NH4+ ion and H2O2. We recently reported that AA aerobic oxidation to MG also takes place immediately upon addition of catalytic amounts of copper and iron ions. Taking into account that (i) MG and H2O2 are reportedly cytotoxic to insulin-producing cell lineages such as RINm5f and that (ii) the metal-catalyzed oxidation of AA is propagated by O-2(center dot-) radical anion, we decided to investigate the possible pro-oxidant action of AA on these cells taken here as a reliable model system for pancreatic beta-cells. Indeed, we show that AA (0.10-5.0 mM) administration to RINm5f cultures induces cell death. Ferrous (50-300 mu M) and Fe3+ ion (100 mu M) addition to the cell cultures had no effect, whereas Cu2+ (5.0-100 mu M) significantly increased cell death. Supplementation of the AA- and Cu2+-containing culture medium with antioxidants, such as catalase (5.0 mu M), superoxide dismutase (SOD, 50 U/mL), and N-acetylcysteine (NAC, 5.0 mM) led to partial protection. mRNA expression of MnSOD, CuZnSOD, glutathione peroxidase, and glutathione reductase, but not of catalase, is higher in cells treated with AA (0.50-1.0 mM) plus Cu2+ ions (10-50 mu M) relative to control cultures. This may imply higher activity of antioxidant enzymes C, in RINm5f AA-treated cells. In addition, we have found that AA (0.50-1.0 mM) Plus Cu2+ (100 mu M) (i) increase RINm5f cytosolic calcium; (ii) promote DNA fragmentation; and (iii) increase the pro-apoptotic (Bax)/antiapoptotic (Bcl-2) ratio at the level of mRNA expression. In conclusion, although both normal and pathological concentrations of AA are probably much lower than those used here, it is tempting to propose that excess AA in diabetic patients may drive oxidative damage and eventually the death of pancreatic beta-cells.

Human endogenous RNAs: Implications for the immunomodulation of Toll-like receptor 3

Toll-like receptors (TLRs), a family of mammalian receptors, are able to recognize nucleic acids. TLR3 recognizes double-stranded (ds)RNA, a product of the replication of certain viruses. Polyinosinic-polycytidylic acid, referred to as poly(I:C), an analog of viral dsRNA, interacts with TLR3 thereby eliciting immunoinflammatory responses characteristic of viral infection or down-regulating the expression of chemokine receptor CXCR4. It is known that dsRNA also directly activates interferon (IFN)-induced enzymes, such as the RNA-dependent protein kinase (PKR). In the present study, the mRNA expression of TLR3, CXCR4, IFN gamma and PKR was investigated in a culture of peripheral blood mononuclear cells (PBMCs) stimulated with poly(I:C) and endogenous RNA from human PBMCs. No cytotoxic effect on the cells or on the proliferation of CD3(+), CD4(+) and CD8(+) cells was observed. TLR3 expression in the PBMCs in the presence of poly(I:C) was up-regulated 9.5-fold, and TLR3 expression in the PBMCs treated with endogenous RNA was down-regulated 1.8-fold (p=0.002). The same trend was observed for IFN gamma where in the presence of poly(I:C) an 8.7-fold increase was noted and in the presence of endogenous RNA a 3.1-fold decrease was observed. In the culture activated with poly(1:C), mRNA expression of CXCR4 increased 8.0-fold and expression of PKR increased 33.0-fold. Expression of these genes decreased in the culture treated with endogenous RNA when compared to the culture without stimulus. Thus, high expression of mRNA for TLR3, IFN gamma, CXCR4 and PKR was observed in the presence of poly(I:C) and low expression was observed in the cells cultured with endogenous RNA. In conclusion, TLR3 may play major physiological roles that are not in the context of viral infection. It is possible that RNA released from cells could contain enough double-stranded structures to regulate cell activation. The involvement of endogenous RNA in endogenous gene expression and its implications in the regulation thereof, are still being studied, and will have significant implications in the future.

Fluoxetine and norfluoxetine analysis by direct injection of human plasma in a column switching liquid chromatographic system

A column switching LC method is presented for the analysis of fluoxetine (FLU) and norfluoxetine (NFLU) by direct injection of human plasma using a lab-made restricted access media (RAM) column. A RAM-BSA-octadecyl silica (C-18) column (40 min x 4.6 mm, 10 mu m) is evaluated in both backflush and foreflush elution modes and coupled with a C-18 lab-made (50 mm x 4.6 mm, 3 pm) analytical column in order to perform online sample preparation. Direct injection of 100 mu L, of plasma samples is possible with the developed approach. In addition, reduction of sample handling is obtained when compared with traditional liquid-liquid extraction (LLE) and SPE. The total analysis time is around 20 min. A LOQ of 15 ng/mL is achieved in a concentration range of 15-500 ng/mL, allowing the therapeutic drug monitoring of clinical samples. The precision values achieved are lower than 15% for all the evaluated points with adequate recovery and accuracy. Furthermore, no matrix interferences are found in the analysis and the proposed method shows to be an adequate alternative for analysis of FLU in plasma.

An Analysis of the Telecommunications Business in China by Linear Regression

In this paper, we study the influence of the National Telecom Business Volume by the data in 2008 that have been published in China Statistical Yearbook of Statistics. We illustrate the procedure of modeling “National Telecom Business Volume” on the following eight variables, GDP, Consumption Levels, Retail Sales of Social Consumer Goods Total Renovation Investment, the Local Telephone Exchange Capacity, Mobile Telephone Exchange Capacity, Mobile Phone End Users, and the Local Telephone End Users. The testing of heteroscedasticity and multicollinearity for model evaluation is included. We also consider AIC and BIC criterion to select independent variables, and conclude the result of the factors which are the optimal regression model for the amount of telecommunications business and the relation between independent variables and dependent variable. Based on the final results, we propose several recommendations about how to improve telecommunication services and promote the economic development.

Testing linear cointegration against smooth-transition cointegration

This paper studies a smooth-transition (ST) type cointegration. The proposed ST cointegration allows for regime switching structure in a cointegrated system. It nests the linear cointegration developed by Engle and Granger (1987) and the threshold cointegration studied by Balke and Fomby (1997). We develop F-type tests to examine linear cointegration against ST cointegration in ST-type cointegrating regression models with or without time trends. The null asymptotic distributions of the tests are derived with stationary transition variables in ST cointegrating regression models. And it is shown that our tests have nonstandard limiting distributions expressed in terms of standard Brownian motion when regressors are pure random walks, while have standard asymptotic distributions when regressors contain random walks with nonzero drift. Finite-sample distributions of those tests are studied by Monto Carlo simulations. The small-sample performance of the tests states that our F-type tests have a better power when the system contains ST cointegration than when the system is linearly cointegrated. An empirical example for the purchasing power parity (PPP) data (monthly US dollar, Italy lira and dollar-lira exchange rate from 1973:01 to 1989:10) is illustrated by applying the testing procedures in this paper. It is found that there is no linear cointegration in the system, but there exits the ST-type cointegration in the PPP data.

A Note on Proxy Variables and Instrument Variable Regression

This is a note about proxy variables and instruments for identification of structural parameters in regression models. We have experienced that in the econometric textbooks these two issues are treated separately, although in practice these two concepts are very often combined. Usually, proxy variables are inserted in instrument variable regressions with the motivation they are exogenous. Implicitly meaning they are exogenous in a reduced form model and not in a structural model. Actually if these variables are exogenous they should be redundant in the structural model, e.g. IQ as a proxy for ability. Valid proxies reduce unexplained variation and increases the efficiency of the estimator of the structural parameter of interest. This is especially important in situations when the instrument is weak. With a simple example we demonstrate what is required of a proxy and an instrument when they are combined. It turns out that when a researcher has a valid instrument the requirements on the proxy variable is weaker than if no such instrument exists

Computerized identification of motor complications in Parkinson's disease

Objective: To investigate whether spirography-based objective measures are able to effectively characterize the severity of unwanted symptom states (Off and dyskinesia) and discriminate them from motor state of healthy elderly subjects. Background: Sixty-five patients with advanced Parkinson’s disease (PD) and 10 healthy elderly (HE) subjects performed repeated assessments of spirography, using a touch screen telemetry device in their home environments. On inclusion, the patients were either treated with levodopa-carbidopa intestinal gel or were candidates for switching to this treatment. On each test occasion, the subjects were asked trace a pre-drawn Archimedes spiral shown on the screen, using an ergonomic pen stylus. The test was repeated three times and was performed using dominant hand. A clinician used a web interface which animated the spiral drawings, allowing him to observe different kinematic features, like accelerations and spatial changes, during the drawing process and to rate different motor impairments. Initially, the motor impairments of drawing speed, irregularity and hesitation were rated on a 0 (normal) to 4 (extremely severe) scales followed by marking the momentary motor state of the patient into 2 categories that is Off and Dyskinesia. A sample of spirals drawn by HE subjects was randomly selected and used in subsequent analysis. Methods: The raw spiral data, consisting of stylus position and timestamp, were processed using time series analysis techniques like discrete wavelet transform, approximate entropy and dynamic time warping in order to extract 13 quantitative measures for representing meaningful motor impairment information. A principal component analysis (PCA) was used to reduce the dimensions of the quantitative measures into 4 principal components (PC). In order to classify the motor states into 3 categories that is Off, HE and dyskinesia, a logistic regression model was used as a classifier to map the 4 PCs to the corresponding clinically assigned motor state categories. A stratified 10-fold cross-validation (also known as rotation estimation) was applied to assess the generalization ability of the logistic regression classifier to future independent data sets. To investigate mean differences of the 4 PCs across the three categories, a one-way ANOVA test followed by Tukey multiple comparisons was used. Results: The agreements between computed and clinician ratings were very good with a weighted area under the receiver operating characteristic curve (AUC) coefficient of 0.91. The mean PC scores were different across the three motor state categories, only at different levels. The first 2 PCs were good at discriminating between the motor states whereas the PC3 was good at discriminating between HE subjects and PD patients. The mean scores of PC4 showed a trend across the three states but without significant differences. The Spearman’s rank correlations between the first 2 PCs and clinically assessed motor impairments were as follows: drawing speed (PC1, 0.34; PC2, 0.83), irregularity (PC1, 0.17; PC2, 0.17), and hesitation (PC1, 0.27; PC2, 0.77). Conclusions: These findings suggest that spirography-based objective measures are valid measures of spatial- and time-dependent deficits and can be used to distinguish drug-related motor dysfunctions between Off and dyskinesia in PD. These measures can be potentially useful during clinical evaluation of individualized drug-related complications such as over- and under-medications thus maximizing the amount of time the patients spend in the On state.