= P k,j = l,,nwhere P ij
|
|
- Patrick Hensley
- 5 years ago
- Views:
Transcription
1 UNIVERSITY OF NORTH CAROLINA Department of Statistics Chapel Hill, N. C. ON A DISTRIBUTION OF SUM OF IDENrICALLY DISTRIBUTED CORRElATED GAMMA VARIABIES by S. Kotz and John W. Adams June 1963 Grant No. AF-AFISR \ ( A distribution of a sum of identically distributed Gamma-variables correlated according to an 'exponential' autocorrelation law P kj = P k,j = l,,nwhere P ij is the correlation coefficient between the i-th and j-th random varia~les /k-j '< ) and 0 < P < 1 is a given number is derived. This is performed by determining the characteristic roots of the appropriate variance - covariance matrix, using a special method, and by applying Robbins' and Pitman's result on a mixtures of distributions to the case of Ga.mm&ovariables. An "approximate" distribution of the sum of these variables under the assumption that the sum itself is a Gamma variable is given. A comparison between an "exact" and an "approximate"distribution for certa.in values of the correlation coefficient, the number of variables in the sum and the values of parameters of the initial distributions are presented. This research was supported by the Mathematics Division of the Air Force Office of Scientific Research. Institute of Statistics Mimeo. Series No. 367
2 ON A DISTRIBUTION OF SUM OF IDENTICALLY DISTRIBUTED CORRELATED GJIMMA VARIABIES I by Samuel Kotz and John W. Adams University of North Carolina ==== == = = = = O. Abstract: = = = = ===== = ======== = = = = == = == = = = A distribution of a sum of identica~.ly distributed Gamma-variables correlated according to an 'exponential' autocorre~tion law Pkj = plk- j (k l j = l~ n) where P ij is the correlation coefficien,t betileen the i-th and j-th random variables and 0 < P < 1 is a given n1jil1per 1& oerived. This is performed by determining the characteristic roots of the appropriate variance - covariance matrix l using a special method, and by applying Robbins f and Pitman's result on mixtures of distributions to the case of Gamma-variables. An "approximate" distribution of the sum of these variables under the assllillption that the sum itself is a Gamma variable is given. A comparison between an "exact" and an "approximate" distributions for certain values of the correlation coefficient, the number of variables in the sum and the values of parameters of the initial distributions are presented. 1. Introduction and general remarks. Distribution of sum of correlated Gamma variables has 8 significant interest and many applications in engineering1 insurance and other areas. Besides the applications mentioned in I)} and f~7 'We should like to note the usefulness of this distribution in problems connected with representation of precipitation amounts. Weekly or monthly, etc., precipitation amounts are regarded as sums of ~his research was supported by the Mathematics Division of ~he Air Force Office of Scientific Research.
3 daily amounts, the latter being well fitted by Gamma distribution L'd.7. A partieular solution for the case of a constant correlation between each ~air of variables in the sum is ~resent.ed in I)}. In this note we shall extend this result for the case of " an exponential autocorrelation scheme" between the variables, where each one of the variables has the marginal density function given by f(x) X: 1 --g = e X r - l, r(r)q"j.. x>o = 0 x<o In meteorological ap~lications it is sojaetimes assumed that the sum random variable is itself a Gamma. variable and a.n "approximate" Gamma distribution whose first two moments are identical with the "exact" distribution of the sum is used L'd.7, (see also, for exam~le, ~~7 for an earlier a~plication in another field). In the case of identically distributed r.v. and of stationary exponential correlation model (namely when the correlation coefficient between the i-th and the j-th r.v. in the sum is given by p i-j for all i and j, p being a given positive constant) it is easy to verify L'd.7 that the appropriate parameters rand Q of the "appr~fjilaten Gamma. variable, representing the sum r.v., n n are given by: and n2. r n = r 2. 1 n Q = n n+...e(n _ -=.e ). l-p l-p n n+ ~( n-- l-p) l..p l-p n when n is the number variables in the sum and rand Q are the parameters Q of initial r.v. distributed according to (1.1). Since the first two moments
4 3 com.:p1etely determine Gamma-distribution the "a:p:proximate" distribution is" unique 4 'rhe :pur:pose of this note is to cozn:pare this "a:p:proximate" distribution (which is quite easy com.:putab1e and a:p:p1ied) with an "exact distribution of the sum (which is much more difficult to be actually determined). We also note that, similar to the (lase treated in ["J}, the obtained ifnatura1" exact distribution may not be unique, since the correlation scheme, in general, does not determine uniquely a multivariate Gamma distribution as it will be noted in some more detail in the following section. 2. JA ;Joint Characteristic P\.mcti. on of I!lx;Ponentia11y Correlated Gamma Variables Consider the characteristic function (t, t,, t ), 1 2 n, t) = n. II - r Q T vi -r where Q and r are :positive constants, I is the nxn identity matrix, T is the nxn nxn diagonal matrix with the elements t jj = t j, and V is an arbitrary :positive definite matrix. This chara.c tel'~stic function leads to a joint :probability density function whose margin?ls are given by (1.1) and whose matrix of second moments is some :positive definite matrix V*, say. In general, a joint d density function with marginal given by (1.1) and a given covariance matrix V* is not uniquely determined; (for exam:p1e.. as it is known, a multivariate Gamma distribution with constant correlation between the components is not unique L27). An explicit expression for a joint distribution with a general matrix is rather complicated.. but can be obtained in terms of a series of Laguerre polynomials by applying the derivations given in LE7. We will not present these expressions here since the main concern of this paper is the validity Elind applicability of approximations to the distribution of the sum of the component random variables. Following Gurland LJ:7 we shall regard (2.1) as the characteristic function of a natural multivariate analogue of the Gamma distribution.
5 In particular, if the elements of V are given by V ij = p~i-jj. 4 i, j =1, n (0 < P < 1 a given constant) it is easy to verify by differentiation of (2.1) that the corresponding elements of V* will be given by 2 2ri-jf V1'j = rq p...' i,j =1,... n. The characteristic function of the distribution of the sum of the random variables whose joint distribution has the characteristic function given in (2.1) is: (t) = JII - i Q t vi -r and (2.3) may be expressed in the form (t) = n--'(l -i Q A. j tr r j:;l where the A. j are the characteristic roots of the matrix V. In section 4 we shall describe a procedure for an explicit determination of these characte~istic roots. 3. The Distribution of the Random Variables ~ose Characteristic Function is (t) Robbins p~c:red in 17 that if {\-.ll is a ~equence of constants such that '1t? 0, k = ~)... and k:o ~ = 1 and if. { F. k ( )J is a sequence of distrib~tien functions with corresponding sequence of characteristic functions { kl'lj. then (i) FC,-) = I: ;c k Fk.~ d is a distribution function, (ii) the characteristic function (.) of F(.) is given by (.) = I: ~k k(.)' and (iii) for any Borel measurable function g( ) the following equality holds: Jg(x) df(x) = I: c k j g(x) dfk(x), where 1.* is the approp~ia.te space. These ~ k ~. results is valid for a general class of spaces and in particular for any finite dimensional Euclidean space. The distribution function of the Gamma variable with positive parameters rand Q denoted here by F () is given by:. r
6 5 ' 1 F ( X») = r l or rem) o r-1 u du I x> 0 x<o The characteristie function of this distribution is *(t) = (1 - i Q t)-r Let Y be the random variable whose characteristic function is given by (2.4). The Y has the same distribution as the random variable X, X = n ~ A X. j=l j J (,.,) where X jl j =1, 2, n, are independent identically distributed random variables, each following a Gamma distribution with parameters ~ and Q as given in (,.1) Let A* = min A j Then j Prey < y) = Pr (~* < t*) (,.4) and the characteristic function of the random variable X/A* is n Aj-r ~ (t) = 1T (1 - i 0 ~ t) j=l Using the method developed by Pitman and Robbins in f l, we express the right hand side of (,.5) in the form n A-r 1T (1 - i Q A~ t) j=l -nr = (1 - i 0 t)
7 6 Since" (3.7) for all real Q and t the binomial theorem may be applied to expand the r.h.s. of (3.6). We thus obtain fr (1 _ i Q 7-. j t)-r =.1=1 7-.* where the coefficients c are determined by the identity k n I *)-r!Xl 17/ ( ~ ) 1 - (1 - r-) z_7, = I: Ok zk.1=1 l j j k=o It f01~ows from (3.6) that all. c k > 0" and putting (1. i Q t) =1 in (3.6) and (3.8) shows that I: c k =1 Thus the results of Robbins 17 apply for c k determined 'by (3.9) and hence the r.h.s. of (3.8) is the characteristic function of the distribution F(.. ) given by!xl F(x) = I: c k F -.Lk (x) (3.10) k=o nol" where Fnr+k( ) (k =0"1,, ) are defined by (3.1). Hence Pr(Y < y) co = I: c k Fnr+k (Y/7-.*) = k=o ẏ u Xi --. (U)nr+k-l Q ~. ' ~ e \,I,U o (3.11) and the corresponding probability density function is given by:
8 x 7 co \ nr+k-1 - G* fiy ) - I 1 1: (~~Q) e ~ -I, x>o k=o f(nr+k1 ta~ Q 0 x<o (The term by term differentiation is justified by the uniform convergence of the series). The upper bound on the error of truncation is given by the following formula: for any integers 0 ~ P1 ~ P2 and for any xl.. P2 o ~Pr(Y < y) - 1: P1 i P1- l P P2 + sup J L Fnr+k(Y/~*) L (1-1:!Ok - 1: 2 ~) ~ 1-1: ok. :J 0 P1 P1 We should like to point out that the expansions (3.8) and (3.9) are closely related to the general results on expansion of the hypergeametric functions of matrix variables into series of Go ~ca.llea."zonal polynomials" flo, 10~7. 4. The Characteristic Roots of the Matrix V. i-j The inverse matrix of V, whose elements are v ij = P,i, j = 1,,n, is V- l given by V- = «l+p )I - P A - P B)(l-p ) (4.1) '\-There I is the nxn identity matrix, A is the nm matrix with elements all = ann = 1 and all other elements 0,and B is the nxn matrix with b ij = 1 for li..j J = 1 and all other elements O. Hence the characteris tic roots of V- are (l-p) Y j I ~ =l,2,,n, where the Y j are those value of Y which satisfy the equation;
9 Let bon{n,p) be the f.h.s. of (4.2). Expanding An(A,p) by its minors we obtain 8 bon(y,p) = (1 - y)2 D n _ 2 (y,p) - 2p2(1-Y)Dn_,(y,p) 4 + P D n _ 4 (y,p) (4.,) where Dn(y,p) is the determinant of the nxn matrix I U, U = (1 + (:,2 - y) - P B (4.4) Expanding D n {YIP) by its minors we get the following difference equation and, is it is well know., (see e.g. 27), this difference equation has the solution where = pn sin(n+l} Q sin Q P - Y = 2p Cos Q (4.6) equation Substituting (4.6) in (4.;) and equating bon(y,p) to zero ~ obtain the ( 1 _ y)2 pn-2 sin(n-l)q _ 2 (1 _ ) n-l sin(n-2)q sin Q.. Y P. sin Q + + n sin(n-,)q P sin Q = o Using the complex representation of the sine function we get after multiplying (4.7) by sin Q (excluding the value of Q for which sin Q = 0), i(n-,}q.~ e... ~. l-ye iq) -p =,:tei(n-,) 2 Q {( l-ye ) -iq -P, ) «) and after a few a.lgebraic manipulations the equation (4.8) reduces to (4.8)
10 9 and Sin (n+l) Q Sin (n-l) Q = p 2 2 Cos (n+l) Cos (n-l) Q Q. 2 = P 2 (~.9) (4.10) If Q j j = 1~2, n are the values of Q which satisfy one or the other of the equations (4.9) and (4.10), then the characteristic roots of the matrix V 1 are given by Yj = (1-2p Cos Q j + P ) (1- p), and the characteristic roots of V are ~j =Yj = (1-2p Cos Q j + p) (1 - p ), j = l,2,,n. Apparently it is not possible to express the solutions of (4.9) and (4.10) as explicit functions of I&' and p. However, for 0 ~ p ~ 1" it can be verified that the solutions of (4.9) are located one in each of the intervals ( 2k-l 2k3t ) n n-l.. ~ 3t, n+l ' k = 1,2,... '2 for T depending on whether n ~s odd; and that the solutions of (4.10) are one in each of the intervals -n 3t, li.'+"l tt" k = 1,2,, '2 or T depending on whether n ( 2k-2 2k-l) n n+l of odd. even or is even Using this observation we can obtain the solutio~ of these equations. Consider the iterative equations: Qjk 2 ( -l( n-l» = ii+i ktt - Sin p Sin 2" Qj-l"k ' for k odd r ~. ln~l (ktt + Sin- l (p Sin n;l Qj-l,k»' for k even (4.11) where Q (2k-13t 2K1t) ok -n' ii+i ' n (n-1) and k =1,2,... '2 T for 7C -1(. n-l ) tt - '2 ~ Sinp Si~ Qj-l"k ~'2' j =0,1,, n even (odd); and
11 where Q e (2k-2 n. 2k-l n). n < C -l( Cos n-l Q ) < ~ ok n ' n+ l ' - "2 - os p :2 j-l,k - t=. ' n (n+l) ( ) j =0,1, and k = 1,2, "2 2 for n even odd. The series of solutions obtained using this iterative procedure actually converges to the solutions of equations (4.9) and (4.10). 10 This can easily be seen by means of the following theorem of differential calculus. (See e.g. 117). o If Q e 1\, where '--0 1\ is an interval containing a root of the equation o Q = g(q), and if g(q) is a diffe;rential function of Q throughoutr and, (...:..' in addition, Ig'(Q) I < M < 1 for all Q e I?, then the sequence of iterative solutions of the equation Q j = g(qj-l)' starting with the point Qo' converges to a root of the equation. In the case of the equations (4.11) and (4.12) the conditions of this theorem are clearly satisfied. In fact, Q 1. belongs, in each case, to an interval 01.\,. containing exactly one solution of (4.11) or (4.12) respectively. absolute values of the derivatives of (4.11) and (4.12) are: ( n-l) n-l! P Cos T Q ii+i1'1_p2 6in2 ( > Q Moreover, the and P Sin (n-l) Q. 2 (4.14) which for all p, ipl < 1 and for all real Q, are, in each case, less than!!:1: < 1 ntl 5. Comparison between the "a;p;proximate lt and "exact lt distributions of sums of identically distributed exponentially correlated Gamma-variables. In order to obtain the distribution of sum of identically distributed Gamma variables correlated according to the stationary exponential law (r (Xi' X ) j = p,i- j /, i,j =1, n) we replace in (2.1) the matrix
12 li-jl/ 11 V = V ij = P f.t-~ by the matrix 1""1 =[w ij } =ip 2 so that the matrix J "{-1* vrhich corresponds to the matrix V* 'Will have its elements.,1~ t'1'j1 =lr Q2 p li-j1' ~here "- ~l is the common variance each of the gammavariables Xi i = l, n. Using the UNIVAC al05 of the Computation Center of the University of North Carolina we computed the parameters { ~J and "-*.of the distribution given by (:5.11) for several different values of n" p, Q" t and compared several percentiles of this distribution 'With the corresponding percentiles of the "approximate" Gamma distribution with the parameters given by (1.2) and (1.'). The results of the computations are presented in Table 1 on page 12. We recall that the first two moments of these two distributions are identical; (this fact has been used for checking the correctness of the numerical calculations). Q'\.ling to computational difficulties and complications especially on computing the sequences t~j we cannot claim higher accuracy than the second significant figure. Comparing these two types of distributions we observe that the functional dependence of parameters on n follows the same pattern in both cases. The parameter r is asymptotically linear with. n while Q is independent of n. The discrepancies between the corresponding distributions are relatively small and may be considered insignificant for most of the practical purposes in the cases of small values of pep So.,) even for values of n as small as 5. The computations also indicate that for higher values of p (p = 0.5 and greater) the relatively large deviations at the lower tail decline sharply with the increase of the values of n. The number of. c k ': necessary to obtain the cumulative sum equal to 1
13 The number of iterative steps necessary to obtain the values of the characteristic roots ~j with accuracy up to the 6-th significant digit also depends strongly on nand p and varies for the cases presented below between 5 and 18. TABLE I Comparison between the approximate a.nd "exact" distributions of sum of identically distributed exponentially correlated Gamma variables 12 IPercentiles of _ct distribution corresponding to the follmdng percentiles of the approximate distribution n p r Q ; ; , ; ;.01 2; ;4 ; ;; * *This case corresponds to a distribution of precipitation amounts at a certain station in..
14 ACKNOWLEDGE:ME:1'T. The authors wish to express their thanks to Professor W. HoefftUng and especially to Professor N. L. Johnson for their valuable ccmm.ents. We are also much indebted to Mr. P. J. Bro'WIl who carried out the programming of all the computations on the UNIVAC computer. Studien M1. 3 (1961)" p , 1'2.7 S. Kotz and J. Neumann" 'ton distribution of precipitation amounts for the periods of increasing lengths,," to be published in J. Geoph. Research (1963).!:.7 Welch" B. L." "The significance of difference between two means when the e population variances are unequal," Biometrika" Vol. 29 (1938)" p J. Gurland" Correction to "Distribution of the Maximum of the Arithmetic Mean of Correlated Random Variables,," Ann. Math. Stat., Vol. 30 (1959), p J A. S. Krishnamoorthy and M. Parthasarathy" "A Multivariate Gamma-Type Distribution," Ann. Math. Stat." Vol. 22 (1951)" p L17 H. Robbins" ''Mixture sf Distributions,," Ann. Math. Stat." Vol. 19 (1948)".L.7 p H. Robbins and E. J. G. Pitman" "Application of the method of mi~ures to quadratic forms in normal variates, n Ann. Math, Stat." Vol. 20 (1949)" 1'.552.L27 F. B. Hildebrand" Introduction to Numerical Analysis" McGraw Hill, J.g7 A. '.diw1l~'t J,~.he, d~s.ri.but:j.ott af~ 1a1teX!.t roots of the covariance matrix," Anp. :t-1ath. Stat., Vol.. 31 (1960)1. p L10!7 A. T. James" "A comparison of Univariate and Multivariate Distributions,," special invited ~aper presented at the 94th Eastern Regional Meeting of the n;.m (May 6, 1963)., ["1]} L. Brand, Advanced Calculus" Wiley and Son, 1955.
Inferences on a Normal Covariance Matrix and Generalized Variance with Monotone Missing Data
Journal of Multivariate Analysis 78, 6282 (2001) doi:10.1006jmva.2000.1939, available online at http:www.idealibrary.com on Inferences on a Normal Covariance Matrix and Generalized Variance with Monotone
More informationA NOTE ON CONFIDENCE BOUNDS CONNECTED WITH ANOVA AND MANOVA FOR BALANCED AND PARTIALLY BALANCED INCOMPLETE BLOCK DESIGNS. v. P.
... --... I. A NOTE ON CONFIDENCE BOUNDS CONNECTED WITH ANOVA AND MANOVA FOR BALANCED AND PARTIALLY BALANCED INCOMPLETE BLOCK DESIGNS by :}I v. P. Bhapkar University of North Carolina. '".". This research
More informationDepartment of Mathematics Temple University~ Philadelphia
lsupported by the U.S. Army Research Office, Durham, under Contract No. DAHC0471C0042. 2Supported by the Air Force Office of Scientific Research under Grant AFOSR722386. EXTENDED OCCUPANCY PROBABILITY
More information~,. :'lr. H ~ j. l' ", ...,~l. 0 '" ~ bl '!; 1'1. :<! f'~.., I,," r: t,... r':l G. t r,. 1'1 [<, ."" f'" 1n. t.1 ~- n I'>' 1:1 , I. <1 ~'..
,, 'l t (.) :;,/.I I n ri' ' r l ' rt ( n :' (I : d! n t, :?rj I),.. fl.),. f!..,,., til, ID f-i... j I. 't' r' t II!:t () (l r El,, (fl lj J4 ([) f., () :. -,,.,.I :i l:'!, :I J.A.. t,.. p, - ' I I I
More informationMultivariate Normal-Laplace Distribution and Processes
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
More informationTesting the homogeneity of variances in a two-way classification
Biomelrika (1982), 69, 2, pp. 411-6 411 Printed in Ortal Britain Testing the homogeneity of variances in a two-way classification BY G. K. SHUKLA Department of Mathematics, Indian Institute of Technology,
More informationJUST THE MATHS UNIT NUMBER DIFFERENTIATION APPLICATIONS 5 (Maclaurin s and Taylor s series) A.J.Hobson
JUST THE MATHS UNIT NUMBER.5 DIFFERENTIATION APPLICATIONS 5 (Maclaurin s and Taylor s series) by A.J.Hobson.5. Maclaurin s series.5. Standard series.5.3 Taylor s series.5.4 Exercises.5.5 Answers to exercises
More informationwhich arises when we compute the orthogonal projection of a vector y in a subspace with an orthogonal basis. Hence assume that P y = A ij = x j, x i
MODULE 6 Topics: Gram-Schmidt orthogonalization process We begin by observing that if the vectors {x j } N are mutually orthogonal in an inner product space V then they are necessarily linearly independent.
More informationON THE DIVERGENCE OF FOURIER SERIES
ON THE DIVERGENCE OF FOURIER SERIES RICHARD P. GOSSELIN1 1. By a well known theorem of Kolmogoroff there is a function whose Fourier series diverges almost everywhere. Actually, Kolmogoroff's proof was
More informationExplicit evaluation of the transmission factor T 1. Part I: For small dead-time ratios. by Jorg W. MUller
Rapport BIPM-87/5 Explicit evaluation of the transmission factor T (8,E) Part I: For small dead-time ratios by Jorg W. MUller Bureau International des Poids et Mesures, F-930 Sevres Abstract By a detailed
More informationSOME RESULTS ON THE MULTIPLE GROUP DISCRIMINANT PROBLEM. Peter A. Lachenbruch
.; SOME RESULTS ON THE MULTIPLE GROUP DISCRIMINANT PROBLEM By Peter A. Lachenbruch Department of Biostatistics University of North Carolina, Chapel Hill, N.C. Institute of Statistics Mimeo Series No. 829
More informationStein s Method and Characteristic Functions
Stein s Method and Characteristic Functions Alexander Tikhomirov Komi Science Center of Ural Division of RAS, Syktyvkar, Russia; Singapore, NUS, 18-29 May 2015 Workshop New Directions in Stein s method
More informationThe final is cumulative, but with more emphasis on chapters 3 and 4. There will be two parts.
Math 141 Review for Final The final is cumulative, but with more emphasis on chapters 3 and 4. There will be two parts. Part 1 (no calculator) graphing (polynomial, rational, linear, exponential, and logarithmic
More informationContinuous Univariate Distributions
Continuous Univariate Distributions Volume 2 Second Edition NORMAN L. JOHNSON University of North Carolina Chapel Hill, North Carolina SAMUEL KOTZ University of Maryland College Park, Maryland N. BALAKRISHNAN
More informationSOLUTION TO RUBEL S QUESTION ABOUT DIFFERENTIALLY ALGEBRAIC DEPENDENCE ON INITIAL VALUES GUY KATRIEL PREPRINT NO /4
SOLUTION TO RUBEL S QUESTION ABOUT DIFFERENTIALLY ALGEBRAIC DEPENDENCE ON INITIAL VALUES GUY KATRIEL PREPRINT NO. 2 2003/4 1 SOLUTION TO RUBEL S QUESTION ABOUT DIFFERENTIALLY ALGEBRAIC DEPENDENCE ON INITIAL
More information1.1 Basic Algebra. 1.2 Equations and Inequalities. 1.3 Systems of Equations
1. Algebra 1.1 Basic Algebra 1.2 Equations and Inequalities 1.3 Systems of Equations 1.1 Basic Algebra 1.1.1 Algebraic Operations 1.1.2 Factoring and Expanding Polynomials 1.1.3 Introduction to Exponentials
More informationChapter R - Basic Algebra Operations (94 topics, no due date)
Course Name: Math 00024 Course Code: N/A ALEKS Course: College Algebra Instructor: Master Templates Course Dates: Begin: 08/15/2014 End: 08/15/2015 Course Content: 207 topics Textbook: Barnett/Ziegler/Byleen/Sobecki:
More informationVALUES FOR THE CUMULATIVE DISTRIBUTION FUNCTION OF THE STANDARD MULTIVARIATE NORMAL DISTRIBUTION. Carol Lindee
VALUES FOR THE CUMULATIVE DISTRIBUTION FUNCTION OF THE STANDARD MULTIVARIATE NORMAL DISTRIBUTION Carol Lindee LindeeEmail@netscape.net (708) 479-3764 Nick Thomopoulos Illinois Institute of Technology Stuart
More informationIDENTIFIABILITY OF THE MULTIVARIATE NORMAL BY THE MAXIMUM AND THE MINIMUM
Surveys in Mathematics and its Applications ISSN 842-6298 (electronic), 843-7265 (print) Volume 5 (200), 3 320 IDENTIFIABILITY OF THE MULTIVARIATE NORMAL BY THE MAXIMUM AND THE MINIMUM Arunava Mukherjea
More informationAsymptotic Multivariate Kriging Using Estimated Parameters with Bayesian Prediction Methods for Non-linear Predictands
Asymptotic Multivariate Kriging Using Estimated Parameters with Bayesian Prediction Methods for Non-linear Predictands Elizabeth C. Mannshardt-Shamseldin Advisor: Richard L. Smith Duke University Department
More informationA Test for Order Restriction of Several Multivariate Normal Mean Vectors against all Alternatives when the Covariance Matrices are Unknown but Common
Journal of Statistical Theory and Applications Volume 11, Number 1, 2012, pp. 23-45 ISSN 1538-7887 A Test for Order Restriction of Several Multivariate Normal Mean Vectors against all Alternatives when
More informationOn Limit Distributions Arising from I terated Random Subdivisions of an Interval
On Limit Distributions Arising from I terated Rom Subdivisions of an Interval Norman L. Johnson University of North Carolina at Chapel Hill Samuel Kotz University of Maryl College Park ABSTRACT A method
More informationTHE LENGTH OF THE CONTINUED FRACTION EXPANSION FOR A CLASS OF RATIONAL FUNCTIONS IN f q {X)
Proceedings of the Edinburgh Mathematical Society (1991) 34, 7-17 THE LENGTH OF THE CONTINUED FRACTION EXPANSION FOR A CLASS OF RATIONAL FUNCTIONS IN f q {X) by ARNOLD KNOPFMACHER (Received 6th December
More information2 Functions of random variables
2 Functions of random variables A basic statistical model for sample data is a collection of random variables X 1,..., X n. The data are summarised in terms of certain sample statistics, calculated as
More informationContinuous Univariate Distributions
Continuous Univariate Distributions Volume 1 Second Edition NORMAN L. JOHNSON University of North Carolina Chapel Hill, North Carolina SAMUEL KOTZ University of Maryland College Park, Maryland N. BALAKRISHNAN
More informationNEW GENERATING FUNCTIONS FOR CLASSICAL POLYNOMIALS1. J. w. brown
NEW GENERATING FUNCTIONS FOR CLASSICAL POLYNOMIALS1 J. w. brown 1. Introduction. As recently as 1951 Brafman [l]-[4] initiated a series of papers in which he gave a variety of new unusual generating functions
More informationA NOTE ON A DISTRIBUTION OF WEIGHTED SUMS OF I.I.D. RAYLEIGH RANDOM VARIABLES
Sankhyā : The Indian Journal of Statistics 1998, Volume 6, Series A, Pt. 2, pp. 171-175 A NOTE ON A DISTRIBUTION OF WEIGHTED SUMS OF I.I.D. RAYLEIGH RANDOM VARIABLES By P. HITCZENKO North Carolina State
More informationThe Free Central Limit Theorem: A Combinatorial Approach
The Free Central Limit Theorem: A Combinatorial Approach by Dennis Stauffer A project submitted to the Department of Mathematical Sciences in conformity with the requirements for Math 4301 (Honour s Seminar)
More informationContents 1. Contents
Contents 1 Contents 6 Distributions of Functions of Random Variables 2 6.1 Transformation of Discrete r.v.s............. 3 6.2 Method of Distribution Functions............. 6 6.3 Method of Transformations................
More informationUnivariate Discrete Distributions
Univariate Discrete Distributions Second Edition NORMAN L. JOHNSON University of North Carolina Chapel Hill, North Carolina SAMUEL KOTZ University of Maryland College Park, Maryland ADRIENNE W. KEMP University
More informationH NT Z N RT L 0 4 n f lt r h v d lt n r n, h p l," "Fl d nd fl d " ( n l d n l tr l t nt r t t n t nt t nt n fr n nl, th t l n r tr t nt. r d n f d rd n t th nd r nt r d t n th t th n r lth h v b n f
More informationMathematics Standards for High School Financial Algebra A and Financial Algebra B
Mathematics Standards for High School Financial Algebra A and Financial Algebra B Financial Algebra A and B are two semester courses that may be taken in either order or one taken without the other; both
More informationn
p l p bl t n t t f Fl r d, D p rt nt f N t r l R r, D v n f nt r r R r, B r f l. n.24 80 T ll h, Fl. : Fl r d D p rt nt f N t r l R r, B r f l, 86. http://hdl.handle.net/2027/mdp.39015007497111 r t v n
More informationPR D NT N n TR T F R 6 pr l 8 Th Pr d nt Th h t H h n t n, D D r r. Pr d nt: n J n r f th r d t r v th tr t d rn z t n pr r f th n t d t t. n
R P RT F TH PR D NT N N TR T F R N V R T F NN T V D 0 0 : R PR P R JT..P.. D 2 PR L 8 8 J PR D NT N n TR T F R 6 pr l 8 Th Pr d nt Th h t H h n t n, D.. 20 00 D r r. Pr d nt: n J n r f th r d t r v th
More informationA Note on Some Wishart Expectations
Metrika, Volume 33, 1986, page 247-251 A Note on Some Wishart Expectations By R. J. Muirhead 2 Summary: In a recent paper Sharma and Krishnamoorthy (1984) used a complicated decisiontheoretic argument
More informationThe Distributions of Sums, Products and Ratios of Inverted Bivariate Beta Distribution 1
Applied Mathematical Sciences, Vol. 2, 28, no. 48, 2377-2391 The Distributions of Sums, Products and Ratios of Inverted Bivariate Beta Distribution 1 A. S. Al-Ruzaiza and Awad El-Gohary 2 Department of
More informationThe Distribution of Mixing Times in Markov Chains
The Distribution of Mixing Times in Markov Chains Jeffrey J. Hunter School of Computing & Mathematical Sciences, Auckland University of Technology, Auckland, New Zealand December 2010 Abstract The distribution
More informationThe Behavior of Multivariate Maxima of Moving Maxima Processes
The Behavior of Multivariate Maxima of Moving Maxima Processes Zhengjun Zhang Department of Mathematics Washington University Saint Louis, MO 6313-4899 USA Richard L. Smith Department of Statistics University
More informationTHE APPROXIMATE SOLUTION OF / = F(x,y)
PACIFIC JOURNAL OF MATHEMATICS Vol. 24, No. 2, 1968 THE APPROXIMATE SOLUTION OF / = F(x,y) R. G. HϋFFSTUTLER AND F. MAX STEIN In general the exact solution of the differential system y' = F(x, y), 7/(0)
More informationBasic Math Review for CS1340
Basic Math Review for CS1340 Dr. Mihail January 15, 2015 (Dr. Mihail) Math Review for CS1340 January 15, 2015 1 / 34 Sets Definition of a set A set is a collection of distinct objects, considered as an
More informationSection 9.7 and 9.10: Taylor Polynomials and Approximations/Taylor and Maclaurin Series
Section 9.7 and 9.10: Taylor Polynomials and Approximations/Taylor and Maclaurin Series Power Series for Functions We can create a Power Series (or polynomial series) that can approximate a function around
More informationApplied Probability and Stochastic Processes
Applied Probability and Stochastic Processes In Engineering and Physical Sciences MICHEL K. OCHI University of Florida A Wiley-Interscience Publication JOHN WILEY & SONS New York - Chichester Brisbane
More informationSolutions for Math 217 Assignment #3
Solutions for Math 217 Assignment #3 (1) Which of the following sets in R n are open? Which are closed? Which are neither open nor closed? (a) {(x, y) R 2 : x 2 y 2 = 1}. (b) {(x, y, z) R 3 : 0 < x + y
More informationNYS Algebra II and Trigonometry Suggested Sequence of Units (P.I's within each unit are NOT in any suggested order)
1 of 6 UNIT P.I. 1 - INTEGERS 1 A2.A.1 Solve absolute value equations and inequalities involving linear expressions in one variable 1 A2.A.4 * Solve quadratic inequalities in one and two variables, algebraically
More informationVariance Function Estimation in Multivariate Nonparametric Regression
Variance Function Estimation in Multivariate Nonparametric Regression T. Tony Cai 1, Michael Levine Lie Wang 1 Abstract Variance function estimation in multivariate nonparametric regression is considered
More informationEcon 508B: Lecture 5
Econ 508B: Lecture 5 Expectation, MGF and CGF Hongyi Liu Washington University in St. Louis July 31, 2017 Hongyi Liu (Washington University in St. Louis) Math Camp 2017 Stats July 31, 2017 1 / 23 Outline
More informationBasic Math Review for CS4830
Basic Math Review for CS4830 Dr. Mihail August 18, 2016 (Dr. Mihail) Math Review for CS4830 August 18, 2016 1 / 35 Sets Definition of a set A set is a collection of distinct objects, considered as an object
More informationMathematics High School Algebra
Mathematics High School Algebra Expressions. An expression is a record of a computation with numbers, symbols that represent numbers, arithmetic operations, exponentiation, and, at more advanced levels,
More informationAlgebra 2 (3 rd Quad Expectations) CCSS covered Key Vocabulary Vertical
Algebra 2 (3 rd Quad Expectations) CCSS covered Key Vocabulary Vertical Chapter (McGraw-Hill Algebra 2) Chapter 7 (Suggested Pacing 14 Days) Lesson 7-1: Graphing Exponential Functions Lesson 7-2: Solving
More informationSchool District of Marshfield Course Syllabus
School District of Marshfield Course Syllabus Course Name: Algebra II Length of Course: 1 Year Credit: 1 Program Goal: The School District of Marshfield Mathematics Program will prepare students for college
More informationOn the exact computation of the density and of the quantiles of linear combinations of t and F random variables
Journal of Statistical Planning and Inference 94 3 www.elsevier.com/locate/jspi On the exact computation of the density and of the quantiles of linear combinations of t and F random variables Viktor Witkovsky
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 informationRoot Finding For NonLinear Equations Bisection Method
Root Finding For NonLinear Equations Bisection Method P. Sam Johnson November 17, 2014 P. Sam Johnson (NITK) Root Finding For NonLinear Equations Bisection MethodNovember 17, 2014 1 / 26 Introduction The
More informationPart 6: Multivariate Normal and Linear Models
Part 6: Multivariate Normal and Linear Models 1 Multiple measurements Up until now all of our statistical models have been univariate models models for a single measurement on each member of a sample of
More informationNegative Multinomial Model and Cancer. Incidence
Generalized Linear Model under the Extended Negative Multinomial Model and Cancer Incidence S. Lahiri & Sunil K. Dhar Department of Mathematical Sciences, CAMS New Jersey Institute of Technology, Newar,
More informationACM 116: Lectures 3 4
1 ACM 116: Lectures 3 4 Joint distributions The multivariate normal distribution Conditional distributions Independent random variables Conditional distributions and Monte Carlo: Rejection sampling Variance
More informationREQUIRED MATHEMATICAL SKILLS FOR ENTERING CADETS
REQUIRED MATHEMATICAL SKILLS FOR ENTERING CADETS The Department of Applied Mathematics administers a Math Placement test to assess fundamental skills in mathematics that are necessary to begin the study
More informationMAXIMUM LIKELIHOOD ESTIMATION FOR THE GROWTH CURVE MODEL WITH UNEQUAL DISPERSION MATRICES. S. R. Chakravorti
Zo. J MAXIMUM LIKELIHOOD ESTIMATION FOR THE GROWTH CURVE MODEL WITH UNEQUAL DISPERSION MATRICES By S. R. Chakravorti Department of Biostatistics University of North Carolina and University of Calcutta,
More informationAlgebra I Number and Quantity The Real Number System (N-RN)
Number and Quantity The Real Number System (N-RN) Use properties of rational and irrational numbers N-RN.3 Explain why the sum or product of two rational numbers is rational; that the sum of a rational
More informationMATH 205C: STATIONARY PHASE LEMMA
MATH 205C: STATIONARY PHASE LEMMA For ω, consider an integral of the form I(ω) = e iωf(x) u(x) dx, where u Cc (R n ) complex valued, with support in a compact set K, and f C (R n ) real valued. Thus, I(ω)
More information- 2 ' a 2 =- + _ y'2 + _ 21/4. a 3 =! +! y'2 +! 21/4 +!. / (!:_ +! y'2) 21/ Y 2 2 c!+!v'2+!21/4). lc!+!v'2) 21/4.
....................................... ~. The Arithmetic Geometric Mean The arithmetic mean of two numbers a and b is defined as the "average" of the numbers, namely, aib while the geometric mean is given
More informationA Library of Functions
LibraryofFunctions.nb 1 A Library of Functions Any study of calculus must start with the study of functions. Functions are fundamental to mathematics. In its everyday use the word function conveys to us
More informationMath 361: Homework 1 Solutions
January 3, 4 Math 36: Homework Solutions. We say that two norms and on a vector space V are equivalent or comparable if the topology they define on V are the same, i.e., for any sequence of vectors {x
More informationAlgebra II Curriculum Crosswalk
Algebra II Curriculum Crosswalk The following document is to be used to compare the 2003 North Carolina Mathematics Course of Study for Algebra II and the State s for Mathematics for Algebra II. As noted
More information828.^ 2 F r, Br n, nd t h. n, v n lth h th n l nd h d n r d t n v l l n th f v r x t p th l ft. n ll n n n f lt ll th t p n nt r f d pp nt nt nd, th t
2Â F b. Th h ph rd l nd r. l X. TH H PH RD L ND R. L X. F r, Br n, nd t h. B th ttr h ph rd. n th l f p t r l l nd, t t d t, n n t n, nt r rl r th n th n r l t f th f th th r l, nd d r b t t f nn r r pr
More informationJ2 e-*= (27T)- 1 / 2 f V* 1 '»' 1 *»
THE NORMAL APPROXIMATION TO THE POISSON DISTRIBUTION AND A PROOF OF A CONJECTURE OF RAMANUJAN 1 TSENG TUNG CHENG 1. Summary. The Poisson distribution with parameter X is given by (1.1) F(x) = 23 p r where
More informationEngg. Math. II (Unit-IV) Numerical Analysis
Dr. Satish Shukla of 33 Engg. Math. II (Unit-IV) Numerical Analysis Syllabus. Interpolation and Curve Fitting: Introduction to Interpolation; Calculus of Finite Differences; Finite Difference and Divided
More information4 4 N v b r t, 20 xpr n f th ll f th p p l t n p pr d. H ndr d nd th nd f t v L th n n f th pr v n f V ln, r dn nd l r thr n nt pr n, h r th ff r d nd
n r t d n 20 20 0 : 0 T P bl D n, l d t z d http:.h th tr t. r pd l 4 4 N v b r t, 20 xpr n f th ll f th p p l t n p pr d. H ndr d nd th nd f t v L th n n f th pr v n f V ln, r dn nd l r thr n nt pr n,
More informationTrimester 2 Expectations. Chapter (McGraw-Hill. CCSS covered Key Vocabulary Vertical. Alignment
Algebra 2 Trimester 2 Expectations Chapter (McGraw-Hill Algebra 2) Chapter 5 (Suggested Pacing 14 Days) Polynomials and Polynomial Functions Lesson 5-1: Operations with Polynomials Lesson 5-2: Dividing
More informationarxiv: v1 [stat.me] 28 Mar 2011
arxiv:1103.5447v1 [stat.me] 28 Mar 2011 On matrix variance inequalities G. Afendras and N. Papadatos Department of Mathematics, Section of Statistics and O.R., University of Athens, Panepistemiopolis,
More informationReview (Probability & Linear Algebra)
Review (Probability & Linear Algebra) CE-725 : Statistical Pattern Recognition Sharif University of Technology Spring 2013 M. Soleymani Outline Axioms of probability theory Conditional probability, Joint
More informationPre-calculus 12 Curriculum Outcomes Framework (110 hours)
Curriculum Outcomes Framework (110 hours) Trigonometry (T) (35 40 hours) General Curriculum Outcome: Students will be expected to develop trigonometric reasoning. T01 Students will be expected to T01.01
More informationOn Bayesian Inference with Conjugate Priors for Scale Mixtures of Normal Distributions
Journal of Applied Probability & Statistics Vol. 5, No. 1, xxx xxx c 2010 Dixie W Publishing Corporation, U. S. A. On Bayesian Inference with Conjugate Priors for Scale Mixtures of Normal Distributions
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 informationn r t d n :4 T P bl D n, l d t z d th tr t. r pd l
n r t d n 20 20 :4 T P bl D n, l d t z d http:.h th tr t. r pd l 2 0 x pt n f t v t, f f d, b th n nd th P r n h h, th r h v n t b n p d f r nt r. Th t v v d pr n, h v r, p n th pl v t r, d b p t r b R
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 informationTAYLOR AND MACLAURIN SERIES
TAYLOR AND MACLAURIN SERIES. Introduction Last time, we were able to represent a certain restricted class of functions as power series. This leads us to the question: can we represent more general functions
More informationP(I -ni < an for all n > in) = 1 - Pm# 1
ITERATED LOGARITHM INEQUALITIES* By D. A. DARLING AND HERBERT ROBBINS UNIVERSITY OF CALIFORNIA, BERKELEY Communicated by J. Neyman, March 10, 1967 1. Introduction.-Let x,x1,x2,... be a sequence of independent,
More informationPascal Eigenspaces and Invariant Sequences of the First or Second Kind
Pascal Eigenspaces and Invariant Sequences of the First or Second Kind I-Pyo Kim a,, Michael J Tsatsomeros b a Department of Mathematics Education, Daegu University, Gyeongbu, 38453, Republic of Korea
More informationAlgebra 2 Standards. Essential Standards:
Benchmark 1: Essential Standards: 1. Alg2.M.F.LE.A.02 (linear): I can create linear functions if provided either a graph, relationship description or input-output tables. - 15 Days 2. Alg2.M.A.APR.B.02a
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 informationStat 451 Lecture Notes Numerical Integration
Stat 451 Lecture Notes 03 12 Numerical Integration Ryan Martin UIC www.math.uic.edu/~rgmartin 1 Based on Chapter 5 in Givens & Hoeting, and Chapters 4 & 18 of Lange 2 Updated: February 11, 2016 1 / 29
More informationTAYLOR SERIES [SST 8.8]
TAYLOR SERIES [SST 8.8] TAYLOR SERIES: Every function f C (c R, c + R) has a unique Taylor series about x = c of the form: f (k) (c) f(x) = (x c) k = f(c) + f (c) (x c) + f (c) (x c) 2 + f (c) (x c) 3
More informationFitting Linear Statistical Models to Data by Least Squares I: Introduction
Fitting Linear Statistical Models to Data by Least Squares I: Introduction Brian R. Hunt and C. David Levermore University of Maryland, College Park Math 420: Mathematical Modeling January 25, 2012 version
More informationMULTISECTION METHOD AND FURTHER FORMULAE FOR π. De-Yin Zheng
Indian J. pure appl. Math., 39(: 37-55, April 008 c Printed in India. MULTISECTION METHOD AND FURTHER FORMULAE FOR π De-Yin Zheng Department of Mathematics, Hangzhou Normal University, Hangzhou 30036,
More informationApplications of Basu's TheorelTI. Dennis D. Boos and Jacqueline M. Hughes-Oliver I Department of Statistics, North Car-;'lina State University
i Applications of Basu's TheorelTI by '. Dennis D. Boos and Jacqueline M. Hughes-Oliver I Department of Statistics, North Car-;'lina State University January 1997 Institute of Statistics ii-limeo Series
More informationAlgebra II Unit Breakdown (Curriculum Map Outline)
QUARTER 1 Unit 1 (Arithmetic and Geometric Sequences) A. Sequences as Functions 1. Identify finite and infinite sequences 2. Identify an arithmetic sequence and its parts 3. Identify a geometric sequence
More informationAlgebra Performance Level Descriptors
Limited A student performing at the Limited Level demonstrates a minimal command of Ohio s Learning Standards for Algebra. A student at this level has an emerging ability to A student whose performance
More informationB. L. Raktoe* and W. T. Federer University of Guelph and Cornell University. Abstract
BALANCED OPTIMAL SA'IURATED MAIN EFFECT PLANS OF 'IHE 2n FACTORIAL AND THEIR RELATION TO (v,k,'x.) CONFIGURATIONS BU-406-M by January, 1S72 B. L. Raktoe* and W. T. Federer University of Guelph and Cornell
More informationZero-Mean Difference Discrimination and the. Absolute Linear Discriminant Function. Peter A. Lachenbruch
Zero-Mean Difference Discrimination and the Absolute Linear Discriminant Function By Peter A. Lachenbruch University of California, Berkeley and University of North Carolina, Chapel Hill Institute of Statistics
More informationGENERATING FUNCTIONS OF CENTRAL VALUES IN GENERALIZED PASCAL TRIANGLES. CLAUDIA SMITH and VERNER E. HOGGATT, JR.
GENERATING FUNCTIONS OF CENTRAL VALUES CLAUDIA SMITH and VERNER E. HOGGATT, JR. San Jose State University, San Jose, CA 95. INTRODUCTION In this paper we shall examine the generating functions of the central
More informationUniversity of Lisbon, Portugal
Development and comparative study of two near-exact approximations to the distribution of the product of an odd number of independent Beta random variables Luís M. Grilo a,, Carlos A. Coelho b, a Dep.
More informationMATHEMATICS. Higher 2 (Syllabus 9740)
MATHEMATICS Higher (Syllabus 9740) CONTENTS Page AIMS ASSESSMENT OBJECTIVES (AO) USE OF GRAPHING CALCULATOR (GC) 3 LIST OF FORMULAE 3 INTEGRATION AND APPLICATION 3 SCHEME OF EXAMINATION PAPERS 3 CONTENT
More informationWA State Common Core Standards - Mathematics
Number & Quantity The Real Number System Extend the properties of exponents to rational exponents. 1. Explain how the definition of the meaning of rational exponents follows from extending the properties
More informationMath 160 Final Exam Info and Review Exercises
Math 160 Final Exam Info and Review Exercises Fall 2018, Prof. Beydler Test Info Will cover almost all sections in this class. This will be a 2-part test. Part 1 will be no calculator. Part 2 will be scientific
More informationON THE MOMENTS OF ITERATED TAIL
ON THE MOMENTS OF ITERATED TAIL RADU PĂLTĂNEA and GHEORGHIŢĂ ZBĂGANU The classical distribution in ruins theory has the property that the sequence of the first moment of the iterated tails is convergent
More informationOn Expected Gaussian Random Determinants
On Expected Gaussian Random Determinants Moo K. Chung 1 Department of Statistics University of Wisconsin-Madison 1210 West Dayton St. Madison, WI 53706 Abstract The expectation of random determinants whose
More informationContinuous RVs. 1. Suppose a random variable X has the following probability density function: π, zero otherwise. f ( x ) = sin x, 0 < x < 2
STAT 4 Exam I Continuous RVs Fall 7 Practice. Suppose a random variable X has the following probability density function: f ( x ) = sin x, < x < π, zero otherwise. a) Find P ( X < 4 π ). b) Find µ = E
More informationAMSCO Algebra 2. Number and Quantity. The Real Number System
AMSCO Algebra 2 Number and Quantity The Real Number System Extend the properties of exponents to rational exponents. N-RN.1 Explain how the definition of the meaning of rational exponents follows from
More informationHANDBOOK OF APPLICABLE MATHEMATICS
HANDBOOK OF APPLICABLE MATHEMATICS Chief Editor: Walter Ledermann Volume II: Probability Emlyn Lloyd University oflancaster A Wiley-Interscience Publication JOHN WILEY & SONS Chichester - New York - Brisbane
More information