Exact Likelihood Equations for Autoregression Models with Multivariate Elliptically Contoured Distributions

Abstract

The multivariate elliptically contoured distributions provide a viable framework for modeling time-series data. It includes the multivariate normal, power exponential, t, and Cauchy distributions as special cases. For multivariate elliptically contoured autoregressive models, we derive the exact likelihood equations for the model parameters. They are closely related to the Yule-Walker equations and involve simple function of the data. The maximum likelihood estimators are obtained by alternately solving two linear systems and illustrated using the simulation data.     

Concomitants of multivariate order statistics from multivariate elliptical distributions

ABSTRACT

In this article, we consider a (k + 1)n-dimensional elliptically contoured random vector (X^T₁, X₂^T, …, X^T_k, Z^T)^T = (X₁₁, …, X_1n, …, X_k1, …, X_kn, Z₁, …, Z_n)^T and derive the distribution of concomitant of multivariate order statistics arising from X₁, X₂, …, X_k. Specially, we derive a mixture representation for concomitant of bivariate order statistics. The joint distribution of the concomitant of bivariate order statistics is also obtained. Finally, the usefulness of our result is illustrated by a real-life data.     

Shrinkage Estimation Under Multivariate Elliptic Models

The estimation of the location vector of a p-variate elliptically contoured distribution (ECD) is considered using independent random samples from two multivariate elliptically contoured populations when it is apriori suspected that the location vectors of the two populations are equal. For the setting where the covariance structure of the populations is the same, we define the maximum likelihood, Stein-type shrinkage and positive-rule shrinkage estimators. The exact expressions for the bias and quadratic risk functions of the estimators are derived. The comparison of the quadratic risk functions reveals the dominance of the Stein-type estimators if p ≥ 3. A graphical illustration of the risk functions under a “typical” member of the elliptically contoured family of distributions is provided to confirm the analytical results.     

Regression model with elliptically contoured errors

For the regression model y=X β+ε where the errors follow the elliptically contoured distribution, we consider the least squares, restricted least squares, preliminary test, Stein-type shrinkage and positive-rule shrinkage estimators for the regression parameters, β.

We compare the quadratic risks of the estimators to determine the relative dominance properties of the five estimators.     

A new family of life distributions based on the elliptically contoured distributions

In this article we propose a new family of life distributions, generated from an elliptically contoured distribution, in which the density and some of its properties are obtained. Explicit expressions for the density are found for a large number of specific elliptical distributions, such as Pearson type VII, t, Cauchy, Kotz type, normal, Bessel, Laplace and logistic.     

Concomitants of order statistics from multivariate elliptical distributions

In this paper, by considering a 2n-dimensional elliptically contoured random vector (X^T,Y^T)^T=(X₁,…,X_n,Y₁,…,Y_n)^T, we derive the exact joint distribution of linear combinations of concomitants of order statistics arising from X. Specifically, we establish a mixture representation for the distribution of the rth concomitant order statistic, and also for the joint distribution of the rth order statistic and its concomitant. We show that these distributions are indeed mixtures of multivariate unified skew-elliptical distributions. The two most important special cases of multivariate normal and multivariate t distributions are then discussed in detail. Finally, an application of the established results in an inferential problem is outlined.     

A note on the orthant probability of the elliptically contoured distribution

When a scale matrix satisfies certain conditions, the orthant probability of the elliptically contoured distribution with the scale matrix is expressed as the same probability of the equicorrelated normal distribution.     

Asymptotic tail behavior of a random sum with conditionally dependent subexponential summands

In this paper, we obtain some results for the asymptotic behavior of the tail probability of a random sum S_τ = ∑^τ_{k = 1}X_k, where the summands X_k, k = 1, 2, …, are conditionally dependent random variables with a common subexponential distribution F, and the random number τ is a non negative integer-valued random variable, independent of {X_k: k ? 1}.     

Bayesian analysis in multivariate regression models with conjugate priors

In this paper, we consider the full rank multivariate regression model with matrix elliptically contoured distributed errors. We formulate a conjugate prior distribution for matrix elliptical models and derive the posterior distributions of mean and scale matrices. In the sequel, some characteristics of regression matrix parameters are also proposed.     

A study of generalized skew-normal distribution

Following the paper by Genton and Loperfido [Generalized skew-elliptical distributions and their quadratic forms, Ann. Inst. Statist. Math. 57 (2005), pp. 389–401], we say that Z has a generalized skew-normal distribution, if its probability density function (p.d.f.) is given by f(z)=2φ _p (z; ξ, Ω)π (z?ξ), z∈? ^p , where φ _p (·; ξ, Ω) is the p-dimensional normal p.d.f. with location vector ξ and scale matrix Ω, ξ∈? ^p , Ω>0, and π is a skewing function from ? ^p to ?, that is 0≤π (z)≤1 and π (?z)=1?π (z), ? z∈? ^p . First the distribution of linear transformations of Z are studied, and some moments of Z and its quadratic forms are derived. Next we obtain the joint moment-generating functions (m.g.f.’s) of linear and quadratic forms of Z and then investigate conditions for their independence. Finally explicit forms for the above distributions, m.g.f.’s and moments are derived when π (z)=κ (α′z), where α∈? ^p and κ is the normal, Laplace, logistic or uniform distribution function.     

A note on predictive inference for multivariate elliptically contoured distributions

The prediction problem is considered for the multivariate regression model with an elliptically contoured error distribution. We show that the predictive distribution under elliptical errors assumption is the same as that obtained under normally distributed error in both the Bayesian approach using an im-proper prior and the classical approach. This gives inference robustness with respect to departures from the reference case of independent sampling from the normal distribution.     

Geometric approach to the skewed normal distribution

The representations of the skewed normal distribution given in Propositions 1–4 in Genton (Ed., 2004) are considered here from a unified geometric point of view and are, based upon this, generalized in two respects. On the one hand, the four concrete representations motivate us for a unified and much more general algebraic–geometric representation of the skewed normal distribution (Theorems 1 and 2 as well as Remarks 2 and 3); on the other hand, the mentioned representations are generalized to the elliptically contoured case (Propositions and Corollaries 1c–4c).     

A bayesian approach to prediction in analysis of variance settings

This paper takes the results of Lindley and smith ( 1972 ) one step further, by finding the predictive distribution of an observation y^* whose distribution is normal, and centred at A^* ₁θ₁ We then apply this distribution to the case of prediction based on data obtained in one and two wau ANOVA situations. For Example, it turns out that for two way ANOVA with interaction, the predictive mean, (which we would use as the predictor) is a weighted combination of sample main effects and interaction effects     

Correlation analysis of ordered symmetrically dependent observations and their concomitants of order statistics

Given two jointly observed random vectors Y and Z of the same dimension, let Y be a reordered version of Y and Z the resulting vector of concomitants of order statistics. When X is a covariate of interest, also jointly observed with Y, the authors obtain the joint covariance structure of (X, y, Z) and the related correlation parameters explicitly, under the assumption that the vector (X, Y, Z) is normal and that its joint covariance structure is permutation symmetric. They also discuss extensions to elliptically contoured distributions.     

Generalized Tobit models: diagnostics and application in econometrics

The standard Tobit model is constructed under the assumption of a normal distribution and has been widely applied in econometrics. Atypical/extreme data have a harmful effect on the maximum likelihood estimates of the standard Tobit model parameters. Then, we need to count with diagnostic tools to evaluate the effect of extreme data. If they are detected, we must have available a Tobit model that is robust to this type of data. The family of elliptically contoured distributions has the Laplace, logistic, normal and Student-t cases as some of its members. This family has been largely used for providing generalizations of models based on the normal distribution, with excellent practical results. In particular, because the Student-t distribution has an additional parameter, we can adjust the kurtosis of the data, providing robust estimates against extreme data. We propose a methodology based on a generalization of the standard Tobit model with errors following elliptical distributions. Diagnostics in the Tobit model with elliptical errors are developed. We derive residuals and global/local influence methods considering several perturbation schemes. This is important because different diagnostic methods can detect different atypical data. We implement the proposed methodology in an R package. We illustrate the methodology with real-world econometrical data by using the R package, which shows its potential applications. The Tobit model based on the Student-t distribution with a small quantity of degrees of freedom displays an excellent performance reducing the influence of extreme cases in the maximum likelihood estimates in the application presented. It provides new empirical evidence on the capabilities of the Student-t distribution for accommodation of atypical data.     

Prediction intervals for growth curves

The growth curve model Y^n×p = A^{n×p ξ mtimes;k}B^k×p+ E^nxp, where Y is an observation matrix, &sigma is a matrix of unknown parameters, A is a known matrix of rank m, B is a known matrix of rank k with 1'= (1, …, 1) as its first row, and the rows of E are independent each distributed as N^p(0,Σ,) is considered. The problem of constructing the prediction intervals for future observations using the above model is considered and approximate intervals assuming different structures on σ are derived. The results are illustrated with several data sets.     

A generalization of the multivariate slash distribution

In this paper, we propose a generalization of the multivariate slash distribution and investigate some of its properties. We show that the new distribution belongs to the elliptically contoured distributions family, and can have heavier tails than the multivariate slash distribution. Therefore, this generalization of the multivariate slash distribution can be considered as an alternative heavy-tailed distribution for modeling data sets in a variety of settings. We apply the generalized multivariate slash distribution to two real data sets to provide some illustrative examples.     

Reliability analysis of systems with components having two dependent subcomponents

In this article, a system that consists of n independent components each having two dependent subcomponents (A_i, B_i), i = 1, …, n is considered. The system is assumed to compose of components that have two correlated subcomponents (A_i, B_i), and functions iff both systems of subcomponents A₁, A₂, …, A_n and B₁, B₂, …, B_n work under certain structural rules. The expressions for reliability and mean time to failure of such systems are obtained. A sufficient condition to compare two systems of bivariate components in terms of stochastic ordering is also presented.     

On the Distribution of Matrix Quadratic Forms

A characterization of the distribution of the multivariate quadratic form given by X A X′, where X is a p × n normally distributed matrix and A is an n × n symmetric real matrix, is presented. We show that the distribution of the quadratic form is the same as the distribution of a weighted sum of non central Wishart distributed matrices. This is applied to derive the distribution of the sample covariance between the rows of X when the expectation is the same for every column and is estimated with the regular mean.