Determining source cumulants in femtoscopy with Gram-Charlier and Edgeworth series

Size: px
Start display at page:

Download "Determining source cumulants in femtoscopy with Gram-Charlier and Edgeworth series"

Transcription

1 Determining source cumulants in femtoscopy with Gram-Charlier and Edgeworth series M.B. de Kock a H.C. Eggers a J. Schmiegel b a University of Stellenbosch, South Africa b Aarhus University, Denmark VI Workshop on Particle Correlations and Femtoscopy, September 010, Kiev

2 Outline 1 Objective 3 4

3 Nomenclature Measured correlation function C(q) = ρ(q) ρ ref (q) Interpret as a probability density function f (q) = C(q) 1 [C(q) 1 dq Assume noninteracting, identical particles C(q) 1 = d 3 xs(x)e iq x C( q) = C(q)

4 Objective: Describe the shape of the source S(x) in terms of the shape of the correlator C(q) using cumulants. C(q) 1 = d 3 xs(x)e iq x

5 Objective: Describe the shape of the source S(x) in terms of the shape of the correlator C(q) using cumulants. C(q) 1 = d 3 xs(x)e iq x κ (q) r Known κ (x) r Unknown

6 Objective: Describe the shape of the source S(x) in terms of the shape of the correlator C(q) using cumulants. C(q) 1 = d 3 xs(x)e iq x κ (q) r Known Gram-Charlier Edgeworth κ (x) r Unknown

7 How cumulants describe shape Mean κ 1 < 0 < κ 1 Variance κ < κ Skewness κ 3 < 0 < κ 3 Kurtosis κ 4 < κ 4

8 Useful properties of cumulants Translation invariant Additive under sums of random variables Simpler than moments, Gaussian cumulants: κ 1 = µ, κ = σ, κ r = 0 for all r 3 Shift and squeeze, q q µ σ Yields standardised cumulants: γ 1 = 0, γ = 1 γ r = κ r σ r Guassian as a baseline, with γ 4 0 the distance from it

9 Useful properties of cumulants Translation invariant Additive under sums of random variables Simpler than moments, Gaussian cumulants: κ 1 = µ, κ = σ, κ r = 0 for all r 3 Shift and squeeze, q q µ σ Yields standardised cumulants: γ 1 = 0, γ = 1 γ r = κ r σ r Guassian as a baseline, with γ 4 0 the distance from it

10 Useful properties of cumulants Translation invariant Additive under sums of random variables Simpler than moments, Gaussian cumulants: κ 1 = µ, κ = σ, κ r = 0 for all r 3 Shift and squeeze, q q µ σ Yields standardised cumulants: γ 1 = 0, γ = 1 γ r = κ r σ r Guassian as a baseline, with γ 4 0 the distance from it

11 Useful properties of cumulants Translation invariant Additive under sums of random variables Simpler than moments, Gaussian cumulants: κ 1 = µ, κ = σ, κ r = 0 for all r 3 Shift and squeeze, q q µ σ Yields standardised cumulants: γ 1 = 0, γ = 1 γ r = κ r σ r Guassian as a baseline, with γ 4 0 the distance from it

12 Measuring cumulants Defining moments µ (q) r = q r f (q)dq Cumulants in moments: κ (q) r κ (q) = µ(q) κ (q) 4 = µ(q) 4 3(µ(q) ) κ (q) 6 = µ(q) 6 15µ(q) 4 µ(q) + 30(µ(q) ) q Π q q Π q q 4 Π

13 Symmetry of cumulants and generating functions Generating function: q-cumulants κ (q) r = ( i) r d r log S(x) dx r Generating function: Source cumulants κ (x) r = ( i) r d r f (q) log dqr f (0) x=0 q= Easy to measure Hard to measure

14 Gram-Charlier expansion Series expansion: Gram-Charlier series f (q) = f 0 (q) c 1 f (1) 0 (q) + c! f () 0 (q) c 3 3! f (3) 0 (q) + Reference function: Gaussian Distribution f 0 (q) = 1 e q π Derivatives: Hermite functions ( ) d r 1 e q 1 = H r (q) e q dq π π

15 Gram-Charlier in pictures Gaussian + Second Hermite + Fourth Hermite + Sixth Hermite

16 Relating coefficients to cumulants Gram-Charlier with Gaussian reference function: f (q) = f 0 (q) c 4 4! H 4(q)f 0 (q) + c 6 6! H 6(q)f 0 (q) + Fourier Transform S(x) f (0) = x [ e 1 + c 4 4! (ix)4 + c 6 6! (ix)6 + c 8 8! (ix)8 + Expand in cumulants [ S(x) f (0) = exp κ (q)! (ix) + κ(q) 4 4! (ix)4 + κ(q) 6 6! (ix)6 +

17 Relating coefficients to cumulants Gram-Charlier with Gaussian reference function: f (q) = f 0 (q) c 4 4! H 4(q)f 0 (q) + c 6 6! H 6(q)f 0 (q) + Fourier Transform S(x) f (0) = x [ e 1 + c 4 4! (ix)4 + c 6 6! (ix)6 + c 8 8! (ix)8 + Expand in cumulants [ S(x) f (0) = exp κ (q)! (ix) + κ(q) 4 4! (ix)4 + κ(q) 6 6! (ix)6 +

18 Relating coefficients to cumulants Gram-Charlier with Gaussian reference function: f (q) = f 0 (q) c 4 4! H 4(q)f 0 (q) + c 6 6! H 6(q)f 0 (q) + Fourier Transform S(x) f (0) = x [ e 1 + c 4 4! (ix)4 + c 6 6! (ix)6 + c 8 8! (ix)8 + Expand in cumulants [ S(x) f (0) = exp κ (q)! (ix) + κ(q) 4 4! (ix)4 + κ(q) 6 6! (ix)6 +

19 Coefficients in terms of cumulants Equate series and compare powers in x m e x = e x [ 1 + c 4 4! (ix)4 + c 6 6! (ix)6 + c 8 8! (ix)8 + [ 1 + κ(q) 4 4! (ix)4 + κ(q) 6 6! (ix)6 + κ(q) 8 +35(κ(q) 4 ) 8! (ix) 8 For symmetrical distributions c 4 = κ (q) 4 c 6 = κ (q) 6 c 8 = κ (q) (κ(q) 4 ) c 10 = κ (q) κ(q) 6 κ(q) 4 c 1 = κ (q) κ(q) 4 κ(q) (κ(q) 6 ) (κ (q) 4 )3

20 Source cumulants in terms of q cumulants Connecting cumulants { κ (x) = 1 f () (q) f (q) κ (x) = ( i)r d r q=0 κ (q) r = 1 κ (q) dq log f (q) r f (0) [ 1+ 5!! 4! γ(q) 4 7!! 6! γ(q) 6 + 9!! 8! (γ(q) 8 +35(γ(q) 4 ) ) !! 4! γ(q) 4 5!! 6! γ(q) 6 + 7!! 8! (γ(q) 8 +35(γ(q) 4 ) )+... } q=0 Ratio of two infinite series [ !! 4! γ(q) 4 κ (q) 1 + 3!! 4! γ(q) 4 Where do we truncate x m, m = 4?

21 Source cumulants in terms of q cumulants Connecting cumulants { κ (x) = 1 f () (q) f (q) κ (x) = ( i)r d r q=0 κ (q) r = 1 κ (q) dq log f (q) r f (0) [ 1+ 5!! 4! γ(q) 4 7!! 6! γ(q) 6 + 9!! 8! (γ(q) 8 +35(γ(q) 4 ) ) !! 4! γ(q) 4 5!! 6! γ(q) 6 + 7!! 8! (γ(q) 8 +35(γ(q) 4 ) )+... } q=0 Ratio of two infinite series [ !! 4! γ(q) 4 7!! 6! γ(q) 6 κ (q) 1 + 3!! 4! γ(q) 4 5!! 6! γ(q) 6 Where do we truncate x m, m = 6?

22 Source cumulants in terms of q cumulants Connecting cumulants { κ (x) = 1 f () (q) f (q) κ (x) = ( i)r d r q=0 κ (q) r = 1 κ (q) dq log f (q) r f (0) [ 1+ 5!! 4! γ(q) 4 7!! 6! γ(q) 6 + 9!! 8! (γ(q) 8 +35(γ(q) 4 ) ) !! 4! γ(q) 4 5!! 6! γ(q) 6 + 7!! 8! (γ(q) 8 +35(γ(q) 4 ) )+... } q=0 Ratio of two infinite series 1 κ (q) [ 1 + 5!! 4! γ(q) 4 7!! 6! γ(q) 6 + 9!! 8! (γ(q) (γ(q) 4 ) ) 1 + 3!! 4! γ(q) 4 5!! 6! γ(q) 6 + 7!! 8! (γ(q) (γ(q) 4 ) ) Where do we truncate x m, m = 8?

23 Test with a Toy Model Test with Toy Model (Symmetrical Normal Inverse Gaussian): ) (q) 3e3/γ 4 K 1 (3γ (q) q f (q γ (q), γ(q) 4 ) = 3γ (q) 4 γ(q) (γ ) (q) π 4 q + 3γ (q) 4 γ(q) Adjustable kurtosis. Exact κ (q) 0.6, κ(q) 4, κ(q) 6, κ(x), κ(x) 4, κ(x)

24 Compare Toy model and Gram-Charlier Keep x m terms in series. Percentage Deviation κ (x) j (κ (x) j ) SNIG 1 Against q-kurtosis Gaussian limit Asymptotic series becomes worse Unexpected disaster 50%-100% error

25 Compare Toy model and Gram-Charlier Keep x m terms in series. Percentage Deviation κ (x) j (κ (x) j ) SNIG 1 Against q-kurtosis Gaussian limit Asymptotic series becomes worse Unexpected disaster 50%-100% error

26 How do we fix this? Cumulants of Toy model γ (q) r ( ) γ r (q) 1 4 Suggests truncation in γ4 r. Assume f (q) is n-divisible q = q 1 + q + + q n Product of n independent components ( x i S(x) = S σ n Resulting cumulants γ (q) r ( ) r 1 1 n ) n

27 Edgeworth series Gram-Charlier series Edgeworth series S(x) = exp S(x) = exp [ [ x + r=3 γ r (q) (ix) r r! x + γ r (q) (ix) r n r/ 1 r! r=3 Central Limit Theorem lim S(x) = exp [ x n

28 Edgeworth and Gram-Charlier S(x) f (0) = [ exp exp x + [ x + r=3 γ (q) r (ix) r r=3 r! γ (q) r (ix) r n r/ 1 r! = x m F m Gram-Charlier = n w C w Edgeworth Gram-Charlier Powers in x m γ 4 H 4 (q) 4! H γ 6 (q) 6 6! [ γ8 + 35γ 4 H8 (q) 8! Edgeworth Powers in n w [ 1 n γ 4 H 4 (q) 4! 1 [35γ H 8 (q) H n 4 8! + γ 6 (q) 6 6! [ 1 H γ 8 (q) H n 3 8 8! + 10γ 6 γ 10 (q) 4 10! γ4 3 H 1 (q) 1!

29 Compare Edgeworth series and Toy model Keep terms up to order n w Percentage Deviation κ (x) j (κ (x) j ) SNIG 1 Against q-kurtosis Edgeworth accuracy 1% Improving approximation Independent of the value of n

30 Summary Edgeworth series can relate correlation cumulants to source cumulants. Gram-Charlier cannot. Gram-Charlier is a orthogonal polynomial expansion. Simple problem and text book method is surprisingly complex. Reordering gives a large decrease in error. Edgeworth groups terms according to their approach to normality. Why this works is still unclear.

Applications of Good s Generalized Diversity Index. A. J. Baczkowski Department of Statistics, University of Leeds Leeds LS2 9JT, UK

Applications of Good s Generalized Diversity Index. A. J. Baczkowski Department of Statistics, University of Leeds Leeds LS2 9JT, UK Applications of Good s Generalized Diversity Index A. J. Baczkowski Department of Statistics, University of Leeds Leeds LS2 9JT, UK Internal Report STAT 98/11 September 1998 Applications of Good s Generalized

More information

arxiv: v1 [math.st] 18 Sep 2018

arxiv: v1 [math.st] 18 Sep 2018 Gram Charlier and Edgeworth expansion for sample variance arxiv:809.06668v [math.st] 8 Sep 08 Eric Benhamou,* A.I. SQUARE CONNECT, 35 Boulevard d Inkermann 900 Neuilly sur Seine, France and LAMSADE, Universit

More information

Application of Moment Expansion Method to Option Square Root Model

Application of Moment Expansion Method to Option Square Root Model Application of Moment Expansion Method to Option Square Root Model Yun Zhou Advisor: Professor Steve Heston University of Maryland May 5, 2009 1 / 19 Motivation Black-Scholes Model successfully explain

More information

The rth moment of a real-valued random variable X with density f(x) is. x r f(x) dx

The rth moment of a real-valued random variable X with density f(x) is. x r f(x) dx 1 Cumulants 1.1 Definition The rth moment of a real-valued random variable X with density f(x) is µ r = E(X r ) = x r f(x) dx for integer r = 0, 1,.... The value is assumed to be finite. Provided that

More information

Improved Inference for Signal Discovery Under Exceptionally Low False Positive Error Rates

Improved Inference for Signal Discovery Under Exceptionally Low False Positive Error Rates Improved Inference for Signal Discovery Under Exceptionally Low False Positive Error Rates (to appear in Journal of Instrumentation) Igor Volobouev & Alex Trindade Dept. of Physics & Astronomy, Texas Tech

More information

Gram-Charlierand Edgeworthexpansions for nongaussiancorrelations in femtoscopy

Gram-Charlierand Edgeworthexpansions for nongaussiancorrelations in femtoscopy Gram-Charlierand Edgeworthepansions or nongaussianorrelations in emtosopy Zimányi9 Winter Shool on eavy Ion Physis Mihiel de Kok University o Stellenbosh South Aria Eperimental Femtosopy Fireball Position

More information

6.2 Normal Distribution. Normal Distributions

6.2 Normal Distribution. Normal Distributions 6.2 Normal Distribution Normal Distributions 1 Homework Read Sec 6-1, and 6-2. Make sure you have a good feel for the normal curve. Do discussion question p302 2 3 Objective Identify Complete normal model

More information

Version A. Problem 1. Let X be the continuous random variable defined by the following pdf: 1 x/2 when 0 x 2, f(x) = 0 otherwise.

Version A. Problem 1. Let X be the continuous random variable defined by the following pdf: 1 x/2 when 0 x 2, f(x) = 0 otherwise. Math 224 Q Exam 3A Fall 217 Tues Dec 12 Version A Problem 1. Let X be the continuous random variable defined by the following pdf: { 1 x/2 when x 2, f(x) otherwise. (a) Compute the mean µ E[X]. E[X] x

More information

2.4 STATISTICAL FOUNDATIONS

2.4 STATISTICAL FOUNDATIONS 2.4 STATISTICAL FOUNDATIONS Characteristics of Return Distributions Moments of Return Distribution Correlation Standard Deviation & Variance Test for Normality of Distributions Time Series Return Volatility

More information

ROM SIMULATION Exact Moment Simulation using Random Orthogonal Matrices

ROM SIMULATION Exact Moment Simulation using Random Orthogonal Matrices ROM SIMULATION Exact Moment Simulation using Random Orthogonal Matrices Bachelier Finance Society Meeting Toronto 2010 Henley Business School at Reading Contact Author : d.ledermann@icmacentre.ac.uk Alexander

More information

Chapter 5. Continuous Random Variables and Probability Distributions. 5.1 Continuous Random Variables

Chapter 5. Continuous Random Variables and Probability Distributions. 5.1 Continuous Random Variables Chapter 5 Continuous Random Variables and Probability Distributions 5.1 Continuous Random Variables 1 2CHAPTER 5. CONTINUOUS RANDOM VARIABLES AND PROBABILITY DISTRIBUTIONS Probability Distributions Probability

More information

The Normal Distribution

The Normal Distribution 5.1 Introduction to Normal Distributions and the Standard Normal Distribution Section Learning objectives: 1. How to interpret graphs of normal probability distributions 2. How to find areas under the

More information

Probability theory: basic notions

Probability theory: basic notions 1 Probability theory: basic notions All epistemologic value of the theory of probability is based on this: that large scale random phenomena in their collective action create strict, non random regularity.

More information

Bivariate Birnbaum-Saunders Distribution

Bivariate Birnbaum-Saunders Distribution Department of Mathematics & Statistics Indian Institute of Technology Kanpur January 2nd. 2013 Outline 1 Collaborators 2 3 Birnbaum-Saunders Distribution: Introduction & Properties 4 5 Outline 1 Collaborators

More information

Edgeworth Binomial Trees

Edgeworth Binomial Trees Mark Rubinstein Paul Stephens Professor of Applied Investment Analysis University of California, Berkeley a version published in the Journal of Derivatives (Spring 1998) Abstract This paper develops a

More information

Reliability and Risk Analysis. Survival and Reliability Function

Reliability and Risk Analysis. Survival and Reliability Function Reliability and Risk Analysis Survival function We consider a non-negative random variable X which indicates the waiting time for the risk event (eg failure of the monitored equipment, etc.). The probability

More information

Testing the significance of the RV coefficient

Testing the significance of the RV coefficient 1 / 19 Testing the significance of the RV coefficient Application to napping data Julie Josse, François Husson and Jérôme Pagès Applied Mathematics Department Agrocampus Rennes, IRMAR CNRS UMR 6625 Agrostat

More information

INDIAN INSTITUTE OF SCIENCE STOCHASTIC HYDROLOGY. Lecture -5 Course Instructor : Prof. P. P. MUJUMDAR Department of Civil Engg., IISc.

INDIAN INSTITUTE OF SCIENCE STOCHASTIC HYDROLOGY. Lecture -5 Course Instructor : Prof. P. P. MUJUMDAR Department of Civil Engg., IISc. INDIAN INSTITUTE OF SCIENCE STOCHASTIC HYDROLOGY Lecture -5 Course Instructor : Prof. P. P. MUJUMDAR Department of Civil Engg., IISc. Summary of the previous lecture Moments of a distribubon Measures of

More information

A continuous random variable is one that can theoretically take on any value on some line interval. We use f ( x)

A continuous random variable is one that can theoretically take on any value on some line interval. We use f ( x) Section 6-2 I. Continuous Probability Distributions A continuous random variable is one that can theoretically take on any value on some line interval. We use f ( x) to represent a probability density

More information

Continuous Probability Distributions & Normal Distribution

Continuous Probability Distributions & Normal Distribution Mathematical Methods Units 3/4 Student Learning Plan Continuous Probability Distributions & Normal Distribution 7 lessons Notes: Students need practice in recognising whether a problem involves a discrete

More information

Prob and Stats, Nov 7

Prob and Stats, Nov 7 Prob and Stats, Nov 7 The Standard Normal Distribution Book Sections: 7.1, 7.2 Essential Questions: What is the standard normal distribution, how is it related to all other normal distributions, and how

More information

The statistics of polynomial roots

The statistics of polynomial roots The statistics of polynomial roots Herbert E. Müller Published in 14 on http://herbert-mueller.info/ Abstract A set of real numbers can be described summarily by its statistics mean value, standard deviation,

More information

Computation of one-sided probability density functions from their cumulants

Computation of one-sided probability density functions from their cumulants Journal of Mathematical Chemistry, Vol. 41, No. 1, January 27 26) DOI: 1.17/s191-6-969-x Computation of one-sided probability density functions from their cumulants Mário N. Berberan-Santos Centro de Química-Física

More information

Continuous Distributions

Continuous Distributions Quantitative Methods 2013 Continuous Distributions 1 The most important probability distribution in statistics is the normal distribution. Carl Friedrich Gauss (1777 1855) Normal curve A normal distribution

More information

Favorite Distributions

Favorite Distributions Favorite Distributions Binomial, Poisson and Normal Here we consider 3 favorite distributions in statistics: Binomial, discovered by James Bernoulli in 1700 Poisson, a limiting form of the Binomial, found

More information

FINITE SAMPLE DISTRIBUTIONS OF RISK-RETURN RATIOS

FINITE SAMPLE DISTRIBUTIONS OF RISK-RETURN RATIOS Available Online at ESci Journals Journal of Business and Finance ISSN: 305-185 (Online), 308-7714 (Print) http://www.escijournals.net/jbf FINITE SAMPLE DISTRIBUTIONS OF RISK-RETURN RATIOS Reza Habibi*

More information

Normal Probability Distributions

Normal Probability Distributions Normal Probability Distributions Properties of Normal Distributions The most important probability distribution in statistics is the normal distribution. Normal curve A normal distribution is a continuous

More information

1/12/2011. Chapter 5: z-scores: Location of Scores and Standardized Distributions. Introduction to z-scores. Introduction to z-scores cont.

1/12/2011. Chapter 5: z-scores: Location of Scores and Standardized Distributions. Introduction to z-scores. Introduction to z-scores cont. Chapter 5: z-scores: Location of Scores and Standardized Distributions Introduction to z-scores In the previous two chapters, we introduced the concepts of the mean and the standard deviation as methods

More information

I. Time Series and Stochastic Processes

I. Time Series and Stochastic Processes I. Time Series and Stochastic Processes Purpose of this Module Introduce time series analysis as a method for understanding real-world dynamic phenomena Define different types of time series Explain the

More information

Portfolios that Contain Risky Assets 12 Growth Rate Mean and Variance Estimators

Portfolios that Contain Risky Assets 12 Growth Rate Mean and Variance Estimators Portfolios that Contain Risky Assets 12 Growth Rate Mean and Variance Estimators C. David Levermore University of Maryland, College Park Math 420: Mathematical Modeling April 11, 2017 version c 2017 Charles

More information

Statistics 511 Supplemental Materials

Statistics 511 Supplemental Materials Gaussian (or Normal) Random Variable In this section we introduce the Gaussian Random Variable, which is more commonly referred to as the Normal Random Variable. This is a random variable that has a bellshaped

More information

Strategies for Improving the Efficiency of Monte-Carlo Methods

Strategies for Improving the Efficiency of Monte-Carlo Methods Strategies for Improving the Efficiency of Monte-Carlo Methods Paul J. Atzberger General comments or corrections should be sent to: paulatz@cims.nyu.edu Introduction The Monte-Carlo method is a useful

More information

Moments and Measures of Skewness and Kurtosis

Moments and Measures of Skewness and Kurtosis Moments and Measures of Skewness and Kurtosis Moments The term moment has been taken from physics. The term moment in statistical use is analogous to moments of forces in physics. In statistics the values

More information

Continuous random variables

Continuous random variables Continuous random variables probability density function (f(x)) the probability distribution function of a continuous random variable (analogous to the probability mass function for a discrete random variable),

More information

Probability and Random Variables A FINANCIAL TIMES COMPANY

Probability and Random Variables A FINANCIAL TIMES COMPANY Probability Basics Probability and Random Variables A FINANCIAL TIMES COMPANY 2 Probability Probability of union P[A [ B] =P[A]+P[B] P[A \ B] Conditional Probability A B P[A B] = Bayes Theorem P[A \ B]

More information

A New Multivariate Kurtosis and Its Asymptotic Distribution

A New Multivariate Kurtosis and Its Asymptotic Distribution A ew Multivariate Kurtosis and Its Asymptotic Distribution Chiaki Miyagawa 1 and Takashi Seo 1 Department of Mathematical Information Science, Graduate School of Science, Tokyo University of Science, Tokyo,

More information

VI. Continuous Probability Distributions

VI. Continuous Probability Distributions VI. Continuous Proaility Distriutions A. An Important Definition (reminder) Continuous Random Variale - a numerical description of the outcome of an experiment whose outcome can assume any numerical value

More information

Cambridge University Press Risk Modelling in General Insurance: From Principles to Practice Roger J. Gray and Susan M.

Cambridge University Press Risk Modelling in General Insurance: From Principles to Practice Roger J. Gray and Susan M. adjustment coefficient, 272 and Cramér Lundberg approximation, 302 existence, 279 and Lundberg s inequality, 272 numerical methods for, 303 properties, 272 and reinsurance (case study), 348 statistical

More information

12 The Bootstrap and why it works

12 The Bootstrap and why it works 12 he Bootstrap and why it works For a review of many applications of bootstrap see Efron and ibshirani (1994). For the theory behind the bootstrap see the books by Hall (1992), van der Waart (2000), Lahiri

More information

Technical Note: An Improved Range Chart for Normal and Long-Tailed Symmetrical Distributions

Technical Note: An Improved Range Chart for Normal and Long-Tailed Symmetrical Distributions Technical Note: An Improved Range Chart for Normal and Long-Tailed Symmetrical Distributions Pandu Tadikamalla, 1 Mihai Banciu, 1 Dana Popescu 2 1 Joseph M. Katz Graduate School of Business, University

More information

Modeling Obesity and S&P500 Using Normal Inverse Gaussian

Modeling Obesity and S&P500 Using Normal Inverse Gaussian Modeling Obesity and S&P500 Using Normal Inverse Gaussian Presented by Keith Resendes and Jorge Fernandes University of Massachusetts, Dartmouth August 16, 2012 Diabetes and Obesity Data Data obtained

More information

Financial Time Series and Their Characteristics

Financial Time Series and Their Characteristics Financial Time Series and Their Characteristics Egon Zakrajšek Division of Monetary Affairs Federal Reserve Board Summer School in Financial Mathematics Faculty of Mathematics & Physics University of Ljubljana

More information

UQ, STAT2201, 2017, Lectures 3 and 4 Unit 3 Probability Distributions.

UQ, STAT2201, 2017, Lectures 3 and 4 Unit 3 Probability Distributions. UQ, STAT2201, 2017, Lectures 3 and 4 Unit 3 Probability Distributions. Random Variables 2 A random variable X is a numerical (integer, real, complex, vector etc.) summary of the outcome of the random experiment.

More information

Section 7.5 The Normal Distribution. Section 7.6 Application of the Normal Distribution

Section 7.5 The Normal Distribution. Section 7.6 Application of the Normal Distribution Section 7.6 Application of the Normal Distribution A random variable that may take on infinitely many values is called a continuous random variable. A continuous probability distribution is defined by

More information

Data Analysis and Statistical Methods Statistics 651

Data Analysis and Statistical Methods Statistics 651 Data Analysis and Statistical Methods Statistics 651 http://www.stat.tamu.edu/~suhasini/teaching.html Suhasini Subba Rao The binomial: mean and variance Recall that the number of successes out of n, denoted

More information

Dependence Structure and Extreme Comovements in International Equity and Bond Markets

Dependence Structure and Extreme Comovements in International Equity and Bond Markets Dependence Structure and Extreme Comovements in International Equity and Bond Markets René Garcia Edhec Business School, Université de Montréal, CIRANO and CIREQ Georges Tsafack Suffolk University Measuring

More information

An Improved Saddlepoint Approximation Based on the Negative Binomial Distribution for the General Birth Process

An Improved Saddlepoint Approximation Based on the Negative Binomial Distribution for the General Birth Process Computational Statistics 17 (March 2002), 17 28. An Improved Saddlepoint Approximation Based on the Negative Binomial Distribution for the General Birth Process Gordon K. Smyth and Heather M. Podlich Department

More information

Saddlepoint Approximation Methods for Pricing. Financial Options on Discrete Realized Variance

Saddlepoint Approximation Methods for Pricing. Financial Options on Discrete Realized Variance Saddlepoint Approximation Methods for Pricing Financial Options on Discrete Realized Variance Yue Kuen KWOK Department of Mathematics Hong Kong University of Science and Technology Hong Kong * This is

More information

Approximating random inequalities with. Edgeworth expansions

Approximating random inequalities with. Edgeworth expansions Approximating random inequalities with Edgeworth expansions John D. Cook November 3, 2012 Random inequalities of the form Abstract Prob(X > Y + δ often appear as part of Bayesian clinical trial methods.

More information

Asymmetric fan chart a graphical representation of the inflation prediction risk

Asymmetric fan chart a graphical representation of the inflation prediction risk Asymmetric fan chart a graphical representation of the inflation prediction ASYMMETRIC DISTRIBUTION OF THE PREDICTION RISK The uncertainty of a prediction is related to the in the input assumptions for

More information

Appendix A (Pornprasertmanit & Little, in press) Mathematical Proof

Appendix A (Pornprasertmanit & Little, in press) Mathematical Proof Appendix A (Pornprasertmanit & Little, in press) Mathematical Proof Definition We begin by defining notations that are needed for later sections. First, we define moment as the mean of a random variable

More information

Financial Econometrics

Financial Econometrics Financial Econometrics Introduction to Financial Econometrics Gerald P. Dwyer Trinity College, Dublin January 2016 Outline 1 Set Notation Notation for returns 2 Summary statistics for distribution of data

More information

Basic Data Analysis. Stephen Turnbull Business Administration and Public Policy Lecture 4: May 2, Abstract

Basic Data Analysis. Stephen Turnbull Business Administration and Public Policy Lecture 4: May 2, Abstract Basic Data Analysis Stephen Turnbull Business Administration and Public Policy Lecture 4: May 2, 2013 Abstract Introduct the normal distribution. Introduce basic notions of uncertainty, probability, events,

More information

Approximation of probability density functions for SPDEs using truncated series expansions 0

Approximation of probability density functions for SPDEs using truncated series expansions 0 Approximation of probability density functions for SPDEs using truncated series expansions 0 Giacomo Capodaglio Max Gunzburger Henry P. Wynn Abstract The probability density function (PDF) of a random

More information

Two hours. To be supplied by the Examinations Office: Mathematical Formula Tables and Statistical Tables THE UNIVERSITY OF MANCHESTER

Two hours. To be supplied by the Examinations Office: Mathematical Formula Tables and Statistical Tables THE UNIVERSITY OF MANCHESTER Two hours MATH20802 To be supplied by the Examinations Office: Mathematical Formula Tables and Statistical Tables THE UNIVERSITY OF MANCHESTER STATISTICAL METHODS Answer any FOUR of the SIX questions.

More information

A New Test for Correlation on Bivariate Nonnormal Distributions

A New Test for Correlation on Bivariate Nonnormal Distributions Journal of Modern Applied Statistical Methods Volume 5 Issue Article 8 --06 A New Test for Correlation on Bivariate Nonnormal Distributions Ping Wang Great Basin College, ping.wang@gbcnv.edu Ping Sa University

More information

Terms & Characteristics

Terms & Characteristics NORMAL CURVE Knowledge that a variable is distributed normally can be helpful in drawing inferences as to how frequently certain observations are likely to occur. NORMAL CURVE A Normal distribution: Distribution

More information

Theoretical Problems in Credit Portfolio Modeling 2

Theoretical Problems in Credit Portfolio Modeling 2 Theoretical Problems in Credit Portfolio Modeling 2 David X. Li Shanghai Advanced Institute of Finance (SAIF) Shanghai Jiaotong University(SJTU) November 3, 2017 Presented at the University of South California

More information

Two Hours. Mathematical formula books and statistical tables are to be provided THE UNIVERSITY OF MANCHESTER. 22 January :00 16:00

Two Hours. Mathematical formula books and statistical tables are to be provided THE UNIVERSITY OF MANCHESTER. 22 January :00 16:00 Two Hours MATH38191 Mathematical formula books and statistical tables are to be provided THE UNIVERSITY OF MANCHESTER STATISTICAL MODELLING IN FINANCE 22 January 2015 14:00 16:00 Answer ALL TWO questions

More information

Quantitative Methods for Economics, Finance and Management (A86050 F86050)

Quantitative Methods for Economics, Finance and Management (A86050 F86050) Quantitative Methods for Economics, Finance and Management (A86050 F86050) Matteo Manera matteo.manera@unimib.it Marzio Galeotti marzio.galeotti@unimi.it 1 This material is taken and adapted from Guy Judge

More information

Convergence of statistical moments of particle density time series in scrape-off layer plasmas

Convergence of statistical moments of particle density time series in scrape-off layer plasmas Convergence of statistical moments of particle density time series in scrape-off layer plasmas R. Kube and O. E. Garcia Particle density fluctuations in the scrape-off layer of magnetically confined plasmas,

More information

Statistics for Business and Economics

Statistics for Business and Economics Statistics for Business and Economics Chapter 5 Continuous Random Variables and Probability Distributions Ch. 5-1 Probability Distributions Probability Distributions Ch. 4 Discrete Continuous Ch. 5 Probability

More information

High Dimensional Bayesian Optimisation and Bandits via Additive Models

High Dimensional Bayesian Optimisation and Bandits via Additive Models 1/20 High Dimensional Bayesian Optimisation and Bandits via Additive Models Kirthevasan Kandasamy, Jeff Schneider, Barnabás Póczos ICML 15 July 8 2015 2/20 Bandits & Optimisation Maximum Likelihood inference

More information

Lecture 5: Fundamentals of Statistical Analysis and Distributions Derived from Normal Distributions

Lecture 5: Fundamentals of Statistical Analysis and Distributions Derived from Normal Distributions Lecture 5: Fundamentals of Statistical Analysis and Distributions Derived from Normal Distributions ELE 525: Random Processes in Information Systems Hisashi Kobayashi Department of Electrical Engineering

More information

Statistical Analysis of Data from the Stock Markets. UiO-STK4510 Autumn 2015

Statistical Analysis of Data from the Stock Markets. UiO-STK4510 Autumn 2015 Statistical Analysis of Data from the Stock Markets UiO-STK4510 Autumn 2015 Sampling Conventions We observe the price process S of some stock (or stock index) at times ft i g i=0,...,n, we denote it by

More information

University of Texas, MD Anderson Cancer Center

University of Texas, MD Anderson Cancer Center University of Texas, MD Anderson Cancer Center UT MD Anderson Cancer Center Department of Biostatistics Working Paper Series Year 2012 Paper 78 Approximating random inequalities with Edgeworth expansions

More information

12. Conditional heteroscedastic models (ARCH) MA6622, Ernesto Mordecki, CityU, HK, 2006.

12. Conditional heteroscedastic models (ARCH) MA6622, Ernesto Mordecki, CityU, HK, 2006. 12. Conditional heteroscedastic models (ARCH) MA6622, Ernesto Mordecki, CityU, HK, 2006. References for this Lecture: Robert F. Engle. Autoregressive Conditional Heteroscedasticity with Estimates of Variance

More information

2.1 Properties of PDFs

2.1 Properties of PDFs 2.1 Properties of PDFs mode median epectation values moments mean variance skewness kurtosis 2.1: 1/13 Mode The mode is the most probable outcome. It is often given the symbol, µ ma. For a continuous random

More information

Business Statistics 41000: Probability 3

Business Statistics 41000: Probability 3 Business Statistics 41000: Probability 3 Drew D. Creal University of Chicago, Booth School of Business February 7 and 8, 2014 1 Class information Drew D. Creal Email: dcreal@chicagobooth.edu Office: 404

More information

Normal Inverse Gaussian (NIG) Process

Normal Inverse Gaussian (NIG) Process With Applications in Mathematical Finance The Mathematical and Computational Finance Laboratory - Lunch at the Lab March 26, 2009 1 Limitations of Gaussian Driven Processes Background and Definition IG

More information

Statistics for Managers Using Microsoft Excel 7 th Edition

Statistics for Managers Using Microsoft Excel 7 th Edition Statistics for Managers Using Microsoft Excel 7 th Edition Chapter 7 Sampling Distributions Statistics for Managers Using Microsoft Excel 7e Copyright 2014 Pearson Education, Inc. Chap 7-1 Learning Objectives

More information

Random Variables Handout. Xavier Vilà

Random Variables Handout. Xavier Vilà Random Variables Handout Xavier Vilà Course 2004-2005 1 Discrete Random Variables. 1.1 Introduction 1.1.1 Definition of Random Variable A random variable X is a function that maps each possible outcome

More information

Dependence Modeling and Credit Risk

Dependence Modeling and Credit Risk Dependence Modeling and Credit Risk Paola Mosconi Banca IMI Bocconi University, 20/04/2015 Paola Mosconi Lecture 6 1 / 53 Disclaimer The opinion expressed here are solely those of the author and do not

More information

CH 5 Normal Probability Distributions Properties of the Normal Distribution

CH 5 Normal Probability Distributions Properties of the Normal Distribution Properties of the Normal Distribution Example A friend that is always late. Let X represent the amount of minutes that pass from the moment you are suppose to meet your friend until the moment your friend

More information

SAQ KONTROLL AB Box 49306, STOCKHOLM, Sweden Tel: ; Fax:

SAQ KONTROLL AB Box 49306, STOCKHOLM, Sweden Tel: ; Fax: ProSINTAP - A Probabilistic Program for Safety Evaluation Peter Dillström SAQ / SINTAP / 09 SAQ KONTROLL AB Box 49306, 100 29 STOCKHOLM, Sweden Tel: +46 8 617 40 00; Fax: +46 8 651 70 43 June 1999 Page

More information

Math 227 Elementary Statistics. Bluman 5 th edition

Math 227 Elementary Statistics. Bluman 5 th edition Math 227 Elementary Statistics Bluman 5 th edition CHAPTER 6 The Normal Distribution 2 Objectives Identify distributions as symmetrical or skewed. Identify the properties of the normal distribution. Find

More information

Box-Cox Transforms for Realized Volatility

Box-Cox Transforms for Realized Volatility Box-Cox Transforms for Realized Volatility Sílvia Gonçalves and Nour Meddahi Université de Montréal and Imperial College London January 1, 8 Abstract The log transformation of realized volatility is often

More information

Since his score is positive, he s above average. Since his score is not close to zero, his score is unusual.

Since his score is positive, he s above average. Since his score is not close to zero, his score is unusual. Chapter 06: The Standard Deviation as a Ruler and the Normal Model This is the worst chapter title ever! This chapter is about the most important random variable distribution of them all the normal distribution.

More information

NORMAL RANDOM VARIABLES (Normal or gaussian distribution)

NORMAL RANDOM VARIABLES (Normal or gaussian distribution) NORMAL RANDOM VARIABLES (Normal or gaussian distribution) Many variables, as pregnancy lengths, foot sizes etc.. exhibit a normal distribution. The shape of the distribution is a symmetric bell shape.

More information

Random variables. Contents

Random variables. Contents Random variables Contents 1 Random Variable 2 1.1 Discrete Random Variable............................ 3 1.2 Continuous Random Variable........................... 5 1.3 Measures of Location...............................

More information

Chapter 3 - Lecture 4 Moments and Moment Generating Funct

Chapter 3 - Lecture 4 Moments and Moment Generating Funct Chapter 3 - Lecture 4 and s October 7th, 2009 Chapter 3 - Lecture 4 and Moment Generating Funct Central Skewness Chapter 3 - Lecture 4 and Moment Generating Funct Central Skewness The expected value of

More information

Asymptotic refinements of bootstrap tests in a linear regression model ; A CHM bootstrap using the first four moments of the residuals

Asymptotic refinements of bootstrap tests in a linear regression model ; A CHM bootstrap using the first four moments of the residuals Asymptotic refinements of bootstrap tests in a linear regression model ; A CHM bootstrap using the first four moments of the residuals Pierre-Eric Treyens To cite this version: Pierre-Eric Treyens. Asymptotic

More information

Standard Normal, Inverse Normal and Sampling Distributions

Standard Normal, Inverse Normal and Sampling Distributions Standard Normal, Inverse Normal and Sampling Distributions Section 5.5 & 6.6 Cathy Poliak, Ph.D. cathy@math.uh.edu Office in Fleming 11c Department of Mathematics University of Houston Lecture 9-3339 Cathy

More information

On multivariate Multi-Resolution Analysis, using generalized (non homogeneous) polyharmonic splines. or: A way for deriving RBF and associated MRA

On multivariate Multi-Resolution Analysis, using generalized (non homogeneous) polyharmonic splines. or: A way for deriving RBF and associated MRA MAIA conference Erice (Italy), September 6, 3 On multivariate Multi-Resolution Analysis, using generalized (non homogeneous) polyharmonic splines or: A way for deriving RBF and associated MRA Christophe

More information

A market risk model for asymmetric distributed series of return

A market risk model for asymmetric distributed series of return University of Wollongong Research Online University of Wollongong in Dubai - Papers University of Wollongong in Dubai 2012 A market risk model for asymmetric distributed series of return Kostas Giannopoulos

More information

Week 1 Quantitative Analysis of Financial Markets Distributions B

Week 1 Quantitative Analysis of Financial Markets Distributions B Week 1 Quantitative Analysis of Financial Markets Distributions B Christopher Ting http://www.mysmu.edu/faculty/christophert/ Christopher Ting : christopherting@smu.edu.sg : 6828 0364 : LKCSB 5036 October

More information

The mean-variance portfolio choice framework and its generalizations

The mean-variance portfolio choice framework and its generalizations The mean-variance portfolio choice framework and its generalizations Prof. Massimo Guidolin 20135 Theory of Finance, Part I (Sept. October) Fall 2014 Outline and objectives The backward, three-step solution

More information

Frequency Distribution Models 1- Probability Density Function (PDF)

Frequency Distribution Models 1- Probability Density Function (PDF) Models 1- Probability Density Function (PDF) What is a PDF model? A mathematical equation that describes the frequency curve or probability distribution of a data set. Why modeling? It represents and summarizes

More information

PORTFOLIO THEORY. Master in Finance INVESTMENTS. Szabolcs Sebestyén

PORTFOLIO THEORY. Master in Finance INVESTMENTS. Szabolcs Sebestyén PORTFOLIO THEORY Szabolcs Sebestyén szabolcs.sebestyen@iscte.pt Master in Finance INVESTMENTS Sebestyén (ISCTE-IUL) Portfolio Theory Investments 1 / 60 Outline 1 Modern Portfolio Theory Introduction Mean-Variance

More information

ECON Introductory Econometrics. Lecture 1: Introduction and Review of Statistics

ECON Introductory Econometrics. Lecture 1: Introduction and Review of Statistics ECON4150 - Introductory Econometrics Lecture 1: Introduction and Review of Statistics Monique de Haan (moniqued@econ.uio.no) Stock and Watson Chapter 1-2 Lecture outline 2 What is econometrics? Course

More information

Hypothesis Tests: One Sample Mean Cal State Northridge Ψ320 Andrew Ainsworth PhD

Hypothesis Tests: One Sample Mean Cal State Northridge Ψ320 Andrew Ainsworth PhD Hypothesis Tests: One Sample Mean Cal State Northridge Ψ320 Andrew Ainsworth PhD MAJOR POINTS Sampling distribution of the mean revisited Testing hypotheses: sigma known An example Testing hypotheses:

More information

Two-term Edgeworth expansions of the distributions of fit indexes under fixed alternatives in covariance structure models

Two-term Edgeworth expansions of the distributions of fit indexes under fixed alternatives in covariance structure models Economic Review (Otaru University of Commerce), Vo.59, No.4, 4-48, March, 009 Two-term Edgeworth expansions of the distributions of fit indexes under fixed alternatives in covariance structure models Haruhiko

More information

Probability Weighted Moments. Andrew Smith

Probability Weighted Moments. Andrew Smith Probability Weighted Moments Andrew Smith andrewdsmith8@deloitte.co.uk 28 November 2014 Introduction If I asked you to summarise a data set, or fit a distribution You d probably calculate the mean and

More information

Chapter Seven. The Normal Distribution

Chapter Seven. The Normal Distribution Chapter Seven The Normal Distribution 7-1 Introduction Many continuous variables have distributions that are bellshaped and are called approximately normally distributed variables, such as the heights

More information

Diploma in Business Administration Part 2. Quantitative Methods. Examiner s Suggested Answers

Diploma in Business Administration Part 2. Quantitative Methods. Examiner s Suggested Answers Cumulative frequency Diploma in Business Administration Part Quantitative Methods Examiner s Suggested Answers Question 1 Cumulative Frequency Curve 1 9 8 7 6 5 4 3 1 5 1 15 5 3 35 4 45 Weeks 1 (b) x f

More information

Dynamic Response of Jackup Units Re-evaluation of SNAME 5-5A Four Methods

Dynamic Response of Jackup Units Re-evaluation of SNAME 5-5A Four Methods ISOPE 2010 Conference Beijing, China 24 June 2010 Dynamic Response of Jackup Units Re-evaluation of SNAME 5-5A Four Methods Xi Ying Zhang, Zhi Ping Cheng, Jer-Fang Wu and Chee Chow Kei ABS 1 Main Contents

More information

2011 Pearson Education, Inc

2011 Pearson Education, Inc Statistics for Business and Economics Chapter 4 Random Variables & Probability Distributions Content 1. Two Types of Random Variables 2. Probability Distributions for Discrete Random Variables 3. The Binomial

More information

Statistical Analysis of Sample-Size Effects in ICA

Statistical Analysis of Sample-Size Effects in ICA Statistical Analysis of Sample-Size Effects in ICA J. Michael Herrmann 1,2 Fabian J. Theis 1,3 1 Bernstein Center for Computational Neuroscience Göttingen 2 Göttingen University, Institute for Nonlinear

More information

SADDLEPOINT APPROXIMATIONS TO OPTION PRICES 1. By L. C. G. Rogers and O. Zane University of Bath and First Chicago NBD

SADDLEPOINT APPROXIMATIONS TO OPTION PRICES 1. By L. C. G. Rogers and O. Zane University of Bath and First Chicago NBD The Annals of Applied Probability 1999, Vol. 9, No. 2, 493 53 SADDLEPOINT APPROXIMATIONS TO OPTION PRICES 1 By L. C. G. Rogers and O. Zane University of Bath and First Chicago NBD The use of saddlepoint

More information

Lecture 8. The Binomial Distribution. Binomial Distribution. Binomial Distribution. Probability Distributions: Normal and Binomial

Lecture 8. The Binomial Distribution. Binomial Distribution. Binomial Distribution. Probability Distributions: Normal and Binomial Lecture 8 The Binomial Distribution Probability Distributions: Normal and Binomial 1 2 Binomial Distribution >A binomial experiment possesses the following properties. The experiment consists of a fixed

More information