Information Matrix for Pareto(IV), Burr, and Related Distributions
|
|
- Madeleine Little
- 6 years ago
- Views:
Transcription
1 MARCEL DEKKER INC. 7 MADISON AVENUE NEW YORK NY 6 3 Marcel Dekker Inc. All rights reserved. This material may not be used or reproduced in any form without the express written permission of Marcel Dekker Inc. COMMUNICATIONS IN STATISTICS Theory and Methods Vol. 3 No. pp Information Matrix for Pareto(IV) Burr and Related Distributions Vytaras Brazauskas* Department of Mathematical Sciences University of Wisconsin-Milwaukee Milwaukee Wisconsin USA ABSTRACT In this paper the exact form of information matrix for Pareto(IV) and related distributions is determined. The Pareto(IV) family being very general includes more specialized families of Pareto(I) Pareto(II) and Pareto(III) and the Burr family of distributions as special cases. These distributions for example arise as tractable parametric models in actuarial science economics finance and telecommunications. Additionally a useful mathematical result with its own domain of importance is obtained. In particular explicit formula for the improper integral R ððlog xþ m =ð þ xþ þb Þ dx with b > and non-negative integer m is derived. *Correspondence: Vytaras Brazauskas Department of Mathematical Sciences University of Wisconsin-Milwaukee P.O. Box 43 Milwaukee WI USA; vytarasuwm.edu. 35 DOI:.8/STA-888 Copyright & 3 by Marcel Dekker Inc (Print); 53-45X (Online)
2 MARCEL DEKKER INC. 7 MADISON AVENUE NEW YORK NY 6 3 Marcel Dekker Inc. All rights reserved. This material may not be used or reproduced in any form without the express written permission of Marcel Dekker Inc. 36 Brazauskas Key Words: Gamma function; Information; Pareto models; Polygamma functions.. INTRODUCTION In this paper the exact form of (Fisher) information matrix for Pareto(IV) and related distributions is determined. It is well-known that information matrix serves as a valuable tool for derivation of covariance matrix in the asymptotic distribution of maximum likelihood estimators (MLE). As discussed in Serfling (98) Sec. 4. under suitable regularity conditions the determinant (divided by the sample size) of the asymptotic covariance matrix of MLE reaches an optimal lower bound for the volume of the spread ellipsoid of joint estimators. In the univariate case this optimality property of MLE is widely used in the robustness versus efficiency studies as a quantitative benchmark for efficiency considerations. See for example Brazauskas and Serfling (a b) Hampel et al. (986) Huber (98) Kimber (983a 983b) and Lehmann (983) Chapter 5. The Pareto(IV) family being very general includes more specialized families of Pareto(I) Pareto(II) and Pareto(III) and the Burr family of distributions as special cases. These distributions are suitable for situations involving relatively high probability in the upper tail. More specifically such models have been formulated in the context of actuarial science economics finance and teletraffic for example for distributions of variables such as sizes of insurance claims sizes of firms incomes in a population of people stock price fluctuations and length of telephone calls. (See Arnold (983) and Johnson et al. (994) Chapter 9 for a broad discussion of Pareto models and their diverse applications.) New application contexts continue to be found. For example Crato et al. (997) have recently discovered Pareto-type tail behavior in the cost distributions of combinatorial search algorithms. As a by-product of the main result a useful mathematical formula (which plays a very important role in our derivations) is obtained. In particular explicit formula for the improper integral ðlog xþ m dx þb ð þ xþ ð:þ with b > and non-negative integer m is derived. Surprisingly however computation of integral (.) has not been addressed in the mathematical
3 MARCEL DEKKER INC. 7 MADISON AVENUE NEW YORK NY 6 3 Marcel Dekker Inc. All rights reserved. This material may not be used or reproduced in any form without the express written permission of Marcel Dekker Inc. Information Matrix for Pareto(IV) 37 literature including such comprehensive handbooks of mathematics as Abramowitz and Stegun (97) and Harris and Stocker (998). A hierarchy of Pareto models as well as their relation to the Burr family is discussed in Sec.. Elements of the information matrix for Pareto(IV) Burr and related distributions are computed in Sec. 3. Derivation of integral (.) along with other intermediate integrals and formulas are presented in the Appendix.. PARETO(IV) AND RELATED DISTRIBUTIONS As discussed in Arnold (983) Chapter 3 a hierarchy of Pareto distributions is established by starting with the classical Pareto(I) distribution and subsequently introducing additional parameters which relate to location scale shape and inequality. Such an approach leads to a very general family of distributions called the Pareto(IV) family with the cumulative distribution function FðxÞ þ x = x > ð:þ where <<þ is the location parameter > is the scale parameter > is the inequality parameter and > is the shape parameter which characterizes the tail of the distribution. We denote this distribution by Pareto(IV) ð Þ. Note that in general statistical science there is no such type of parameters like inequality. Parameter is called the inequality parameter because of its interpretation in the economics context. That is if we choose and in expression (.) then parameter ( ) is precisely the Gini index of inequality. Clearly the other three types of Pareto distributions can be identified as special cases of the Pareto(IV) family by appropriately choosing parameters in Eq. (.) (see Arnold (983) pp ): Pareto(I) ð Þ Pareto(IV) ð Þ Pareto(II) ð Þ Pareto(IV) ð Þ Pareto(III) ð Þ Pareto(IV) ð Þ: The Burr family of distributions is also sufficiently flexible and enjoys long popularity in the actuarial science literature (see for example
4 MARCEL DEKKER INC. 7 MADISON AVENUE NEW YORK NY 6 3 Marcel Dekker Inc. All rights reserved. This material may not be used or reproduced in any form without the express written permission of Marcel Dekker Inc. 38 Brazauskas Daykin et al. (994) and Klugman et al. (998)). However even this family can be treated as a special case of Pareto(IV): Burr ð Þ Pareto(IV) ð = Þ (see Klugman et al. (998) p. 574). In order to make the Pareto(IV) distribution a regular family we assume that parameter is known and without loss of generality equal to. (For regularity conditions on F and their interpretation see for example Serfling (98) pp ) Also we note that this assumption is not too restrictive for modeling purposes because in typical applications the lower limit of variables of interest is known. In insurance and reinsurance context for example the lower limit is pre-defined by a contract and can be represented as a deductible or a retention level. (This perhaps is one of the main reasons why generality of the Burr distribution remains sufficient for actuarial modeling.) 3. INFORMATION MATRIX FOR PARETO(IV) Suppose X is a random variable with the probability density function f ðþ where ð... k Þ. Then the information matrix IðÞ is the k k matrix with elements log f I ij ðþ E ðxþ log f ðxþ : i j For the Pareto(IV) ð Þ distribution we have ð 3 Þ ð Þ and the density function ðx=þ = f ðxþ ð þðx=þ = Þ þ : ð3:þ Thus the required partial derivatives are log f ðxþ log f ðxþ log f ðxþ 3 log f ðxþ log f ðxþ log f ðxþ þ þðx=þ = log x þ logðx=þ þðx=þ = log þ x = þ :
5 MARCEL DEKKER INC. 7 MADISON AVENUE NEW YORK NY 6 3 Marcel Dekker Inc. All rights reserved. This material may not be used or reproduced in any form without the express written permission of Marcel Dekker Inc. Information Matrix for Pareto(IV) 39 Note that for this problem the first two partial derivatives = i j of f ðþ exist and therefore an alternative (based on the second-order derivatives) formula for computation of information matrix elements may be used. We found however that except for a couple of cases this approach does not simplify our derivations. Therefore it is not pursued here. Since the information matrix IðÞ is symmetric it is enough to find elements I ðþ I ðþ I 3 ðþ I ðþ I 3 ðþ and I 33 ðþ. Derivation of these elements is based on the following strategy: first we express each I ij ðþ i j 3 in terms of integrals A A5 B B4 C C D D which are defined (and their explicit formulas are presented) in Appendix A.; then tedious algebraic simplifications yield the following formulas. log f ðxþ I ðþ f ðxþ dx ð þ Þ ðþ A þ þ A ðþ ð þ Þ log f ðxþ log f ðxþ I ðþ f ðxþ dx ð þ Þ 3 B 3 ð þ Þ I 3 ðþ log f ðxþ B ð þ Þ 3 ðþ ðþ ðþ þ log f ðxþ f ðxþ dx B3 þ þ þ A3 A þ þ A5 ð þ Þ log f ðxþ I ðþ f ðxþ dx ð þ Þ B4 þ 4 4 C þ ð þ Þ 4 ð þ Þ C ðþ ðþ þ ðþþ ð þ Þ ðþ ðþ ðþ ð þ Þ B þ 3 3 B ð Þ þ ð þ Þ A ðþ ðþ ðþ
6 MARCEL DEKKER INC. 7 MADISON AVENUE NEW YORK NY 6 3 Marcel Dekker Inc. All rights reserved. This material may not be used or reproduced in any form without the express written permission of Marcel Dekker Inc. 3 Brazauskas I 3 ðþ log f ðxþ log f ðxþ f ðxþ dx þ B B D þ þ D þ A3 ðþ ð þ Þ ðþ ðþ log f ðxþ I 33 ðþ f ðxþ dx A3 þ A4 : Finally elements I ðþ I ðþandi 3 ðþ can be written in a simpler form by using the polygamma functions ðnþ ðaþ ðd n =da n Þð ðaþ=ðaþþ for a > and integer n. Specifically we use digamma ðaþ ðþ ðaþ and trigamma ðaþ functions. Thus the information matrix I P ðþ for the Pareto(IV) ð Þ distribution is given by: ðþ ð þ Þ ½ ðþ ðþþš ð þ Þ B ð þ Þ ½ ðþ ðþþš ð þ Þ ð ðþ ðþþ þ ðþþ ðþ þ ð ðþ ðþþ ð þ Þ ðþ ðþ ð þ Þ ð þ Þ ðþ ðþ ð þ Þ C A 3.. Special Cases 3... Burr( ) Distribution Since the Burr distribution is a reparametrization of Pareto(IV) ð Þ it follows from Lehmann (983) Sec..7 that its information matrix I B ðþ can be derived from J I P ðþ J where J is the Jacobian matrix of the transformation of variables. Thus I B ðþ is then given by: B ðþ ð þ Þ ½ ðþ ðþšþ ð þ Þ ð þ Þ ½ ðþ ðþšþ ð þ Þ ð ðþ ðþþ þ ðþþ ðþ þ ð ðþ ðþþ þ ½ ðþ ðþþš þ ð þ Þ ½ ðþ ðþþš þ C A
7 MARCEL DEKKER INC. 7 MADISON AVENUE NEW YORK NY 6 3 Marcel Dekker Inc. All rights reserved. This material may not be used or reproduced in any form without the express written permission of Marcel Dekker Inc. Information Matrix for Pareto(IV) Pareto(III) ( ) Distribution This is a special case of Pareto(IV) with. Therefore third row and third column (these represent information about parameter ) in I P ðþ vanish. And into expressions of the remaining elements we substitute. This yields I 3ðÞ ð; Þ B þ C ðþ A 3ðÞ B C þ 3 A Pareto(II) ( ) Distribution This is a special case of Pareto(IV) with. Therefore second row and second column (these represent information about parameter ) in I P ðþ vanish. And into the remaining formulas we substitute. This yields I ð þ Þ ð;þ B ð þ Þ ð þ Þ C A It should be noted here that information matrices I ð Þ and I ð Þ are readily available in Arnold (983) p.. APPENDIX A.. Three Lemmas Lemma. For b > and non-negative integer m ðlog xþ m ð þ xþ þb dx ð þ bþ X m i m ðþ mi ðiþ ðþ ðmiþ ðbþ i ða:þ
8 MARCEL DEKKER INC. 7 MADISON AVENUE NEW YORK NY 6 3 Marcel Dekker Inc. All rights reserved. This material may not be used or reproduced in any form without the express written permission of Marcel Dekker Inc. 3 Brazauskas where ðaþ R x a e x dx and ðnþ ðaþ R ðlog xþn x a e x dx are gamma function and its n-th order derivative respectively. Proof. We start with the following fact (Abramowitz and Stegun (97) p. 55): x a ðaþ t a e xt dt x > a > : Next after applying the above formula to the integral in Eq. (A.) and interchanging the order of integration we have ðlog xþ m ð þ xþ þb dx ð þ bþ t b e t e tx ðlog xþ m dx dt: Further substitution of variables z tx and application of the binomial formula to ðlog z log tþ m lead to the following expression Z Z " # X t b e t m m ðlog zþ i ð log tþ mi e z dz dt: ð þ bþ i i Finally simple reorganization of terms yields the result. v Lemma. For b > and integer n ðnþ ðbþ ðb Þ ðnþ ðb Þþn ðnþ ðb Þ ða:þ where ðþ and ðnþ ðþ denote gamma function and its n-th order derivative respectively. Proof. Integration by parts of the first derivative of the gamma function ðbþ leads to ðbþ ðb Þ ðb Þþðb Þ: Differentiation of this equation n times yields Eq. (A.). v Lemma 3. For b > logð þ xþ log x ð þ xþ þb dx b ðbþ ðbþ ðbþ ðbþ b! ðbþ ðbþ ðþ ða:3þ
9 MARCEL DEKKER INC. 7 MADISON AVENUE NEW YORK NY 6 3 Marcel Dekker Inc. All rights reserved. This material may not be used or reproduced in any form without the express written permission of Marcel Dekker Inc. Information Matrix for Pareto(IV) 33 where ðþ and ðnþ ðþ denote gamma function and its n-th order derivative respectively. Proof. Note that logð þ xþ log x ð þ xþ þb dx d db log x ð þ xþ þb dx : Formula (A.) and differentiation of the right-hand side yield (A.3). v A.. Useful Integrals In all expressions below function f ðxþ is the Pareto(IV) ð Þ density function given by Eq. (3.). Integrals AA5 are derived via straightforward integration. A A A3 A4 A5 f ðxþ dx = þðx=þ þ f ðxþ þðx=þ = dx þ log þðx=þ = f ðxþ dx log þðx=þ = f ðxþ dx log þðx=þ = þðx=þ = f ðxþ dx ð þ Þ : Integrals B B4 are computed by making use of formula (A.). B B B3 B4 logðx=þ f ðxþ dx ðþ ðþ ðþ logðx=þ f ðxþ dx = þðx=þ þ ðþ ðþ ðþ logðx=þ ð þðx=þ = f ðxþ dx Þ þ log ðx=þ f ðxþ dx ðþ ðþ ðþ ðþ þ ðþ ðþ þ : ðþ ðþ ðþ
10 MARCEL DEKKER INC. 7 MADISON AVENUE NEW YORK NY 6 3 Marcel Dekker Inc. All rights reserved. This material may not be used or reproduced in any form without the express written permission of Marcel Dekker Inc. 34 Brazauskas Lemmas and are applied in computation of integrals C and C. log ðx=þ C f ðxþ dx = þðx=þ þ ðþ ðþ ðþ ðþ þ ðþ ðþ ðþ log ðx=þ C ð þðx=þ = f ðxþ dx Þ þ ðþ ðþ þ þ þ ðþ ðþ þ ð þ Þ ðþ ðþ þ! þ ðþ : Integrals D and D are derived by employing Lemmas 3. D D þ logðx=þ logð þðx=þ = Þ f ðxþ dx ðþ ðþ ðþ ðþ! ðþ ðþ ðþ logðx=þ logð þðx=þ = Þ þðx=þ = f ðxþ dx ðþ ðþ ðþ ðþ þ ðþ ðþ ðþ! þ : ð þ Þ REFERENCES Abramowitz M. Stegun I. A. (97). Handbook of Mathematical Functions. National Bureau of Standards Applied Mathematics Series No. 55. Arnold B. C. (983). Pareto Distributions. Fairland Maryland: International Cooperative Publishing House. Brazauskas V. Serfling R. (a). Robust estimation of tail parameters for two-parameter Pareto and exponential models via generalized quantile statistics. Extremes 3(3):3 49. Brazauskas V. Serfling R. (b). Robust and efficient estimation of the tail index of a single-parameter Pareto distribution. North American Actuarial Journal 4(4): 7.
11 MARCEL DEKKER INC. 7 MADISON AVENUE NEW YORK NY 6 3 Marcel Dekker Inc. All rights reserved. This material may not be used or reproduced in any form without the express written permission of Marcel Dekker Inc. Information Matrix for Pareto(IV) 35 Daykin C. D. Pentika inen T. Pesonen M. (994). Practical Risk Theory for Actuaries. London: Chapman and Hall. Gomes C. P. Selman B. Crato N. (997). Heavy-tailed distributions in combinatorial search. In: Smolka G. ed. Principles and Practice of Constraint Programming CP-97. Lecture Notes in Computer Science. Vol Hampel F. R. Ronchetti E. M. Rousseeuw P. J. Stahel W. A. (986). Robust Statistics: The Approach Based On Influence Functions. New York: Wiley. Harris J. W. Stocker H. (998). Handbook of Mathematics and Computational Science. New York: Springer. Huber P. J. (98). Robust Statistics. New York: Wiley. Johnson N. L. Kotz S. Balakrishnan N. (994). Continuous Univariate Distributions. Vol.. nd ed. New York: Wiley. Kimber A. C. (983a). Trimming in gamma samples. Applied Statistics 3():7 4. Kimber A. C. (983b). Comparison of some robust estimators of scale in gamma samples with known shape. Journal of Statistical Computation and Simulation 8: Klugman S. A. Panjer H. H. Willmot G. E. (998). Loss Models: From Data to Decisions. New York: Wiley. Lehmann E. L. (983). Theory of Point Estimation. New York: Wiley. Serfling R. J. (98). Approximation Theorems of Mathematical Statistics. New York: Wiley.
A Statistical Application of the Quantile Mechanics Approach: MTM Estimators for the Parameters of t and Gamma Distributions
Under consideration for publication in Euro. Jnl of Applied Mathematics A Statistical Application of the Quantile Mechanics Approach: MTM Estimators for the Parameters of t and Gamma Distributions A. KLEEFELD
More informationShannon entropy in generalized order statistics from Pareto-type distributions
Int. J. Nonlinear Anal. Appl. 4 (203 No., 79-9 ISSN: 2008-6822 (electronic http://www.ijnaa.semnan.ac.ir Shannon entropy in generalized order statistics from Pareto-type distributions B. Afhami a, M. Madadi
More informationROBUST AND EFFICIENT ESTIMATION OF THE TAIL INDEX OF A SINGLE-PARAMETER PARETO DISTRIBUTION
ROBUST AND EFFICIENT ESTIMATION OF THE TAIL INDEX OF A SINGLE-PARAMETER PARETO DISTRIBUTION Vytaras Brazausas* and Robert Serfling ABSTRACT Estimation of the tail index parameter of a single-parameter
More informationOPTIMAL B-ROBUST ESTIMATORS FOR THE PARAMETERS OF THE GENERALIZED HALF-NORMAL DISTRIBUTION
REVSTAT Statistical Journal Volume 15, Number 3, July 2017, 455 471 OPTIMAL B-ROBUST ESTIMATORS FOR THE PARAMETERS OF THE GENERALIZED HALF-NORMAL DISTRIBUTION Authors: Fatma Zehra Doğru Department of Econometrics,
More informationRobust and Efficient Fitting of Loss Models: Diagnostic Tools and Insights
Robust and Efficient Fitting of Loss Models: Diagnostic Tools and Insights Vytaras Brazauskas 1 University of Wisconsin-Milwaukee Abstract We consider robust and efficient fitting of claim severity models
More informationCONSISTENT ESTIMATION THROUGH WEIGHTED HARMONIC MEAN OF INCONSISTENT ESTIMATORS IN REPLICATED MEASUREMENT ERROR MODELS
ECONOMETRIC REVIEWS, 20(4), 507 50 (200) CONSISTENT ESTIMATION THROUGH WEIGHTED HARMONIC MEAN OF INCONSISTENT ESTIMATORS IN RELICATED MEASUREMENT ERROR MODELS Shalabh Department of Statistics, anjab University,
More informationRobust and Efficient Fitting of the Generalized Pareto Distribution with Actuarial Applications in View
Robust and Efficient Fitting of the Generalized Pareto Distribution with Actuarial Applications in View Vytaras Brazauskas 1 University of Wisconsin-Milwaukee Andreas Kleefeld 2 University of Wisconsin-Milwaukee
More informationRobust estimation of the Pareto index: A Monte Carlo Analysis
Working Papers No. 32/2013 (117) MICHAŁ BRZEZIŃSKI Robust estimation of the Pareto index: A Monte Carlo Analysis Warsaw 2013 Robust estimation of the Pareto index: A Monte Carlo Analysis MICHAL BRZEZINSKI
More informationAsymptotic efficiency and small sample power of a locally most powerful linear rank test for the log-logistic distribution
Math Sci (2014) 8:109 115 DOI 10.1007/s40096-014-0135-4 ORIGINAL RESEARCH Asymptotic efficiency and small sample power of a locally most powerful linear rank test for the log-logistic distribution Hidetoshi
More informationAPPROXIMATING THE GENERALIZED BURR-GAMMA WITH A GENERALIZED PARETO-TYPE OF DISTRIBUTION A. VERSTER AND D.J. DE WAAL ABSTRACT
APPROXIMATING THE GENERALIZED BURR-GAMMA WITH A GENERALIZED PARETO-TYPE OF DISTRIBUTION A. VERSTER AND D.J. DE WAAL ABSTRACT In this paper the Generalized Burr-Gamma (GBG) distribution is considered to
More informationOptimal Reinsurance Strategy with Bivariate Pareto Risks
University of Wisconsin Milwaukee UWM Digital Commons Theses and Dissertations May 2014 Optimal Reinsurance Strategy with Bivariate Pareto Risks Evelyn Susanne Gaus University of Wisconsin-Milwaukee Follow
More informationTail Dependence of Multivariate Pareto Distributions
!#"%$ & ' ") * +!-,#. /10 243537698:6 ;=@?A BCDBFEHGIBJEHKLB MONQP RS?UTV=XW>YZ=eda gihjlknmcoqprj stmfovuxw yy z {} ~ ƒ }ˆŠ ~Œ~Ž f ˆ ` š œžÿ~ ~Ÿ œ } ƒ œ ˆŠ~ œ
More informationTail Index Estimation of Heavy-tailed Distributions
CHAPTER 2 Tail Index Estimation of Heavy-tailed Distributions 2.1 Introduction In many diverse fields such as meteriology, finance, hydrology, climatology, environmental sciences, telecommunication, insurance
More informationON THE TRUNCATED COMPOSITE WEIBULL-PARETO MODEL
ON THE TRUNCATED COMPOSITE WEIBULL-PARETO MODEL SANDRA TEODORESCU and EUGENIA PANAITESCU The composite Weibull-Pareto model 2] was introduced as an alternative to the composite Lognormal-Pareto ] used
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 informationPossible numbers of ones in 0 1 matrices with a given rank
Linear and Multilinear Algebra, Vol, No, 00, Possible numbers of ones in 0 1 matrices with a given rank QI HU, YAQIN LI and XINGZHI ZHAN* Department of Mathematics, East China Normal University, Shanghai
More informationStat 5101 Lecture Notes
Stat 5101 Lecture Notes Charles J. Geyer Copyright 1998, 1999, 2000, 2001 by Charles J. Geyer May 7, 2001 ii Stat 5101 (Geyer) Course Notes Contents 1 Random Variables and Change of Variables 1 1.1 Random
More informationSAMPLE SIZE AND OPTIMAL DESIGNS IN STRATIFIED COMPARATIVE TRIALS TO ESTABLISH THE EQUIVALENCE OF TREATMENT EFFECTS AMONG TWO ETHNIC GROUPS
MARCEL DEKKER, INC. 70 MADISON AVENUE NEW YORK, NY 006 JOURNAL OF BIOPHARMACEUTICAL STATISTICS Vol., No. 4, pp. 553 566, 00 SAMPLE SIZE AND OPTIMAL DESIGNS IN STRATIFIED COMPARATIVE TRIALS TO ESTABLISH
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 informationGB2 Regression with Insurance Claim Severities
GB2 Regression with Insurance Claim Severities Mitchell Wills, University of New South Wales Emiliano A. Valdez, University of New South Wales Edward W. (Jed) Frees, University of Wisconsin - Madison UNSW
More informationAsymptotic Variance Formulas, Gamma Functions, and Order Statistics
Statistical Models and Methods for Lifetime Data, Second Edition by Jerald F. Lawless Copyright 2003 John Wiley & Sons, Inc. APPENDIX Β Asymptotic Variance Formulas, Gamma Functions, and Order Statistics
More informationInverse-Type Sampling Procedure for Estimating the Number of Multinomial Classes
MARCEL DEKKER, INC. 27 MADISON AVENUE NEW YORK, NY 116 23 Marcel Dekker, Inc. All rights reserved. This material may not be used or reproduced in any form without the express written permission of Marcel
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 informationModeling Severity and Measuring Tail Risk of Norwegian Fire Claims
Modeling Severity and Measuring Tail Risk of Norwegian Fire Claims Vytaras Brazauskas 1 University of Wisconsin-Milwaukee Andreas Kleefeld 2 Brandenburg University of Technology To appear in North American
More informationA New Class of Positively Quadrant Dependent Bivariate Distributions with Pareto
International Mathematical Forum, 2, 27, no. 26, 1259-1273 A New Class of Positively Quadrant Dependent Bivariate Distributions with Pareto A. S. Al-Ruzaiza and Awad El-Gohary 1 Department of Statistics
More informationUnconditional Distributions Obtained from Conditional Specification Models with Applications in Risk Theory
Unconditional Distributions Obtained from Conditional Specification Models with Applications in Risk Theory E. Gómez-Déniz a and E. Calderín Ojeda b Abstract Bivariate distributions, specified in terms
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 informationOn Five Parameter Beta Lomax Distribution
ISSN 1684-840 Journal of Statistics Volume 0, 01. pp. 10-118 On Five Parameter Beta Lomax Distribution Muhammad Rajab 1, Muhammad Aleem, Tahir Nawaz and Muhammad Daniyal 4 Abstract Lomax (1954) developed
More informationStat260: Bayesian Modeling and Inference Lecture Date: February 10th, Jeffreys priors. exp 1 ) p 2
Stat260: Bayesian Modeling and Inference Lecture Date: February 10th, 2010 Jeffreys priors Lecturer: Michael I. Jordan Scribe: Timothy Hunter 1 Priors for the multivariate Gaussian Consider a multivariate
More informationRobust estimation of the Pareto tail index: a Monte Carlo analysis
Empir Econ (2016) 51:1 30 DOI 10.1007/s00181-015-0989-9 Robust estimation of the Pareto tail index: a Monte Carlo analysis Michal Brzezinski 1 Received: 19 November 2013 / Accepted: 4 May 2015 / Published
More informationApplying the proportional hazard premium calculation principle
Applying the proportional hazard premium calculation principle Maria de Lourdes Centeno and João Andrade e Silva CEMAPRE, ISEG, Technical University of Lisbon, Rua do Quelhas, 2, 12 781 Lisbon, Portugal
More informationResearch Article Strong Convergence Bound of the Pareto Index Estimator under Right Censoring
Hindawi Publishing Corporation Journal of Inequalities and Applications Volume 200, Article ID 20956, 8 pages doi:0.55/200/20956 Research Article Strong Convergence Bound of the Pareto Index Estimator
More informationTail Approximation of Value-at-Risk under Multivariate Regular Variation
Tail Approximation of Value-at-Risk under Multivariate Regular Variation Yannan Sun Haijun Li July 00 Abstract This paper presents a general tail approximation method for evaluating the Valueat-Risk of
More informationJoint modelling of daily maximum wind strengths through the Multivariate Burr Gamma distribution
Journal of Wind Engineering and Industrial Aerodynamics 92 (2004) 1025 1037 www.elsevier.com/locate/jweia Joint modelling of daily maximum wind strengths through the Multivariate Burr Gamma distribution
More informationSpecification Tests for Families of Discrete Distributions with Applications to Insurance Claims Data
Journal of Data Science 18(2018), 129-146 Specification Tests for Families of Discrete Distributions with Applications to Insurance Claims Data Yue Fang China Europe International Business School, Shanghai,
More informationParameters Estimation Methods for the Negative Binomial-Crack Distribution and Its Application
Original Parameters Estimation Methods for the Negative Binomial-Crack Distribution and Its Application Pornpop Saengthong 1*, Winai Bodhisuwan 2 Received: 29 March 2013 Accepted: 15 May 2013 Abstract
More informationParameter Estimation for Partially Complete Time and Type of Failure Data
Biometrical Journal 46 (004), 65 79 DOI 0.00/bimj.00004 arameter Estimation for artially Complete Time and Type of Failure Data Debasis Kundu Department of Mathematics, Indian Institute of Technology Kanpur,
More informationAsymptotic Relative Efficiency in Estimation
Asymptotic Relative Efficiency in Estimation Robert Serfling University of Texas at Dallas October 2009 Prepared for forthcoming INTERNATIONAL ENCYCLOPEDIA OF STATISTICAL SCIENCES, to be published by Springer
More informationA Comparison: Some Approximations for the Aggregate Claims Distribution
A Comparison: ome Approximations for the Aggregate Claims Distribution K Ranee Thiagarajah ABTRACT everal approximations for the distribution of aggregate claims have been proposed in the literature. In
More informationA New Generalization of the Pareto Distribution and Its Application to Insurance Data
Journal of Risk and Financial Management Article A New Generalization of the Pareto Distribution and Its Application to Insurance Data Mohamed E. Ghitany 1, *, Emilio Gómez-Déniz 2 and Saralees Nadarajah
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 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 informationASTIN Colloquium 1-4 October 2012, Mexico City
ASTIN Colloquium 1-4 October 2012, Mexico City Modelling and calibration for non-life underwriting risk: from empirical data to risk capital evaluation Gian Paolo Clemente Dept. Mathematics, Finance Mathematics
More informationIntroduction to Algorithmic Trading Strategies Lecture 10
Introduction to Algorithmic Trading Strategies Lecture 10 Risk Management Haksun Li haksun.li@numericalmethod.com www.numericalmethod.com Outline Value at Risk (VaR) Extreme Value Theory (EVT) References
More informationREFERENCES AND FURTHER STUDIES
REFERENCES AND FURTHER STUDIES by..0. on /0/. For personal use only.. Afifi, A. A., and Azen, S. P. (), Statistical Analysis A Computer Oriented Approach, Academic Press, New York.. Alvarez, A. R., Welter,
More informationModified Kolmogorov-Smirnov Test of Goodness of Fit. Catalonia-BarcelonaTECH, Spain
152/304 CoDaWork 2017 Abbadia San Salvatore (IT) Modified Kolmogorov-Smirnov Test of Goodness of Fit G.S. Monti 1, G. Mateu-Figueras 2, M. I. Ortego 3, V. Pawlowsky-Glahn 2 and J. J. Egozcue 3 1 Department
More information(ii) Given that gðxþ ¼ð4x 1Þðax 2 þ bx þ cþ, find the values of the integers a, b and c. (3 marks)
Centre Number Candidate Number For Examiner s Use Surname Other Names Candidate Signature Examiner s Initials Mathematics Unit Pure Core 4 Tuesday 15 June 2010 General Certificate of Education Advanced
More informationShooting Methods for Numerical Solution of Stochastic Boundary-Value Problems
STOCHASTIC ANALYSIS AND APPLICATIONS Vol. 22, No. 5, pp. 1295 1314, 24 Shooting Methods for Numerical Solution of Stochastic Boundary-Value Problems Armando Arciniega and Edward Allen* Department of Mathematics
More informationApproximating the Conway-Maxwell-Poisson normalizing constant
Filomat 30:4 016, 953 960 DOI 10.98/FIL1604953S Published by Faculty of Sciences and Mathematics, University of Niš, Serbia Available at: http://www.pmf.ni.ac.rs/filomat Approximating the Conway-Maxwell-Poisson
More informationEstimation of Truncated Data Samples in Operational Risk Modeling
2015 The Journal of Risk and Insurance. Vol. 83, No. 3, 613 640 2016). DOI: 10.1111/jori.12062 Estimation of Truncated Data Samples in Operational Risk Modeling Bakhodir Ergashev Konstantin Pavlikov Stan
More informationOverview of Extreme Value Theory. Dr. Sawsan Hilal space
Overview of Extreme Value Theory Dr. Sawsan Hilal space Maths Department - University of Bahrain space November 2010 Outline Part-1: Univariate Extremes Motivation Threshold Exceedances Part-2: Bivariate
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 informationROBUST TESTS BASED ON MINIMUM DENSITY POWER DIVERGENCE ESTIMATORS AND SADDLEPOINT APPROXIMATIONS
ROBUST TESTS BASED ON MINIMUM DENSITY POWER DIVERGENCE ESTIMATORS AND SADDLEPOINT APPROXIMATIONS AIDA TOMA The nonrobustness of classical tests for parametric models is a well known problem and various
More informationApproximations to the distribution of sum of independent non-identically gamma random variables
Math Sci (205) 9:205 23 DOI 0.007/s40096-05-069-2 ORIGINAL RESEARCH Approximations to the distribution of sum of independent non-identically gamma random variables H. Murakami Received: 3 May 205 / Accepted:
More informationPolynomial approximation of mutivariate aggregate claim amounts distribution
Polynomial approximation of mutivariate aggregate claim amounts distribution Applications to reinsurance P.O. Goffard Axa France - Mathematics Institute of Marseille I2M Aix-Marseille University 19 th
More informationIMPROVING THE SMALL-SAMPLE EFFICIENCY OF A ROBUST CORRELATION MATRIX: A NOTE
IMPROVING THE SMALL-SAMPLE EFFICIENCY OF A ROBUST CORRELATION MATRIX: A NOTE Eric Blankmeyer Department of Finance and Economics McCoy College of Business Administration Texas State University San Marcos
More informationABC methods for phase-type distributions with applications in insurance risk problems
ABC methods for phase-type with applications problems Concepcion Ausin, Department of Statistics, Universidad Carlos III de Madrid Joint work with: Pedro Galeano, Universidad Carlos III de Madrid Simon
More informationMathematics (JAN12MPC401) General Certificate of Education Advanced Level Examination January Unit Pure Core TOTAL
Centre Number Candidate Number For Examiner s Use Surname Other Names Candidate Signature Examiner s Initials Mathematics Unit Pure Core 4 Monday 23 January 2012 General Certificate of Education Advanced
More informationFisher information for generalised linear mixed models
Journal of Multivariate Analysis 98 2007 1412 1416 www.elsevier.com/locate/jmva Fisher information for generalised linear mixed models M.P. Wand Department of Statistics, School of Mathematics and Statistics,
More informationexclusive prepublication prepublication discount 25 FREE reprints Order now!
Dear Contributor Please take advantage of the exclusive prepublication offer to all Dekker authors: Order your article reprints now to receive a special prepublication discount and FREE reprints when you
More informationDA Freedman Notes on the MLE Fall 2003
DA Freedman Notes on the MLE Fall 2003 The object here is to provide a sketch of the theory of the MLE. Rigorous presentations can be found in the references cited below. Calculus. Let f be a smooth, scalar
More informationHeriot-Watt University
Heriot-Watt University Heriot-Watt University Research Gateway Prediction of settlement delay in critical illness insurance claims by using the generalized beta of the second kind distribution Dodd, Erengul;
More informationThe Instability of Correlations: Measurement and the Implications for Market Risk
The Instability of Correlations: Measurement and the Implications for Market Risk Prof. Massimo Guidolin 20254 Advanced Quantitative Methods for Asset Pricing and Structuring Winter/Spring 2018 Threshold
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 informationAnalysis of Gamma and Weibull Lifetime Data under a General Censoring Scheme and in the presence of Covariates
Communications in Statistics - Theory and Methods ISSN: 0361-0926 (Print) 1532-415X (Online) Journal homepage: http://www.tandfonline.com/loi/lsta20 Analysis of Gamma and Weibull Lifetime Data under a
More informationModelling Dropouts by Conditional Distribution, a Copula-Based Approach
The 8th Tartu Conference on MULTIVARIATE STATISTICS, The 6th Conference on MULTIVARIATE DISTRIBUTIONS with Fixed Marginals Modelling Dropouts by Conditional Distribution, a Copula-Based Approach Ene Käärik
More informationCHAPTER 2 A Function for Size Distribution of Incomes
CHAPTER 2 A Function for Size Distribution of Incomes S. K. Singh and G. S. Maddala Abstract The paper derives a function that describes the size distribution of incomes. The two functions most often used
More informationPANJER CLASS UNITED One formula for the probabilities of the Poisson, Binomial, and Negative Binomial distribution.
PANJER CLASS UNITED One formula for the probabilities of the Poisson, Binomial, and Negative Binomial distribution Michael Facler 1 Abstract. This paper gives a formula representing all discrete loss distributions
More informationACTEX CAS EXAM 3 STUDY GUIDE FOR MATHEMATICAL STATISTICS
ACTEX CAS EXAM 3 STUDY GUIDE FOR MATHEMATICAL STATISTICS TABLE OF CONTENTS INTRODUCTORY NOTE NOTES AND PROBLEM SETS Section 1 - Point Estimation 1 Problem Set 1 15 Section 2 - Confidence Intervals and
More informationSequences, Series, and the Binomial Formula
CHAPTER Sequences, Series, nd the Binomil Formul. SEQUENCES. ; ; ; ; 6 ; 6 6. ðþ ; ðþ ; ð Þ 6; ðþ ; ðþ 6; 6 ð6þ. ðþ ; ðþ : ðþ ; ðþ ; ðþ ; 6 ðþ 6 6 6. ; ; ; ; ; 6 6 6. 0 ; ; ; 8 ; 6 8 ; 6. 0; ; 6 ; ; 6
More informationBayesian Sensitivity Analysis for VaR and TVaR
Bayesian Sensitivity Analysis for VaR and TVaR Edgar Anguiano ceganguiano@ciencias.unam.mx ABSTRACT It is known that VaR and TVaR are useful to measure huge quantities at low probabilities i.e. measures
More informationSome Theoretical Properties and Parameter Estimation for the Two-Sided Length Biased Inverse Gaussian Distribution
Journal of Probability and Statistical Science 14(), 11-4, Aug 016 Some Theoretical Properties and Parameter Estimation for the Two-Sided Length Biased Inverse Gaussian Distribution Teerawat Simmachan
More informationBurr Type X Distribution: Revisited
Burr Type X Distribution: Revisited Mohammad Z. Raqab 1 Debasis Kundu Abstract In this paper, we consider the two-parameter Burr-Type X distribution. We observe several interesting properties of this distribution.
More informationStochastic Comparisons of Weighted Sums of Arrangement Increasing Random Variables
Portland State University PDXScholar Mathematics and Statistics Faculty Publications and Presentations Fariborz Maseeh Department of Mathematics and Statistics 4-7-2015 Stochastic Comparisons of Weighted
More informationEstimation of the Bivariate Generalized. Lomax Distribution Parameters. Based on Censored Samples
Int. J. Contemp. Math. Sciences, Vol. 9, 2014, no. 6, 257-267 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ijcms.2014.4329 Estimation of the Bivariate Generalized Lomax Distribution Parameters
More informationBayesian Inference for the Multivariate Normal
Bayesian Inference for the Multivariate Normal Will Penny Wellcome Trust Centre for Neuroimaging, University College, London WC1N 3BG, UK. November 28, 2014 Abstract Bayesian inference for the multivariate
More informationRisk Aggregation with Dependence Uncertainty
Introduction Extreme Scenarios Asymptotic Behavior Challenges Risk Aggregation with Dependence Uncertainty Department of Statistics and Actuarial Science University of Waterloo, Canada Seminar at ETH Zurich
More informationSTA2603/205/1/2014 /2014. ry II. Tutorial letter 205/1/
STA263/25//24 Tutorial letter 25// /24 Distribution Theor ry II STA263 Semester Department of Statistics CONTENTS: Examination preparation tutorial letterr Solutions to Assignment 6 2 Dear Student, This
More informationA Hypothesis-Free Multiple Scan Statistic with Variable Window
Biometrical Journal 50 (2008) 2, 299 310 DOI: 10.1002/bimj.200710412 299 A Hypothesis-Free Multiple Scan Statistic with Variable Window L. Cucala * Institut de Math matiques de Toulouse, Universit Paul
More informationIntroduction to Probability and Statistics (Continued)
Introduction to Probability and Statistics (Continued) Prof. icholas Zabaras Center for Informatics and Computational Science https://cics.nd.edu/ University of otre Dame otre Dame, Indiana, USA Email:
More informationFisher Information & Efficiency
Fisher Information & Efficiency Robert L. Wolpert Department of Statistical Science Duke University, Durham, NC, USA 1 Introduction Let f(x θ) be the pdf of X for θ Θ; at times we will also consider a
More informationOn the existence of maximum likelihood estimators in a Poissongamma HGLM and a negative binomial regression model
Int. Statistical Inst.: Proc. 58th World Statistical Congress, 2011, Dublin Session CPS027) p.6433 On the existence of maximum likelihood estimators in a Poissongamma HGLM and a negative binomial regression
More informationON STATISTICAL INFERENCE UNDER ASYMMETRIC LOSS. Abstract. We introduce a wide class of asymmetric loss functions and show how to obtain
ON STATISTICAL INFERENCE UNDER ASYMMETRIC LOSS FUNCTIONS Michael Baron Received: Abstract We introduce a wide class of asymmetric loss functions and show how to obtain asymmetric-type optimal decision
More informationRuin Theory. A thesis submitted in partial fulfillment of the requirements for the degree of Master of Science at George Mason University
Ruin Theory A thesis submitted in partial fulfillment of the requirements for the degree of Master of Science at George Mason University by Ashley Fehr Bachelor of Science West Virginia University, Spring
More informationEstimation for Mean and Standard Deviation of Normal Distribution under Type II Censoring
Communications for Statistical Applications and Methods 2014, Vol. 21, No. 6, 529 538 DOI: http://dx.doi.org/10.5351/csam.2014.21.6.529 Print ISSN 2287-7843 / Online ISSN 2383-4757 Estimation for Mean
More informationSOME INEQUALITIES FOR THE q-digamma FUNCTION
Volume 10 (009), Issue 1, Article 1, 8 pp SOME INEQUALITIES FOR THE -DIGAMMA FUNCTION TOUFIK MANSOUR AND ARMEND SH SHABANI DEPARTMENT OF MATHEMATICS UNIVERSITY OF HAIFA 31905 HAIFA, ISRAEL toufik@mathhaifaacil
More informationCHAIN LADDER FORECAST EFFICIENCY
CHAIN LADDER FORECAST EFFICIENCY Greg Taylor Taylor Fry Consulting Actuaries Level 8, 30 Clarence Street Sydney NSW 2000 Australia Professorial Associate, Centre for Actuarial Studies Faculty of Economics
More informationNotes on a skew-symmetric inverse double Weibull distribution
Journal of the Korean Data & Information Science Society 2009, 20(2), 459 465 한국데이터정보과학회지 Notes on a skew-symmetric inverse double Weibull distribution Jungsoo Woo 1 Department of Statistics, Yeungnam
More informationSPRING 2007 EXAM C SOLUTIONS
SPRING 007 EXAM C SOLUTIONS Question #1 The data are already shifted (have had the policy limit and the deductible of 50 applied). The two 350 payments are censored. Thus the likelihood function is L =
More informationExact Inference for the Two-Parameter Exponential Distribution Under Type-II Hybrid Censoring
Exact Inference for the Two-Parameter Exponential Distribution Under Type-II Hybrid Censoring A. Ganguly, S. Mitra, D. Samanta, D. Kundu,2 Abstract Epstein [9] introduced the Type-I hybrid censoring scheme
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 a Normal Covariance Matrix for Small Samples with Monotone Missing Data
Applied Mathematical Sciences, Vol 3, 009, no 54, 695-70 Testing a Normal Covariance Matrix for Small Samples with Monotone Missing Data Evelina Veleva Rousse University A Kanchev Department of Numerical
More informationTHE RUNNING INTERVAL SMOOTHER: A CONFIDENCE BAND HAVING SOME SPECIFIED SIMULTANEOUS PROBABILITY COVERAGE
International Journal of Statistics: Advances in Theory and Applications Vol. 1, Issue 1, 2017, Pages 21-43 Published Online on April 12, 2017 2017 Jyoti Academic Press http://jyotiacademicpress.org THE
More informationIntroduction to Rare Event Simulation
Introduction to Rare Event Simulation Brown University: Summer School on Rare Event Simulation Jose Blanchet Columbia University. Department of Statistics, Department of IEOR. Blanchet (Columbia) 1 / 31
More informationStatistical Mechanics of Money, Income, and Wealth
Statistical Mechanics of Money, Income, and Wealth Victor M. Yakovenko Adrian A. Dragulescu and A. Christian Silva Department of Physics, University of Maryland, College Park, USA http://www2.physics.umd.edu/~yakovenk/econophysics.html
More informationClosed-Form Estimators for the Gamma Distribution Derived from Likelihood Equations
The American Statistician ISSN: 3-135 (Print) 1537-2731 (Online) Journal homepage: http://www.tandfonline.com/loi/utas2 Closed-Form Estimators for the Gamma Distribution Derived from Likelihood Equations
More informationAn algorithm for robust fitting of autoregressive models Dimitris N. Politis
An algorithm for robust fitting of autoregressive models Dimitris N. Politis Abstract: An algorithm for robust fitting of AR models is given, based on a linear regression idea. The new method appears to
More informationINFORMATION THEORY AND STATISTICS
INFORMATION THEORY AND STATISTICS Solomon Kullback DOVER PUBLICATIONS, INC. Mineola, New York Contents 1 DEFINITION OF INFORMATION 1 Introduction 1 2 Definition 3 3 Divergence 6 4 Examples 7 5 Problems...''.
More informationEstimating VaR in credit risk: Aggregate vs single loss distribution
Estimating VaR in credit risk: Aggregate vs single loss distribution M. Assadsolimani and D. Chetalova arxiv:172.4388v1 [q-fin.cp] 14 Feb 217 Abstract Using Monte Carlo simulation to calculate the Value
More informationExercises Chapter 4 Statistical Hypothesis Testing
Exercises Chapter 4 Statistical Hypothesis Testing Advanced Econometrics - HEC Lausanne Christophe Hurlin University of Orléans December 5, 013 Christophe Hurlin (University of Orléans) Advanced Econometrics
More information