arxiv: v1 [stat.me] 6 Jun 2016
|
|
- Noreen Cooper
- 5 years ago
- Views:
Transcription
1 A gamma approximation to the Bayesian posterior distribution of a discrete parameter of the Generalized Poisson model arxiv: v1 [stat.me] 6 Jun 2016 Tsung Fei Khang Institute of Mathematical Sciences, Faculty of Science, University of Malaya, Kuala Lumpur, Malaysia. Abstract Let X have a Generalized Poisson distribution with mean kb, where b is a known constant in the unit interval and k is a discrete, non-negative parameter. We show that if an uninformative uniform prior for k is assumed, then the posterior distribution of k can be approximated using the gamma distribution when b is small. Keywords: Generalized Poisson distribution posterior distribution gamma distribution approximation 2000 MSC: 62E17, 62F15 1. Introduction The family of Generalized Poisson distributions (GP) (Consul and Jain, 1973) has been used for more than 40 years to model count data that may be overdispersed or underdispersed. Some of its interesting theoretical properties include a Poisson mixture interpretation (Joe and Zhu, 2005), and a heavier tail compared to the negative binomial distribution (Joe and Zhu, 2005; Nikoloulopoulos and Karlis, 2008). Various chance mechanisms have been found to generate the GP distribution (Shoukri and Consul, 1987). Numerous applications are given in Consul (1989), and in recent years, it has gained increasing popularity in bioinformatics for modelling RNA- Seqcount data(srivastava and Chen,2010; Li and Jiang,2012; Zhang et al., address: tfkhang@um.edu.my (Tsung Fei Khang) Preprint submitted to a future journal June 7, 2016
2 2015; Wang et al., 2015). In this work, we study the posterior distribution of a discrete parameter of the GP model under a particular parametrization. Let X be a random variable following GP(λ 1,λ 2 ). Its probability mass function (pmf) is given by P(X = x λ 1,λ 2 ) = λ 1(λ 1 +xλ 2 ) x 1 e (λ 1+xλ 2 ), (1) x! where x = 0,1,2,..., λ 1 > 0 and λ 2 < 1. Negative values of λ 2 correspond to overdispersion, positive values to underdispersion, and λ 2 = 0 reduces eq.(1) to the Poisson distribution with mean λ 1. Consider the following parametrization: λ 2 = 1 m, λ 1 = kb m, with m = e ab+c. We assume that a,c R and 0 < b < 1 are known constants, and focus our interest on the discrete non-negative parameter k. Under this parametrization, The mean and the variance of the GP model are given by E(X k) = λ 1 /(1 λ 2 ) = kb, Var(X k) = E(X k)/(1 λ 2 ) 2 = kb/m, We are concerned with the posterior probability of the k parameter. In the absence of any prior information, a Bayesian formulation using an improper uniform prior P(k) = 1,k = 0,1,2,... on k yields the posterior distribution of k as P(k X = x) = = = P(X = x k) j=xp(x = x j) k(bk m+x(1 m)) x 1 e bk m j=x j(bj m+x(1 m)) x 1 e bj m k(k +g(x)) x 1 e bk m j=x j(j +g(x))x 1 e b mj, (2) for k x, where g(x) = {(1 m)/(b m)}x. It is easy to check that P(k X = x) is proper even though an improper uniform prior distribution is used. The aim of this paper is to derive a continuous approximation to the posterior distribution eq.(2), so that the posterior mean and variance can be determined directly from the theoretical properties of the approximating distribution. 2
3 2. Results We show that under certain conditions, the gamma distribution approximates the posterior distribution of k. Theorem 1. If P(k X = x) is treated as a density function, then P(k X = x) is approximately equal to the probability density function of the gamma distribution with mean (x+1)/(b m) and variance (x+1)/(b 2 m) for some k l where l > x. Proof. First, we note that the denominator in eq.(2) can be written as (j +g(x)) x e b mj g(x) (j +g(x)) x 1 e b mj. j=x j=x The Lerch transcendent Φ(z,s,a) is given by Φ(z,s,a) = k=0 z k (a+k) s, where z < 1, a 0, 1, 2,..., and s 1,2,... Representing the denominator using the Lerch transcendent, we get e b mx Φ(e b m, x,wx) (w 1)xe b mx Φ(e b m, (x 1),wx), (3) where w = 1+(1 m)/(b m). The following identity (eq.1.11(11) in Erdélyi (1953)) relates the Lerch transcendent with negative argument for s to the Bernoulli polynomials: Φ(z, h,v) = h! ( log 1 ) (h+1) 1 z v z z v B h+r+1 (v)(logz) r, (4) r!(h+r +1) where logz < 2π, v 0, 1, 2,..., h 1, 2,..., and B n (x) is the n-th Bernoulli polynomial with argument x. The Bernoulli polynomial is defined as n ( ) n B n (x) = b n j x j, j j=0 3
4 where b i is the ith Bernoulli number. If we now substitute z = e b m, h = x, v = wx into the identity eq.(4), then we obtain Φ(e b m, x,wx) = e b m(wx) (b m) (x+1) 1 e b m(wx) B x+r+1 (wx)( b m) r r!(x+r +1) For small values of b (0 < b < 1), we note that the sum of Bernoulli polynomials is dominated by the zero-th term. Hence, for the first term in eq.(3), e bw mx Φ(e b m, x,wx) = (b m) B x+1(wx) +O(b m), (5) x+1 x+1 as b 0. An identity involving the Bernoulli polynomials and the sum of nth powers (eq. 1.13(10) in Erdélyi (1953)) gives B n+1 (x) b n+1 n+1 Since B x+1 (x) dominates b x+1, eq.(5) becomes x 1 = r n. e bw mx Φ(e b m, x,wx) [wx] 1 (b m) r x, x+1 where [wx] is the integer part of wx. Here, provided that as x. Now, e bw mx Φ(e b m, x,wx) (b m) x+1, [wx] 1 ( ) r x = o (b, m) x+1 (b m) x+1 [wx] 1 r x < b m{(bw mx) x +(bw mx) x + +(bw mx) x } x! = (bw mx) x+1 x! 4
5 Let the upper bound on the right-hand-side be bounded by some constant ǫ > 0, which can be made arbitrarily close to 0. Thus, (bw mx) x+1 /x! < ǫ, and multiplying both sides by x! and then taking logarithm, we get [(1 m+b m)x] < exp { 1 x+1 ( x )} logi+logǫ. (6) This implies that for some fixed ǫ, the approximation should be reasonably good as long as a,b,c and x are such that the inequality (6) is satisfied. For a,b,c such that m 1, eq.(3) can be approximated as (b m) x+1 ( 1 (w 1)x b m x i=1 ) (b m) x+1. Itfollows thatthedensity functionp(k X = x) (eq.(2)) canbeapproximated using the probability density function of the gamma distribution with mean (x+1)/(b m) and variance (x+1)/(b 2 m). Thus, mk P(k X = x) (b m) x+1 k(k +g(x)) x 1 e b k = k +g(x) (b m) x+1 (k +g(x)) x+1 1 e b mk (b m) x+1 k x+1 1 e b for k l, since g(x) = {(1 m)/(b m)}x and there exists some l > x, such that k +g(x) = k +o(k) when k l. mk, Corollary 1. For k l, where l > x, the probability mass function P(k X = x) can be approximated as P(k X = x) = k+0.5 (b m) x+1 t x+1 1 e b mt dt k {γ(x+1,b m(k +1/2)) γ(x+1,b m(k 1/2))}, where γ(u, v) is the lower incomplete gamma function: γ(u,v) = v 0 t u 1 e t dt. 5
6 2.1. Computational validation The preceding results establish that the gamma family of distribution is a valid approximation to eq.(2). Computational validation result suggests that an improvement to the fit can be obtained by replacing the shape and scale parameters of the gamma distribution in Theorem 1 as follows. Let µ post and σ 2 post be the posterior mean and the posterior variance, respectively. For a gamma distribution with shape parameter α and scale parameter β, its mean and variance are given by αβ and αβ 2, respectively. Given a,b,c and x, we can compute the exact mean and variance of the posterior distribution of k: µ post = σ 2 post = kp(k X = x), k=x k 2 P(k X = x) µ 2 post. k=x Simple algebra then yields α = µ 2 post /σ2 post and β = σ2 post /µ post. The gamma approximation with α and β thus computed fits eq.(2) better than α = x+1, β = 1/(b m) given in Theorem 1. Figure 1 shows the how well the gamma approximation fits the posterior distribution of k with and without adjustment to α and β parameters. In both cases the quality of the gamma approximation deteriorates when x becomes relatively large if the α and β parameters are not adjusted. 3. Acknowledgment I thank Martti Tammi and Joel Zi-Bin Low for stimulating the current work by discussing an applied problem in RNA-Seq count data analysis with me. Vanamamalai Seshadri read and provided feedback on the initial draft. 6
7 a) Density Observed count (X) x b) Density Observed count (X) x Figure 1: Examples of fitting the gamma distribution to the posterior distribution of k for two cases: a) relatively small b parameter: a = 1.5,b = 0.1,c = 0.05 (m = 1.1); b) relatively large b parameters: a = 1.5,b = 0.5,c = 0.05 (m = 2.0). Solid lines indicate gamma approximation with shape and scale parameters determined by matching moments of the posterior distribution of k and the gamma distribution. Broken lines indicate gamma approximation with shape and scale parameters given in Theorem 1. 7
8 References P. C. Consul, G. C. Jain, A generalization of the Poisson distribution., Technometrics 15 (1973) H. Joe, R. Zhu, Generalized Poisson distribution: the property of mixture of Poisson and comparison with negative binomial distribution, Biometrical Journal 47 (2005) A. K. Nikoloulopoulos, D. Karlis, On modelling count data: a comparison of some well-known discrete distributions, Journal of Statistical Computation and Simulation 78 (2008) M. M. Shoukri, P. C. Consul, Some chance mechanisms generating the Generalized Poisson probability models, McGraw-Hill, New York, P. C. Consul, Generalized Poisson Distribution: Properties and Applications, Marcel Dekker, New York, S. Srivastava, L. Chen, A two-parameter generalized Poisson model to improve the analysis of RNA-Seq data, Nucleic Acids Research 47 (2010) W. Li, T. Jiang, Transcriptome assembly and isoform expression level estimation from biased RNA-Seq reads, Bioinformatics 28 (2012) J. Zhang, C. C. J. Kuo, L. Chen, WemIQ: an accurate and robust isoform quantification method for RNA-seq data, Bioinformatics 31 (2015) Z. Wang, J. Wang, C. Wu, M. Deng, Estimation of isoform expression in RNA-seq data using a hierarchical Bayesian model, Journal of Bioinformatics and Computational Biology 13 (2015) A. Erdélyi, Higher Transcendental Function Vol. I, McGraw-Hill, New York,
Discrete Distributions Chapter 6
Discrete Distributions Chapter 6 Negative Binomial Distribution section 6.3 Consider k r, r +,... independent Bernoulli trials with probability of success in one trial being p. Let the random variable
More informationEVALUATING THE BIVARIATE COMPOUND GENERALIZED POISSON DISTRIBUTION
An. Şt. Univ. Ovidius Constanţa Vol. 92), 2001, 181 192 EVALUATING THE BIVARIATE COMPOUND GENERALIZED POISSON DISTRIBUTION Raluca Vernic Dedicated to Professor Mirela Ştefănescu on the occasion of her
More informationExponential Families
Exponential Families David M. Blei 1 Introduction We discuss the exponential family, a very flexible family of distributions. Most distributions that you have heard of are in the exponential family. Bernoulli,
More informationMAS113 Introduction to Probability and Statistics. Proofs of theorems
MAS113 Introduction to Probability and Statistics Proofs of theorems Theorem 1 De Morgan s Laws) See MAS110 Theorem 2 M1 By definition, B and A \ B are disjoint, and their union is A So, because m is a
More informationBell-shaped curves, variance
November 7, 2017 Pop-in lunch on Wednesday Pop-in lunch tomorrow, November 8, at high noon. Please join our group at the Faculty Club for lunch. Means If X is a random variable with PDF equal to f (x),
More informationProbability and Statistics
Kristel Van Steen, PhD 2 Montefiore Institute - Systems and Modeling GIGA - Bioinformatics ULg kristel.vansteen@ulg.ac.be Chapter 3: Parametric families of univariate distributions CHAPTER 3: PARAMETRIC
More informationarxiv: v1 [stat.me] 28 Mar 2011
arxiv:1103.5447v1 [stat.me] 28 Mar 2011 On matrix variance inequalities G. Afendras and N. Papadatos Department of Mathematics, Section of Statistics and O.R., University of Athens, Panepistemiopolis,
More informationMA/ST 810 Mathematical-Statistical Modeling and Analysis of Complex Systems
MA/ST 810 Mathematical-Statistical Modeling and Analysis of Complex Systems Review of Basic Probability The fundamentals, random variables, probability distributions Probability mass/density functions
More informationSTAT/MATH 395 A - PROBABILITY II UW Winter Quarter Moment functions. x r p X (x) (1) E[X r ] = x r f X (x) dx (2) (x E[X]) r p X (x) (3)
STAT/MATH 395 A - PROBABILITY II UW Winter Quarter 07 Néhémy Lim Moment functions Moments of a random variable Definition.. Let X be a rrv on probability space (Ω, A, P). For a given r N, E[X r ], if it
More informationContinuous Random Variables and Continuous Distributions
Continuous Random Variables and Continuous Distributions Continuous Random Variables and Continuous Distributions Expectation & Variance of Continuous Random Variables ( 5.2) The Uniform Random Variable
More informationSystem Simulation Part II: Mathematical and Statistical Models Chapter 5: Statistical Models
System Simulation Part II: Mathematical and Statistical Models Chapter 5: Statistical Models Fatih Cavdur fatihcavdur@uludag.edu.tr March 20, 2012 Introduction Introduction The world of the model-builder
More informationFINAL EXAM: Monday 8-10am
ECE 30: Probabilistic Methods in Electrical and Computer Engineering Fall 016 Instructor: Prof. A. R. Reibman FINAL EXAM: Monday 8-10am Fall 016, TTh 3-4:15pm (December 1, 016) This is a closed book exam.
More informationNon Uniform Bounds on Geometric Approximation Via Stein s Method and w-functions
Communications in Statistics Theory and Methods, 40: 45 58, 20 Copyright Taylor & Francis Group, LLC ISSN: 036-0926 print/532-45x online DOI: 0.080/036092090337778 Non Uniform Bounds on Geometric Approximation
More informationLecture 2: Repetition of probability theory and statistics
Algorithms for Uncertainty Quantification SS8, IN2345 Tobias Neckel Scientific Computing in Computer Science TUM Lecture 2: Repetition of probability theory and statistics Concept of Building Block: Prerequisites:
More informationSystem Simulation Part II: Mathematical and Statistical Models Chapter 5: Statistical Models
System Simulation Part II: Mathematical and Statistical Models Chapter 5: Statistical Models Fatih Cavdur fatihcavdur@uludag.edu.tr March 29, 2014 Introduction Introduction The world of the model-builder
More informationMath 3339 Homework 6 (Sections )
Math 3339 Homework 6 (Sections 5. 5.4) Name: Key PeopleSoft ID: Instructions: Homework will NOT be accepted through email or in person. Homework must be submitted through CourseWare BEFORE the deadline.
More informationCommon Discrete Distributions
Common Discrete Distributions Statistics 104 Autumn 2004 Taken from Statistics 110 Lecture Notes Copyright c 2004 by Mark E. Irwin Common Discrete Distributions There are a wide range of popular discrete
More informationarxiv: v2 [math.pr] 8 Feb 2016
Noname manuscript No will be inserted by the editor Bounds on Tail Probabilities in Exponential families Peter Harremoës arxiv:600579v [mathpr] 8 Feb 06 Received: date / Accepted: date Abstract In this
More informationMAS113 Introduction to Probability and Statistics. Proofs of theorems
MAS113 Introduction to Probability and Statistics Proofs of theorems Theorem 1 De Morgan s Laws) See MAS110 Theorem 2 M1 By definition, B and A \ B are disjoint, and their union is A So, because m is a
More informationPATTERN RECOGNITION AND MACHINE LEARNING CHAPTER 2: PROBABILITY DISTRIBUTIONS
PATTERN RECOGNITION AND MACHINE LEARNING CHAPTER 2: PROBABILITY DISTRIBUTIONS Parametric Distributions Basic building blocks: Need to determine given Representation: or? Recall Curve Fitting Binary Variables
More informationMixtures of Negative Binomial distributions for modelling overdispersion in RNA-Seq data
Mixtures of Negative Binomial distributions for modelling overdispersion in RNA-Seq data Cinzia Viroli 1 joint with E. Bonafede 1, S. Robin 2 & F. Picard 3 1 Department of Statistical Sciences, University
More informationClinical Trial Design Using A Stopped Negative Binomial Distribution. Michelle DeVeaux and Michael J. Kane and Daniel Zelterman
Statistics and Its Interface Volume 0 (0000) 1 10 Clinical Trial Design Using A Stopped Negative Binomial Distribution Michelle DeVeaux and Michael J. Kane and Daniel Zelterman arxiv:1508.01264v8 [math.st]
More informationIntroduction to Probability and Statistics (Continued)
Introduction to Probability and Statistics (Continued) Prof. icholas Zabaras Center for Informatics and Computational Science https://cics.nd.edu/ University of otre Dame otre Dame, Indiana, USA Email:
More informationIn-house germination methods validation studies: analysis
1 In-house germination methods validation studies: analysis Jean-Louis Laffont - ISTA Statistics Committee Design assumptions for the validation study Based on peer validation guidelines From Table 1 in
More informationTopic 3: The Expectation of a Random Variable
Topic 3: The Expectation of a Random Variable Course 003, 2017 Page 0 Expectation of a discrete random variable Definition (Expectation of a discrete r.v.): The expected value (also called the expectation
More informationSlides 8: Statistical Models in Simulation
Slides 8: Statistical Models in Simulation Purpose and Overview The world the model-builder sees is probabilistic rather than deterministic: Some statistical model might well describe the variations. An
More informationDispersion modeling for RNAseq differential analysis
Dispersion modeling for RNAseq differential analysis E. Bonafede 1, F. Picard 2, S. Robin 3, C. Viroli 1 ( 1 ) univ. Bologna, ( 3 ) CNRS/univ. Lyon I, ( 3 ) INRA/AgroParisTech, Paris IBC, Victoria, July
More informationLimits at Infinity. Horizontal Asymptotes. Definition (Limits at Infinity) Horizontal Asymptotes
Limits at Infinity If a function f has a domain that is unbounded, that is, one of the endpoints of its domain is ±, we can determine the long term behavior of the function using a it at infinity. Definition
More informationStat 5101 Lecture Notes
Stat 5101 Lecture Notes Charles J. Geyer Copyright 1998, 1999, 2000, 2001 by Charles J. Geyer May 7, 2001 ii Stat 5101 (Geyer) Course Notes Contents 1 Random Variables and Change of Variables 1 1.1 Random
More informationADDITIONAL MATHEMATICS
ADDITIONAL MATHEMATICS GCE NORMAL ACADEMIC LEVEL (016) (Syllabus 4044) CONTENTS Page INTRODUCTION AIMS ASSESSMENT OBJECTIVES SCHEME OF ASSESSMENT 3 USE OF CALCULATORS 3 SUBJECT CONTENT 4 MATHEMATICAL FORMULAE
More informationAlgebra I. Book 2. Powered by...
Algebra I Book 2 Powered by... ALGEBRA I Units 4-7 by The Algebra I Development Team ALGEBRA I UNIT 4 POWERS AND POLYNOMIALS......... 1 4.0 Review................ 2 4.1 Properties of Exponents..........
More informationChapter 3: Random Variables 1
Chapter 3: Random Variables 1 Yunghsiang S. Han Graduate Institute of Communication Engineering, National Taipei University Taiwan E-mail: yshan@mail.ntpu.edu.tw 1 Modified from the lecture notes by Prof.
More informationExam C Solutions Spring 2005
Exam C Solutions Spring 005 Question # The CDF is F( x) = 4 ( + x) Observation (x) F(x) compare to: Maximum difference 0. 0.58 0, 0. 0.58 0.7 0.880 0., 0.4 0.680 0.9 0.93 0.4, 0.6 0.53. 0.949 0.6, 0.8
More informationSeverity Models - Special Families of Distributions
Severity Models - Special Families of Distributions Sections 5.3-5.4 Stat 477 - Loss Models Sections 5.3-5.4 (Stat 477) Claim Severity Models Brian Hartman - BYU 1 / 1 Introduction Introduction Given that
More informationBinomial and Poisson Probability Distributions
Binomial and Poisson Probability Distributions Esra Akdeniz March 3, 2016 Bernoulli Random Variable Any random variable whose only possible values are 0 or 1 is called a Bernoulli random variable. What
More informationDiscrete Probability Refresher
ECE 1502 Information Theory Discrete Probability Refresher F. R. Kschischang Dept. of Electrical and Computer Engineering University of Toronto January 13, 1999 revised January 11, 2006 Probability theory
More informationExpectation. DS GA 1002 Statistical and Mathematical Models. Carlos Fernandez-Granda
Expectation DS GA 1002 Statistical and Mathematical Models http://www.cims.nyu.edu/~cfgranda/pages/dsga1002_fall16 Carlos Fernandez-Granda Aim Describe random variables with a few numbers: mean, variance,
More information3.4. ZEROS OF POLYNOMIAL FUNCTIONS
3.4. ZEROS OF POLYNOMIAL FUNCTIONS What You Should Learn Use the Fundamental Theorem of Algebra to determine the number of zeros of polynomial functions. Find rational zeros of polynomial functions. Find
More informationCLASS NOTES Models, Algorithms and Data: Introduction to computing 2018
CLASS NOTES Models, Algorithms and Data: Introduction to computing 208 Petros Koumoutsakos, Jens Honore Walther (Last update: June, 208) IMPORTANT DISCLAIMERS. REFERENCES: Much of the material (ideas,
More informationAlgorithms for Uncertainty Quantification
Algorithms for Uncertainty Quantification Tobias Neckel, Ionuț-Gabriel Farcaș Lehrstuhl Informatik V Summer Semester 2017 Lecture 2: Repetition of probability theory and statistics Example: coin flip Example
More informationThree hours. To be supplied by the Examinations Office: Mathematical Formula Tables and Statistical Tables THE UNIVERSITY OF MANCHESTER.
Three hours To be supplied by the Examinations Office: Mathematical Formula Tables and Statistical Tables THE UNIVERSITY OF MANCHESTER EXTREME VALUES AND FINANCIAL RISK Examiner: Answer QUESTION 1, QUESTION
More informationADDITIONAL MATHEMATICS
ADDITIONAL MATHEMATICS GCE Ordinary Level (06) (Syllabus 4047) CONTENTS Page INTRODUCTION AIMS ASSESSMENT OBJECTIVES SCHEME OF ASSESSMENT 3 USE OF CALCULATORS 3 SUBJECT CONTENT 4 MATHEMATICAL FORMULAE
More information3 Modeling Process Quality
3 Modeling Process Quality 3.1 Introduction Section 3.1 contains basic numerical and graphical methods. familiar with these methods. It is assumed the student is Goal: Review several discrete and continuous
More informationSTAT 414: Introduction to Probability Theory
STAT 414: Introduction to Probability Theory Spring 2016; Homework Assignments Latest updated on April 29, 2016 HW1 (Due on Jan. 21) Chapter 1 Problems 1, 8, 9, 10, 11, 18, 19, 26, 28, 30 Theoretical Exercises
More informationGenerating and characteristic functions. Generating and Characteristic Functions. Probability generating function. Probability generating function
Generating and characteristic functions Generating and Characteristic Functions September 3, 03 Probability generating function Moment generating function Power series expansion Characteristic function
More informationChapter 2. Random Variable. Define single random variables in terms of their PDF and CDF, and calculate moments such as the mean and variance.
Chapter 2 Random Variable CLO2 Define single random variables in terms of their PDF and CDF, and calculate moments such as the mean and variance. 1 1. Introduction In Chapter 1, we introduced the concept
More informationSTAT J535: Chapter 5: Classes of Bayesian Priors
STAT J535: Chapter 5: Classes of Bayesian Priors David B. Hitchcock E-Mail: hitchcock@stat.sc.edu Spring 2012 The Bayesian Prior A prior distribution must be specified in a Bayesian analysis. The choice
More informationBayesian Statistics Part III: Building Bayes Theorem Part IV: Prior Specification
Bayesian Statistics Part III: Building Bayes Theorem Part IV: Prior Specification Michael Anderson, PhD Hélène Carabin, DVM, PhD Department of Biostatistics and Epidemiology The University of Oklahoma
More informationExpectation. DS GA 1002 Probability and Statistics for Data Science. Carlos Fernandez-Granda
Expectation DS GA 1002 Probability and Statistics for Data Science http://www.cims.nyu.edu/~cfgranda/pages/dsga1002_fall17 Carlos Fernandez-Granda Aim Describe random variables with a few numbers: mean,
More informationPolynomial Functions
Polynomial Functions Polynomials A Polynomial in one variable, x, is an expression of the form a n x 0 a 1 x n 1... a n 2 x 2 a n 1 x a n The coefficients represent complex numbers (real or imaginary),
More informationMassachusetts Institute of Technology
Problem. (0 points) Massachusetts Institute of Technology Final Solutions: December 15, 009 (a) (5 points) We re given that the joint PDF is constant in the shaded region, and since the PDF must integrate
More informationThings to remember when learning probability distributions:
SPECIAL DISTRIBUTIONS Some distributions are special because they are useful They include: Poisson, exponential, Normal (Gaussian), Gamma, geometric, negative binomial, Binomial and hypergeometric distributions
More informationMathematical Induction Assignments
1 Mathematical Induction Assignments Prove the Following using Principle of Mathematical induction 1) Prove that for any positive integer number n, n 3 + 2 n is divisible by 3 2) Prove that 1 3 + 2 3 +
More information18.175: Lecture 8 Weak laws and moment-generating/characteristic functions
18.175: Lecture 8 Weak laws and moment-generating/characteristic functions Scott Sheffield MIT 18.175 Lecture 8 1 Outline Moment generating functions Weak law of large numbers: Markov/Chebyshev approach
More information4 Sums of Independent Random Variables
4 Sums of Independent Random Variables Standing Assumptions: Assume throughout this section that (,F,P) is a fixed probability space and that X 1, X 2, X 3,... are independent real-valued random variables
More informationPolyGamma Functions of Negative Order
Carnegie Mellon University Research Showcase @ CMU Computer Science Department School of Computer Science -998 PolyGamma Functions of Negative Order Victor S. Adamchik Carnegie Mellon University Follow
More informationDiscrete Distributions
Chapter 2 Discrete Distributions 2.1 Random Variables of the Discrete Type An outcome space S is difficult to study if the elements of S are not numbers. However, we can associate each element/outcome
More informationWeek 9 The Central Limit Theorem and Estimation Concepts
Week 9 and Estimation Concepts Week 9 and Estimation Concepts Week 9 Objectives 1 The Law of Large Numbers and the concept of consistency of averages are introduced. The condition of existence of the population
More informationA Note on Bootstraps and Robustness. Tony Lancaster, Brown University, December 2003.
A Note on Bootstraps and Robustness Tony Lancaster, Brown University, December 2003. In this note we consider several versions of the bootstrap and argue that it is helpful in explaining and thinking about
More informationEE/CpE 345. Modeling and Simulation. Fall Class 5 September 30, 2002
EE/CpE 345 Modeling and Simulation Class 5 September 30, 2002 Statistical Models in Simulation Real World phenomena of interest Sample phenomena select distribution Probabilistic, not deterministic Model
More informationChapters 3.2 Discrete distributions
Chapters 3.2 Discrete distributions In this section we study several discrete distributions and their properties. Here are a few, classified by their support S X. There are of course many, many more. For
More informationStat 135, Fall 2006 A. Adhikari HOMEWORK 6 SOLUTIONS
Stat 135, Fall 2006 A. Adhikari HOMEWORK 6 SOLUTIONS 1a. Under the null hypothesis X has the binomial (100,.5) distribution with E(X) = 50 and SE(X) = 5. So P ( X 50 > 10) is (approximately) two tails
More informationarxiv: v1 [stat.me] 29 Nov 2017
Extended Poisson INAR(1) processes with equidispersion, underdispersion and overdispersion arxiv:1711.10940v1 [stat.me] 29 Nov 2017 Marcelo Bourguignon 1, Josemar Rodrigues 2,, and Manoel Santos-Neto 3,4,
More informationMathematics Textbook Correlation to the 2016 Algebra I Standards of Learning and Curriculum Framework
and Curriculum Framework Publisher: McGraw-Hill School Education Text: Algebra 1 Copyright date 2018 A.1 The student will a) represent verbal quantitative situations algebraically; and TE: 5-9, 23-29,
More informationChapter 6: Rational Expr., Eq., and Functions Lecture notes Math 1010
Section 6.1: Rational Expressions and Functions Definition of a rational expression Let u and v be polynomials. The algebraic expression u v is a rational expression. The domain of this rational expression
More informationOutline. Binomial, Multinomial, Normal, Beta, Dirichlet. Posterior mean, MAP, credible interval, posterior distribution
Outline A short review on Bayesian analysis. Binomial, Multinomial, Normal, Beta, Dirichlet Posterior mean, MAP, credible interval, posterior distribution Gibbs sampling Revisit the Gaussian mixture model
More informationPetr Volf. Model for Difference of Two Series of Poisson-like Count Data
Petr Volf Institute of Information Theory and Automation Academy of Sciences of the Czech Republic Pod vodárenskou věží 4, 182 8 Praha 8 e-mail: volf@utia.cas.cz Model for Difference of Two Series of Poisson-like
More informationChapter 8.8.1: A factorization theorem
LECTURE 14 Chapter 8.8.1: A factorization theorem The characterization of a sufficient statistic in terms of the conditional distribution of the data given the statistic can be difficult to work with.
More informationMATHEMATICS. Higher 2 (Syllabus 9740)
MATHEMATICS Higher (Syllabus 9740) CONTENTS Page AIMS ASSESSMENT OBJECTIVES (AO) USE OF GRAPHING CALCULATOR (GC) 3 LIST OF FORMULAE 3 INTEGRATION AND APPLICATION 3 SCHEME OF EXAMINATION PAPERS 3 CONTENT
More informationMath Bootcamp 2012 Miscellaneous
Math Bootcamp 202 Miscellaneous Factorial, combination and permutation The factorial of a positive integer n denoted by n!, is the product of all positive integers less than or equal to n. Define 0! =.
More information1 Bernoulli Distribution: Single Coin Flip
STAT 350 - An Introduction to Statistics Named Discrete Distributions Jeremy Troisi Bernoulli Distribution: Single Coin Flip trial of an experiment that yields either a success or failure. X Bern(p),X
More informationCMPSCI 240: Reasoning Under Uncertainty
CMPSCI 240: Reasoning Under Uncertainty Lecture 7 Prof. Hanna Wallach wallach@cs.umass.edu February 14, 2012 Reminders Check the course website: http://www.cs.umass.edu/ ~wallach/courses/s12/cmpsci240/
More informationINFINITE SEQUENCES AND SERIES
11 INFINITE SEQUENCES AND SERIES INFINITE SEQUENCES AND SERIES In section 11.9, we were able to find power series representations for a certain restricted class of functions. INFINITE SEQUENCES AND SERIES
More informationMFM Practitioner Module: Quantitative Risk Management. John Dodson. September 23, 2015
MFM Practitioner Module: Quantitative Risk Management September 23, 2015 Mixtures Mixtures Mixtures Definitions For our purposes, A random variable is a quantity whose value is not known to us right now
More informationPUTNAM TRAINING POLYNOMIALS. Exercises 1. Find a polynomial with integral coefficients whose zeros include
PUTNAM TRAINING POLYNOMIALS (Last updated: December 11, 2017) Remark. This is a list of exercises on polynomials. Miguel A. Lerma Exercises 1. Find a polynomial with integral coefficients whose zeros include
More informationReview session Midterm 1
AS.110.109: Calculus II (Eng) Review session Midterm 1 Yi Wang, Johns Hopkins University Fall 2018 7.1: Integration by parts Basic integration method: u-sub, integration table Integration By Parts formula
More informationMathematical Statistics
Mathematical Statistics Chapter Three. Point Estimation 3.4 Uniformly Minimum Variance Unbiased Estimator(UMVUE) Criteria for Best Estimators MSE Criterion Let F = {p(x; θ) : θ Θ} be a parametric distribution
More informationStat410 Probability and Statistics II (F16)
Stat4 Probability and Statistics II (F6 Exponential, Poisson and Gamma Suppose on average every /λ hours, a Stochastic train arrives at the Random station. Further we assume the waiting time between two
More informationFinal Exam C Name i D) 2. Solve the equation by factoring. 4) x2 = x + 72 A) {1, 72} B) {-8, 9} C) {-8, -9} D) {8, 9} 9 ± i
Final Exam C Name First, write the value(s) that make the denominator(s) zero. Then solve the equation. 7 ) x + + 3 x - = 6 (x + )(x - ) ) A) No restrictions; {} B) x -, ; C) x -; {} D) x -, ; {2} Add
More informationMMSE Dimension. snr. 1 We use the following asymptotic notation: f(x) = O (g(x)) if and only
MMSE Dimension Yihong Wu Department of Electrical Engineering Princeton University Princeton, NJ 08544, USA Email: yihongwu@princeton.edu Sergio Verdú Department of Electrical Engineering Princeton University
More informationProbability Distributions Columns (a) through (d)
Discrete Probability Distributions Columns (a) through (d) Probability Mass Distribution Description Notes Notation or Density Function --------------------(PMF or PDF)-------------------- (a) (b) (c)
More informationON THE BIVARIATE GENERALIZED POISSON DISTRIBUTION. University "Ovidius" Constanta, Romania
ON THE BIVARIATE GENERALIZED POISSON DISTRIBUTION RALUCA VERNIC University "Ovidius" Constanta, Romania ABSTRACT This paper deals with the bivariate generalized Poisson distribution. The distribution is
More informationChapter 3, 4 Random Variables ENCS Probability and Stochastic Processes. Concordia University
Chapter 3, 4 Random Variables ENCS6161 - Probability and Stochastic Processes Concordia University ENCS6161 p.1/47 The Notion of a Random Variable A random variable X is a function that assigns a real
More informationPROBABILITY DISTRIBUTIONS. J. Elder CSE 6390/PSYC 6225 Computational Modeling of Visual Perception
PROBABILITY DISTRIBUTIONS Credits 2 These slides were sourced and/or modified from: Christopher Bishop, Microsoft UK Parametric Distributions 3 Basic building blocks: Need to determine given Representation:
More informationMarch Algebra 2 Question 1. March Algebra 2 Question 1
March Algebra 2 Question 1 If the statement is always true for the domain, assign that part a 3. If it is sometimes true, assign it a 2. If it is never true, assign it a 1. Your answer for this question
More informationPROBABILITY THEORY LECTURE 3
PROBABILITY THEORY LECTURE 3 Per Sidén Division of Statistics Dept. of Computer and Information Science Linköping University PER SIDÉN (STATISTICS, LIU) PROBABILITY THEORY - L3 1 / 15 OVERVIEW LECTURE
More informationWrite your Registration Number, Test Centre, Test Code and the Number of this booklet in the appropriate places on the answersheet.
2016 Booklet No. Test Code : PSA Forenoon Questions : 30 Time : 2 hours Write your Registration Number, Test Centre, Test Code and the Number of this booklet in the appropriate places on the answersheet.
More information13: Variational inference II
10-708: Probabilistic Graphical Models, Spring 2015 13: Variational inference II Lecturer: Eric P. Xing Scribes: Ronghuo Zheng, Zhiting Hu, Yuntian Deng 1 Introduction We started to talk about variational
More informationEXAMINATIONS OF THE ROYAL STATISTICAL SOCIETY
EXAMINATIONS OF THE ROYAL STATISTICAL SOCIETY GRADUATE DIPLOMA, 2016 MODULE 1 : Probability distributions Time allowed: Three hours Candidates should answer FIVE questions. All questions carry equal marks.
More informationMathematical Statistics 1 Math A 6330
Mathematical Statistics 1 Math A 6330 Chapter 3 Common Families of Distributions Mohamed I. Riffi Department of Mathematics Islamic University of Gaza September 28, 2015 Outline 1 Subjects of Lecture 04
More informationSome Statistics. V. Lindberg. May 16, 2007
Some Statistics V. Lindberg May 16, 2007 1 Go here for full details An excellent reference written by physicists with sample programs available is Data Reduction and Error Analysis for the Physical Sciences,
More informationChapter R - Basic Algebra Operations (94 topics, no due date)
Course Name: Math 00024 Course Code: N/A ALEKS Course: College Algebra Instructor: Master Templates Course Dates: Begin: 08/15/2014 End: 08/15/2015 Course Content: 207 topics Textbook: Barnett/Ziegler/Byleen/Sobecki:
More informationarxiv: v1 [cs.sc] 29 Jul 2009
An Explicit Construction of Gauss-Jordan Elimination Matrix arxiv:09075038v1 [cssc] 29 Jul 2009 Yi Li Laboratory of Computer Reasoning and Trustworthy Computing, University of Electronic Science and Technology
More informationIntroduction to Statistical Data Analysis Lecture 3: Probability Distributions
Introduction to Statistical Data Analysis Lecture 3: Probability Distributions James V. Lambers Department of Mathematics The University of Southern Mississippi James V. Lambers Statistical Data Analysis
More informationPage Max. Possible Points Total 100
Math 3215 Exam 2 Summer 2014 Instructor: Sal Barone Name: GT username: 1. No books or notes are allowed. 2. You may use ONLY NON-GRAPHING and NON-PROGRAMABLE scientific calculators. All other electronic
More informationExponential Distribution and Poisson Process
Exponential Distribution and Poisson Process Stochastic Processes - Lecture Notes Fatih Cavdur to accompany Introduction to Probability Models by Sheldon M. Ross Fall 215 Outline Introduction Exponential
More informationExplicit formulas for computing Bernoulli numbers of the second kind and Stirling numbers of the first kind
Filomat 28:2 (24), 39 327 DOI.2298/FIL4239O Published by Faculty of Sciences and Mathematics, University of Niš, Serbia Available at: http://www.pmf.ni.ac.rs/filomat Explicit formulas for computing Bernoulli
More informationTaylor and Maclaurin Series. Copyright Cengage Learning. All rights reserved.
11.10 Taylor and Maclaurin Series Copyright Cengage Learning. All rights reserved. We start by supposing that f is any function that can be represented by a power series f(x)= c 0 +c 1 (x a)+c 2 (x a)
More informationON THE TAYLOR COEFFICIENTS OF THE HURWITZ ZETA FUNCTION
ON THE TAYLOR COEFFICIENTS OF THE HURWITZ ZETA FUNCTION Khristo N. Boyadzhiev Department of Mathematics, Ohio Northern University, Ada, Ohio, 45810 k-boyadzhiev@onu.edu Abstract. We find a representation
More informationStatistics, Data Analysis, and Simulation SS 2013
Statistics, Data Analysis, and Simulation SS 213 8.128.73 Statistik, Datenanalyse und Simulation Dr. Michael O. Distler Mainz, 23. April 213 What we ve learned so far Fundamental
More information