Multivariate Normal-Laplace Distribution and Processes
|
|
- Raymond Baldwin
- 5 years ago
- Views:
Transcription
1 CHAPTER 4 Multivariate Normal-Laplace Distribution and Processes The normal-laplace distribution, which results from the convolution of independent normal and Laplace random variables is introduced by Reed and Jorgensen (2004. It is the distribution of the stopped state of a Brownian motion with normally distributed starting value if the stopping hazard rate is constant. Normal-Laplace distribution is a new distribution which (in its symmetric form behaves somewhat like the normal distribution in the middle of its range, and like the Laplace distribution in its tails. Reed and Jorgensen (2004 also introduced a generalized normal-laplace distribution, which is useful in financial applications for obtaining an alternative stochastic process model to Brownian motion for logarithmic prices, in which the increments exhibit fatter tails than the normal distribution. Reed (2007 developed Brownian-Laplace motion for modelling financial asset price returns. Multivariate normal distribution, a generalization of univariate normal distribution is 74
2 studied by various authors. They have applications in statistical inference, image processing etc. Ernst (998 introduced multivariate extension of symmetric Laplace distributions via an elliptic contouring. Many properties in the univariate laws can be extended to this class of distributions. With an appropriate limit of the parameters of multivariate hyperbolic distributions, one can obtain a multivariate and asymmetric extension of the Laplace laws, see Blaesid (98. In the present chapter, we consider the univariate normal-laplace distribution, introduce multivariate normal-laplace distribution and study its properties. First order autoregressive processes with multivariate normal-laplace marginals is developed. We introduce multivariate geometric normal-laplace distribution. Multivariate generalized normal- Laplace distribution is introduced. We introduce the geometric generalized normal-laplace distribution and study its properties. We consider the estimation of parameters. Some applications are also discussed. 4. Normal-Laplace Distribution and its Properties The normal-laplace distribution introduced by Reed and Jorgensen (2004, arises as the convolution of independent normal and asymmetric Laplace densities. A normal-laplace random variable X with parameters µ, σ 2, α and β can be represented as X d = Z + W, (4.. where Z and W are independent random variables with Z following normal distribution with mean µ and variance σ 2 and W following an asymmetric Laplace distribution with parameters α, β. The corresponding normal-laplace distribution shall be denoted by NL(µ, σ 2, α, β. Various results on normal-laplace distribution are available in Reed (2006. The characteristic function (c.f. of NL (µ, σ 2, α, β can be obtained as the product of the c.f. s of its normal and Laplace components and is given by Reed and Jorgensen (
3 as, φ X (t = ] [ [exp(iµt σ2 2 t2 αβ (α it(β + it ]. (4..2 The normal-laplace distribution is infinitely divisible and is closed under linear transformation. The mean, variance and cumulants exist for the distribution. 4.2 Multivariate Normal Laplace Distribution A multivariate extension of the normal-laplace distribution of Reed and Jorgensen (2004, namely the multivariate normal-laplace distribution can be obtained as the convolution of multivariate normal and multivariate symmetric Laplace random vectors. The c.f. of multivariate normal-laplace distribution is given by ϕ X (t = exp it µ ( + 2 t V t, t R p, Σ > 0, V > 0. (4.2. Some properties A p-variate normal-laplace distribution with parameters µ, Σ and V can be denoted by NL p (µ, Σ, V. Let X NL p (µ, Σ, V, then X can be expressed as X d = Z + Y, (4.2.2 where Z and Y are independent random vectors with Z following a p-variate normal distribution with mean vector (µ and dispersion matrix Σ (N p (µ, Σ and Y following a p-variate symmetric Laplace distribution with parameter V (L p (V. Another representation is X d = Z + W Y, (4.2.3 where Z follows a p-variate normal distribution with mean vector (µ and dispersion matrix Σ, W is a standard exponential variable and Y follows p-variate normal distribution with mean vector 0 and dispersion matrix Σ. 76
4 Moments. The k th moment M k (x, µ of a random vector x and µ is defined by M k (x, µ = E[(x µ (x µ (x µ (x µ], if k = 2m and M k (x, µ = E[(x µ (x µ (x µ (x µ ], if k = 2m. The Kronecker product of matrices A = (a ij : m n and B : p q is a mp nq matrix A B = [a ij B]. When µ = 0, we write only M k (x. By the differentiation of the c.f. of a random vector, we can obtain, M k = k ϕ X (t i k t t... The mean vector and variance-covariance matrix of NL p (µ, Σ, V can be obtained as t=0. E(X = µ Cov(X = Σ + V. (4.2.4 Cumulants. The k th cumulant of random vector x is denoted by C k (x. The cumulants are obtained as the matrix derivatives of the function φ X (t = ln ϕ X (t by C k = k φ X (t i k t t... For multivariate normal-laplace distribution, t=0. C = µ C 2 = Σ + V. 77
5 Skewness and Kurtosis. Let x be a random vector with mean vector θ and covariance matrix Σ and y = Σ /2 (x θ with mean vector 0, covariance matrix I p and with third and fourth moments M 3 (y, M 4 (y and cumulants C 3 (y, C 4 (y. Then skewness measure β p (x and kurtosis characteristic β 2p (x are defined by β p (x = tr[c 3(yC 3 (y] = tr[m 3(yM 3 (y] and β 2p (x = tr[m 4 (y] = tr[c 4 (y] + p 2 + 2p The class of elliptical distributions. The class of elliptical distributions, introduced by Kelker (970, is a generalization of multivariate normal distributions. These distributions are symmetric and may not adequately represent the data when some asymmetry is present. Definition The random vector X has a multivariate elliptical distribution, if its c.f. can be expressed as ϕ X (t = exp(it µψ( (4.2.5 for some column vector µ, n n positive matrix Σ and for some function ψ(t ψ n, which is called the characteristic generator. The multivariate normal and symmetric Laplace distributions belong to elliptical family, since the c.f. of multivariate normal and symmetric Laplace distributions can be factorized as ( The multivariate normal-laplace distribution belongs to the class of elliptical distributions. It can be easily followed from the result that the sum of elliptical distributions is elliptical, see Fang et al. (987. This property is very important when we deal with portfolio of assets represented by sum. Marginal distributions of elliptical distributions are also elliptical distributions. So the 78
6 marginal distributions of multivariate normal-laplace distributions are also elliptical. Also all odd order moments of an elliptical distribution are zero and hence the result holds for multivariate normal-laplace distribution. Infinite divisibility. The multivariate normal-laplace distribution is infinitely divisible. Since the c.f. of NL p (µ, Σ, V can be written as ϕ X (t = [e (it µn 2 t Σn t ( ] n n +, 2 t V t for any integer n > 0. The term in brackets is the c.f. of a random vector expressed as Z + Y, where Z N p ( µ, Σ and Y following a multivariate generalized symmetric Laplace n n distribution with parameters V, n (L p(v, n. Theorem If X NL p (µ, Σ, V and Y = AX + b, where A is a p p matrix and b R p, then Y NL p (Aµ + b, AΣA, AV A. Proof. The c.f. of Y is ϕ Y (t = E Hence Y NL p (Aµ + b, AΣA, AV A. ( e it (AX+b = e it b ϕ X (A t (Aµ+b = eit t AΣA t AV A. t Theorem If X NL p (µ, Σ, V, then a Xa R p follows a univariate normal- Laplace distribution. Proof. This theorem can also be easily be proved using the c.f. 79
7 Theorem Let X NL p (µ, Σ, V and partition X, µ, Σ and V as X = X X 2 µ = µ µ 2 Σ = Σ Σ 2 Σ 2 Σ 22 V = V V 2 V 2 V 22 where X and µ are k vectors and Σ and V are k k matrices. Then X NL p (µ, Σ, V and X NL p (µ 2, Σ 22, V 22. When Σ 2 = 0, X and X 2 are independently distributed. Proof. By using the joint c.f., we can easily prove the result. Definition A random variable X with c.f. φ is said to be semi-self decomposable, if for some 0 < a <, there exists a c.f. φ a such that φ(t = φ X (at φ a (t, t R. If this relation holds for every 0 < a <, then φ is self-decomposable or the corresponding distribution belongs to L class. Theorem Multivariate normal-laplace c.f. ϕ is self decomposable or the corresponding distribution belongs to the class L. Proof. The c.f. of multivariate normal-laplace distribution is ϕ(t = = exp it µ ( exp it µ + 2 t V t ( + 2 t V t exp it ( aµ 2 t ( a 2 Σt = ϕ X (atϕ a (t. [ a 2 + ( a t V t, ]. The second term in the expression is also a c.f., which can be seen in (
8 4.3 Multivariate normal-laplace processes Consider the usual linear, additive first order autoregressive model given by X n = ax n + ε n, 0 < a, n = 0, ±, ±2,..., (4.3. where X n and innovations ε n are independent p- variate random vectors. In terms of c.f., we have The c.f. of {ε n } can be obtained as ϕ ɛ (t = ϕ X (t = ϕ X (atϕ ε (t. exp it ( aµ [ 2 t ( a 2 Σt a 2 + ( a 2 From this, we can obtain the distribution of the innovation sequence as + 2 t V t ]. (4.3.2 ε d = Z + L. where Z N p (( aµ, ( a 2 Σ and L can be treated as a sequence of random vectors of the form 0, with probability a 2 L = L p, with probability ( a 2, where L p s are independently and identically distributed (i.i.d. Laplace random vectors. symmetric multivariate Remark The process is stationary with NL p (µ, Σ, V marginals. Proof. We can prove this by the method of induction. We assume that X n NL p (µ, Σ, V. 8
9 Then ϕ Xn (t = ϕ Xn (atϕ εn (t = exp iat µ ( 2 a2 t Σt + 2 a2 t V t exp it ( aµ [ ] 2 t ( a 2 Σt a 2 + ( a t V t = exp it µ ( + 2 t V t which is the c.f. of NL p (µ, Σ, V. Therefore {X n } is strictly stationary with NL p (µ, Σ, V marginals. Remark If X 0 is distributed arbitrary, then also the process is asymptotically Markovian with multivariate normal-laplace marginal distribution. Proof. We have In terms of c.f., we get n X n = ax n + ε n = a n X 0 + a k ε n k. n ϕ Xn (t = ϕ X0 (a n t ϕ ε (a k t k=0 n = ϕ X0 (a n exp ( it a k 2 t a2k t Σt + 2 a2(k+ t V t + k=0 2 a2k t V t exp ( it a k+ 2 a2(k+ t Σt exp it µ ( + as n. 2 t V t Hence even if X 0 is arbitrarily distributed, the process is asymptotically stationary Markovian with multivariate normal-laplace marginals. k=0 82
10 4.3. Distribution of sums and joint distribution of (X n, X n+ Consider a stationary sequence {X n } satisfying (4.3.. Then we have X n+j = a j X n + a j ε n+ + a j 2 ε n ε n+j. Hence The c.f. of T r is given by T r = X n + X n+ + + X n+r r = [a j X n + a j ε n+ + + ε n+j ] j=0 = X n ( a r a r + j= ( a r j ε n+j. a ( ϕ Tr (t = ϕ Xn t r ( ar ϕ ε t ar j a a j= [ a r = exp it µ ( ( a r 2 t Σt] ( a 2 a + a r 2 2 a t V t r { exp it ( a r j µ 2 ( ar j ( + at Σt j= [ a 2 + ( a ( a r j a 2 t V t The distribution of T r can be obtained by inverting the above expression. The joint distribution of contiguous observation vectors (X n, X n+ can be given in ]}. 83
11 terms of c.f. as ϕ Xn,X n+ (t, t 2 = E[exp(it X n + it 2X n+ ] = E[exp(it X n + it 2(aX n + ε n+ ] = E[exp(i(t + at 2X n + it 2ε n+ ] = ϕ Xn (t + at 2 ϕ εn+ (t 2 = exp i(t + at 2 µ 2 (t + at 2 Σ(t + at 2 ( + 2 (t + at 2 V (t + at 2 [ ] a 2 + ( a t 2V t 2 Here ϕ Xn,X n+ (t, t 2 ϕ Xn,X n+ (t 2, t. Therefore the process is not time-reversible. 4.4 Multivariate geometric normal-laplace distribution The c.f. of the Multivariate Geometric Normal-Laplace distribution denoted by GNL p (µ, Σ, V, is given by ψ(t = it µ + + ln( + 2 t V t. The above c.f. can be written in the form ϕ(t = exp (, where ϕ(t is the c.f. of ψ(t NL p (µ, Σ, V. Hence the multivariate geometric normal-laplace distribution is geometrically infinitely divisible. 4.5 Multivariate generalized normal-laplace distribution A multivariate generalized normal-laplace distribution can be defined by introducing an additional parameter. Multivariate generalized normal-laplace distribution can be obtained as the convolution of multivariate normal and multivariate generalized symmetric Laplace random vectors. The c.f. of multivariate generalized normal-laplace distribution is given 84
12 by ϕ X (t = [ exp it µ ( ] ν +, t R p, Σ > 0, V > 0, ν > 0. 2 t V t (4.5. A p-variate generalized normal-laplace distribution with parameters µ, Σ, V and ν can be denoted by NL p (µ, Σ, V, ν. Let X NL p (µ, Σ, V, ν, then X can be expressed as X d = Z + Y (4.5.2 where Z and Y are independent random vectors with Z following a p-variate normal distribution with mean vector (νµ and dispersion matrix νσ (N p (νµ, νσ and Y following a p-variate generalized symmetric Laplace distribution with parameters V, ν (L p (V, ν. NL p (µ, Σ, V, ν is infinitely divisible and self-decomposable. 4.6 Multivariate geometric generalized normal-laplace distribution The c.f. of the Multivariate Geometric Generalized Normal-Laplace distribution denoted by GGNL p (µ, Σ, V, ν, is given by ψ(t = it νµ + 2 νt Σt + ν ln( + 2 t V t. Theorem Let X, X 2, are i.i.d. GGNL p (µ, Σ, V, ν random vectors and N be geometric with mean p, such that P [N = k] = p( pk, k =, 2,..., 0 < p <. This establishes that Y GGNL p (µ, Σ, V, ν/p Proof. The c.f. of Y is 85
13 ψ Y (t = = = [ψ X (t] k p( p k k= p/( it νµ + 2 νt Σt + ν ln( + 2 t V t (( p/( it νµ + 2 νt Σt + ν ln( + 2 t V t it ν µ + ν p 2 p t Σt + ν ln( + p 2 t V t. Hence Y is distributed as GGNL p (µ, Σ, V, ν/p. Theorem Multivariate geometric generalized normal-laplace distribution GGNL p (µ, Σ, V, ν is the limiting distribution of multivariate generalized normal-laplace distribution. Proof. We have [ exp it µ ( + 2 t V t = { + ] ν/n [ exp it µ ( ] ν/n + } 2 t V t is the c.f. of a probability distribution since Multivariate generalized normal-laplace distribution is infinitely divisible. Hence by Lemma 3.2 of Pillai (990. ψ n (t = { + n [ exp it µ ( ] ν/n + } 2 t V t is the c.f. of geometric sum of independently and identically distributed multivariate gener- 86
14 alized normal-laplace random vectors. Taking limit as n. ψ(t = lim ψ n (t n { { = = + lim n n [ exp it µ ( ] } ν/n + 2 t V t it νµ + 2 νt Σt + ν ln( + 2 t V t. Hence GGNL p (µ, Σ, V, ν is the limiting distribution of multivariate generalized normal- Laplace distribution. 4.7 Estimation of parameters Kollo and Srivastava (2004 discussed the estimation of Multivariate Laplace distribution. We can use the method of moments to estimate the mean, covariance matrix and the skewness and kurtosis measures of multivariate normal-laplace distribution. Let x, x 2,, x n are i.i.d. as NL p (µ, Σ, V. Then the estimates of mean vector and covariance matrix are obtained as and S = n x = n n i= x i n (x i x(x i x i= respectively. Skewness and kurtosis measures can be estimated using the random vector defined by y i = S /2 (x i x, i =,..., n, where S /2 is any square root of S such that S = S /2 ( S /2. Then the estimates of third and fourth moments are given by } M 3 = n n (y i y i y i i= 87
15 and M 4 = n + n n n (y i y i y i y i. i= These estimates are unbiased estimators, see Mardia (970. The estimates of skewness and kurtosis measures are given by β p = tr[ M 3 M 3 ] and β 2p = tr[ M 4 ]. 4.8 Multivariate slash normal-laplace distribution Now let us define the slash version of the multivariate Normal-Laplace distribution. Definition A random vector X R p has a p-variate slash Normal-Laplace (SNL p distribution with location parameter µ, positive definite scale matrix parameter Σ and tail parameter q > 0, denoted by X SNL p (µ, Σ, q, if X = Y U /q + µ, where Y is Normal-Laplace random vector with characteristic function given by ϕ X (t = exp it µ (, t R p, Σ > 0, V > 0 and U U(0,, which is + 2 t V t independent of Y. The stochastic representations of the random vector Y can be considered for simulation study. Let Y NL p (µ, Σ, V, then Y can be expressed as Y d = Z + L (4.8. where Z and L are independent random vectors with Z following a p-variate normal distribution with mean vector (µ and dispersion matrix Σ (N p (µ, Σ and L following a p-variate symmetric Laplace distribution with parameter V (L p (V. Another representation is Y d = Z + EZ (
16 where Z follows a p-variate normal distribution with mean vector µ and dispersion matrix Σ, E is a standard exponential variable and Z follows p-variate normal distribution with mean vector 0 and dispersion matrix Σ. The characteristic function of multivariate slash Normal-Laplace random vector X can be obtained as ϕ X (t = E (exp(it X = 0 ϕ Y (tu /q du where ϕ Y (t is the characteristic function of multivariate Normal-Laplace distribution given by ϕ Y (t = exp it µ ( + 2 t V t. (4.8.3 As q, the multivariate slash Normal-Laplace distribution tends to a multivariate Normal-Laplace distribution. 4.9 Applications In insurance and financial markets, there is significant need for the development of a standard framework for the risk measurement. Landsman and Valdez (2003 derives explicit formulas for computing tail conditional expectations for elliptical distributions and extends to multivariate case. Multivariate elliptical distributions are useful to model combinations of correlated risks. Since multivariate normal-laplace distribution comes in the class of elliptical distributions and it also convolutes both Gaussian and non-gaussian distributions, it also has a very important role in risk analysis. It seems to provide an attractive tool for actuarial and financial risk management because it allows a multivariate portfolio of risks to have the property of regular variation in the marginal tails. The importance of normal-laplace model lies in the fact that it is the first attempt to combine Gaussian and non-gaussian marginals to model time series data, see Jose et al (2008 and Lishamol and Jose (2009,200,20. Applications of normal-laplace distribution are wide spread in areas like financial modelling, Levy process, Brownian motion etc. 89
17 Reed (2006 showed that it is the distribution of the stopped state of a Brownian motion with normally distributed starting value if the stopping hazard rate is constant. In financial modelling, a normal-laplace model is a more realistic alternative for Gaussian models as logarithmic price returns do not follow exactly a normal distribution. But it is more realistic to consider multivariate data where several variables are discussed. The models developed in this chapter can be used for modelling multivariate time series data. 4.0 Conclusion In this chapter the univariate normal-laplace distribution is considered. We introduced a multivariate normal-laplace distribution and studied its properties. First order autoregressive processes with multivariate normal-laplace marginals is developed. We introduced a multivariate geometric normal-laplace distribution. Multivariate generalized normal- Laplace distribution is also introduced. We introduced the geometric generalized normal- Laplace distribution and studied its properties. We considered the estimation of parameters. Some applications are also discussed. References Blaesid, P. (98. The two-dimensional hyperbolic distribution and related distributions, with an application to Johannsen s bean data, Biometrika 68, Ernst, M.D. (998. A multivariate generalized Laplace distribution, Computational Statistics 3, Fang, Kai-Tai, Kotz, S., Kai Wang Ng. (987. Symmetric multivariate and related distributions, Chapman & Hall.London. Jose, K.K, Lishamol, T., Sreekumar, J. (2008. Autoregressive processes with normal Laplace marginals, Statistics and Probability Letters 78, Kelker, D. (970. Distribution theory of spherical distributions and a location scale parameter generalization, Sankhya, Series A, 32,
18 Kollo, T., Srivastava, M.S. (2004. Estimation and testing of parameters in multivariate Laplace distribution, Communications in Statistics - Theory and Methods 33(0, Landsman, Z.M., Valdez, E.A. (2003. Tail conditional expectations for elliptical distributions, North American Actuarial Journal 7(4, Lishamol, T., Jose, K.K. (2009. Generalized normal-laplace AR process, Statistics and Probability Letters 79, Lishamol, T., Jose, K.K. (200. A unified framework for Gaussian and non-gaussian AR( modelling, Journal of Probability and Statistical Science 8(, Lishamol, T., Jose, K.K. (20. Geometric Normal-Laplace Distributions and Autoregressive Processes, Journal of Applied Statistical Science 8(, Mardia, K.V. (970. Measures of multivariate skewness and kurtosis with applications, Biometrical Journal 3, Pillai, R.N. (990. Harmonic mixtures and geometric infinite divisibility, Journal of Indian Statistical Association 28, Reed, W.J., Jorgensen, M. (2004. The double Pareto lognormal distribution-a new parametric model for size distributions, Communications in Statistics-Theory and Methods 33(8, Reed, W.J. (2006. The normal-laplace distribution and its relatives. Advances in Distribution Theory, Order Statistics, and Inference Statistics for Industry and Technology Part II, Reed, W.J. (2007. Brownian Laplace motion and its use in financial modelling, Communications in Statistics-Theory and Methods 36(3,
Generalized Laplacian Distributions and Autoregressive Processes
CHAPTER Generalized Laplacian Distributions and Autoregressive Processes. Introduction Generalized Laplacian distribution corresponds to the distribution of differences of independently and identically
More informationMultivariate Distributions
IEOR E4602: Quantitative Risk Management Spring 2016 c 2016 by Martin Haugh Multivariate Distributions We will study multivariate distributions in these notes, focusing 1 in particular on multivariate
More informationA Note on Certain Stability and Limiting Properties of ν-infinitely divisible distributions
Int. J. Contemp. Math. Sci., Vol. 1, 2006, no. 4, 155-161 A Note on Certain Stability and Limiting Properties of ν-infinitely divisible distributions Tomasz J. Kozubowski 1 Department of Mathematics &
More informationKatz Family of Distributions and Processes
CHAPTER 7 Katz Family of Distributions and Processes 7. Introduction The Poisson distribution and the Negative binomial distribution are the most widely used discrete probability distributions for the
More informationWeighted Exponential Distribution and Process
Weighted Exponential Distribution and Process Jilesh V Some generalizations of exponential distribution and related time series models Thesis. Department of Statistics, University of Calicut, 200 Chapter
More informationTAMS39 Lecture 2 Multivariate normal distribution
TAMS39 Lecture 2 Multivariate normal distribution Martin Singull Department of Mathematics Mathematical Statistics Linköping University, Sweden Content Lecture Random vectors Multivariate normal distribution
More informationx. Figure 1: Examples of univariate Gaussian pdfs N (x; µ, σ 2 ).
.8.6 µ =, σ = 1 µ = 1, σ = 1 / µ =, σ =.. 3 1 1 3 x Figure 1: Examples of univariate Gaussian pdfs N (x; µ, σ ). The Gaussian distribution Probably the most-important distribution in all of statistics
More informationResearch Article The Laplace Likelihood Ratio Test for Heteroscedasticity
International Mathematics and Mathematical Sciences Volume 2011, Article ID 249564, 7 pages doi:10.1155/2011/249564 Research Article The Laplace Likelihood Ratio Test for Heteroscedasticity J. Martin van
More information3. Probability and Statistics
FE661 - Statistical Methods for Financial Engineering 3. Probability and Statistics Jitkomut Songsiri definitions, probability measures conditional expectations correlation and covariance some important
More informationVolatility. Gerald P. Dwyer. February Clemson University
Volatility Gerald P. Dwyer Clemson University February 2016 Outline 1 Volatility Characteristics of Time Series Heteroskedasticity Simpler Estimation Strategies Exponentially Weighted Moving Average Use
More informationSome results on Denault s capital allocation rule
Faculty of Economics and Applied Economics Some results on Denault s capital allocation rule Steven Vanduffel and Jan Dhaene DEPARTMENT OF ACCOUNTANCY, FINANCE AND INSURANCE (AFI) AFI 0601 Some Results
More informationGARCH Models Estimation and Inference
GARCH Models Estimation and Inference Eduardo Rossi University of Pavia December 013 Rossi GARCH Financial Econometrics - 013 1 / 1 Likelihood function The procedure most often used in estimating θ 0 in
More informationSupermodular ordering of Poisson arrays
Supermodular ordering of Poisson arrays Bünyamin Kızıldemir Nicolas Privault Division of Mathematical Sciences School of Physical and Mathematical Sciences Nanyang Technological University 637371 Singapore
More informationEfficient rare-event simulation for sums of dependent random varia
Efficient rare-event simulation for sums of dependent random variables Leonardo Rojas-Nandayapa joint work with José Blanchet February 13, 2012 MCQMC UNSW, Sydney, Australia Contents Introduction 1 Introduction
More informationMultivariate Statistics
Multivariate Statistics Chapter 2: Multivariate distributions and inference Pedro Galeano Departamento de Estadística Universidad Carlos III de Madrid pedro.galeano@uc3m.es Course 2016/2017 Master in Mathematical
More informationSome Generalizations of Weibull Distribution and Related Processes
Journal of Statistical Theory and Applications, Vol. 4, No. 4 (December 205), 425-434 Some Generalizations of Weibull Distribution and Related Processes K. Jayakumar Department of Statistics, University
More informationSTA 2201/442 Assignment 2
STA 2201/442 Assignment 2 1. This is about how to simulate from a continuous univariate distribution. Let the random variable X have a continuous distribution with density f X (x) and cumulative distribution
More informationMarshall-Olkin Bivariate Exponential Distribution: Generalisations and Applications
CHAPTER 6 Marshall-Olkin Bivariate Exponential Distribution: Generalisations and Applications 6.1 Introduction Exponential distributions have been introduced as a simple model for statistical analysis
More informationGaussian distributions and processes have long been accepted as useful tools for stochastic
Chapter 3 Alpha-Stable Random Variables and Processes Gaussian distributions and processes have long been accepted as useful tools for stochastic modeling. In this section, we introduce a statistical model
More information1 Simulating normal (Gaussian) rvs with applications to simulating Brownian motion and geometric Brownian motion in one and two dimensions
Copyright c 2007 by Karl Sigman 1 Simulating normal Gaussian rvs with applications to simulating Brownian motion and geometric Brownian motion in one and two dimensions Fundamental to many applications
More informationModule 3. Function of a Random Variable and its distribution
Module 3 Function of a Random Variable and its distribution 1. Function of a Random Variable Let Ω, F, be a probability space and let be random variable defined on Ω, F,. Further let h: R R be a given
More informationProduct Autoregressive Models with Log-Laplace and Double Pareto Log-normal Marginal Distributions
CHAPTER 6 Product Autoregressive Models with Log-Laplace and Double Pareto Log-normal Marginal Distributions 6.1 Introduction The relationships between Laplace and log-laplace distributions and double
More information1 Appendix A: Matrix Algebra
Appendix A: Matrix Algebra. Definitions Matrix A =[ ]=[A] Symmetric matrix: = for all and Diagonal matrix: 6=0if = but =0if 6= Scalar matrix: the diagonal matrix of = Identity matrix: the scalar matrix
More informationPractical conditions on Markov chains for weak convergence of tail empirical processes
Practical conditions on Markov chains for weak convergence of tail empirical processes Olivier Wintenberger University of Copenhagen and Paris VI Joint work with Rafa l Kulik and Philippe Soulier Toronto,
More informationDoes k-th Moment Exist?
Does k-th Moment Exist? Hitomi, K. 1 and Y. Nishiyama 2 1 Kyoto Institute of Technology, Japan 2 Institute of Economic Research, Kyoto University, Japan Email: hitomi@kit.ac.jp Keywords: Existence of moments,
More informationEEL 5544 Noise in Linear Systems Lecture 30. X (s) = E [ e sx] f X (x)e sx dx. Moments can be found from the Laplace transform as
L30-1 EEL 5544 Noise in Linear Systems Lecture 30 OTHER TRANSFORMS For a continuous, nonnegative RV X, the Laplace transform of X is X (s) = E [ e sx] = 0 f X (x)e sx dx. For a nonnegative RV, the Laplace
More informationMarshall-Olkin Univariate and Bivariate Logistic Processes
CHAPTER 5 Marshall-Olkin Univariate and Bivariate Logistic Processes 5. Introduction The logistic distribution is the most commonly used probability model for modelling data on population growth, bioassay,
More information(Multivariate) Gaussian (Normal) Probability Densities
(Multivariate) Gaussian (Normal) Probability Densities Carl Edward Rasmussen, José Miguel Hernández-Lobato & Richard Turner April 20th, 2018 Rasmussen, Hernàndez-Lobato & Turner Gaussian Densities April
More informationGumbel Distribution: Generalizations and Applications
CHAPTER 3 Gumbel Distribution: Generalizations and Applications 31 Introduction Extreme Value Theory is widely used by many researchers in applied sciences when faced with modeling extreme values of certain
More informationThe Use of Copulas to Model Conditional Expectation for Multivariate Data
The Use of Copulas to Model Conditional Expectation for Multivariate Data Kääri, Meelis; Selart, Anne; Kääri, Ene University of Tartu, Institute of Mathematical Statistics J. Liivi 2 50409 Tartu, Estonia
More informationMinimization of the root of a quadratic functional under a system of affine equality constraints with application in portfolio management
Minimization of the root of a quadratic functional under a system of affine equality constraints with application in portfolio management Zinoviy Landsman Department of Statistics, University of Haifa.
More informationLecture 11. Multivariate Normal theory
10. Lecture 11. Multivariate Normal theory Lecture 11. Multivariate Normal theory 1 (1 1) 11. Multivariate Normal theory 11.1. Properties of means and covariances of vectors Properties of means and covariances
More informationModelling and Estimation of Stochastic Dependence
Modelling and Estimation of Stochastic Dependence Uwe Schmock Based on joint work with Dr. Barbara Dengler Financial and Actuarial Mathematics and Christian Doppler Laboratory for Portfolio Risk Management
More informationChapter 5 continued. Chapter 5 sections
Chapter 5 sections Discrete univariate distributions: 5.2 Bernoulli and Binomial distributions Just skim 5.3 Hypergeometric distributions 5.4 Poisson distributions Just skim 5.5 Negative Binomial distributions
More informationSlash Distributions and Applications
CHAPTER 2 Slash Distributions and Alications 2.1 Introduction The concet of slash distributions was introduced by Kafadar (1988) as a heavy tailed alternative to the normal distribution. Further literature
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 information5.1 Consistency of least squares estimates. We begin with a few consistency results that stand on their own and do not depend on normality.
88 Chapter 5 Distribution Theory In this chapter, we summarize the distributions related to the normal distribution that occur in linear models. Before turning to this general problem that assumes normal
More informationA Marshall-Olkin Gamma Distribution and Process
CHAPTER 3 A Marshall-Olkin Gamma Distribution and Process 3.1 Introduction Gamma distribution is a widely used distribution in many fields such as lifetime data analysis, reliability, hydrology, medicine,
More informationThe Multivariate Normal Distribution 1
The Multivariate Normal Distribution 1 STA 302 Fall 2017 1 See last slide for copyright information. 1 / 40 Overview 1 Moment-generating Functions 2 Definition 3 Properties 4 χ 2 and t distributions 2
More informationGARCH Models Estimation and Inference. Eduardo Rossi University of Pavia
GARCH Models Estimation and Inference Eduardo Rossi University of Pavia Likelihood function The procedure most often used in estimating θ 0 in ARCH models involves the maximization of a likelihood function
More informationTesting Some Covariance Structures under a Growth Curve Model in High Dimension
Department of Mathematics Testing Some Covariance Structures under a Growth Curve Model in High Dimension Muni S. Srivastava and Martin Singull LiTH-MAT-R--2015/03--SE Department of Mathematics Linköping
More informationTail dependence coefficient of generalized hyperbolic distribution
Tail dependence coefficient of generalized hyperbolic distribution Mohalilou Aleiyouka Laboratoire de mathématiques appliquées du Havre Université du Havre Normandie Le Havre France mouhaliloune@gmail.com
More informationEcon 424 Time Series Concepts
Econ 424 Time Series Concepts Eric Zivot January 20 2015 Time Series Processes Stochastic (Random) Process { 1 2 +1 } = { } = sequence of random variables indexed by time Observed time series of length
More informationFinancial Econometrics and Volatility Models Copulas
Financial Econometrics and Volatility Models Copulas Eric Zivot Updated: May 10, 2010 Reading MFTS, chapter 19 FMUND, chapters 6 and 7 Introduction Capturing co-movement between financial asset returns
More informationConditional Tail Expectations for Multivariate Phase Type Distributions
Conditional Tail Expectations for Multivariate Phase Type Distributions Jun Cai Department of Statistics and Actuarial Science University of Waterloo Waterloo, ON N2L 3G1, Canada Telphone: 1-519-8884567,
More informationRandom Vectors and Multivariate Normal Distributions
Chapter 3 Random Vectors and Multivariate Normal Distributions 3.1 Random vectors Definition 3.1.1. Random vector. Random vectors are vectors of random 75 variables. For instance, X = X 1 X 2., where each
More informationGaussian Processes. Le Song. Machine Learning II: Advanced Topics CSE 8803ML, Spring 2012
Gaussian Processes Le Song Machine Learning II: Advanced Topics CSE 8803ML, Spring 01 Pictorial view of embedding distribution Transform the entire distribution to expected features Feature space Feature
More informationAsymptotic Statistics-III. Changliang Zou
Asymptotic Statistics-III Changliang Zou The multivariate central limit theorem Theorem (Multivariate CLT for iid case) Let X i be iid random p-vectors with mean µ and and covariance matrix Σ. Then n (
More informationMAS223 Statistical Inference and Modelling Exercises
MAS223 Statistical Inference and Modelling Exercises The exercises are grouped into sections, corresponding to chapters of the lecture notes Within each section exercises are divided into warm-up questions,
More informationSYSM 6303: Quantitative Introduction to Risk and Uncertainty in Business Lecture 4: Fitting Data to Distributions
SYSM 6303: Quantitative Introduction to Risk and Uncertainty in Business Lecture 4: Fitting Data to Distributions M. Vidyasagar Cecil & Ida Green Chair The University of Texas at Dallas Email: M.Vidyasagar@utdallas.edu
More informationGARCH Models. Eduardo Rossi University of Pavia. December Rossi GARCH Financial Econometrics / 50
GARCH Models Eduardo Rossi University of Pavia December 013 Rossi GARCH Financial Econometrics - 013 1 / 50 Outline 1 Stylized Facts ARCH model: definition 3 GARCH model 4 EGARCH 5 Asymmetric Models 6
More informationCHAPTER 1 DISTRIBUTION THEORY 1 CHAPTER 1: DISTRIBUTION THEORY
CHAPTER 1 DISTRIBUTION THEORY 1 CHAPTER 1: DISTRIBUTION THEORY CHAPTER 1 DISTRIBUTION THEORY 2 Basic Concepts CHAPTER 1 DISTRIBUTION THEORY 3 Random Variables (R.V.) discrete random variable: probability
More informationMFM Practitioner Module: Risk & Asset Allocation. John Dodson. February 3, 2010
MFM Practitioner Module: Risk & Asset Allocation Estimator February 3, 2010 Estimator Estimator In estimation we do not endow the sample with a characterization; rather, we endow the parameters with a
More informationChapter 5. Chapter 5 sections
1 / 43 sections Discrete univariate distributions: 5.2 Bernoulli and Binomial distributions Just skim 5.3 Hypergeometric distributions 5.4 Poisson distributions Just skim 5.5 Negative Binomial distributions
More informationMultiple Random Variables
Multiple Random Variables Joint Probability Density Let X and Y be two random variables. Their joint distribution function is F ( XY x, y) P X x Y y. F XY ( ) 1, < x
More informationDiscrete 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 information5. Random Vectors. probabilities. characteristic function. cross correlation, cross covariance. Gaussian random vectors. functions of random vectors
EE401 (Semester 1) 5. Random Vectors Jitkomut Songsiri probabilities characteristic function cross correlation, cross covariance Gaussian random vectors functions of random vectors 5-1 Random vectors we
More informationLifetime Dependence Modelling using a Generalized Multivariate Pareto Distribution
Lifetime Dependence Modelling using a Generalized Multivariate Pareto Distribution Daniel Alai Zinoviy Landsman Centre of Excellence in Population Ageing Research (CEPAR) School of Mathematics, Statistics
More informationThe Laplace driven moving average a non-gaussian stationary process
The Laplace driven moving average a non-gaussian stationary process 1, Krzysztof Podgórski 2, Igor Rychlik 1 1 Mathematical Sciences, Mathematical Statistics, Chalmers 2 Centre for Mathematical Sciences,
More informationUnit Roots in White Noise?!
Unit Roots in White Noise?! A.Onatski and H. Uhlig September 26, 2008 Abstract We show that the empirical distribution of the roots of the vector auto-regression of order n fitted to T observations of
More informationTail negative dependence and its applications for aggregate loss modeling
Tail negative dependence and its applications for aggregate loss modeling Lei Hua Division of Statistics Oct 20, 2014, ISU L. Hua (NIU) 1/35 1 Motivation 2 Tail order Elliptical copula Extreme value copula
More informationExtreme Value Analysis and Spatial Extremes
Extreme Value Analysis and Department of Statistics Purdue University 11/07/2013 Outline Motivation 1 Motivation 2 Extreme Value Theorem and 3 Bayesian Hierarchical Models Copula Models Max-stable Models
More informationHIGHER ORDER CUMULANTS OF RANDOM VECTORS, DIFFERENTIAL OPERATORS, AND APPLICATIONS TO STATISTICAL INFERENCE AND TIME SERIES
HIGHER ORDER CUMULANTS OF RANDOM VECTORS DIFFERENTIAL OPERATORS AND APPLICATIONS TO STATISTICAL INFERENCE AND TIME SERIES S RAO JAMMALAMADAKA T SUBBA RAO AND GYÖRGY TERDIK Abstract This paper provides
More informationNonlinear Time Series Modeling
Nonlinear Time Series Modeling Part II: Time Series Models in Finance Richard A. Davis Colorado State University (http://www.stat.colostate.edu/~rdavis/lectures) MaPhySto Workshop Copenhagen September
More informationThe Multivariate Normal Distribution 1
The Multivariate Normal Distribution 1 STA 302 Fall 2014 1 See last slide for copyright information. 1 / 37 Overview 1 Moment-generating Functions 2 Definition 3 Properties 4 χ 2 and t distributions 2
More informationCONTAGION VERSUS FLIGHT TO QUALITY IN FINANCIAL MARKETS
EVA IV, CONTAGION VERSUS FLIGHT TO QUALITY IN FINANCIAL MARKETS Jose Olmo Department of Economics City University, London (joint work with Jesús Gonzalo, Universidad Carlos III de Madrid) 4th Conference
More information1: PROBABILITY REVIEW
1: PROBABILITY REVIEW Marek Rutkowski School of Mathematics and Statistics University of Sydney Semester 2, 2016 M. Rutkowski (USydney) Slides 1: Probability Review 1 / 56 Outline We will review the following
More informationSummary of basic probability theory Math 218, Mathematical Statistics D Joyce, Spring 2016
8. For any two events E and F, P (E) = P (E F ) + P (E F c ). Summary of basic probability theory Math 218, Mathematical Statistics D Joyce, Spring 2016 Sample space. A sample space consists of a underlying
More informationStat 5101 Notes: Brand Name Distributions
Stat 5101 Notes: Brand Name Distributions Charles J. Geyer September 5, 2012 Contents 1 Discrete Uniform Distribution 2 2 General Discrete Uniform Distribution 2 3 Uniform Distribution 3 4 General Uniform
More information3 Continuous Random Variables
Jinguo Lian Math437 Notes January 15, 016 3 Continuous Random Variables Remember that discrete random variables can take only a countable number of possible values. On the other hand, a continuous random
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 informationdiscrete random variable: probability mass function continuous random variable: probability density function
CHAPTER 1 DISTRIBUTION THEORY 1 Basic Concepts Random Variables discrete random variable: probability mass function continuous random variable: probability density function CHAPTER 1 DISTRIBUTION THEORY
More informationQuantile prediction of a random eld extending the gaussian setting
Quantile prediction of a random eld extending the gaussian setting 1 Joint work with : Véronique Maume-Deschamps 1 and Didier Rullière 2 1 Institut Camille Jordan Université Lyon 1 2 Laboratoire des Sciences
More informationLecture 3. Probability - Part 2. Luigi Freda. ALCOR Lab DIAG University of Rome La Sapienza. October 19, 2016
Lecture 3 Probability - Part 2 Luigi Freda ALCOR Lab DIAG University of Rome La Sapienza October 19, 2016 Luigi Freda ( La Sapienza University) Lecture 3 October 19, 2016 1 / 46 Outline 1 Common Continuous
More informationGeneralized normal mean variance mixtures and subordinated Brownian motions in Finance. Elisa Luciano and Patrizia Semeraro
Generalized normal mean variance mixtures and subordinated Brownian motions in Finance. Elisa Luciano and Patrizia Semeraro Università degli Studi di Torino Dipartimento di Statistica e Matematica Applicata
More information4. Distributions of Functions of Random Variables
4. Distributions of Functions of Random Variables Setup: Consider as given the joint distribution of X 1,..., X n (i.e. consider as given f X1,...,X n and F X1,...,X n ) Consider k functions g 1 : R n
More informationExam 2. Jeremy Morris. March 23, 2006
Exam Jeremy Morris March 3, 006 4. Consider a bivariate normal population with µ 0, µ, σ, σ and ρ.5. a Write out the bivariate normal density. The multivariate normal density is defined by the following
More informationOn q-gamma Distributions, Marshall-Olkin q-gamma Distributions and Minification Processes
CHAPTER 4 On q-gamma Distributions, Marshall-Olkin q-gamma Distributions and Minification Processes 4.1 Introduction Several skewed distributions such as logistic, Weibull, gamma and beta distributions
More informationDiscrete time processes
Discrete time processes Predictions are difficult. Especially about the future Mark Twain. Florian Herzog 2013 Modeling observed data When we model observed (realized) data, we encounter usually the following
More information. Find E(V ) and var(v ).
Math 6382/6383: Probability Models and Mathematical Statistics Sample Preliminary Exam Questions 1. A person tosses a fair coin until she obtains 2 heads in a row. She then tosses a fair die the same number
More informationJump-type Levy Processes
Jump-type Levy Processes Ernst Eberlein Handbook of Financial Time Series Outline Table of contents Probabilistic Structure of Levy Processes Levy process Levy-Ito decomposition Jump part Probabilistic
More informationTest Code: STA/STB (Short Answer Type) 2013 Junior Research Fellowship for Research Course in Statistics
Test Code: STA/STB (Short Answer Type) 2013 Junior Research Fellowship for Research Course in Statistics The candidates for the research course in Statistics will have to take two shortanswer type tests
More informationStochastic Processes
Stochastic Processes Stochastic Process Non Formal Definition: Non formal: A stochastic process (random process) is the opposite of a deterministic process such as one defined by a differential equation.
More informationLecture 2: Univariate Time Series
Lecture 2: Univariate Time Series Analysis: Conditional and Unconditional Densities, Stationarity, ARMA Processes Prof. Massimo Guidolin 20192 Financial Econometrics Spring/Winter 2017 Overview Motivation:
More informationCopula Regression RAHUL A. PARSA DRAKE UNIVERSITY & STUART A. KLUGMAN SOCIETY OF ACTUARIES CASUALTY ACTUARIAL SOCIETY MAY 18,2011
Copula Regression RAHUL A. PARSA DRAKE UNIVERSITY & STUART A. KLUGMAN SOCIETY OF ACTUARIES CASUALTY ACTUARIAL SOCIETY MAY 18,2011 Outline Ordinary Least Squares (OLS) Regression Generalized Linear Models
More informationNotes on Random Vectors and Multivariate Normal
MATH 590 Spring 06 Notes on Random Vectors and Multivariate Normal Properties of Random Vectors If X,, X n are random variables, then X = X,, X n ) is a random vector, with the cumulative distribution
More informationTwo hours. To be supplied by the Examinations Office: Mathematical Formula Tables and Statistical Tables THE UNIVERSITY OF MANCHESTER.
Two hours MATH38181 To be supplied by the Examinations Office: Mathematical Formula Tables and Statistical Tables THE UNIVERSITY OF MANCHESTER EXTREME VALUES AND FINANCIAL RISK Examiner: Answer any FOUR
More informationWeek 1 Quantitative Analysis of Financial Markets Distributions A
Week 1 Quantitative Analysis of Financial Markets Distributions A Christopher Ting http://www.mysmu.edu/faculty/christophert/ Christopher Ting : christopherting@smu.edu.sg : 6828 0364 : LKCSB 5036 October
More informationRESEARCH REPORT. Estimation of sample spacing in stochastic processes. Anders Rønn-Nielsen, Jon Sporring and Eva B.
CENTRE FOR STOCHASTIC GEOMETRY AND ADVANCED BIOIMAGING www.csgb.dk RESEARCH REPORT 6 Anders Rønn-Nielsen, Jon Sporring and Eva B. Vedel Jensen Estimation of sample spacing in stochastic processes No. 7,
More informationMonotonicity and Aging Properties of Random Sums
Monotonicity and Aging Properties of Random Sums Jun Cai and Gordon E. Willmot Department of Statistics and Actuarial Science University of Waterloo Waterloo, Ontario Canada N2L 3G1 E-mail: jcai@uwaterloo.ca,
More informationCalculation of Bayes Premium for Conditional Elliptical Risks
1 Calculation of Bayes Premium for Conditional Elliptical Risks Alfred Kume 1 and Enkelejd Hashorva University of Kent & University of Lausanne February 1, 13 Abstract: In this paper we discuss the calculation
More informationRandom Vectors 1. STA442/2101 Fall See last slide for copyright information. 1 / 30
Random Vectors 1 STA442/2101 Fall 2017 1 See last slide for copyright information. 1 / 30 Background Reading: Renscher and Schaalje s Linear models in statistics Chapter 3 on Random Vectors and Matrices
More informationIf we want to analyze experimental or simulated data we might encounter the following tasks:
Chapter 1 Introduction If we want to analyze experimental or simulated data we might encounter the following tasks: Characterization of the source of the signal and diagnosis Studying dependencies Prediction
More informationAsymptotic Tail Probabilities of Sums of Dependent Subexponential Random Variables
Asymptotic Tail Probabilities of Sums of Dependent Subexponential Random Variables Jaap Geluk 1 and Qihe Tang 2 1 Department of Mathematics The Petroleum Institute P.O. Box 2533, Abu Dhabi, United Arab
More informationLeast Squares Estimators for Stochastic Differential Equations Driven by Small Lévy Noises
Least Squares Estimators for Stochastic Differential Equations Driven by Small Lévy Noises Hongwei Long* Department of Mathematical Sciences, Florida Atlantic University, Boca Raton Florida 33431-991,
More informationFinancial Econometrics
Financial Econometrics Estimation and Inference Gerald P. Dwyer Trinity College, Dublin January 2013 Who am I? Visiting Professor and BB&T Scholar at Clemson University Federal Reserve Bank of Atlanta
More informationMultivariate Distribution Models
Multivariate Distribution Models Model Description While the probability distribution for an individual random variable is called marginal, the probability distribution for multiple random variables is
More informationFrontier estimation based on extreme risk measures
Frontier estimation based on extreme risk measures by Jonathan EL METHNI in collaboration with Ste phane GIRARD & Laurent GARDES CMStatistics 2016 University of Seville December 2016 1 Risk measures 2
More informationVaR bounds in models with partial dependence information on subgroups
VaR bounds in models with partial dependence information on subgroups L. Rüschendorf J. Witting February 23, 2017 Abstract We derive improved estimates for the model risk of risk portfolios when additional
More informationLARGE DEVIATION PROBABILITIES FOR SUMS OF HEAVY-TAILED DEPENDENT RANDOM VECTORS*
LARGE EVIATION PROBABILITIES FOR SUMS OF HEAVY-TAILE EPENENT RANOM VECTORS* Adam Jakubowski Alexander V. Nagaev Alexander Zaigraev Nicholas Copernicus University Faculty of Mathematics and Computer Science
More information