A TEST OF FIT FOR THE GENERALIZED PARETO DISTRIBUTION BASED ON TRANSFORMS
|
|
- Cora Baker
- 6 years ago
- Views:
Transcription
1 A TEST OF FIT FOR THE GENERALIZED PARETO DISTRIBUTION BASED ON TRANSFORMS Dimitrios Konstantinides, Simos G. Meintanis Department of Statistics and Acturial Science, University of the Aegean, Karlovassi, Samos, Greece and Department of Economics, National and Kapodistrian University of Athens, 8 Pesmazoglou Street, Athens, Greece Abstract. The generalized Pareto distribution is a very popular two parameter model for extreme events. In this article we develop a class of goodness of fit tests for the composite null hypothesis that the data at hand come from a Generalized Pareto distribution with unspecified parameter values. In doing so, the parameters are estimated, and the resulting estimates are employed in transforming the data, so that under the null hypothesis the transformed data approximately follow a unit exponential distribution. Then the null hypothesis of exponentiality is tested instead, by utilizing the empirical Laplace transform. The method is shown to be consistent, and the asymptotic null distribution of the test statistic is derived. The results of a Monte Carlo study are presented which include, apart from the proposed test, the standard methods based on the empirical distribution function, implemented by employing moment and maximum likelihood estimates, as well as probability weighted moments estimates. Keywords. Exterme events, Goodness of fit test, Empirical Laplace transform. AMS 2000 classification numbers: 62G10, 62G20 1
2 1 Introduction The analysis of series of observations consisting of the largest (or smallest) values was tradinionally based on the family of generalized extreme value distributions. However this approach has recently been critisized mainly because of its inefficient use of available data. One approach towards incorporating this loss of information is to consider several of the larger order statistics. Then the resulting exceedances (differences between these order statistics and a given high threshold) are typically modelled by the generalized Pareto distribution (GPD). The GPD is a two parameter model with a shape parameter α and a scale parameter c. Its distribution function is F (x; α, c) = 1 (1 (αx/c)) 1/α, α IR, c > 0, and has support x > 0 ( resp. 0 < x < c/α), if a 0 (resp. a > 0). We write GP(α, c) to denote the GPD with parameters α and c. For α < 0, the distribution is related to the P areto distribution of the second kind (see Johnson et al., 1994, p. 575), whereas for α = 0 and α = 1, respectively, the exponential with scale c and the uniform in (0,c) result. For applications of the GPD to such diverse areas of applied research as metereology, economics, ecology and reliability, the reader is referred to Choulakian and Stephens (2001), Singh and Ahmad (2004), and references therein. Recently, Choulakian and Stephens (2001) proposed classical goodness of fit tests for the null hypothesis H 0 : The law of X is GP(α, c) for some α IR and c > 0, based on independent copies {X j } n of the random variable X 0. These tests, namely the Kolmogorov Smirnov, the Crámer von Mises and the Anderson Darling test, utilize the empirical distribution function (EDF). In this article we develop a class of goodness of fit tests for the GPD based on the empirical Laplace transform. To this end consider, instead of X, the random variable Y = (1/α) log(1 (αx/c)), which under H 0 follows a unit exponential distribution. Consequently the Laplace transform (LT) L(t) = E(e ty ) of Y satisfies (1.1) (1 + t)l(t) 1 = 0, for all t > 0. 2
3 Hence we may test H 0, by employing the data (1.2) Y j = (1/ˆα) log(1 (ˆαX j /ĉ)), j = 1, 2,..., n, where ˆα (resp. ĉ) denotes a consistent estimator of the parameter α (resp. c). Under H 0 and for large n, {Y j } n will approximately follow a unit exponential distribution. Therefore it would be natural (in view of (1.1)) to base a test for H 0 on a measure of deviation from zero of the random function D n (t) = (1 + t)l n (t) 1 on [0, ], where L n (t) = 1 n exp( ty j ), is the empirical LT of the transformed data {Y j } n. Such a test for exponentiality was developed by Henze and Meintanis (2002) and takes the form (1.3) T n,a = n 0 D 2 n(t)e at dt, a > 0. Since it was shown that the test that rejects exponentiality for large values of T n,a, is more powerful than the classical tests referred to above, it is hoped that this will be so in the present, more general, situation. It should be emphasized that the empirical LT has proved a valuable tool for statistical inference, not only in the case of testing exponentiality, but also in goodness of fit tests for the inverse Gaussian (Henze and Klar, 2002) and the gamma distribution (Henze and Meintanis, 2004), as well as in the estimation context (Csörgő and Teugels, 1990, Feuerverger, 1989), to name a few. The paper is organized as follows. In Section 2 we derive the limit null distribution of T n,a and establish the consistency of the test based on T n,a against general alternatives. The results of a Monte Carlo study are presented in Section 3, where the new tests are compared to the EDF goodness of fit tests for the GPD law. The paper concludes with real data examples given in Section 4. 2 Theoretical results We briefly review the setting for asymptotic distribution theory. For more details the reader is referred to Henze and Meintanis (2004). A convenient setting is the separable Hilbert space H = L 2 (IR +, B, e at dt) of (equivalence classes of) measurable functions 3
4 f : IR + IR + satisfying 0 f 2 (t)e at dt <. The inner product and the norm in H are denoted by < f, g > and f, respectively. Here and in what follows, the notation D P means convergence in distribution of random elements and random variables, means convergence in probability, o P (1) stands for convergence in probability to 0, O P (1) denotes boundedness in probability, and i.i.d. means independent and identically distributed. In view of the parametric bootstrap procedure, which will be carried out to obtain critical values for T n,a, we study the null asymptotic distribution of T n,a under a triangular array X n1, X n2,..., X nn, n 1, of rowwise i.i.d. random variables, where X n1 GP(α n, 1), α n IR, and lim n α n = α. For the bootstrap procedure, α n = ˆα n is chosen. We further consider only regular estimators ˆα n and ĉ n of α and c, respectively. Namely we assume that under the triangular array referred to above, n(ˆαn α n ) = 1 n n(ĉn 1) = 1 n ψ 1 (X nj ; α n ) + ɛ n,1, ψ 2 (X nj ; α n ) + ɛ n,2, where for k = 1, 2, ɛ n,k = o p (1) and ψ k has zero mean, and finite second moment which satisfies lim n Eψk 2(X n1; α n ) = Eψk 2 (X; α), with X GP(α, 1). 3 Finite sample comparisons This section presents the results of a Monte Carlo study conducted at a 10% nominal level with replications to assess the performance of the new tests. In order to avoid reliance on asymptotic critical values and since the null distribution of all test statistics considered depends on the (unknown) value of the shape parameter α we performed a parametric bootstrap to obtain the critical point p n of the test as follows: Conditionally on the observed value of ˆα n, generate 100 bootstrap samples from GP(ˆα n, 1). Calculate the value of the test statistic, say T j, (j = 1, 2,..., 100), for each bootstrap sample. Obtain p n as T(90), where T (j), j = 1, 2,..., 100 denotes the ordered values T j values. We have adapted the choice in Gürtler and Henze (2000) 4
5 and used the modified critical point p n = T(90) (T (91) T (90) ), which leads to a more accurate empirical level of the test. i) Methods of estimation: We consider the moment (MO) estimators where ˆα n = 1 2 [ ( ) 2 [ ( Xn 1], ĉ n = 1 ) 2 S n 2 X Xn n + 1], S n X n = 1 n X j, S 2 n = 1 n (X j X n ) 2. The MO estimators are regular, in the sense of (2.3), for α > 1/4 (Hosking and Wallis, 1987). Hosking and Wallis (1987) also consider a related class of estimators, which are regular for α > 1/2. These are termed probability weighted moments (PWM) estimators and may be written as ˆα n = X n X n 2w 2, ĉ n = 2w X n X n 2w, where w = 1 n (1 π j )X (j), π j = j 0.35 n, j = 1, 2,..., n, and X (1) X (2)... X (n) are the order statistics of {X j } n. The maximum likelihood (ML) estimators are regular for α < 1/2 (Smith, 1984). We have employed the routine of Grimshaw (1993) in obtaining the ML estimators, supplemented by a bisection routine to locate the maximum of the profile log likelihood, [ ] L(ϑ) = n log(1 ϑx j ) n log 1 log(1 ϑx j ), nϑ with respect to ϑ = α/c. Let ˆϑ n denote the maximizing value. Then the ML estimators of α and c, are ˆα n = 1 n log(1 ˆϑ n X j ), ĉ n = ˆα n ˆϑ n. 5
6 Comparative evaluation for the performance of estimators of the GPD, including the estimators considered here, may be found in Peng and Welsh (2001) and Singh and Ahmad (2004). ii) Test statistics: From (1.3) it follows by straightforward algebra that the new test statistic may conviniently be written as T n,a = 1 n j,k=1,n 1 + (Y j + Y k + a + 1) 2 (Y j + Y k + a) 3 2 with Y j, j = 1, 2,..., n defined in (1.2). 1 + Y j + a (Y j + a) 2, The new test is compared with the EDF procedures developed by Choulakian and Stephens (2001). Let ˆF (x) = F (x; ˆα n, ĉ n ). Then the Kolmogorov Smirnov (KS) statistic is where KS = max{d +, D }, D + = max { j,2,...,n n ˆF (X (j) )}, D = max { ˆF (X (j) j 1,2,...,n n )}. The Crámer von Mises (CM) statistic is CM = 1 ( 12n + ˆF (X (j) ) 2j 1 ) 2, 2n and the Anderson Darling (AD) statistic is AD = n 1 n [ (2j 1) log ˆF (X (j) ) + (2(n j) + 1) log(1 ˆF ] (X (j) )). iii) Simulation results: All calculations were done at the Department of Economics, University of Athens, using double precision arithmetic in FORTRAN and routines from the IMSL library, whenever available. Apart from the GPD, the gamma with density Γ(θ) 1 x θ 1 exp ( x), the Weibull with density θx θ 1 exp( x θ ), and the log normal distribution with density (θx) 1 (2π) 1/2 exp[ (log x) 2 /(2θ 2 )], are considered. These distributions will be denoted by Γ(θ), W (θ) and LN(θ), respectively, while the GPD with unit scale will simply be denoted by GP (α). 6
7 Table 1: Percentage of rejection for Monte Carlo samples of size n = 50 (left part) n = 75 (middle part) and n = 100 (right part). Estimation: MO (top), PWM (middle) and ML (bottom). T 0.25 T 0.5 T 0.75 T 0.25 T 0.50 T 0.75 T 0.25 T 0.50 T 0.75 GP ( 0.20) GP ( 0.10) GP (0.0) GP (0.20) GP (0.50) GP (1.0) Γ(2.0) W (0.75) W (1.25) W (1.50) LN(1.0) GP ( 0.40) GP ( 0.20) GP (0.0) GP (0.20) GP (0.50) GP (1.0) Γ(2.0) W (0.75) W (1.25) W (1.50) LN(1.0) GP ( 1.0) GP ( 0.75) GP ( 0.50) GP ( 0.25) GP (0.0) GP (0.20) GP (0.40) Γ(2.0) W (0.75) W (1.25) W (1.50) LN(1.0)
8 Table 2: Percentage of rejection for Monte Carlo samples of size n = 50 (left part) n = 75 (middle part) and n = 100 (right part). Estimation: MO (top), PWM (middle) and ML (bottom). CM KS AD CM KS AD CM KS AD GP ( 0.20) GP ( 0.10) GP (0.0) GP (0.20) GP (0.50) GP (1.0) Γ(2.0) W (0.75) W (1.25) W (1.50) LN(1.0) GP ( 0.40) GP ( 0.20) GP (0.0) GP (0.20) GP (0.50) GP (1.0) Γ(2.0) W (0.75) W (1.25) W (1.50) LN(1.0) GP ( 1.0) GP ( 0.75) GP ( 0.50) GP ( 0.25) GP (0.0) GP (0.20) GP (0.40) Γ(2.0) W (0.75) W (1.25) W (1.50) LN(1.0)
9 Table 1 shows results (percentage of rejection rounded to the nearest integer) for T n,a, a = 0.25, 0.50, In the table we simply write T a instead of T n,a. Corresponding results for the KS, CM, and AD tests are shown in Table 2. 9
Fitting the generalized Pareto distribution to data using maximum goodness-of-fit estimators
Computational Statistics & Data Analysis 51 (26) 94 917 www.elsevier.com/locate/csda Fitting the generalized Pareto distribution to data using maximum goodness-of-fit estimators Alberto Luceño E.T.S. de
More informationGOODNESS OF FIT TESTS FOR PARAMETRIC REGRESSION MODELS BASED ON EMPIRICAL CHARACTERISTIC FUNCTIONS
K Y B E R N E T I K A V O L U M E 4 5 2 0 0 9, N U M B E R 6, P A G E S 9 6 0 9 7 1 GOODNESS OF FIT TESTS FOR PARAMETRIC REGRESSION MODELS BASED ON EMPIRICAL CHARACTERISTIC FUNCTIONS Marie Hušková and
More informationStatistic Distribution Models for Some Nonparametric Goodness-of-Fit Tests in Testing Composite Hypotheses
Communications in Statistics - Theory and Methods ISSN: 36-926 (Print) 532-45X (Online) Journal homepage: http://www.tandfonline.com/loi/lsta2 Statistic Distribution Models for Some Nonparametric Goodness-of-Fit
More informationENTROPY-BASED GOODNESS OF FIT TEST FOR A COMPOSITE HYPOTHESIS
Bull. Korean Math. Soc. 53 (2016), No. 2, pp. 351 363 http://dx.doi.org/10.4134/bkms.2016.53.2.351 ENTROPY-BASED GOODNESS OF FIT TEST FOR A COMPOSITE HYPOTHESIS Sangyeol Lee Abstract. In this paper, we
More informationTESTING FOR EQUAL DISTRIBUTIONS IN HIGH DIMENSION
TESTING FOR EQUAL DISTRIBUTIONS IN HIGH DIMENSION Gábor J. Székely Bowling Green State University Maria L. Rizzo Ohio University October 30, 2004 Abstract We propose a new nonparametric test for equality
More informationGoodness-of-fit tests for randomly censored Weibull distributions with estimated parameters
Communications for Statistical Applications and Methods 2017, Vol. 24, No. 5, 519 531 https://doi.org/10.5351/csam.2017.24.5.519 Print ISSN 2287-7843 / Online ISSN 2383-4757 Goodness-of-fit tests for randomly
More informationTesting Exponentiality by comparing the Empirical Distribution Function of the Normalized Spacings with that of the Original Data
Testing Exponentiality by comparing the Empirical Distribution Function of the Normalized Spacings with that of the Original Data S.Rao Jammalamadaka Department of Statistics and Applied Probability University
More informationJournal of Environmental Statistics
jes Journal of Environmental Statistics February 2010, Volume 1, Issue 3. http://www.jenvstat.org Exponentiated Gumbel Distribution for Estimation of Return Levels of Significant Wave Height Klara Persson
More informationarxiv: v1 [stat.me] 25 May 2018
Body and Tail arxiv:1805.10040v1 [stat.me] 5 May 018 Separating the distribution function by an efficient tail-detecting procedure in risk management Abstract Ingo Hoffmann a,, Christoph J. Börner a a
More informationRecall the Basics of Hypothesis Testing
Recall the Basics of Hypothesis Testing The level of significance α, (size of test) is defined as the probability of X falling in w (rejecting H 0 ) when H 0 is true: P(X w H 0 ) = α. H 0 TRUE H 1 TRUE
More informationTest of fit for a Laplace distribution against heavier tailed alternatives
DEPARTMENT OF STATISTICS AND ACTUARIAL SCIENCE University of Waterloo, 200 University Avenue West Waterloo, Ontario, Canada, N2L 3G1 519-888-4567, ext. 00000 Fax: 519-746-1875 www.stats.uwaterloo.ca UNIVERSITY
More informationInvestigation of goodness-of-fit test statistic distributions by random censored samples
d samples Investigation of goodness-of-fit test statistic distributions by random censored samples Novosibirsk State Technical University November 22, 2010 d samples Outline 1 Nonparametric goodness-of-fit
More informationEstimation of Quantiles
9 Estimation of Quantiles The notion of quantiles was introduced in Section 3.2: recall that a quantile x α for an r.v. X is a constant such that P(X x α )=1 α. (9.1) In this chapter we examine quantiles
More informationTesting Goodness-of-Fit for Exponential Distribution Based on Cumulative Residual Entropy
This article was downloaded by: [Ferdowsi University] On: 16 April 212, At: 4:53 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 172954 Registered office: Mortimer
More informationPackage EWGoF. October 5, 2017
Type Package Package EWGoF October 5, 2017 Title Goodness-of-Fit Tests for the Eponential and Two-Parameter Weibull Distributions Version 2.2.1 Author Meryam Krit Maintainer Meryam Krit
More informationAPPROXIMATING THE GENERALIZED BURR-GAMMA WITH A GENERALIZED PARETO-TYPE OF DISTRIBUTION A. VERSTER AND D.J. DE WAAL ABSTRACT
APPROXIMATING THE GENERALIZED BURR-GAMMA WITH A GENERALIZED PARETO-TYPE OF DISTRIBUTION A. VERSTER AND D.J. DE WAAL ABSTRACT In this paper the Generalized Burr-Gamma (GBG) distribution is considered to
More informationOn robust and efficient estimation of the center of. Symmetry.
On robust and efficient estimation of the center of symmetry Howard D. Bondell Department of Statistics, North Carolina State University Raleigh, NC 27695-8203, U.S.A (email: bondell@stat.ncsu.edu) Abstract
More informationResearch Article Tests of Fit for the Logarithmic Distribution
Hindawi Publishing Corporation Journal of Applied Mathematics and Decision Sciences Volume 2008 Article ID 463781 8 pages doi:10.1155/2008/463781 Research Article Tests of Fit for the Logarithmic Distribution
More informationAsymptotic Statistics-VI. Changliang Zou
Asymptotic Statistics-VI Changliang Zou Kolmogorov-Smirnov distance Example (Kolmogorov-Smirnov confidence intervals) We know given α (0, 1), there is a well-defined d = d α,n such that, for any continuous
More informationRecall that in order to prove Theorem 8.8, we argued that under certain regularity conditions, the following facts are true under H 0 : 1 n
Chapter 9 Hypothesis Testing 9.1 Wald, Rao, and Likelihood Ratio Tests Suppose we wish to test H 0 : θ = θ 0 against H 1 : θ θ 0. The likelihood-based results of Chapter 8 give rise to several possible
More informationOne-Sample Numerical Data
One-Sample Numerical Data quantiles, boxplot, histogram, bootstrap confidence intervals, goodness-of-fit tests University of California, San Diego Instructor: Ery Arias-Castro http://math.ucsd.edu/~eariasca/teaching.html
More informationFinancial Econometrics and Volatility Models Extreme Value Theory
Financial Econometrics and Volatility Models Extreme Value Theory Eric Zivot May 3, 2010 1 Lecture Outline Modeling Maxima and Worst Cases The Generalized Extreme Value Distribution Modeling Extremes Over
More informationPackage homtest. February 20, 2015
Version 1.0-5 Date 2009-03-26 Package homtest February 20, 2015 Title Homogeneity tests for Regional Frequency Analysis Author Alberto Viglione Maintainer Alberto Viglione
More informationComputer simulation on homogeneity testing for weighted data sets used in HEP
Computer simulation on homogeneity testing for weighted data sets used in HEP Petr Bouř and Václav Kůs Department of Mathematics, Faculty of Nuclear Sciences and Physical Engineering, Czech Technical University
More informationDoes k-th Moment Exist?
Does k-th Moment Exist? Hitomi, K. 1 and Y. Nishiyama 2 1 Kyoto Institute of Technology, Japan 2 Institute of Economic Research, Kyoto University, Japan Email: hitomi@kit.ac.jp Keywords: Existence of moments,
More informationThe Goodness-of-fit Test for Gumbel Distribution: A Comparative Study
MATEMATIKA, 2012, Volume 28, Number 1, 35 48 c Department of Mathematics, UTM. The Goodness-of-fit Test for Gumbel Distribution: A Comparative Study 1 Nahdiya Zainal Abidin, 2 Mohd Bakri Adam and 3 Habshah
More informationFall 2012 Analysis of Experimental Measurements B. Eisenstein/rev. S. Errede
Hypothesis Testing: Suppose we have two or (in general) more simple hypotheses which can describe a set of data Simple means explicitly defined, so if parameters have to be fitted, that has already been
More informationExact Statistical Inference in. Parametric Models
Exact Statistical Inference in Parametric Models Audun Sektnan December 2016 Specialization Project Department of Mathematical Sciences Norwegian University of Science and Technology Supervisor: Professor
More informationA Very Brief Summary of Statistical Inference, and Examples
A Very Brief Summary of Statistical Inference, and Examples Trinity Term 2009 Prof. Gesine Reinert Our standard situation is that we have data x = x 1, x 2,..., x n, which we view as realisations of random
More informationA Very Brief Summary of Statistical Inference, and Examples
A Very Brief Summary of Statistical Inference, and Examples Trinity Term 2008 Prof. Gesine Reinert 1 Data x = x 1, x 2,..., x n, realisations of random variables X 1, X 2,..., X n with distribution (model)
More informationTESTING PROCEDURES BASED ON THE EMPIRICAL CHARACTERISTIC FUNCTIONS I: GOODNESS-OF-FIT, TESTING FOR SYMMETRY AND INDEPENDENCE. 1.
Tatra Mt. Math. Publ. 39 (2008), 225 233 t m Mathematical Publications TESTING PROCEDURES BASED ON THE EMPIRICAL CHARACTERISTIC FUNCTIONS I: GOODNESS-OF-FIT, TESTING FOR SYMMETRY AND INDEPENDENCE Marie
More informationApplication of Homogeneity Tests: Problems and Solution
Application of Homogeneity Tests: Problems and Solution Boris Yu. Lemeshko (B), Irina V. Veretelnikova, Stanislav B. Lemeshko, and Alena Yu. Novikova Novosibirsk State Technical University, Novosibirsk,
More informationParameter Estimation
Parameter Estimation Chapters 13-15 Stat 477 - Loss Models Chapters 13-15 (Stat 477) Parameter Estimation Brian Hartman - BYU 1 / 23 Methods for parameter estimation Methods for parameter estimation Methods
More informationLecture 8 Inequality Testing and Moment Inequality Models
Lecture 8 Inequality Testing and Moment Inequality Models Inequality Testing In the previous lecture, we discussed how to test the nonlinear hypothesis H 0 : h(θ 0 ) 0 when the sample information comes
More informationSpline Density Estimation and Inference with Model-Based Penalities
Spline Density Estimation and Inference with Model-Based Penalities December 7, 016 Abstract In this paper we propose model-based penalties for smoothing spline density estimation and inference. These
More informationASTIN Colloquium 1-4 October 2012, Mexico City
ASTIN Colloquium 1-4 October 2012, Mexico City Modelling and calibration for non-life underwriting risk: from empirical data to risk capital evaluation Gian Paolo Clemente Dept. Mathematics, Finance Mathematics
More informationHypothesis testing: theory and methods
Statistical Methods Warsaw School of Economics November 3, 2017 Statistical hypothesis is the name of any conjecture about unknown parameters of a population distribution. The hypothesis should be verifiable
More informationLecture 13: Subsampling vs Bootstrap. Dimitris N. Politis, Joseph P. Romano, Michael Wolf
Lecture 13: 2011 Bootstrap ) R n x n, θ P)) = τ n ˆθn θ P) Example: ˆθn = X n, τ n = n, θ = EX = µ P) ˆθ = min X n, τ n = n, θ P) = sup{x : F x) 0} ) Define: J n P), the distribution of τ n ˆθ n θ P) under
More informationA Conditional Approach to Modeling Multivariate Extremes
A Approach to ing Multivariate Extremes By Heffernan & Tawn Department of Statistics Purdue University s April 30, 2014 Outline s s Multivariate Extremes s A central aim of multivariate extremes is trying
More informationFall 2017 STAT 532 Homework Peter Hoff. 1. Let P be a probability measure on a collection of sets A.
1. Let P be a probability measure on a collection of sets A. (a) For each n N, let H n be a set in A such that H n H n+1. Show that P (H n ) monotonically converges to P ( k=1 H k) as n. (b) For each n
More informationUNIVERSITY OF CALGARY. A New Hybrid Estimation Method for the. Generalized Exponential Distribution. Shan Zhu A THESIS
UNIVERSITY OF CALGARY A New Hybrid Estimation Method for the Generalized Exponential Distribution by Shan Zhu A THESIS SUBMITTED TO THE FACULTY OF GRADUATE STUDIES IN PARTIAL FULFILLMENT OF THE REQUIREMENTS
More informationDepartment of Statistics, School of Mathematical Sciences, Ferdowsi University of Mashhad, Iran.
JIRSS (2012) Vol. 11, No. 2, pp 191-202 A Goodness of Fit Test For Exponentiality Based on Lin-Wong Information M. Abbasnejad, N. R. Arghami, M. Tavakoli Department of Statistics, School of Mathematical
More informationA goodness of fit test for the Pareto distribution
Chilean Journal of Statistics Vol. 9, No. 1, April 2018, 33-46 A goodness of fit test for the Pareto distribution Javier Suárez-Espinosa, José A. Villaseñor-Alva, Annel Hurtado-Jaramillo and Paulino Pérez-Rodríguez
More informationModified Kolmogorov-Smirnov Test of Goodness of Fit. Catalonia-BarcelonaTECH, Spain
152/304 CoDaWork 2017 Abbadia San Salvatore (IT) Modified Kolmogorov-Smirnov Test of Goodness of Fit G.S. Monti 1, G. Mateu-Figueras 2, M. I. Ortego 3, V. Pawlowsky-Glahn 2 and J. J. Egozcue 3 1 Department
More informationA comparison of inverse transform and composition methods of data simulation from the Lindley distribution
Communications for Statistical Applications and Methods 2016, Vol. 23, No. 6, 517 529 http://dx.doi.org/10.5351/csam.2016.23.6.517 Print ISSN 2287-7843 / Online ISSN 2383-4757 A comparison of inverse transform
More informationA New Two Sample Type-II Progressive Censoring Scheme
A New Two Sample Type-II Progressive Censoring Scheme arxiv:609.05805v [stat.me] 9 Sep 206 Shuvashree Mondal, Debasis Kundu Abstract Progressive censoring scheme has received considerable attention in
More informationIrr. Statistical Methods in Experimental Physics. 2nd Edition. Frederick James. World Scientific. CERN, Switzerland
Frederick James CERN, Switzerland Statistical Methods in Experimental Physics 2nd Edition r i Irr 1- r ri Ibn World Scientific NEW JERSEY LONDON SINGAPORE BEIJING SHANGHAI HONG KONG TAIPEI CHENNAI CONTENTS
More informationExtreme Value Analysis and Spatial Extremes
Extreme Value Analysis and Department of Statistics Purdue University 11/07/2013 Outline Motivation 1 Motivation 2 Extreme Value Theorem and 3 Bayesian Hierarchical Models Copula Models Max-stable Models
More informationSTAT 135 Lab 5 Bootstrapping and Hypothesis Testing
STAT 135 Lab 5 Bootstrapping and Hypothesis Testing Rebecca Barter March 2, 2015 The Bootstrap Bootstrap Suppose that we are interested in estimating a parameter θ from some population with members x 1,...,
More informationChapter 1: A Brief Review of Maximum Likelihood, GMM, and Numerical Tools. Joan Llull. Microeconometrics IDEA PhD Program
Chapter 1: A Brief Review of Maximum Likelihood, GMM, and Numerical Tools Joan Llull Microeconometrics IDEA PhD Program Maximum Likelihood Chapter 1. A Brief Review of Maximum Likelihood, GMM, and Numerical
More informationClassical Extreme Value Theory - An Introduction
Chapter 1 Classical Extreme Value Theory - An Introduction 1.1 Introduction Asymptotic theory of functions of random variables plays a very important role in modern statistics. The objective of the asymptotic
More informationInferences on stress-strength reliability from weighted Lindley distributions
Inferences on stress-strength reliability from weighted Lindley distributions D.K. Al-Mutairi, M.E. Ghitany & Debasis Kundu Abstract This paper deals with the estimation of the stress-strength parameter
More informationDistribution Fitting (Censored Data)
Distribution Fitting (Censored Data) Summary... 1 Data Input... 2 Analysis Summary... 3 Analysis Options... 4 Goodness-of-Fit Tests... 6 Frequency Histogram... 8 Comparison of Alternative Distributions...
More informationA comparison study of the nonparametric tests based on the empirical distributions
통계연구 (2015), 제 20 권제 3 호, 1-12 A comparison study of the nonparametric tests based on the empirical distributions Hyo-Il Park 1) Abstract In this study, we propose a nonparametric test based on the empirical
More informationA note on vector-valued goodness-of-fit tests
A note on vector-valued goodness-of-fit tests Vassilly Voinov and Natalie Pya Kazakhstan Institute of Management, Economics and Strategic Research 050010 Almaty, Kazakhstan (e-mail: voinovv@kimep.kz, pya@kimep.kz)
More informationTests of fit for symmetric variance gamma distributions
Tests of fit for symmetri variane gamma distributions Fragiadakis Kostas UADPhilEon, National and Kapodistrian University of Athens, 4 Euripidou Street, 05 59 Athens, Greee. Keywords: Variane Gamma Distribution,
More informationA TEST OF RANDOMNESS BASED ON THE DISTANCE BETWEEN CONSECUTIVE RANDOM NUMBER PAIRS. Matthew J. Duggan John H. Drew Lawrence M.
Proceedings of the 2005 Winter Simulation Conference M. E. Kuhl, N. M. Steiger, F. B. Armstrong, and J. A. Joines, eds. A TEST OF RANDOMNESS BASED ON THE DISTANCE BETWEEN CONSECUTIVE RANDOM NUMBER PAIRS
More informationAsymptotic Statistics-III. Changliang Zou
Asymptotic Statistics-III Changliang Zou The multivariate central limit theorem Theorem (Multivariate CLT for iid case) Let X i be iid random p-vectors with mean µ and and covariance matrix Σ. Then n (
More informationStatistics 3858 : Maximum Likelihood Estimators
Statistics 3858 : Maximum Likelihood Estimators 1 Method of Maximum Likelihood In this method we construct the so called likelihood function, that is L(θ) = L(θ; X 1, X 2,..., X n ) = f n (X 1, X 2,...,
More informationGoodness of Fit Tests for Rayleigh Distribution Based on Phi-Divergence
Revista Colombiana de Estadística July 2017, Volume 40, Issue 2, pp. 279 to 290 DOI: http://dx.doi.org/10.15446/rce.v40n2.60375 Goodness of Fit Tests for Rayleigh Distribution Based on Phi-Divergence Pruebas
More informationarxiv:math/ v1 [math.pr] 9 Sep 2003
arxiv:math/0309164v1 [math.pr] 9 Sep 003 A NEW TEST FOR THE MULTIVARIATE TWO-SAMPLE PROBLEM BASED ON THE CONCEPT OF MINIMUM ENERGY G. Zech and B. Aslan University of Siegen, Germany August 8, 018 Abstract
More informationBivariate Rainfall and Runoff Analysis Using Entropy and Copula Theories
Entropy 2012, 14, 1784-1812; doi:10.3390/e14091784 Article OPEN ACCESS entropy ISSN 1099-4300 www.mdpi.com/journal/entropy Bivariate Rainfall and Runoff Analysis Using Entropy and Copula Theories Lan Zhang
More informationStep-Stress Models and Associated Inference
Department of Mathematics & Statistics Indian Institute of Technology Kanpur August 19, 2014 Outline Accelerated Life Test 1 Accelerated Life Test 2 3 4 5 6 7 Outline Accelerated Life Test 1 Accelerated
More informationExtreme Value Theory as a Theoretical Background for Power Law Behavior
Extreme Value Theory as a Theoretical Background for Power Law Behavior Simone Alfarano 1 and Thomas Lux 2 1 Department of Economics, University of Kiel, alfarano@bwl.uni-kiel.de 2 Department of Economics,
More informationGoodness of Fit: an axiomatic approach
Goodness of Fit: an axiomatic approach by Frank A. Cowell STICERD London School of Economics Houghton Street London, WC2A 2AE, UK email: f.cowell@lse.ac.uk Russell Davidson AMSE-GREQAM Department of Economics
More informationChapter 31 Application of Nonparametric Goodness-of-Fit Tests for Composite Hypotheses in Case of Unknown Distributions of Statistics
Chapter Application of Nonparametric Goodness-of-Fit Tests for Composite Hypotheses in Case of Unknown Distributions of Statistics Boris Yu. Lemeshko, Alisa A. Gorbunova, Stanislav B. Lemeshko, and Andrey
More informationModel Fitting, Bootstrap, & Model Selection
Model Fitting, Bootstrap, & Model Selection G. Jogesh Babu Penn State University http://www.stat.psu.edu/ babu http://astrostatistics.psu.edu Model Fitting Non-linear regression Density (shape) estimation
More informationAdvanced Statistics II: Non Parametric Tests
Advanced Statistics II: Non Parametric Tests Aurélien Garivier ParisTech February 27, 2011 Outline Fitting a distribution Rank Tests for the comparison of two samples Two unrelated samples: Mann-Whitney
More informationAbstract: In this short note, I comment on the research of Pisarenko et al. (2014) regarding the
Comment on Pisarenko et al. Characterization of the Tail of the Distribution of Earthquake Magnitudes by Combining the GEV and GPD Descriptions of Extreme Value Theory Mathias Raschke Institution: freelancer
More information1.0 I] MICROCOPY RESOLUTION TEST CHART ..: -,.,-..., ,,..,e ' ilhi~am &32. Il W tle p 10 A..
AD-R124 835 A NEW GOODNESS OF FIT TEST FOR THE UNIFORM DISTRIBUTION i WITH UNSPECIFIED PRRRMETERS(U) AIR FORCE INST OF TECH WRIGHT-PRTTERSON RFB OH SCHOOL OF ENGI. L B WOODBURY UNLSIIDDEC 82 RFIT/GOR/NA/82D-6
More informationA copula goodness-of-t approach. conditional probability integral transform. Daniel Berg 1 Henrik Bakken 2
based on the conditional probability integral transform Daniel Berg 1 Henrik Bakken 2 1 Norwegian Computing Center (NR) & University of Oslo (UiO) 2 Norwegian University of Science and Technology (NTNU)
More informationUse of mean residual life in testing departures from exponentiality
Nonparametric Statistics Vol. 18, No. 3, April 2006, 277 292 Use of mean residual life in testing departures from exponentiality S. RAO JAMMALAMADAKA and EMANUELE TAUFER* Department of Statistics and Applied
More informationExact goodness-of-fit tests for censored data
Exact goodness-of-fit tests for censored data Aurea Grané Statistics Department. Universidad Carlos III de Madrid. Abstract The statistic introduced in Fortiana and Grané (23, Journal of the Royal Statistical
More informationTESTING PROCEDURES BASED ON THE EMPIRICAL CHARACTERISTIC FUNCTIONS II: k-sample PROBLEM, CHANGE POINT PROBLEM. 1. Introduction
Tatra Mt. Math. Publ. 39 (2008), 235 243 t m Mathematical Publications TESTING PROCEDURES BASED ON THE EMPIRICAL CHARACTERISTIC FUNCTIONS II: k-sample PROBLEM, CHANGE POINT PROBLEM Marie Hušková Simos
More informationTESTS BASED ON EMPIRICAL DISTRIBUTION FUNCTION. Submitted in partial fulfillment of the requirements for the award of the degree of
TESTS BASED ON EMPIRICAL DISTRIBUTION FUNCTION Submitted in partial fulfillment of the requirements for the award of the degree of MASTER OF SCIENCE IN MATHEMATICS AND COMPUTING Submitted by Gurpreet Kaur
More informationRobustness and Distribution Assumptions
Chapter 1 Robustness and Distribution Assumptions 1.1 Introduction In statistics, one often works with model assumptions, i.e., one assumes that data follow a certain model. Then one makes use of methodology
More informationLecture 2: CDF and EDF
STAT 425: Introduction to Nonparametric Statistics Winter 2018 Instructor: Yen-Chi Chen Lecture 2: CDF and EDF 2.1 CDF: Cumulative Distribution Function For a random variable X, its CDF F () contains all
More informationChapter 6. Order Statistics and Quantiles. 6.1 Extreme Order Statistics
Chapter 6 Order Statistics and Quantiles 61 Extreme Order Statistics Suppose we have a finite sample X 1,, X n Conditional on this sample, we define the values X 1),, X n) to be a permutation of X 1,,
More informationAnalysis of Gamma and Weibull Lifetime Data under a General Censoring Scheme and in the presence of Covariates
Communications in Statistics - Theory and Methods ISSN: 0361-0926 (Print) 1532-415X (Online) Journal homepage: http://www.tandfonline.com/loi/lsta20 Analysis of Gamma and Weibull Lifetime Data under a
More informationSome New Aspects of Dose-Response Models with Applications to Multistage Models Having Parameters on the Boundary
Some New Aspects of Dose-Response Models with Applications to Multistage Models Having Parameters on the Boundary Bimal Sinha Department of Mathematics & Statistics University of Maryland, Baltimore County,
More informationIE 303 Discrete-Event Simulation L E C T U R E 6 : R A N D O M N U M B E R G E N E R A T I O N
IE 303 Discrete-Event Simulation L E C T U R E 6 : R A N D O M N U M B E R G E N E R A T I O N Review of the Last Lecture Continuous Distributions Uniform distributions Exponential distributions and memoryless
More informationIntroduction to Rare Event Simulation
Introduction to Rare Event Simulation Brown University: Summer School on Rare Event Simulation Jose Blanchet Columbia University. Department of Statistics, Department of IEOR. Blanchet (Columbia) 1 / 31
More informationChange Point Analysis of Extreme Values
Change Point Analysis of Extreme Values TIES 2008 p. 1/? Change Point Analysis of Extreme Values Goedele Dierckx Economische Hogeschool Sint Aloysius, Brussels, Belgium e-mail: goedele.dierckx@hubrussel.be
More informationACTEX CAS EXAM 3 STUDY GUIDE FOR MATHEMATICAL STATISTICS
ACTEX CAS EXAM 3 STUDY GUIDE FOR MATHEMATICAL STATISTICS TABLE OF CONTENTS INTRODUCTORY NOTE NOTES AND PROBLEM SETS Section 1 - Point Estimation 1 Problem Set 1 15 Section 2 - Confidence Intervals and
More informationNAG Library Chapter Introduction. G08 Nonparametric Statistics
NAG Library Chapter Introduction G08 Nonparametric Statistics Contents 1 Scope of the Chapter.... 2 2 Background to the Problems... 2 2.1 Parametric and Nonparametric Hypothesis Testing... 2 2.2 Types
More informationinferences on stress-strength reliability from lindley distributions
inferences on stress-strength reliability from lindley distributions D.K. Al-Mutairi, M.E. Ghitany & Debasis Kundu Abstract This paper deals with the estimation of the stress-strength parameter R = P (Y
More information1 Mixed effect models and longitudinal data analysis
1 Mixed effect models and longitudinal data analysis Mixed effects models provide a flexible approach to any situation where data have a grouping structure which introduces some kind of correlation between
More informationEstimation of parametric functions in Downton s bivariate exponential distribution
Estimation of parametric functions in Downton s bivariate exponential distribution George Iliopoulos Department of Mathematics University of the Aegean 83200 Karlovasi, Samos, Greece e-mail: geh@aegean.gr
More informationParameter Estimation
Parameter Estimation Consider a sample of observations on a random variable Y. his generates random variables: (y 1, y 2,, y ). A random sample is a sample (y 1, y 2,, y ) where the random variables y
More informationVariable inspection plans for continuous populations with unknown short tail distributions
Variable inspection plans for continuous populations with unknown short tail distributions Wolfgang Kössler Abstract The ordinary variable inspection plans are sensitive to deviations from the normality
More informationLecture 7 Introduction to Statistical Decision Theory
Lecture 7 Introduction to Statistical Decision Theory I-Hsiang Wang Department of Electrical Engineering National Taiwan University ihwang@ntu.edu.tw December 20, 2016 1 / 55 I-Hsiang Wang IT Lecture 7
More informationNew Consistent Integral-Type Tests for Stochastic Dominance
New Consistent Integral-Type Tests for Stochastic Dominance Job Market Paper Chris J. Bennett January 5, 28 Abstract This paper proposes and examines several new statistics for testing stochastic dominance
More informationHypothesis Test. The opposite of the null hypothesis, called an alternative hypothesis, becomes
Neyman-Pearson paradigm. Suppose that a researcher is interested in whether the new drug works. The process of determining whether the outcome of the experiment points to yes or no is called hypothesis
More informationTESTS OF FIT FOR THE RAYLEIGH DISTRIBUTION BASED ON THE EMPIRICAL LAPLACE TRANSFORM
Ann. Inst. Statist. Math. Vol. 55, No. 1, 137-151 (2003) (~)2003 The Institute of Statistical Mathematics TESTS OF FIT FOR THE RAYLEIGH DISTRIBUTION BASED ON THE EMPIRICAL LAPLACE TRANSFORM SIMOS MEINTANIS
More informationDifferent methods of estimation for generalized inverse Lindley distribution
Different methods of estimation for generalized inverse Lindley distribution Arbër Qoshja & Fatmir Hoxha Department of Applied Mathematics, Faculty of Natural Science, University of Tirana, Albania,, e-mail:
More informationCherry Blossom run (1) The credit union Cherry Blossom Run is a 10 mile race that takes place every year in D.C. In 2009 there were participants
18.650 Statistics for Applications Chapter 5: Parametric hypothesis testing 1/37 Cherry Blossom run (1) The credit union Cherry Blossom Run is a 10 mile race that takes place every year in D.C. In 2009
More informationReview. DS GA 1002 Statistical and Mathematical Models. Carlos Fernandez-Granda
Review DS GA 1002 Statistical and Mathematical Models http://www.cims.nyu.edu/~cfgranda/pages/dsga1002_fall16 Carlos Fernandez-Granda Probability and statistics Probability: Framework for dealing with
More informationNon-parametric Inference and Resampling
Non-parametric Inference and Resampling Exercises by David Wozabal (Last update. Juni 010) 1 Basic Facts about Rank and Order Statistics 1.1 10 students were asked about the amount of time they spend surfing
More informationarxiv: v2 [math.st] 14 Sep 2015
Two-Sample Smooth Tests for the Equality of Distributions Wen-Xin Zhou, Chao Zheng, and Zhen Zhang arxiv:59.3459v2 [math.st] 4 Sep 25 Abstract This paper considers the problem of testing the equality of
More informationIntroduction to Empirical Processes and Semiparametric Inference Lecture 09: Stochastic Convergence, Continued
Introduction to Empirical Processes and Semiparametric Inference Lecture 09: Stochastic Convergence, Continued Michael R. Kosorok, Ph.D. Professor and Chair of Biostatistics Professor of Statistics and
More information