1988: Logarithmic Series Distribution and Its Use In

Total Page:16

File Type:pdf, Size:1020Kb

1988: Logarithmic Series Distribution and Its Use In LOGARITHMIC SERIES DISTRIBUTION AND ITS USE IN ANALYZING DISCRETE DATA Jeffrey R, Wilson, Arizona State University Tempe, Arizona 85287 1. Introduction is assumed to be 7. for any unit and any trial. In a previous paper, Wilson and Koehler Furthermore, for ~ach unit the m trials are (1988) used the generalized-Dirichlet identical and independent so upon summing across multinomial model to account for extra trials X. = (XI~. X 2 ...... X~.)', the vector of variation. The model allows for a second order counts ~r the j th ~nit has ~3multinomial of pairwise correlation among units, a type of distribution with probability vector ~ = (~1' assumption found reasonable in some biological 72 ..... 71) ' and sample size m. Howe~er, data. In that paper, the two-way crossed responses given by the J units at a particular generalized Dirichlet Multinomial model was used trial may be correlated, producing a set of J to analyze repeated measure on the categorical correlated multinomial random vectors, X I, X 2, preferences of insurance customers. The number • o.~ X • of respondents was assumed to be fixed and Ta~is (1962) developed a model for this known. situation, which he called the generalized- In this paper a generalization of the model multinomial distribution in which a single is made allowing the number of respondents m, to parameter P, is used to reflect the common be random. Thus both the number of units m, and dependency between any two of the dependent the underlying probability vector are allowed to multinomial random vectors. The distribution of vary. The model presented here uses the m logarithmic series distribution to account for the category total Xi3. = Z Xij k is binomial the variation among number of units and the k=l Dirichlet distribution to model the with sample size m and parameter ~., for each probabilities. In particular the Dirichlet- unit. Tallis formalized the dependencies among Multinomial distribution is used to incorporate unit totals by specifying the joint moment the two types of variation, and the logarithmic generating function as series distribution is used to account for m J variation among the number of units within a Gj(~) = P { Z pi ( H eUj) i} + given time period• Ignoring either level of i=0 j=1 variation leads to underestimation of the true (2.1) standard errors of estimated proportions. J Tallis (1962) proposed the use of the (l-p) { H p(eUj)} , 1,2 ..... I, generalized-multinomial model for dependent j=l multinomials. Wilson and Koehler (1988) u I iu extended the model to allow for a second random r r and = (~ , ~ , where P(e ) = l pi e ~ I 2 component• The extended model considered can be i=l viewed as multivariate extensions of the beta- .... ~j )'. The parameter P appearing in (2.1) binomial and correlated binomial models is the correlation coefficient between X.. and considered by Kupper and Haseman (1978) and for any j ~ j'. When O ~ 0, GT(U) ll~ a Xi., o aj Crowder (1978) for binary data. Paul (1987) weighted mean of moment generating function for considered a modification of the beta-correlated a distribution with perfect correlation and one binomial as, a means of analyzing affected with complete independence, the weights being p foetuses in litters of live foetuses. Section 2 and I- 0 respectively. Altham (1978) proposed a outlines the generalized-multinomial model. similar model for a joint moment generating Section 3 discusses the Dirichlet-Multinomial function for correlated binary variables. model. Section 4 presents the generalized Consider the overall vector of category Dirichlet-Multinomial model. The logarithmic J distribution is presented in Section 5. An totals X = Z X.. From the moment generating extended generalized Dirichlet-Multinomial model j=l~3 is developed in Section 6. Tests for certain function in (2.1) it can be shown that hypotheses and fit of the model are developed in E(X) = Jm~ and V(X) = Jm{l+(J-l)p}M for the Section 7. Parameter estimates are obtained in ~ 7 Section 8. generalized-multinomial model, wher~ 2. Generalized-Multinomial Model M = diag(~) - ~7' and diag(7) is a diagonal 7 ~ One way to view the generalized-multinomial m~trix with diagonal elements provided by the model is to consider mJ vectors of outcomes for (Jm)--X is an vector ~. Consequently, ~ = Ta~lis (1962) a set of J units that are simultaneously unbiase~ estimator for 7. subjected to a series of m trials. At each proposed estimators forgo, but he did not trial, each unit is classified as being in discuss techniques for making inferences about exactly one of I mutually exclusive states• Let 7. We consider here a technique for making such the random variable X.. k take the value i if the l~nferences. k-th trial of the j-t~3unit is observed to be in One approach is to use the limiting normal the i-th state, and zero otherwise. distribution of X as m÷~. At trial k consider The probability that Xij k takes the value i an IJ-dimensiona~ response vector ~(j) = (X~k, 275 ~2k'7''~$ ~) where = (X k' X ...... ' ~ ' e ~ Xljk)' respons $~ctor ~r thSJ~-th unit • .. ~n t)'"the generalized-multinomial model, the m X. t vectors (j=l,2 ..... J) are identically at the k-th trial. Define ~(j)X = k=iX X 4k(J) ~stributed and are not independent. The observations taken at time t on the J Since the X . vectors are independent and the ~k J individuals are equally pairwise correlated as first and second moments of l~._. are finite, . ~ ~ measured by the parameter P- The vector of the mul~ivariate central llm~ ~eorem Implies n J that m-~(Z... - mu) ÷ N__(0, E), as m÷~, where total counts ~X = X X~t, where Xt= E Xjt for u = i_@ ~, ~ = M @ Q, ~. is a J-dimensional ~J ~ ~ ~j t=l j=l vector of ones, @~denotes direct product between the generalized Dirichlet-Multinomial model, has and Q is a square matrix of dimension J with mean vector E(X) = N~ and covariance matrix ones on the diagonal and p as each off diagonal V(X) = NC{I+ 0(~-I)}M ~, where N=nS is the total element. Now X = GX._. where G = i! @ E . and number of observations, S is the to~al number of E 1 is the identity m~ix of dimension I. IThen, un'its at time t and C = (S+o)(l+o) . Using an by Rao (1973 page 124) the limiting distribution argument similar to the one in section 2 it can of ~ is specified by be shown that n -½ (X-N~) ÷ N I(O, SC{I+(J-{)o}M ) ^ m-½(~- ~ ÷ NI( ~, (mj)-l{l+(J-1)p}M) (2 2) and tests of hypotheses about ~ or vector "kJ functions g(~), where g is a continuous function Given a consistent estimator for 0, with second partial derivatives, can be obtained asymptotic chi-square tests involving using the large sample chi-square distributions sufficiently smooth functions of ~ can be for the Wald statistic obtained as Wald statistics, N{C{I+(J-I)o}}-I(g(~) _ (4.1) 4 = mJ{ l+(J-1) p }-I [g(~)_^ (2.3) g(~))'(DM D')-(g(~) - g(~)), g(~) ]' [DMSD' ]-[g($)-g(~) ] where [DM D'] denotes the generalized inverse where D is the matrix of first partial of DM D', with degrees of freedom equal to rank derivatives of g evaluated at ~, and [DMSD]- is aj of DM D'. The greater imprecision in the a generalized inverse of DM^D. The degrees of estimation for ~ due to variation in vectors of freedom correspond to the r~nk of DM D'. In nu proportion among individuals is accounted for by some applicatioDs it may be necessar~ ~o replace the factor C which cannot be less than one. The replace D with D (i.e. D evaluated at ~). consequence of ignoring this extra variation is 3. Dirichlet-Multinomial Model an inflation of the type I error levels for such For each of N units abserve a multinomial tests. vector of responses, with parameters ~ = (PI' 5. The Univariate Logarithmic P2 .... ~h~l )' and sample size S. Furthermore Series Distribution assume probability vector ~ has a Dirichlet The logarithmic series distribution was distribution! with mean vector ~ = (~i' ~2 , oo., introduced by Fisher, Corbett and Williams ~i ) and scaling parameter ~. (1943) to investigate the distribution of For this model the sum of the vector of butterflies in the Malayan Peninsula. It has counts has a Dirichlet-Multinomial distribution been used in the sampling of quadrants for plant and the vector of proportions has first moment species, the distribution of animal species, and covariance matrix N -I(s+~) (l+a)-iM . The population growth and in economic applications. Chatfield et al (1966) used the logarithmic Dirichlet distribution provides a convenient series distribution to represent the model for describing variation among vectors of distribution of numbers of items of a product proportions since it has relatively simple purchased by a buyer in a specified time period. mathematical properties. The Dirichlet- They point out that the logarithmic series has Multinomial model has been studied by Mosimann the advantage of dependency on only one (1962) and Good (1965). Brier (1980) used the parameter O. model to analyze sample proportions obtained The random variable M has a logarithmic from a single two-stage cluster sample. Koehler series distribution if the probability function and Wilson (1986) extended some of Brier's P(M=k) = a0~/k (k=l,2 ..... ; 0<0<i) results to analyze vectors of proportions where a=-[log(l-0)]-l.
Recommended publications
  • Benford's Law As a Logarithmic Transformation
    Benford’s Law as a Logarithmic Transformation J. M. Pimbley Maxwell Consulting, LLC Abstract We review the history and description of Benford’s Law - the observation that many naturally occurring numerical data collections exhibit a logarithmic distribution of digits. By creating numerous hypothetical examples, we identify several cases that satisfy Benford’s Law and several that do not. Building on prior work, we create and demonstrate a “Benford Test” to determine from the functional form of a probability density function whether the resulting distribution of digits will find close agreement with Benford’s Law. We also discover that one may generalize the first-digit distribution algorithm to include non-logarithmic transformations. This generalization shows that Benford’s Law rests on the tendency of the transformed and shifted data to exhibit a uniform distribution. Introduction What the world calls “Benford’s Law” is a marvelous mélange of mathematics, data, philosophy, empiricism, theory, mysticism, and fraud. Even with all these qualifiers, one easily describes this law and then asks the simple questions that have challenged investigators for more than a century. Take a large collection of positive numbers which may be integers or real numbers or both. Focus only on the first (non-zero) digit of each number. Count how many numbers of the large collection have each of the nine possibilities (1-9) as the first digit. For typical J. M. Pimbley, “Benford’s Law as a Logarithmic Transformation,” Maxwell Consulting Archives, 2014. number collections – which we’ll generally call “datasets” – the first digit is not equally distributed among the values 1-9.
    [Show full text]
  • Notes on Scale-Invariance and Base-Invariance for Benford's
    NOTES ON SCALE-INVARIANCE AND BASE-INVARIANCE FOR BENFORD’S LAW MICHAŁ RYSZARD WÓJCIK Abstract. It is known that if X is uniformly distributed modulo 1 and Y is an arbitrary random variable independent of X then Y + X is also uniformly distributed modulo 1. We prove a converse for any continuous random variable Y (or a reasonable approximation to a continuous random variable) so that if X and Y +X are equally distributed modulo 1 and Y is independent of X then X is uniformly distributed modulo 1 (or approximates the uniform distribution equally reasonably). This translates into a characterization of Benford’s law through a generalization of scale-invariance: from multiplication by a constant to multiplication by an independent random variable. We also show a base-invariance characterization: if a positive continuous random variable has the same significand distribution for two bases then it is Benford for both bases. The set of bases for which a random variable is Benford is characterized through characteristic functions. 1. Introduction Before the early 1970s, handheld electronic calculators were not yet in widespread use and scientists routinely used in their calculations books with tables containing the decimal logarithms of numbers between 1 and 10 spaced evenly with small increments like 0.01 or 0.001. For example, the first page would be filled with numbers 1.01, 1.02, 1.03, . , 1.99 in the left column and their decimal logarithms in the right column, while the second page with 2.00, 2.01, 2.02, 2.03, ..., 2.99, and so on till the ninth page with 9.00, 9.01, 9.02, 9.03, .
    [Show full text]
  • Package 'Distributional'
    Package ‘distributional’ February 2, 2021 Title Vectorised Probability Distributions Version 0.2.2 Description Vectorised distribution objects with tools for manipulating, visualising, and using probability distributions. Designed to allow model prediction outputs to return distributions rather than their parameters, allowing users to directly interact with predictive distributions in a data-oriented workflow. In addition to providing generic replacements for p/d/q/r functions, other useful statistics can be computed including means, variances, intervals, and highest density regions. License GPL-3 Imports vctrs (>= 0.3.0), rlang (>= 0.4.5), generics, ellipsis, stats, numDeriv, ggplot2, scales, farver, digest, utils, lifecycle Suggests testthat (>= 2.1.0), covr, mvtnorm, actuar, ggdist RdMacros lifecycle URL https://pkg.mitchelloharawild.com/distributional/, https: //github.com/mitchelloharawild/distributional BugReports https://github.com/mitchelloharawild/distributional/issues Encoding UTF-8 Language en-GB LazyData true Roxygen list(markdown = TRUE, roclets=c('rd', 'collate', 'namespace')) RoxygenNote 7.1.1 1 2 R topics documented: R topics documented: autoplot.distribution . .3 cdf..............................................4 density.distribution . .4 dist_bernoulli . .5 dist_beta . .6 dist_binomial . .7 dist_burr . .8 dist_cauchy . .9 dist_chisq . 10 dist_degenerate . 11 dist_exponential . 12 dist_f . 13 dist_gamma . 14 dist_geometric . 16 dist_gumbel . 17 dist_hypergeometric . 18 dist_inflated . 20 dist_inverse_exponential . 20 dist_inverse_gamma
    [Show full text]
  • Package 'Univrng'
    Package ‘UnivRNG’ March 5, 2021 Type Package Title Univariate Pseudo-Random Number Generation Version 1.2.3 Date 2021-03-05 Author Hakan Demirtas, Rawan Allozi, Ran Gao Maintainer Ran Gao <[email protected]> Description Pseudo-random number generation of 17 univariate distributions proposed by Demir- tas. (2005) <DOI:10.22237/jmasm/1114907220>. License GPL-2 | GPL-3 NeedsCompilation no Repository CRAN Date/Publication 2021-03-05 18:10:02 UTC R topics documented: UnivRNG-package . .2 draw.beta.alphabeta.less.than.one . .3 draw.beta.binomial . .4 draw.gamma.alpha.greater.than.one . .5 draw.gamma.alpha.less.than.one . .6 draw.inverse.gaussian . .7 draw.laplace . .8 draw.left.truncated.gamma . .8 draw.logarithmic . .9 draw.noncentral.chisquared . 10 draw.noncentral.F . 11 draw.noncentral.t . 12 draw.pareto . 12 draw.rayleigh . 13 draw.t ............................................ 14 draw.von.mises . 15 draw.weibull . 16 draw.zeta . 16 1 2 UnivRNG-package Index 18 UnivRNG-package Univariate Pseudo-Random Number Generation Description This package implements the algorithms described in Demirtas (2005) for pseudo-random number generation of 17 univariate distributions. The following distributions are available: Left Truncated Gamma, Laplace, Inverse Gaussian, Von Mises, Zeta (Zipf), Logarithmic, Beta-Binomial, Rayleigh, Pareto, Non-central t, Non-central Chi-squared, Doubly non-central F , Standard t, Weibull, Gamma with α<1, Gamma with α>1, and Beta with α<1 and β<1. For some distributions, functions that have similar capabilities exist in the base package; the functions herein should be regarded as com- plementary tools. The methodology for each random-number generation procedure varies and each distribution has its own function.
    [Show full text]
  • Hand-Book on STATISTICAL DISTRIBUTIONS for Experimentalists
    Internal Report SUF–PFY/96–01 Stockholm, 11 December 1996 1st revision, 31 October 1998 last modification 10 September 2007 Hand-book on STATISTICAL DISTRIBUTIONS for experimentalists by Christian Walck Particle Physics Group Fysikum University of Stockholm (e-mail: [email protected]) Contents 1 Introduction 1 1.1 Random Number Generation .............................. 1 2 Probability Density Functions 3 2.1 Introduction ........................................ 3 2.2 Moments ......................................... 3 2.2.1 Errors of Moments ................................ 4 2.3 Characteristic Function ................................. 4 2.4 Probability Generating Function ............................ 5 2.5 Cumulants ......................................... 6 2.6 Random Number Generation .............................. 7 2.6.1 Cumulative Technique .............................. 7 2.6.2 Accept-Reject technique ............................. 7 2.6.3 Composition Techniques ............................. 8 2.7 Multivariate Distributions ................................ 9 2.7.1 Multivariate Moments .............................. 9 2.7.2 Errors of Bivariate Moments .......................... 9 2.7.3 Joint Characteristic Function .......................... 10 2.7.4 Random Number Generation .......................... 11 3 Bernoulli Distribution 12 3.1 Introduction ........................................ 12 3.2 Relation to Other Distributions ............................. 12 4 Beta distribution 13 4.1 Introduction .......................................
    [Show full text]
  • A Class of Generalized Beta Distributions, Pareto Power Series and Weibull Power Series
    A CLASS OF GENERALIZED BETA DISTRIBUTIONS, PARETO POWER SERIES AND WEIBULL POWER SERIES ALICE LEMOS DE MORAIS Primary advisor: Prof. Audrey Helen M. A. Cysneiros Secondary advisor: Prof. Gauss Moutinho Cordeiro Concentration area: Probability Disserta¸c˜ao submetida como requerimento parcial para obten¸c˜aodo grau de Mestre em Estat´ısticapela Universidade Federal de Pernambuco. Recife, fevereiro de 2009 Morais, Alice Lemos de A class of generalized beta distributions, Pareto power series and Weibull power series / Alice Lemos de Morais - Recife : O Autor, 2009. 103 folhas : il., fig., tab. Dissertação (mestrado) – Universidade Federal de Pernambuco. CCEN. Estatística, 2009. Inclui bibliografia e apêndice. 1. Probabilidade. I. Título. 519.2 CDD (22.ed.) MEI2009-034 Agradecimentos A` minha m˜ae,Marcia, e ao meu pai, Marcos, por serem meus melhores amigos. Agrade¸co pelo apoio `aminha vinda para Recife e pelo apoio financeiro. Agrade¸coaos meus irm˜aos, Daniel, Gabriel e Isadora, bem como a meus primos, Danilo e Vanessa, por divertirem minhas f´erias. A` minha tia Gilma pelo carinho e por pedir `asua amiga, Mirlana, que me acolhesse nos meus primeiros dias no Recife. A` Mirlana por amenizar minha mudan¸cade ambiente fazendo-me sentir em casa durante meus primeiros dias nesta cidade. Agrade¸comuito a ela e a toda sua fam´ılia. A` toda minha fam´ılia e amigos pela despedida de Belo Horizonte e `aS˜aozinha por preparar meu prato predileto, frango ao molho pardo com angu. A` vov´oLuzia e ao vovˆoTunico por toda amizade, carinho e preocupa¸c˜ao, por todo apoio nos meus dias mais dif´ıceis em Recife e pelo apoio financeiro nos momentos que mais precisei.
    [Show full text]
  • Concentration Inequalities from Likelihood Ratio Method
    Concentration Inequalities from Likelihood Ratio Method ∗ Xinjia Chen September 2014 Abstract We explore the applications of our previously established likelihood-ratio method for deriving con- centration inequalities for a wide variety of univariate and multivariate distributions. New concentration inequalities for various distributions are developed without the idea of minimizing moment generating functions. Contents 1 Introduction 3 2 Likelihood Ratio Method 4 2.1 GeneralPrinciple.................................... ...... 4 2.2 Construction of Parameterized Distributions . ............. 5 2.2.1 WeightFunction .................................... .. 5 2.2.2 ParameterRestriction .............................. ..... 6 3 Concentration Inequalities for Univariate Distributions 7 3.1 BetaDistribution.................................... ...... 7 3.2 Beta Negative Binomial Distribution . ........ 7 3.3 Beta-Prime Distribution . ....... 8 arXiv:1409.6276v1 [math.ST] 1 Sep 2014 3.4 BorelDistribution ................................... ...... 8 3.5 ConsulDistribution .................................. ...... 8 3.6 GeetaDistribution ................................... ...... 9 3.7 GumbelDistribution.................................. ...... 9 3.8 InverseGammaDistribution. ........ 9 3.9 Inverse Gaussian Distribution . ......... 10 3.10 Lagrangian Logarithmic Distribution . .......... 10 3.11 Lagrangian Negative Binomial Distribution . .......... 10 3.12 Laplace Distribution . ....... 11 ∗The author is afflicted with the Department of Electrical
    [Show full text]
  • Ebookdistributions.Pdf
    DOWNLOAD YOUR FREE MODELRISK TRIAL Adapted from Risk Analysis: a quantitative guide by David Vose. Published by John Wiley and Sons (2008). All Rights Reserved. No part of this publication may be reproduced, stored in a retrieval system or transmitted in any form or by any means, electronic, mechanical, photocopying, recording, scanning or otherwise, except under the terms of the Copyright, Designs and Patents Act 1988 or under the terms of a licence issued by the Copyright Licensing Agency Ltd, 90 Tottenham Court Road, London W1T 4LP, UK, without the permission in writing of the Publisher. If you notice any errors or omissions, please contact [email protected] Referencing this document Please use the following reference, following the Harvard system of referencing: Van Hauwermeiren M, Vose D and Vanden Bossche S (2012). A Compendium of Distributions (second edition). [ebook]. Vose Software, Ghent, Belgium. Available from www.vosesoftware.com . Accessed dd/mm/yy. © Vose Software BVBA (2012) www.vosesoftware.com Updated 17 January, 2012. Page 2 Table of Contents Introduction .................................................................................................................... 7 DISCRETE AND CONTINUOUS DISTRIBUTIONS.......................................................................... 7 Discrete Distributions .............................................................................................. 7 Continuous Distributions ........................................................................................
    [Show full text]
  • Compendium of Common Probability Distributions
    Compendium of Common Probability Distributions Michael P. McLaughlin December, 2016 Compendium of Common Probability Distributions Second Edition, v2.7 Copyright c 2016 by Michael P. McLaughlin. All rights reserved. Second printing, with corrections. DISCLAIMER: This document is intended solely as a synopsis for the convenience of users. It makes no guarantee or warranty of any kind that its contents are appropriate for any purpose. All such decisions are at the discretion of the user. Any brand names and product names included in this work are trademarks, registered trademarks, or trade names of their respective holders. This document created using LATEX and TeXShop. Figures produced using Mathematica. Contents Preface iii Legend v I Continuous: Symmetric1 II Continuous: Skewed 21 III Continuous: Mixtures 77 IV Discrete: Standard 101 V Discrete: Mixtures 115 ii Preface This Compendium is part of the documentation for the software package, Regress+. The latter is a utility intended for univariate mathematical modeling and addresses both deteministic models (equations) as well as stochastic models (distributions). Its focus is on the modeling of empirical data so the models it contains are fully-parametrized variants of commonly used formulas. This document describes the distributions available in Regress+ (v2.7). Most of these are well known but some are not described explicitly in the literature. This Compendium supplies the formulas and parametrization as utilized in the software plus additional formulas, notes, etc. plus one or more
    [Show full text]
  • A Generalization of the Exponential-Logarithmic Distribution for Reliability and Life Data Analysis
    A Generalization of the Exponential-Logarithmic Distribution for Reliability and Life Data Analysis Mohieddine Rahmouni,∗ Ayman Orabiy January 5,2018 Abstract In this paper, we introduce a new two-parameter lifetime distribution, called the exponential-generalized truncated logarithmic (EGTL) distribution, by compounding the exponential and generalized truncated logarithmic distributions. Our procedure general- izes the exponential-logarithmic (EL) distribution modelling the reliability of systems by the use of first-order concepts, where the minimum lifetime is considered (Tahmasbi & Rezaei 2008). In our approach, we assume that a system fails if a given number k of the components fails and then, we consider the kth-smallest value of lifetime instead of the minimum lifetime. The reliability and failure rate functions as well as their properties are presented for some special cases. The estimation of the parameters is attained by the maximum likelihood, the expectation maximization algorithm, the method of moments and the Bayesian approach, with a simulation study performed to illustrate the different methods of estimation. The application study is illustrated based on two real data sets used in many applications of reliability. Keywords: Lifetime distributions, reliability, failure rate, order statistics, exponential distribution, truncated logarithmic distribution. 1 Introduction arXiv:1803.01156v1 [math.ST] 3 Mar 2018 Lifetime distributions are often used in reliability theory and survival analysis for modelling real data. They play a fundamental role in reliability in diverse disciplines such as finance, manufacture, biological sciences, physics and engineering. The exponential distribution is a ba- sic model in reliability theory and survival analysis. It is often used to model system reliability at a component level, assuming the failure rate is constant (Balakrishnan & Basu 1995, Barlow & Proschan 1975, Sinha & Kale 1980).
    [Show full text]
  • Stat::Fit® 3 Manual
    Stat::Fit® Statistically Fit® Software Copyright © 2016 Geer Mountain Software Corp. Table of Contents Getting Started . 3 Trees And Notes . 6 Fitting Distributions . 8 Autofit . 9 Manual Setup . 11 Moment Estimates . 13 Maximum Likelihood . 14 Goodness Of Fit . 16 Chi Squared Test . 18 Kolmogorov Smirnov Test. 2 0 Anderson Darling Test . 22 Empirical . 24 No Data . 25 Unbounded . 26 Bounded Above a Minimum . 28 Bounded Between a Minimum and a Maximum . 31 Distributions Continuous Beta . 37 Cauchy . 39 Chi Squared . 40 Erlang . 42 Exponential. 4 4 Extreme Value IA. 4 5 Extreme Value IB. 4 6 Gamma . 47 Inverse Gaussian . 49 Inverse Weibull. 5 0 Johnson SB . ..
    [Show full text]
  • Department of Econometrics and Business Statistics
    ISSN 1440-771X Australia Department of Econometrics and Business Statistics http://www.buseco.monash.edu.au/depts/ebs/pubs/wpapers/ Intermittent demand forecasting for inventory control: A multi-series approach Ralph Snyder, Adrian Beaumont and J. Keith Ord July 2012 Working Paper 15/12 Intermittent demand forecasting for inventory control: A multi-series approach Ralph Snyder, Adrian Beaumont1 and J. Keith Ord2 Abstract This paper is concerned with identifying an effective method for forecasting the lead time demand of slow-moving inventories. Particular emphasis is placed on prediction distributions instead of point predictions alone. It is also placed on methods which work with small samples as well as large samples in recognition of the fact that the typical range of items has a mix of vintages due to different commissioning and decommissioning dates over time. Various forecasting methods are compared using monthly demand data for more than one thousand car parts. It is found that a multi-series version of exponential smoothing coupled with a Pólya (negative binomial) distribution works better than the other twenty-four methods considered, including the Croston method. Author Keywords: demand forecasting; inventory control; shifted Poisson distribution JEL classification: C22 1 Department of Statistics, The University of Melbourne, Parkville, Victoria, Australia. [email protected] 2 McDonough School of Business, Georgetown University, Washington, DC 20057, USA [email protected] 1 1 Introduction Plausible state space models for count time series are developed and benchmarked against standard methods of forecasting (Croston 1972; Harvey and Fernandes 1989). These models permit a rational approach to estimation based on the principle of maximum likelihood and they provide a systematic framework for the derivation of discrete prediction distributions.
    [Show full text]