On elliptical multilevel models

Multilevel models have been widely applied to analyze data sets which present some hierarchical structure. In this paper we propose a generalization of the normal multilevel models, named elliptical multilevel models. This proposal suggests the use of distributions in the elliptical class, thus involving all symmetric continuous distributions, including the normal distribution as a particular case. Elliptical distributions may have lighter or heavier tails than the normal ones. In the case of normal error models with the presence of outlying observations, heavy-tailed error models may be applied to accommodate such observations. In particular, we discuss some aspects of the elliptical multilevel models, such as maximum likelihood estimation and residual analysis to assess features related to the fitting and the model assumptions. Finally, two motivating examples analyzed under normal multilevel models are reanalyzed under Student-t and power exponential multilevel models. Comparisons with the normal multilevel model are performed by using residual analysis.     

Ultrastructural elliptical models

Dolby's (1976) ultrastructural model with no replications is investigated within the class of the elliptical distributions. General asymptotic results are given for the sample covariance matrix S in the presence of incidental parameters. These results are used to study the asymptotic behaviour of some estimators of the slope parameter, unifying and extending existing results in the literature. In particular, under some regularity conditions they are shown to be consistent and asymptotically normal. For the special case of the structural model, some asymptotic relative efficiencies are also reported which show that generalized least squares and the method of moment estimators can be highly inefficient under nonnormality.     

BL-GARCH models with elliptical distributed innovations

In this work, we discuss the class of bilinear GARCH (BL-GARCH) models that are capable of capturing simultaneously two key properties of non-linear time series: volatility clustering and leverage effects. It has often been observed that the marginal distributions of such time series have heavy tails; thus we examine the BL-GARCH model in a general setting under some non-normal distributions. We investigate some probabilistic properties of this model and we conduct a Monte Carlo experiment to evaluate the small-sample performance of the maximum likelihood estimation (MLE) methodology for various models. Finally, within-sample estimation properties were studied using S&P 500 daily returns, when the features of interest manifest as volatility clustering and leverage effects. The main results suggest that the Student-t BL-GARCH seems highly appropriate to describe the S&P 500 daily returns.     

Diagnostics for elliptical linear mixed models with first-order autoregressive errors

For longitudinal time series data, linear mixed models that contain both random effects across individuals and first-order autoregressive errors within individuals may be appropriate. Some statistical diagnostics based on the models under a proposed elliptical error structure are developed in this work. It is well known that the class of elliptical distributions offers a more flexible framework for modelling since it contains both light- and heavy-tailed distributions. Iterative procedures for the maximum-likelihood estimates of the model parameters are presented. Score tests for the presence of autocorrelation and the homogeneity of autocorrelation coefficients among individuals are constructed. The properties of test statistics are investigated through Monte Carlo simulations. The local influence method for the models is also given. The analysed results of a real data set illustrate the values of the models and diagnostic statistics.     

A note on classical Stein-type estimators in elliptically contoured models

In this paper the conditions under which a broad class of Stein-type estimators dominates the best invariant unbiased estimator of the mean of an elliptically contoured population have been established. The superiority conditions are derived for both known and unknown scale structures. Also an example is given when the general scale matrix is assumed to be known in linear regression.     

On influence diagnostics in elliptical multivariate regression models with equicorrelated random errors

In this paper we discuss estimation and diagnostic procedures in elliptical multivariate regression models with equicorrelated random errors. Two procedures are proposed for the parameter estimation and the local influence curvatures are derived under some usual perturbation schemes to assess the sensitivity of the maximum likelihood estimates (MLEs). Two motivating examples preliminarily analyzed under normal errors are reanalyzed considering appropriate elliptical distributions. The local influence approach is used to compare the sensitivity of the model estimates.     

Inference in affine shape theory under elliptical models

This paper studies the elliptical statistical affine shape theory under certain particular conditions on the evenness or oddness of the number of landmarks. In such a case, the related distributions are polynomials, and the inference is easily performed; as an example, a landmark data is studied, and the performance of the polynomial density versus the usual series density is compared.     

Random-effect models with singular precision

We show that smoothing spline, intrinsic autoregression (IAR) and state-space model can be formulated as partially specified random-effect model with singular precision (SP). Various fitting methods have been suggested for the aforementioned models and this paper investigates the relationships among them, once the models have been placed under a single framework. Some methods have been previously shown to give the best linear unbiased predictors (BLUPs) under some random-effect models and here we show that they are in fact uniformly BLUPs (UBLUPs) under a class of models that are generated by the SP of random effects. We offer some new interpretations of the UBLUPs under models of SP and define BLUE and BLUP in these partially specified models without having to specify the covariance. We also show how the full likelihood inferences for random-effect models can be made for these models, so that the maximum likelihood (ML) and restricted maximum likelihood (REML) estimators can be used for the smoothing parameters in splines, etc.     

On diagnostic devices for proposing half-logistic and inverse half-logistic models using generalized (k) record values

In this work we consider the generalized upper (k) record values (GURV’s) and generalized lower (k) record values (GLRV’s) arising from half-logistic distribution (HLD) and inverse half-logistic distribution (IHLD). We derive some characterization results of HLD based on some moment relations of generalized upper (k) record values and those of generalized lower (k) record values and accordingly devised some diagnostic tools to identify HLD as a model to the distribution of a population. Similar characterization theorems and diagnostic tools are developed for IHLD as well. Simulation studies are conducted to validate the diagnostic tools devised for both HLD and IHLD.     

Bayesian meta-analysis using skewed elliptical distributions

Meta-analysis refers to a quantitative method for combining results from independent studies in order to draw overall conclusions. We consider hierarchical models including selection models under a skewed heavy tailed error distribution proposed originally by Chen, Dey, and Shao [M. H. Chen, D. K. Dey, Q. M. Shao, A new skewed link model for dichotomous quantal response data, J. Amer. Statist. Assoc. 94 (1983), pp. 1172–1186.] and Branco and Dey [D. Branco and D.K. Dey, A general class of multivariate skew-elliptical distributions, J. Multivariate Anal. 79, pp. 99–113.]. These rich classes of models combine the information of independent studies, allowing investigation of variability both between and within studies and incorporating weight functions. We constructed a detailed computational scheme under skewed normal and skewed Student's t distribution using the MCMC method. Bayesian model selection was conducted by Bayes factor under a different skewed error. Finally, we illustrated our methodology using a real data example taken from Johnson [M.F. Johnson, Comparative efficacy of Naf and SMFP dentifrices in caries prevention: a meta-analysis overview, J Eur. Organ. Caries Res. 27 (1993), pp. 328–336.].     

Testing Homogeneity of Inverse Gaussian Scale Parameters Based on Generalized Likelihood Ratio

This article presents a new procedure for testing homogeneity of scale parameters from k independent inverse Gaussian populations. Based on the idea of generalized likelihood ratio method, a new generalized p-value is derived. Some simulation results are presented to compare the performance of the proposed method and existing methods. Numerical results show that the proposed test has good size and power performance.     

Empirical transform estimation for indexed stochastic models

We present a method for estimating the parameters in indexed stochastic models via a least squares approach based on empirical transforms. Asymptotic approximations are derived for the distribution of the resulting estimators. An explicit expression for the mean-squared error provides a natural way of selecting the transform variable, and a numerical example illustrates the performance of the resulting method. A common finding, which we term 'diagonal optimization', occurs when multiparameter models are fitted by using transforms. Diagonal optimization arises when optimal performance results from equating the elements of the transform vector, and we provide a heuristic explanation of why this occurs.     

Quantile Regression via the EM Algorithm

The three-parameter asymmetric Laplace distribution (ALD) has received increasing attention in the field of quantile regression due to an important feature between its location and asymmetric parameters. On the basis of the representation of the ALD as a normal-variance–mean mixture with an exponential mixing distribution, this article develops EM and generalized EM algorithms, respectively, for computing regression quantiles of linear and nonlinear regression models. It is interesting to show that the proposed EM algorithm and the MM (Majorization–Minimization) algorithm for quantile regressions are really the same in terms of computation, since the updating formula of them are the same. This provides a good example that connects the EM and MM algorithms. Simulation studies show that the EM algorithm can successfully recover the true parameters in quantile regressions.     

Seemingly unrelated regression on the autoregressive (Ar(p)) singular equation system

This paper studies the estimation of seemingly unrelated regressions (SUR) of singular equation systems with an autoregressive error process (AR(p)) for each equation.Parameter estimates of the autoregressive singular equation system are not generally invariant to the equation deleted. Under the model specification restriction on the autoregressive parameters, the invariance property is preserved, and this paper shows that a single equation generalized least squares (GLS) estimation for a general autoregressive error process is equivalent to the SURGLS estimation of the AR(p) singular equation system.     

On local influence for elliptical linear models

The local influence method plays an important role in regression diagnostics and sensitivity analysis. To implement it, we need the Delta matrix for the underlying scheme of perturbations, in addition to the observed information matrix under the postulated model. Galea, Paula and Bolfarine (1997) has recently given the observed information matrix and the Delta matrix for a scheme of scale perturbations and has assessed of local influence for elliptical linear regression models. In the present paper, we consider the same elliptical linear regression models. We study the schemes of scale, predictor and response perturbations, and obtain their corresponding Delta matrices, respectively. To illustrate the methodology for assessment of local influence for these schemes and the implementation of the obtained results, we give an example.