LOCAL INFLUENCE IN NONLINEAR REPRODUCTIVE DISPERSION MODELS

Nonlinear reproductive dispersion models (NRDM, Jorgensen 1997) include a wider range of distributions and nonlinear models such as the possibility of correlated errors and nonlinear hypotheses dropping the exponential family assumption. Based on the generalized Cook distance and the conformal normal curvature of Poon & Poon (1999), local influence of minor perturbations on the data set is investigated for NRDM. Two examples are used to illustrate our results.     

LEVERAGE ADJUSTMENTS FOR DISPERSION MODELLING IN GENERALIZED NONLINEAR MODELS

For normal linear models, it is generally accepted that residual maximum likelihood estimation is appropriate when covariance components require estimation. This paper considers generalized linear models in which both the mean and the dispersion are allowed to depend on unknown parameters and on covariates. For these models there is no closed form equivalent to residual maximum likelihood except in very special cases. Using a modified profile likelihood for the dispersion parameters, an adjusted score vector and adjusted information matrix are found under an asymptotic development that holds as the leverages in the mean model become small. Subsequently, the expectation of the fitted deviances is obtained directly to show that the adjusted score vector is unbiased at least to O(1/n) . Exact results are obtained in the single‐sample case. The results reduce to residual maximum likelihood estimation in the normal linear case.     

EXPONENTIAL DISPERSION MODELS AND THE GAUSS-NEWTON ALGORITHM

It is known that the Fisher scoring iteration for generalized linear models has the same form as the Gauss-Newton algorithm for normal regression. This note shows that exponential dispersion models are the most general families to preserve this form for the scoring iteration. Therefore exponential dispersion models are the most general extension of generalized linear models for which the analogy with normal regression is preserved. The multinomial distribution is used as an example.     

Maximum-likelihood estimation and influence analysis in multivariate skew-normal reproductive dispersion mixed models for longitudinal data

Various mixed models were developed to capture the features of between- and within-individual variation for longitudinal data under the normality assumption of the random effect and the within-individual random error. However, the normality assumption may be violated in some applications. To this end, this article assumes that the random effect follows a skew-normal distribution and the within-individual error is distributed as a reproductive dispersion model. An expectation conditional maximization (ECME) algorithm together with the Metropolis-Hastings (MH) algorithm within the Gibbs sampler is presented to simultaneously obtain estimates of parameters and random effects. Several diagnostic measures are developed to identify the potentially influential cases and assess the effect of minor perturbation to model assumptions via the case-deletion method and local influence analysis. To reduce the computational burden, we derive the first-order approximations to case-deletion diagnostics. Several simulation studies and a real data example are presented to illustrate the newly developed methodologies.     

IMPROVED SCORE TESTS FOR VON MISES REGRESSION MODELS

We present, in matrix notation, a finite-sample correction formula to improve score tests in von Mises regression models with concentration covariates. The formula only requires simple operations on matrices and can be used to obtain analytically closed-form corrections for score test statistics in a variety of special von Mises models. The paper also provides a numerical comparison of the size of two score test statistics with bootstrap-based critical values.     

INFLUENCE FUNCTIONS IN FUNCTIONAL MEASUREMENT ERROR MODELS WITH REPLICATED DATA

We consider the construction and properties of influence functions in the context of functional measurement error models with replicated data. In these models estimates of the parameters can be affected both by the individual observations and the means of replicated observations. We show that influence function of the means of replicates on the estimate of regression coefficients can be only derived under the assumption that the variances of the errors are known, while one for the individual observations can be only derived simultaneously with their influence function on the estimators of the variances of the errors.     

LOG-LINEAR MODELS FOR MEAN AND DISPERSION IN MIXED POISSON REGRESSION MODELS

This paper is concerned with the analysis of repeated measures count data overdispersed relative to a Poisson distribution, with the overdispersion possibly heterogeneous. To accommodate the overdispersion, the Poisson random variable is compounded with a gamma random variable, and both the mean of the Poisson and the variance of the gamma are modelled using log linear models. Maximum likelihood estimates (MLE) are then obtained. The paper also gives extended quasi-likelihood estimates for a more general class of compounding distributions which are shown to be approximations to the MLEs obtained for the gamma case. The theory is illustrated by modelling the determination of asbestos fibre intensity on membrane filters mounted on microscope slides.     

ADJUSTING RESIDUALS IN ESTIMATION OF DISPERSION EFFECTS

This paper studies a method of adjusting the ordinary least squares residuals, when estimating and comparing dispersions, at various levels of factors in a replicated factorial experiment. Using a general dispersion model, theoretical results demonstrate the benefits of the method of adjusting residuals. An illustrative example is included.     

Influence Diagnostics in log-Birnbaum-Saunders Regression Models

In this paper we present various diagnostic methods for a linear regression model under a logarithmic Birnbaum-Saunders distribution for the errors, which may be applied for accelerated life testing or to compare the median lives of several populations. Some influence methods, such as the local influence, total local influence of an individual and generalized leverage are derived, analysed and discussed. We also present a connection between the local influence and generalized leverage methods. A discussion of the computation of the likelihood displacement as well as the normal curvature in the local influence method are presented. Finally, an example with real data is given for illustration.     

Stepwise local influence in generalized autoregressive conditional heteroskedasticity models

Detection of outliers or influential observations is an important work in statistical modeling, especially for the correlated time series data. In this paper we propose a new procedure to detect patch of influential observations in the generalized autoregressive conditional heteroskedasticity (GARCH) model. Firstly we compare the performance of innovative perturbation scheme, additive perturbation scheme and data perturbation scheme in local influence analysis. We find that the innovative perturbation scheme give better result than other two schemes although this perturbation scheme may suffer from masking effects. Then we use the stepwise local influence method under innovative perturbation scheme to detect patch of influential observations and uncover the masking effects. The simulated studies show that the new technique can successfully detect a patch of influential observations or outliers under innovative perturbation scheme. The analysis based on simulation studies and two real data sets show that the stepwise local influence method under innovative perturbation scheme is efficient for detecting multiple influential observations and dealing with masking effects in the GARCH model.     

MAXIMUM LIKELIHOOD ESTIMATION IN LINEAR MODELS WITH EQUI-CORRELATED RANDOM ERRORS

Necessary and sufficient conditions for the existence of maximum likelihood estimators of unknown parameters in linear models with equi‐correlated random errors are presented. The basic technique we use is that these models are, first, orthogonally transformed into linear models with two variances, and then the maximum likelihood estimation problem is solved in the environment of transformed models. Our results generalize a result of Arnold, S. F. (1981) [The theory of linear models and multivariate analysis. Wiley, New York]. In addition, we give necessary and sufficient conditions for the existence of restricted maximum likelihood estimators of the parameters. The results of Birkes, D. & Wulff, S. (2003) [Existence of maximum likelihood estimates in normal variance‐components models. J Statist Plann. Inference. 113 , 35–47] are compared with our results and differences are pointed out.     

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 Galea, M., Leiva, V. and Paula, G. A. 2004. Influence diagnostics in log-Birnbaum–Saunders regression models. J. Appl. Stat., 31: 1049–1064. [Taylor & Francis Online], [Web of Science ®] [Google Scholar]] 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.     

Perturbation diagnostics of autocorrelation coefficients in non linear mixed-effects models with AR(1) errors based on M-estimation

In this work we propose and analyze non linear mixed-effects models for longitudinal data, which are widely used in the fields of economics, biopharmaceuticals, agriculture, and so on. A robust method to obtain maximum likelihood estimates for the parameters is presented, as well as perturbation diagnostics of autocorrelation coefficient in non linear models based on robust estimates and influence curvature. The obtained results are illustrated by plasma concentrations data presented in Davidian and Giltinan, which was analyzed under the non robust situation.     

Analysis of Circular Longitudinal Data Based on Generalized Estimating Equations

This paper derives estimating equations for modelling circular data with longitudinal structure for a family of circular distributions with two parameters. Estimating equations for modelling the circular mean and the resultant length are given separately. Estimating equations are then derived for a mixed model. This paper shows that the estimators that follow from these equations are consistent and asymptotically normal. The results are illustrated by an example about the direction taken by homing pigeons.     

Detection of influential observations in nonlinear reproductive dispersion models

A robust case-deletion diagnostic which could avoid masking is presented for nonlinear reproductive dispersion models (NRDM) (Jorgensen, 1997) that include generalized linear models and exponential family nonlinear models as special cases. Based on the second-order approximation of log-likelihood displacement and Poon & Poon's (1999) conformal normal curvature, a measure of local influence ranging from 0 to 1 is constructed. An example illustrates application of the techniques.     

Diagnostics for generalized Poisson regression models with errors in variables

In this paper, we develop diagnostic methods for generalized Poisson regression (GPR) models with errors in variables based on the corrected likelihood. The one-step approximations of the estimates in the case-deletion model are given and case-deletion and local influence measures are presented. Meanwhile, based on a corrected score function, the testing statistics for the significance of dispersion parameters in GPR models with measurement errors are investigated. Finally, illustration of our methodology is given through numerical examples.