Evaluating the CDF of the Skew Normal Distribution

Total Page:16

File Type:pdf, Size:1020Kb

Evaluating the CDF of the Skew Normal Distribution Evaluating the CDF of the Skew Normal Distribution Christine Amsler, Michigan State University Alecos Papadopoulos, Athens University of Economics and Business ∗ Peter Schmidt, Michigan State University April 3, 2020 Final version Abstract In this paper we consider various methods of evaluating the cdf of the Skew Nor- mal distribution. This distribution arises in the stochastic frontier model because it is the distribution of the composed error, which is the sum (or difference) of a Normal and a Half Normal random variable. The cdf must be evaluated in models in which the composed error is linked to other errors using a Copula, in some methods of goodness of fit testing, or in the likelihood of models with sample selection bias. We investigate the accuracy of the evaluation of the cdf using expressions based on the bivariate Normal distribution, and also using simulation methods and some approx- imations. We find that the expressions based on the bivariate Normal distribution are quite accurate in the central portion of the distribution, and we propose several new approximations that are accurate in the extreme tails. By a simulated example we show that the use of approximations instead of the theoretical exact expressions may be critical in obtaining meaningful and valid estimation results. Keywords: skew normal distribution, bivariate normal distribution, stochastic frontier, sim- ulation, computational software. JEL classification: C46, C87, C88, C13. ∗Corresponding author, e-mail: [email protected], ORCID ID: 0000-0003-2441-4550. 1 1 Introduction This article deals with the evaluation of the cumulative distribution function of the Skew Normal distribution (hereafter SN-cdf). This is the distribution of the Normal / Half- Normal composed error in stochastic frontier (SF) models, and that is our main interest. However, the Skew Normal distribution also arises in a number of other contexts. For example, the maximum and the minimum of two correlated bivariate standard Normal random variables follow Skew Normal distributions (Loperfido, 2002). As another example, if X,Y are correlated bivariate Normal variables, then the distribution of X Y > 0 is Skew | Normal (Azzalini and Capitanio, 2014, p. 28). As a third example, in the Normal / Half Normal specification of the two-tier stochastic frontier model, the composed error is made up of three components and its density includes the difference of two Skew Normal cdfs (Papadopoulos, 2015). As a final more general example, the Skew Normal distribution is a useful way to model a moderately skewed regression error (which is often the case), and if one wants to account for sample selection bias, the SN-cdf will appear in the likelihood. Our work is an extension of Amsler et al. (2019) –hereafter AST– who addressed the same problem but who considered a much more limited set of methods of evaluation of the SN-cdf than we do in this paper. They proposed a simulation-based method that is accurate in the lower tail, and an approximation due to Tsay et al. (2013) that is accurate in the upper tail. For some uses, such as testing goodness of fit using the Kolmogorov-Smirnov test, the tails are not important. However, for other uses they are important. For example, suppose that a Skew Normal random variable is linked to other random variables using a Copula such as the Gaussian Copula. Then the SN-cdf is an argument of the inverse of the standard Normal cdf. This will explode to infinity (or minus infinity) if the argument is exactly one (or zero), causing the estimation algorithm to break down. So it is important that the evaluations of the SN-cdf are not exactly zero or one, even unreasonably far into the tails. The methods considered by AST achieve this, but it is not clear how accurate they are in the central portion of the distribution, which is presumably where most of the observations will lie. We will consider two other simulation-based methods, two other upper tail and two other lower tail approximations, and a variety of methods of evaluation of closed-form expressions 2 based on the cdf of the bivariate Normal distribution. The numerically evaluated closed- form expressions turn out to be quite accurate in the central portion of the distribution, where “central portion” is defined rather liberally as the range where the cdf value is greater than 10−20 and less than 1 10−20. Our lower tail approximations are also quite accurate − in the extreme lower tail (cdf value less than 10−20 but large enough to calculate, which, given standard floating-point double-precision arithmetic, means larger than about 10−310). Similarly, our upper tail approximations are quite accurate in the extreme upper tail (cdf value greater than 1 10−20, but small enough to calculate, which means smaller than − about 1 10−310). − The plan of the paper is as follows. Section 2 will define notation and present some preliminary results. Sections 3, 4 and 5 will present the various methods of evaluation of the SN-cdf that we will consider. Section 6 will present and discuss the results of our evaluations. In section 7 we provide a simulated example to show that accurate evaluation deep into the tails of the SN-cdf and the use of approximate expressions may be necessary to obtain meaningful and valid estimation results. Finally, Section 8 will give our conclusions. An Appendix and an on-line supplementary file complete the work. 2 Preliminaries We will follow the AST notation. We consider the Skew Normal distribution with location parameter equal to zero; scale parameter σ > 0; and shape parameter λ R. This ∈ distribution is characterized by the Skew Normal density snλ,σ (ε)=(2/σ)φ(ε/σ)Φ(λε/σ), where φ and Φ are the standard Normal pdf and cdf respectively. We want to calculate the Skew Normal cdf P (Q)= P (ε Q)= Q sn (ε) dε. λ,σ ≤ −∞ λ,σ More specifically, we will consider the caseR of σ = 1 and therefore investigate the accuracy of the evaluation of Pλ,1 (Q), which we will write more simply as Pλ(Q). This does not entail any loss of generality, because for σ = 1, one can use the fact that P (Q)= 6 λ,σ P (Q/σ). Similarly, we will consider only the case of λ 0. For λ < 0 we can use the λ ≥ so-called “reflection property” of the distribution, that P− (Q)=1 P ( Q). λ − λ − The connection of the Skew Normal distribution to the SF model is as follows. The “composed error” is ε = v + u, where v N (0, σ2), u N +(0, σ2), and v and u are ∼ v ∼ u 3 2 2 2 independent. Let σ = σu + σv and λ = σu/σv. Then ε has the Skew Normal density snλ,σ (ε). This discussion is for the case of ε = v + u, which would be natural in a cost frontier, and follows the discussion in Tsay et al. (2013) and AST. In the case of a production frontier, as in the original papers of Aigner, Lovell, and Schmidt (1977) and Meeusen and van den Broeck (1977), we would want to consider ε∗ = v u. Then ε∗ has − density sn (ε∗)=(2/σ)φ(ε∗/σ)Φ( λε∗/σ); that is, we just change the sign of λ. λ,σ − AST evaluated the cdf Pλ(Q) by a simulation based on the representation P (Q)= E (Φ[(Q u)/σ ]) . (1) λ u − v Here σ2 =1/ (1 + λ2), σ2 =1 σ2 = λ2/ (1 + λ2), u N(0, σ2)+, and E represents the v u − v ∼ u u expected value over the Half-Normal distribution of u. They estimated this by averaging Φ[(Q u)/σ ] over a large number of draws from the distribution of u. Specifically, they − v used 10,000,000 draws and got results that they regarded as reliable for negative values and small positive values of Q, but not for large positive values of Q. They created a large tabulation of the cdf as a function of λ and Q and suggested interpolation in this table. However, many people would probably prefer a simple calculation that is easier to program than an interpolation, which we will provide. AST also considered an approximation due to Tsay et al. (2013), p. 261, equations (12) and (13). AST concluded that this approximation is accurate for large positive values of Q but not for negative values.1 For the special case of λ = 1, AST also provided an “exact” result, namely P1 (Q) = [Φ (Q)]2.2 They were unaware that this result was already known.3 This result can be regarded as exact because it is not an approximation and it can be calculated very accurately over an extremely wide range of Q. It is useful as a check on the accuracy of the other methods of evaluating the SN-cdf.4 1Ashour and Abdel-hameed (2010) presented also a closed-form approximation to the Skew Normal density and CDF, but they were functions with branches while the partition of the support depended on the value of λ. They considered only the [-4,4] interval. 2This is the cdf of the maximum of two i.i.d. standard Normal variables, and can also be obtained as the limiting case of the results in Loperfido (2002) mentioned earlier. 3Azzalini (1985), p.174. This paper by Adelchi Azzalini has become by all accounts a classic of the statistical literature, and so, the Mark Twain remark on classic works applies. 4Gupta and Chen (2001) used it also as a validation tool for their SN-cdf tables.
Recommended publications
  • A Skew Extension of the T-Distribution, with Applications
    J. R. Statist. Soc. B (2003) 65, Part 1, pp. 159–174 A skew extension of the t-distribution, with applications M. C. Jones The Open University, Milton Keynes, UK and M. J. Faddy University of Birmingham, UK [Received March 2000. Final revision July 2002] Summary. A tractable skew t-distribution on the real line is proposed.This includes as a special case the symmetric t-distribution, and otherwise provides skew extensions thereof.The distribu- tion is potentially useful both for modelling data and in robustness studies. Properties of the new distribution are presented. Likelihood inference for the parameters of this skew t-distribution is developed. Application is made to two data modelling examples. Keywords: Beta distribution; Likelihood inference; Robustness; Skewness; Student’s t-distribution 1. Introduction Student’s t-distribution occurs frequently in statistics. Its usual derivation and use is as the sam- pling distribution of certain test statistics under normality, but increasingly the t-distribution is being used in both frequentist and Bayesian statistics as a heavy-tailed alternative to the nor- mal distribution when robustness to possible outliers is a concern. See Lange et al. (1989) and Gelman et al. (1995) and references therein. It will often be useful to consider a further alternative to the normal or t-distribution which is both heavy tailed and skew. To this end, we propose a family of distributions which includes the symmetric t-distributions as special cases, and also includes extensions of the t-distribution, still taking values on the whole real line, with non-zero skewness. Let a>0 and b>0be parameters.
    [Show full text]
  • A Study of Non-Central Skew T Distributions and Their Applications in Data Analysis and Change Point Detection
    A STUDY OF NON-CENTRAL SKEW T DISTRIBUTIONS AND THEIR APPLICATIONS IN DATA ANALYSIS AND CHANGE POINT DETECTION Abeer M. Hasan A Dissertation Submitted to the Graduate College of Bowling Green State University in partial fulfillment of the requirements for the degree of DOCTOR OF PHILOSOPHY August 2013 Committee: Arjun K. Gupta, Co-advisor Wei Ning, Advisor Mark Earley, Graduate Faculty Representative Junfeng Shang. Copyright c August 2013 Abeer M. Hasan All rights reserved iii ABSTRACT Arjun K. Gupta, Co-advisor Wei Ning, Advisor Over the past three decades there has been a growing interest in searching for distribution families that are suitable to analyze skewed data with excess kurtosis. The search started by numerous papers on the skew normal distribution. Multivariate t distributions started to catch attention shortly after the development of the multivariate skew normal distribution. Many researchers proposed alternative methods to generalize the univariate t distribution to the multivariate case. Recently, skew t distribution started to become popular in research. Skew t distributions provide more flexibility and better ability to accommodate long-tailed data than skew normal distributions. In this dissertation, a new non-central skew t distribution is studied and its theoretical properties are explored. Applications of the proposed non-central skew t distribution in data analysis and model comparisons are studied. An extension of our distribution to the multivariate case is presented and properties of the multivariate non-central skew t distri- bution are discussed. We also discuss the distribution of quadratic forms of the non-central skew t distribution. In the last chapter, the change point problem of the non-central skew t distribution is discussed under different settings.
    [Show full text]
  • A Family of Skew-Normal Distributions for Modeling Proportions and Rates with Zeros/Ones Excess
    S S symmetry Article A Family of Skew-Normal Distributions for Modeling Proportions and Rates with Zeros/Ones Excess Guillermo Martínez-Flórez 1, Víctor Leiva 2,* , Emilio Gómez-Déniz 3 and Carolina Marchant 4 1 Departamento de Matemáticas y Estadística, Facultad de Ciencias Básicas, Universidad de Córdoba, Montería 14014, Colombia; [email protected] 2 Escuela de Ingeniería Industrial, Pontificia Universidad Católica de Valparaíso, 2362807 Valparaíso, Chile 3 Facultad de Economía, Empresa y Turismo, Universidad de Las Palmas de Gran Canaria and TIDES Institute, 35001 Canarias, Spain; [email protected] 4 Facultad de Ciencias Básicas, Universidad Católica del Maule, 3466706 Talca, Chile; [email protected] * Correspondence: [email protected] or [email protected] Received: 30 June 2020; Accepted: 19 August 2020; Published: 1 September 2020 Abstract: In this paper, we consider skew-normal distributions for constructing new a distribution which allows us to model proportions and rates with zero/one inflation as an alternative to the inflated beta distributions. The new distribution is a mixture between a Bernoulli distribution for explaining the zero/one excess and a censored skew-normal distribution for the continuous variable. The maximum likelihood method is used for parameter estimation. Observed and expected Fisher information matrices are derived to conduct likelihood-based inference in this new type skew-normal distribution. Given the flexibility of the new distributions, we are able to show, in real data scenarios, the good performance of our proposal. Keywords: beta distribution; centered skew-normal distribution; maximum-likelihood methods; Monte Carlo simulations; proportions; R software; rates; zero/one inflated data 1.
    [Show full text]
  • Approximating the Distribution of the Product of Two Normally Distributed Random Variables
    S S symmetry Article Approximating the Distribution of the Product of Two Normally Distributed Random Variables Antonio Seijas-Macías 1,2 , Amílcar Oliveira 2,3 , Teresa A. Oliveira 2,3 and Víctor Leiva 4,* 1 Departamento de Economía, Universidade da Coruña, 15071 A Coruña, Spain; [email protected] 2 CEAUL, Faculdade de Ciências, Universidade de Lisboa, 1649-014 Lisboa, Portugal; [email protected] (A.O.); [email protected] (T.A.O.) 3 Departamento de Ciências e Tecnologia, Universidade Aberta, 1269-001 Lisboa, Portugal 4 Escuela de Ingeniería Industrial, Pontificia Universidad Católica de Valparaíso, Valparaíso 2362807, Chile * Correspondence: [email protected] or [email protected] Received: 21 June 2020; Accepted: 18 July 2020; Published: 22 July 2020 Abstract: The distribution of the product of two normally distributed random variables has been an open problem from the early years in the XXth century. First approaches tried to determinate the mathematical and statistical properties of the distribution of such a product using different types of functions. Recently, an improvement in computational techniques has performed new approaches for calculating related integrals by using numerical integration. Another approach is to adopt any other distribution to approximate the probability density function of this product. The skew-normal distribution is a generalization of the normal distribution which considers skewness making it flexible. In this work, we approximate the distribution of the product of two normally distributed random variables using a type of skew-normal distribution. The influence of the parameters of the two normal distributions on the approximation is explored. When one of the normally distributed variables has an inverse coefficient of variation greater than one, our approximation performs better than when both normally distributed variables have inverse coefficients of variation less than one.
    [Show full text]
  • Copulaesque Versions of the Skew-Normal Andskew-Student
    S S symmetry Article Copulaesque Versions of the Skew-Normal and Skew-Student Distributions Christopher Adcock Sheffield University Management School, University College Dublin, Sheffield S10 1FL, UK; c.j.adcock@sheffield.ac.uk Abstract: A recent paper presents an extension of the skew-normal distribution which is a copula. Under this model, the standardized marginal distributions are standard normal. The copula itself depends on the familiar skewing construction based on the normal distribution function. This paper is concerned with two topics. First, the paper presents a number of extensions of the skew-normal copula. Notably these include a case in which the standardized marginal distributions are Student’s t, with different degrees of freedom allowed for each margin. In this case the skewing function need not be the distribution function for Student’s t, but can depend on certain of the special functions. Secondly, several multivariate versions of the skew-normal copula model are presented. The paper contains several illustrative examples. Keywords: copula; multivariate distributions; skew-normal distributions; skew-Student distributions JEL Classification: C18; G01; G10; G12 Citation: Adcock, C.J. Copulaesque 1. Introduction Versions of the Skew-Normal and The recent paper by [1] presents an extension of the skew-normal distribution which Skew-Student Distributions. has subsequently been referred to as a copula. Under this model, the standardized marginal Symmetry 2021, 13, 815. https:// distributions are standard normal and some of the conditional distributions are skew- doi.org/10.3390/sym13050815 normal.The skew-normal distribution itself was introduced in two landmark papers [2] and [3]. These papers have led to a very substantial research effort by numerous authors Academic Editor: Nicola Maria over the last thirty five years.
    [Show full text]
  • Alpha-Skew Generalized T Distribution
    Revista Colombiana de Estadística July 2015, Volume 38, Issue 2, pp. 353 to 370 DOI: http://dx.doi.org/10.15446/rce.v38n2.51665 Alpha-Skew Generalized t Distribution Distribución t generalizada alfa sesgada Sukru Acitas1;a, Birdal Senoglu2;b, Olcay Arslan2;c 1Department of Statistics, Faculty of Science, Anadolu University, Eskisehir, Turkey 2Department of Statistics, Faculty of Science, Ankara University, Ankara, Turkey Abstract The alpha-skew normal (ASN) distribution has been proposed recently in the literature by using standard normal distribution and a skewing ap- proach. Although ASN distribution is able to model both skew and bimodal data, it is shortcoming when data has thinner or thicker tails than normal. Therefore, we propose an alpha-skew generalized t (ASGT) by using the gen- eralized t (GT) distribution and a new skewing procedure. From this point of view, ASGT can be seen as an alternative skew version of GT distribution. However, ASGT differs from the previous skew versions of GT distribution since it is able to model bimodal data sest as well as it nests most commonly used density functions. In this paper, moments and maximum likelihood estimation of the parameters of ASGT distribution are given. Skewness and kurtosis measures are derived based on the first four noncentral moments. The cumulative distribution function (cdf) of ASGT distribution is also ob- tained. In the application part of the study, two real life problems taken from the literature are modeled by using ASGT distribution. Key words: Bimodality, Kurtosis, Maximum Likelihood Estimation, Mod- eling, Skewness. Resumen La distribución normal alfa-sesgada (ASN por sus siglas en inglés) ha sido propuesta recientemente en la literatura mediante el uso de una distribución normal estándar y procedimientos de sesgo.
    [Show full text]
  • Alpha-Skew Generalized Normal Distribution and Its Applications
    Available at http://pvamu.edu/aam Applications and Applied Appl. Appl. Math. Mathematics: ISSN: 1932-9466 An International Journal (AAM) Vol. 14, Issue 2 (December 2019), pp. 784 - 804 Alpha-Skew Generalized Normal Distribution and its Applications 1*Eisa Mahmoudi, 2Hamideh Jafari and 3RahmatSadat Meshkat Department of Statistics Yazd University Yazd, Iran [email protected]; [email protected]; [email protected] * Corresponding Author Received: July 8, 2018; Accepted: November 11, 2019 Abstract The main object of this paper is to introduce a new family of distributions, which is quite flexible to fit both unimodal and bimodal shapes. This new family is entitled alpha-skew generalized normal (ASGN), that skews the symmetric distributions, especially generalized normal distribution through this paper. Here, some properties of this new distribution including cumulative distribution function, survival function, hazard rate function and moments are derived. To estimate the model parameters, the maximum likelihood estimators and the asymptotic distribution of the estimators are discussed. The observed information matrix is derived. Finally, the flexibility of the new distribution, as well as its effectiveness in comparison with other distributions, are shown via an application. Keywords: Alpha-skew distribution; Generalized normal distribution; Maximum likelihood estimator; Skew-normal distribution MSC 2010 No.: 60E05, 62F10 1. Introduction Recently, the skew-symmetric distributions have received considerable amount of attention, for example, the alpha skew normal (ASN) distribution is studied by Elal-Olivero (2010). First, he introduced a new symmetric bimodal-normal distribution with the pdf 784 AAM: Intern. J., Vol. 14, Issue 2 (December 2019) 785 f( y ) y2 ( y ), where is standard normal density function, and then defined alpha-skew-normal distribution as (1 y )2 1 f( y ; ) ( y ), y ¡¡, , 2 2 where is skew parameter.
    [Show full text]
  • Characterization of the Skew-Normal Distribution
    Ann. Inst. Statist. Math. Vol. 56, No. 2, 351-360 (2004) (~)2004 The Institute of Statistical Mathematics CHARACTERIZATION OF THE SKEW-NORMAL DISTRIBUTION ARJUN K. GUPTA 1, TRUC T. NGUYEN I AND JOSE ALMER T. SANQUl 2. 1Department of Mathematics and Statistics, Bowling Green State University, Bowling Green, OH ~3~03-0221, U.S.A. 2Department of Mathematical Sciences, Appalachian State University, ASU Box 3209P, Boone, NC 28608-2092, U.S.A. (Received January 25, 2002; revised July 23, 2003) Abstract. Two characterization results for the skew-normal distribution based on quadratic statistics have been obtained. The results specialize to known character- izations of the standard normal distribution and generalize to the characterizations of members of a larger family of distributions. Results on the decomposition of the family of distributions of random variables whose square is distributed as X~ are obtained. Key words and phrases: Non-normal distribution, chi-square distribution, half- normal distribution, skew-symmetric distribution, sequence of moments, induction, decomposition, characteristic function. 1. Introduction A random variable Z has a skew-normal distribution with parameter A, denoted by Z ~-- SN(A), if its density is given by f(z, A) = 20(Az)r where and r are the standard normal cumulative distribution function and the standard normal probability density function, respectively, and z and A are real numbers (Azzalini (1985)). Some basic properties of the SN(A) distribution given in Azzalini (1985) are: 1. SN(O) = N(0, 1); 2. If Z ~ SN(A) then -Z ~ SN(-A); 3. As )~ -~ +c~, SN(A) tends to the half-normM distribution, i.e., the distribution of • I when X ~ g(0, 1); and 4.
    [Show full text]
  • Truncated Skew-Normal Distributions: Moments, Estimation by Weighted Moments and Application to Climatic Data
    Truncated skew-normal distributions: moments, estimation by weighted moments and application to climatic data by C. Flecher1;2, D. Allard1 and P. Naveau2 1 Biostatistics and Spatial Processes, INRA, Agroparc 84914 Avignon, France 2 Laboratoire des Sciences du Climat et de l'Environnement, CNRS Gif-sur-Yvette, France Running title: Truncated skew-normal distributions Abstract In this paper we derive an expression of the mth order moments and some weighted mo- ments of truncated skew-normal distributions. We link these formulas to previous results for trun- cated distributions and non truncated skew-normal distributions. Methods to estimate skew-normal parameters using weighted moments are proposed and compared to other classical techniques. In a second step we propose to model the distribution of relative humidity with a truncated skew-normal distribution. Key words: Truncated skew-normal distributions, inference methods, classical moments, weighted moments, relative humidity. 1 1 Introduction In many applications, the probability distribution function (pdf) of some observed variables can be simultaneously skewed and restricted to a xed interval. For example, variables such as pH, grades, and humidity in environmental studies, have upper and lower physical bounds and their pdfs are not necessarily symmetric within these bounds. To illustrate such behaviors, Figure 1 shows the histograms per season and the estimated pdfs of daily relative humidity measurements made in Toulouse (France) from 1972 to 1999. Dec−Jan−Feb Mar−Apr−May Jun−Jul−Aug Sep−Oct−Nov 0.06 0.06 0.06 0.06 0.05 0.05 0.05 0.05 0.04 0.04 0.04 0.04 0.03 0.03 0.03 0.03 Density Density Density Density 0.02 0.02 0.02 0.02 0.01 0.01 0.01 0.01 0.00 0.00 0.00 0.00 0 20 40 60 80 100 0 20 40 60 80 100 0 20 40 60 80 100 0 20 40 60 80 100 H (%) H (%) H (%) H (%) Figure 1: Histograms and estimated densities of daily relative humidity measurements made in Toulouse (France) from 1972 to 1999.
    [Show full text]
  • The Skew-Normal Distribution in Spc
    REVSTAT – Statistical Journal Volume 11, Number 1, March 2013, 83–104 THE SKEW-NORMAL DISTRIBUTION IN SPC Authors: Fernanda Figueiredo – CEAUL and Faculdade de Economia da Universidade do Porto, Portugal [email protected] M. Ivette Gomes – Universidade de Lisboa, FCUL, DEIO and CEAUL, Portugal [email protected] Abstract: Modeling real data sets, even when we have some potential (as)symmetric models for • the underlying data distribution, is always a very difficult task due to some uncon- trollable perturbation factors. The analysis of different data sets from diverse areas of application, and in particular from statistical process control (SPC), leads us to notice that they usually exhibit moderate to strong asymmetry as well as light to heavy tails, which leads us to conclude that in most of the cases, fitting a normal distribution to the data is not the best option, despite of the simplicity and popu- larity of the Gaussian distribution. In this paper we consider a class of skew-normal models that include the normal distribution as a particular member. Some properties of the distributions belonging to this class are enhanced in order to motivate their use in applications. To monitor industrial processes some control charts for skew-normal and bivariate normal processes are developed, and their performance analyzed. An application with a real data set from a cork stopper’s process production is presented. Key-Words: bootstrap control charts; false alarm rate; heavy-tails; Monte Carlo simulations; • probability limits; run-length; shewhart control charts; skewness; skew-normal dis- tribution; statistical process control. AMS Subject Classification: 62G05, 62G35, 62P30, 65C05.
    [Show full text]
  • Study of Unified Multivariate Skew Normal Distribution with Applications in Finance and Actuarial Science
    STUDY OF UNIFIED MULTIVARIATE SKEW NORMAL DISTRIBUTION WITH APPLICATIONS IN FINANCE AND ACTUARIAL SCIENCE Mohammad Abdus Samad Aziz A Dissertation Submitted to the Graduate College of Bowling Green State University in partial fulfillment of the requirements for the degree of DOCTOR OF PHILOSOPHY August 2011 Committee: Arjun K. Gupta, Advisor Arthur Yeh, Graduate Faculty Representative John T. Chen Wei Ning ii ABSTRACT Arjun K. Gupta, Advisor The classical work horse in finance and insurance for modeling asset returns is the Gaus- sian model. However, when modeling complex random phenomena, more flexible distribu- tions are needed which are beyond the normal distribution. This is because most of the financial and economic data are skewed and have \fat tails" due to the presence of outliers. Hence symmetric distributions like normal or others may not be good choices while modeling these kinds of data. Flexible distributions like skew normal distribution allow robust model- ing of high-dimensional multimodal and asymmetric data. In this dissertation, we consider a very flexible financial model to construct robust comonotonic lower convex order bounds in approximating the distribution of the sums of dependent log skew normal random vari- ables. The dependence structure of these random variables is based on a recently developed multivariate skew normal distribution, called unified skew normal distribution. In order to accommodate the distribution to the model so considered, we first study inherent properties of this class of skew normal distribution. These properties along with the convex order and comonotonicity of random variables are then used to approximate the distribution function of terminal wealth. By calculating lower bounds in the convex order sense, we consider approximations that reduce the multivariate randomness to univariate randomness.
    [Show full text]
  • The Generalized-Alpha-Beta-Skew-Normal Distribution: Properties and Applications
    The Generalized-Alpha-Beta-Skew-Normal Distribution: Properties and Applications Sricharan Shah1, Subrata Chakraborty1, Partha Jyoti Hazarika1 and M. MasoomAli2 1Department of Statistics, Dibrugarh University, Dibrugarh, Assam 786004 India 2Department of Mathematical Sciences, Ball State University, Muncie, IN 47306 USA Abstract In this paper we have introduced a generalized version of alpha beta skew normal distribution in the same line of Sharafi et al. (2017) and investigated some of its basic properties. The extensions of the proposed distribution have also been studied. The appropriateness of the proposed distribution has been tested by comparing the values of Akaike Information Criterion (AIC) and Bayesian Information Criterion (BIC) with the values of some other known related distributions for better model selection. Likelihood ratio testhas been used for discriminating between nested models. Keywords: Skew Distribution, Alpha-Skew Distribution, Bimodal Distribution, AIC, BIC. Math Subject Classification: 60E05, 62H10, 62H12 1. Introduction In many real life situations the data exhibit many modes as well as asymmetry, e.g., in the fields of demography, insurance, medical sciences, physics, etc. (for details see Hammel et al. (1983), Dimitrov et al. (1997), Sinha (2012), Chakraborty and Hazarika (2011), Chakraborty et al. (2015), among others). Azzalini (1985) first introduced the skew normal distribution and the density function of this distribution, denoted by SN () , is given by fZ (z;) 2(z)(z); z , (1) where, (.) is the density function of standard normal distribution , (.) is the cumulative distribution function (cdf) of standard normal distribution and is the asymmetry parameter. The generalization of the skew normal distribution was proposed by Balakrishnan (2002) as a discussant in Arnold and Beaver (2002) and studied its properties and the density function is given by n f Z (z;, n) (z)[(z)] Cn () (2) 1 n where, n , is a positive integer and Cn () E[ (U)], U ~ N(0,1) .
    [Show full text]