A note on the multivariate generalized asymmetric Laplace motion

Abstract

In this note, we use multivariate subordination to introduce a multivariate extension of the generalized asymmetric Laplace motion. The class introduced provides a unified framework for several multivariate extensions of the popular variance gamma process. We also show that the associated time one distribution extends the multivariate generalized asymmetric Laplace distributions proposed in the statistical literature.     

A Three-Parameter Asymmetric Laplace Distribution and Its Extension

ABSTRACT

In this article, a new three-parameter asymmetric Laplace distribution and its extension are introduced. This includes as special case the symmetric Laplace double-exponential distribution. The distribution has established a direct link to estimation of quantile and quantile regression. Properties of the new distribution are presented. Application is made to a flood data modeling example.     

Influence diagnostics for Grubbs’s model with asymmetric heavy-tailed distributions

Grubbs’s model (Grubbs, Encycl Stat Sci 3:42–549, 1983) is used for comparing several measuring devices, and it is common to assume that the random terms have a normal (or symmetric) distribution. In this paper, we discuss the extension of this model to the class of scale mixtures of skew-normal distributions. Our results provide a useful generalization of the symmetric Grubbs’s model (Osorio et al., Comput Stat Data Anal, 53:1249–1263, 2009) and the asymmetric skew-normal model (Montenegro et al., Stat Pap 51:701–715, 2010). We discuss the EM algorithm for parameter estimation and the local influence method (Cook, J Royal Stat Soc Ser B, 48:133–169, 1986) for assessing the robustness of these parameter estimates under some usual perturbation schemes. The results and methods developed in this paper are illustrated with a numerical example.     

A flexible generalization of the skew normal distribution based on a weighted normal distribution

The skew normal distribution of Azzalini (Scand J Stat 12:171–178, 1985) has been found suitable for unimodal density but with some skewness present. Through this article, we introduce a flexible extension of the Azzalini (Scand J Stat 12:171–178, 1985) skew normal distribution based on a symmetric component normal distribution (Gui et al. in J Stat Theory Appl 12(1):55–66, 2013). The proposed model can efficiently capture the bimodality, skewness and kurtosis criteria and heavy-tail property. The paper presents various basic properties of this family of distributions and provides two stochastic representations which are useful for obtaining theoretical properties and to simulate from the distribution. Further, maximum likelihood estimation of the parameters is studied numerically by simulation and the distribution is investigated by carrying out comparative fitting of three real datasets.     

Tail dependence for skew Laplace distribution and skew Cauchy distribution

ABSTRACT

Coefficient of tail dependence measures the strength of dependence in the tail of a bivariate distribution and it has been found useful in the risk management. In this paper, we derive the upper tail dependence coefficient for a random vector following the skew Laplace distribution and the skew Cauchy distribution, respectively. The result shows that skew Laplace distribution is asymptotically independent in upper tail, however, skew Cauchy distribution has asymptotic upper tail dependence.     

A new family of multivariate skew slash distribution

In this article, the new family of multivariate skew slash distribution is defined. According to the definition, a stochastic representation of the multivariate skew slash distribution is derived. The first four moments and measures of skewness and kurtosis of a random vector with the multivariate skew slash distribution are obtained. The distribution of quadratic forms for the multivariate skew slash distribution and the non central skew slash χ² distribution are studied. Maximum likelihood inference and real data illustration are discussed. In the end, the potential extension of multivariate skew slash distribution is discussed.     

A multivariate two-factor skew model

It is well known that many data, such as the financial or demographic data, exhibit asymmetric distributions. In recent years, researchers have concentrated their efforts to model this asymmetry. Skew normal model is one of such models that are skew and yet possess many properties of the normal model. In this paper, a new multivariate skew model is proposed, along with its statistical properties. It includes the multivariate normal distribution and multivariate skew normal distribution as special cases. The quadratic form of this random vector follows a χ² distribution. The roles of the parameters in the model are investigated using contour plots of bivariate densities.     

On diagnostics in multivariate measurement error models under asymmetric heavy-tailed distributions

In this paper, we discuss the extension of some diagnostic procedures to multivariate measurement error models with scale mixtures of skew-normal distributions (Lachos et?al., Statistics 44:541?C556, 2010c). This class provides a useful generalization of normal (and skew-normal) measurement error models since the random term distributions cover symmetric, asymmetric and heavy-tailed distributions, such as skew-t, skew-slash and skew-contaminated normal, among others. Inspired by the EM algorithm proposed by Lachos et?al. (Statistics 44:541?C556, 2010c), we develop a local influence analysis for measurement error models, following Zhu and Lee??s (J R Stat Soc B 63:111?C126, 2001) approach. This is because the observed data log-likelihood function associated with the proposed model is somewhat complex and Cook??s well-known approach can be very difficult to apply to achieve local influence measures. Some useful perturbation schemes are also discussed. In addition, a score test for assessing the homogeneity of the skewness parameter vector is presented. Finally, the methodology is exemplified through a real data set, illustrating the usefulness of the proposed methodology.     

A robust multivariate Birnbaum–Saunders distribution: EM estimation

We propose here a robust multivariate extension of the bivariate Birnbaum–Saunders (BS) distribution derived by Kundu et al. [Bivariate Birnbaum–Saunders distribution and associated inference. J Multivariate Anal. 2010;101:113–125], based on scale mixtures of normal (SMN) distributions that are used for modelling symmetric data. This resulting multivariate BS-type distribution is an absolutely continuous distribution whose marginal and conditional distributions are of BS-type distribution of Balakrishnan et al. [Estimation in the Birnbaum–Saunders distribution based on scalemixture of normals and the EM algorithm. Stat Oper Res Trans. 2009;33:171–192]. Due to the complexity of the likelihood function, parameter estimation by direct maximization is very difficult to achieve. For this reason, we exploit the nice hierarchical representation of the proposed distribution to propose a fast and accurate EM algorithm for computing the maximum likelihood (ML) estimates of the model parameters. We then evaluate the finite-sample performance of the developed EM algorithm and the asymptotic properties of the ML estimates through empirical experiments. Finally, we illustrate the obtained results with a real data and display the robustness feature of the estimation procedure developed here.     

Comparison of spatial interpolation methods in the first order stationary multiplicative spatial autoregressive models

Three linear prediction methods of a single missing value for a stationary first order multiplicative spatial autoregressive model are proposed based on the quarter observations, observations in the first neighborhood, and observations in the nearest neighborhood. Three different types of innovations including Gaussian (symmetric and thin tailed), exponential (skew to right), and asymmetric Laplace (skew and heavy tailed) are considered. In each case, the proposed predictors are compared based on the two well-known criteria: mean square prediction and Pitman's measure of closeness. Parameter estimation is performed by maximum likelihood, least square errors, and Markov chain Monte Carlo (MCMC).     

Asymmetric generalizations of symmetric univariate probability distributions obtained through quantile splicing

Abstract

Balakrishnan et al. proposed a two-piece skew logistic distribution by making use of the cumulative distribution function (CDF) of half distributions as the building block, to give rise to an asymmetric family of two-piece distributions, through the inclusion of a single shape parameter. This paper proposes the construction of asymmetric families of two-piece distributions by making use of quantile functions of symmetric distributions as building blocks. This proposition will enable the derivation of a general formula for the L-moments of two-piece distributions. Examples will be presented, where the logistic, normal, Student’s t(2) and hyperbolic secant distributions are considered.     

Random number generation and estimation with the bimodal asymmetric power-normal distribution

Recently, Bolfarine et al. [Bimodal symmetric-asymmetric power-normal families. Commun Statist Theory Methods. Forthcoming. doi:10.1080/03610926.2013.765475] introduced a bimodal asymmetric model having the normal and skew normal as special cases. Here, we prove a stochastic representation for their bimodal asymmetric model and use it to generate random numbers from that model. It is shown how the resulting algorithm can be seen as an improvement over the rejection method. We also discuss practical and numerical aspects regarding the estimation of the model parameters by maximum likelihood under simple random sampling. We show that a unique stationary point of the likelihood equations exists except when all observations have the same sign. However, the location-scale extension of the model usually presents two or more roots and this fact is illustrated here. The standard maximization routines available in the R system (Broyden–Fletcher–Goldfarb–Shanno (BFGS), Trust, Nelder–Mead) were considered in our implementations but exhibited similar performance. We show the usefulness of inspecting profile loglikelihoods as a method to obtain starting values for maximization and illustrate data analysis with the location-scale model in the presence of multiple roots. A simple Bayesian model is discussed in the context of a data set which presents a flat likelihood in the direction of the skewness parameter.     

An asymmetric multivariate weibull distribution

Abstract

A class of multivariate laws as an extension of univariate Weibull distribution is presented. A well known representation of the asymmetric univariate Laplace distribution is used as the starting point. This new family of distributions exhibits some similarities to the multivariate normal distribution. Properties of this class of distributions are explored including moments, correlations, densities and simulation algorithms. The distribution is applied to model bivariate exchange rate data. The fit of the proposed model seems remarkably good. Parameters are estimated and a bootstrap study performed to assess the accuracy of the estimators.     

Sampling distribution after splitting experimental units

It has been a common practice to recommend one distribution for a large class of problems. An example is the use of the logarithmic distribution for count data - the point of concern in this paper is its recommendation without regard to the size of the experimental unit. References for this particular distribution go back to Fisher et al. [1943]. Look up Douglas [1980], Patil et al. [1984] and Kotz and Johnson [1985] to pick up additional references in this area -especially through sorting the data base of the American Mathematical Society. We will show a fallacy in this structure; provide a computer algorithm find the actual distributions; and then to check on the divisibility. The language used is APL2. But the users can make up their own programs.     

Robust mixture modeling using the skew <Emphasis Type="Italic">t</Emphasis> distribution

A finite mixture model using the Student's t distribution has been recognized as a robust extension of normal mixtures. Recently, a mixture of skew normal distributions has been found to be effective in the treatment of heterogeneous data involving asymmetric behaviors across subclasses. In this article, we propose a robust mixture framework based on the skew t distribution to efficiently deal with heavy-tailedness, extra skewness and multimodality in a wide range of settings. Statistical mixture modeling based on normal, Student's t and skew normal distributions can be viewed as special cases of the skew t mixture model. We present analytically simple EM-type algorithms for iteratively computing maximum likelihood estimates. The proposed methodology is illustrated by analyzing a real data example.     

The skew generalized t distribution as the scale mixture of a skew exponential power distribution and its applications in robust estimation

In this paper, we consider the family of skew generalized t (SGT) distributions originally introduced by Theodossiou [P. Theodossiou, Financial data and the skewed generalized t distribution, Manage. Sci. Part 1 44 (12) ( 1998), pp. 1650–1661] as a skew extension of the generalized t (GT) distribution. The SGT distribution family warrants special attention, because it encompasses distributions having both heavy tails and skewness, and many of the widely used distributions such as Student's t, normal, Hansen's skew t, exponential power, and skew exponential power (SEP) distributions are included as limiting or special cases in the SGT family. We show that the SGT distribution can be obtained as the scale mixture of the SEP and generalized gamma distributions. We investigate several properties of the SGT distribution and consider the maximum likelihood estimation of the location, scale, and skewness parameters under the assumption that the shape parameters are known. We show that if the shape parameters are estimated along with the location, scale, and skewness parameters, the influence function for the maximum likelihood estimators becomes unbounded. We obtain the necessary conditions to ensure the uniqueness of the maximum likelihood estimators for the location, scale, and skewness parameters, with known shape parameters. We provide a simple iterative re-weighting algorithm to compute the maximum likelihood estimates for the location, scale, and skewness parameters and show that this simple algorithm can be identified as an EM-type algorithm. We finally present two applications of the SGT distributions in robust estimation.     

Testing symmetry under a skew Laplace model

We develop tests of hypothesis about symmetry based on samples from possibly asymmetric Laplace distributions and present exact and limiting distribution of the test statistics. We postulate that the test statistic derived under the Laplace model is a rational choice as a measure of skewness and can be used in testing symmetry for other, quite general classes of skew distributions. Our results are applied to foreign exchange rates for 15 currencies.     

The multivariate hazard gradient and moments of the truncated multinormal distribution

It is well known that the expectation and variance of a truncated normal distribution can be simply expressed in terms of the hazard rate function. This paper shows that it is possible to express the expectation and covariance matrices of a truncated multinormal distribution with similarly simple expressions in which the hazard rate function is generalized to thevector multivariate hazard rate(also: hazard gradient) of Johnson and Kotz. This provides a concise computational form for the mutivariate moments and lends support to the contention that the hazard gradient is the appropriate generalization of the univariate hazard rate.     

Do maternal health problems influence child's worrying status? Evidence from the British Cohort Study

Conventional methods apply symmetric prior distributions such as a normal distribution or a Laplace distribution for regression coefficients, which may be suitable for median regression and exhibit no robustness to outliers. This work develops a quantile regression on linear panel data model without heterogeneity from a Bayesian point of view, i.e. upon a location-scale mixture representation of the asymmetric Laplace error distribution, and provides how the posterior distribution is summarized using Markov chain Monte Carlo methods. Applying this approach to the 1970 British Cohort Study (BCS) data, it finds that a different maternal health problem has different influence on child's worrying status at different quantiles. In addition, applying stochastic search variable selection for maternal health problems to the 1970 BCS data, it finds that maternal nervous breakdown, among the 25 maternal health problems, contributes most to influence the child's worrying status.     

Shannon Entropy and Mutual Information for Multivariate Skew‐Elliptical Distributions

Abstract. The entropy and mutual information index are important concepts developed by Shannon in the context of information theory. They have been widely studied in the case of the multivariate normal distribution. We first extend these tools to the full symmetric class of multivariate elliptical distributions and then to the more flexible families of multivariate skew‐elliptical distributions. We study in detail the cases of the multivariate skew‐normal and skew‐t distributions. We implement our findings to the application of the optimal design of an ozone monitoring station network in Santiago de Chile.