Can we use kernel smoothing to estimate Value at Risk and Tail Value at Risk?
|
|
- Darcy Webster
- 6 years ago
- Views:
Transcription
1 Can we use kernel smoothing to estimate Value at Risk and Tail Value at Risk? Ramon Alemany, Catalina Bolancé and Montserrat Guillén Riskcenter - IREA Universitat de Barcelona ASTIN Colloquium The Hague 22nd-24th May 2013 Alemany, Bolancé & Guillén (Riskcenter UB) ASTIN rd May / 35
2 1 Motivation 2 Estimation alternatives: Emp, CKE and TKE 3 Double transformed kernel estimation 4 Simulation study 5 Data study 6 Conclusions Alemany, Bolancé & Guillén (Riskcenter UB) ASTIN rd May / 35
3 Motivation Analyze the properties of Transformed kernel estimations of Valueat-Risk (VaR) and Tail-Value-at-Risk (TVaR). For VaR estimation the Classical kernel estimator (CKE) and the Transformed kernel estimator (TKE) have been shown to be better than the empirical approach (Chen et al., 2005, and Alemany et al., 2013, respectively). Double transformed kernel density estimation (DTKE) 1 of VaR has a lower mean square error (MSE) than the Classical kernel estimation. Double transformed kernel density estimation can be applied to very extreme risk measurement (in Alemany et al., 2013, the authors reach the confidence level) when we have a loss distribution with heavy tail and is specially suitable when the sample size is large. We assess if this result holds when we estimate TVaR. 1 A particular case of transformed kernel estimation using double transformation Alemany, Bolancé & Guillén (Riskcenter UB) ASTIN rd May / 35
4 Motivation The aim is to estimate: and VaR α (X ) = inf {x, F X (x) α} (1) TVaR α (X ) = E (X X VaR α (X )). (2) where X is a random variable in R + representing a loss with probability distribution function (pdf) f X, and cumulative distribution function (cdf) F X. We compare empirical approach (Emp), classical kernel estimation (CKE) and double transformed kernel estimation (DTKE) using a simulation study and an application to the estimation of the risk of loss in an auto insurance portfolio. Alemany, Bolancé & Guillén (Riskcenter UB) ASTIN rd May / 35
5 Motivation Both, VaR and TVaR, are based on the cdf, F X, associated to the loss variable. Therefore, the theoretical properties of the nonparametric estimations of VaR and TVaR are strongly related to the theoretical properties of the nonparametric estimation of cdf. Alemany, Bolancé & Guillén (Riskcenter UB) ASTIN rd May / 35
6 Estimation alternatives: Emp, CKE and DTKE Empirical approach of cdf Let us assume that X i, i = 1,..., n denotes data observations from the loss random variable X. A natural nonparametric method to estimate cdf is the empirical distribution, F n (x) = 1 n n I (X i x), (3) i=1 where I ( ) is the indicator function, so I ( ) = 1 if condition between parentheses is true. Empirical distribution is an unbiased estimation of cdf F X but it has larger variance than alternative nonparametric approaches. Alemany, Bolancé & Guillén (Riskcenter UB) ASTIN rd May / 35
7 Estimation alternatives: Emp, CKE and TKE Empirical approach of VaR and TVaR Then, the empirical estimate of value-at-risk is: { VaR α (X ) n = inf x, F } n (x) α, (4) and the empirical estimator of tail value-at-risk is: n i=1 TVaR X ii (X i > VaR ) α (X ) n α (X ) n = n i=1 (X I i > VaR ), (5) α (X ) n where VaR α (X ) n is defined in (4). Alemany, Bolancé & Guillén (Riskcenter UB) ASTIN rd May / 35
8 Estimation alternatives: Emp, CKE and TKE Empirical approach of VaR and TVaR Risk estimation from the empirical distribution is very simple, but it cannot extrapolate beyond the maximum observed data point. This is especially troublesome if the sample is not too large, and one may suspect that the probability of a loss larger than the maximum observed loss in the data sample is not zero. This can cause the empirical estimations of VaR and TVaR to underestimate the risk. Alemany, Bolancé & Guillén (Riskcenter UB) ASTIN rd May / 35
9 Estimation alternatives: Emp, CKE and TKE Classical kernel estimation of cdf Classical kernel estimation (CKE) of a cdf, F X, is obtained integrating the classical kernel estimation of its pdf, f X. The usual expression for the kernel estimator of a cdf is: F X (x) = x f X (u)du = 1 n ( ) n i=1 K x Xi b, (6) where K( ) is the cdf of k( ) which is known as the kernel function (usually a symmetric pdf). Some examples of very common kernel functions are the Epanechnikov and the Gaussian kernel (Silverman, 1986). Parameter b is the bandwidth or the smoothing parameter. It controls the smoothness of the estimated cdf. The larger b is, the smoother the resulting cdf. Alemany, Bolancé & Guillén (Riskcenter UB) ASTIN rd May / 35
10 Estimation alternatives: Emp, CKE and TKE Classical kernel estimation of cdf The Classical kernel cdf estimator has less variance than the empirical distribution estimator, but it has some bias which tends to zero if the sample size is large. Azzalini (1981) showed that the MSE of the CKE is less than the MSE of the empirical estimation. Alemany, Bolancé & Guillén (Riskcenter UB) ASTIN rd May / 35
11 Estimation alternatives: Emp, CKE and TKE Classical kernel estimation of VaR and TVaR To estimate VaR α from CKE F X ( ), we use the Newton-Raphson method to solve the equation: F X (x) = α. (7) The kernel estimator of TVaR α is based on the kernel estimation of the survival function, ŜX (x) = 1 F X (x), so: 1 TVaR α (X ) = n(1 α) n ( VaRα (X ) X i X i (1 K b i=1 )), (8) Alemany, Bolancé & Guillén (Riskcenter UB) ASTIN rd May / 35
12 Estimation alternatives: Emp, CKE and TKE Classical kernel estimation of VaR and TVaR: Theoretical properties The properties of CKE of cdf are transferred to the CKE of VaR and TVaR, i.e. it has lower MSE than empirical approach, Emp. However, Chen (2008) showed that in practice these better properties cannot be transferred to the CKE of TVaR. DIFFICULTY: In expression 1 TVaR α (X ) = n(1 α) n ( VaRα (X ) X i X i (1 K b i=1 )), (9) we need to replace theoretical VaR α (X ) by its estimation, then CKE of TVaR has larger MSE than empirical approach. Alemany, Bolancé & Guillén (Riskcenter UB) ASTIN rd May / 35
13 Estimation alternatives: Emp, CKE and TKE Transformed kernel estimation of cdf Let T ( ) be a concave transformation where Y = T (X ) and Y i = T (X i ), i = 1... n are the transformed observed losses. Then the classical kernel estimator of the transformed variable is: F T (X ) (T (x)) = 1 n n ( ) T (x) T (Xi ) K. (10) b i=1 The transformed kernel estimation (TKE) of F X (x) is: F X (x) = F T (X ) (T (x)). Alemany, Bolancé & Guillén (Riskcenter UB) ASTIN rd May / 35
14 Estimation alternatives: Emp, CKE and TKE Transformed kernel estimation of VaR and TVaR In order to estimate VaR α the Newton-Raphson method solves the equation F T (X ) (T (x)) = α and, once the result is obtained, the inverse of the transformation is applied. The transformed kernel estimation of TVaR α is based on the transformed kernel estimator of survival function, i.e. Ŝ X (x) = ŜT (X ) (T (x)) = 1 F T (X ) (T (x)), so, TVaR α (X ) = 1 n(1 α) n ( )) T (VaRα (X )) T (X i ) X i (1 K, b i=1 (11) Alemany, Bolancé & Guillén (Riskcenter UB) ASTIN rd May / 35
15 Estimation alternatives: Emp, CKE and TKE Transformed kernel estimation of VaR and TVaR: Theoretical properties Our last paper (Alemany et al., 2013) showed that the use of a suitable transformation can reduce the variance at the expense of increasing the bias of the estimation, resulting in a lower MSE. These properties are transferred to the TKE of the VaR but, as occurs with CKE, are not transferred to the TVaR estimation, given that we need to replace in TVaR α (X ) the true VaR by the estimated VaR. Can the double transformed kernel estimation (DTKE) of the authors be a good alternative? ANSWER: In some specific cases. Alemany, Bolancé & Guillén (Riskcenter UB) ASTIN rd May / 35
16 Double transformed kernel estimation As is well known, to obtain a smoothing parameter that is asymptotically optimal, it is sufficient to minimize the A-MISE (Asymptotic mean integrated squared error): ( 1 [f 4 b4 Y (y)] 2 dy ) 2 t 2 k (t) dt 1 n b K (t) [1 K (t)] dt. Given b and the kernel k( ), the A-MISE is minimum when functional [f Y (y)]2 dy is minimum. Therefore, the proposed method is based on the transformation of the variable in order to achieve a distribution that minimizes the A-MISE, i.e. that minimizes [f Y (y)]2 dy. Terrell (1990) showed that the density of a Beta (3, 3) distribution defined on the domain [ 1, 1] minimizes [f Y (y)]2 dy, in the set of all densities with known variance. Alemany, Bolancé & Guillén (Riskcenter UB) ASTIN rd May / 35
17 Double transformed kernel estimation Let M( ) be the cdf of the cited Beta (3, 3), the DTKE method requires an initial transformation of the data 2, T (X i ) = Z i, where we get a transformed variable distribution that is close to a Uniform (0, 1) (Buch- Larsen et al, 2005). Afterwards, the data are transformed again using the inverse of the distribution function of a Beta (3, 3), M 1 (Z i ) = Y i. The DTKE is: F M 1 (T (x)) (x) = 1 n n ( M 1 (T (x)) M 1 ) (T (X i )) K.(12) b i=1 2 With a generalized Champernowne cdf: TX (x) = ((x + c) γ c γ ) / ((x + c) γ + (M + c) γ 2c γ ) Alemany, Bolancé & Guillén (Riskcenter UB) ASTIN rd May / 35
18 Double transformed kernel estimation Double transformed kernel estimation of VaR and TVaR Double transformed kernel estimation of VaR is obtained similarly to transformed kernel estimation and has an important advantage: it allows us to calculate the optimal bandwidth (b) that minimizes MSE, given that we know that transformed data have Beta (3, 3) distribution: b = ( m (q) 1 ) 1 1 K 2 (t) dt [ 1 2 m (q) 1 1 t2 k (t) dt ] n 1 3, (13) where m is the pdf of Beta (3, 3), m its first derivative and q is the VaR α of the Beta distribution. Alemany, Bolancé & Guillén (Riskcenter UB) ASTIN rd May / 35
19 Double transformed kernel estimation Double transformed kernel estimation of VaR and TVaR The double transformed kernel estimation of TVaR α is based on the double transformed kernel estimator of the survival function, i.e. so, Ŝ X (x) = ŜM 1 (T (X )) ( M 1 (T (x)) ) = 1 F M 1 (T (X )) ( M 1 (T (x)) ), TVaR α (X ) = 1 n(1 α) n ( M X i (1 1 (T (VaR α (X ))) M 1 )) (T (X i )) K. b i=1 (14) Alemany, Bolancé & Guillén (Riskcenter UB) ASTIN rd May / 35
20 Double transformed kernel estimation Double transformed kernel estimation of VaR and TVaR ( We know that q = M (T 1 VaRα (X ))) is a quantile of Beta (3, 3) distribution. Let p = M(q) be the confidence level of q. We propose to estimate TVaR based on the theoretical quantile of Beta (3, 3), so then the DTKE of TVaR is: TVaR p (X ) = 1 n(1 p) The optimal bandwidth is: b = n ( q M X i (1 1 )) (T (X i )) K. (15) b i=1 ( qm (q) 1 ) 1 1 K 2 (t) dt [ (1 p) 1 (qm (q) + m (q)) 1 1 t2 k (t) dt ] n 1 3, (16) Alemany, Bolancé & Guillén (Riskcenter UB) ASTIN rd May / 35
21 Double transformed kernel estimation Double transformed kernel estimation of VaR and TVaR: Theoretical properties In Alemany et al., 2013, we showed that DTKE of VaR is better than Emp and CKE, specially when we have large sample size and the loss distribution has heavy tail. Now we show that DTKE of TVaR increases bias and reduces the variance of the empirical estimate. The difficulty is that, in some cases, reducing the variance does not compensate the increase in the bias, resulting in greater MSE of DTKE if we compare with the empirical approach. Alemany, Bolancé & Guillén (Riskcenter UB) ASTIN rd May / 35
22 Simulation study A wide simulation study is carried out. We obtain 2000 samples of size 500 and 5000 from different distribution functions with different parameters and different tails shapes: Lognormal, Weibull and mixtures Lognormal-Pareto. To ensure that mixtures Lognormal-Pareto have finite first moment we use parameter of Pareto less than unity. Results show that our DTKE method overestimates TVaR in Lognormal and Weibull distribution. We obtain a larger positive bias which is not compensated by variance reduction. In these cases MSE of DTKE is larger than MSE of the empirical approach. We obtain very good results in mixtures Lognormal-Pareto with large sample size. In these cases MSE of DTKE is much smaller than the MSE of empirical approach when confidence level is 0.95 and Alemany, Bolancé & Guillén (Riskcenter UB) ASTIN rd May / 35
23 Simulation study Results of MSE of TVaR α for a mixture of Lognormal-Pareto n=500 n= % Lognormal-30% Pareto Method α = 0.95 α = α = 0.95 α = Emp DTKEx % Lognormal - 70% Pareto α = 0.95 α = α = 0.95 α = Emp DTKEx Alemany, Bolancé & Guillén (Riskcenter UB) ASTIN rd May / 35
24 Data study We analyze a data set containing a sample with 5,122 automobile claim costs from a Spanish insurance company for traffic accidents in This is an standard insurance data set with observations on the cost of accident claims, i.e. a large sample that looks heavy-tailed with lots of small values and a few large extremes. The original data are divided into two groups: claims from policyholders who were less than 30 years old when the accident took place, and claims from policyholders who were 30 years old or older when they had the accident that originated the claim for compensation. The first group consists of 1,061 observations in the interval of cost of claim from 1 to 126,000 and the second group contains 4,061 observations in the interval from 1 to 17,000. Costs are expressed in monetary units. Alemany, Bolancé & Guillén (Riskcenter UB) ASTIN rd May / 35
25 Data study Summary of the younger and older policyholders claims cost data Data n Mean Median Std. Deviation Coeff. Variation All 5, , Younger 1, , Older 4, Cost of claims are expressed in monetary units. Alemany, Bolancé & Guillén (Riskcenter UB) ASTIN rd May / 35
26 Data study VaR results VaR α results for automobile claim cost data. α=0.95 α=0.995 Method Younger Older All Younger Older All Emp CKE DTKEx Alemany, Bolancé & Guillén (Riskcenter UB) ASTIN rd May / 35
27 Data study VaR results Estimated VaR α for tolerance levels (x-axis) above 99%. Comparison of three methods for All insured. Solid, dashed and dotted lines correspond to the empirical, the classical kernel and the transformed kernel estimation method, respectively. Alemany, Bolancé & Guillén (Riskcenter UB) ASTIN rd May / 35
28 Data study VaR results Estimated VaR α for tolerance levels (x-axis) between 99% and 99.9%. Comparison of three methods for All insured. Solid, dashed and dotted lines correspond to the empirical, the classical kernel and the transformed kernel estimation method, respectively. Alemany, Bolancé & Guillén (Riskcenter UB) ASTIN rd May / 35
29 Data study VaR results Estimated VaR α for tolerance levels (x-axis) above 99%. Value-at-Risk estimated with double transformed kernel estimation given the tolerance level. Solid line and dashed line correspond to older and younger policyholders, respectively. Alemany, Bolancé & Guillén (Riskcenter UB) ASTIN rd May / 35
30 Data study TVaR results TVaR α results for automobile claim cost data. α=0.95 α=0.995 Method Younger Older All Younger Older All Emp CKE DTKEx Alemany, Bolancé & Guillén (Riskcenter UB) ASTIN rd May / 35
31 Data study TVaR results Estimated TVaR α for tolerance levels (x-axis) above 99%. Comparison of two methods for All insured. Solid and dashed lines correspond to the empirical and the transformed kernel estimation method, respectively. Alemany, Bolancé & Guillén (Riskcenter UB) ASTIN rd May / 35
32 Data study TVaR results Estimated TVaR α for tolerance levels (x-axis) above 99%. Tail Value-at-Risk estimated with double transformed kernel estimation given the tolerance level. Solid line and dashed line correspond to older and younger policyholders, respectively. Alemany, Bolancé & Guillén (Riskcenter UB) ASTIN rd May / 35
33 Conclusions Double transformed kernel estimation of TVaR improves empirical approach when we have heavy tailed loss distribution. DTKE allows us extrapolate quantiles of the distribution above the maximum observed with relative efficiency. We believe our approach to be a practical method for risk analysts and pricing departments. Our method can establish a distance between risk classes in terms of differences in the risk of extreme severities. Additionally, a surcharge to the a priori premium can be linked to the loss distribution of severities with the advantage of needing no distributional assumptions and of being easy to implement. Alemany, Bolancé & Guillén (Riskcenter UB) ASTIN rd May / 35
34 References Alemany, R., Bolancé, C. and Guillén, M., 2013, A nonparametric approach to calculating value-at-risk, Insurance: Mathematics and Economics, 52, Azzalini, A., 1981, A note on the estimation of a distribution function and quantiles by a kernel method, Biometrika, 68, Bolancé, C., 2010, Optimal inverse Beta(3,3) transformation in kernel density estimation, SORT-Statistics and Operations Research Transactions, 34, Bolancé, C., Guillén, M. and Nielsen, J.P., 2003, Kernel density estimation of actuarial loss functions, Insurance: Mathematics and Economics, 32, Bolancé, C., Guillén, M. and Nielsen, J.P., 2008, Inverse Beta transformation in kernel density estimation, Statistics& Probability Letters, 78, Bolancé, C., Guillén, M., Gustafsson, J. and Nielsen, J.P., 2012, Quantitative Operational Risk Models, Chapman & Hall/CRC Finance Series, London. Bowman, A.a and Hall, P. and Prvan, T., 1998, Bandwidth selection for smoothing of distribution function, Biometrika, 85, Buch-Larsen, T., Guillén, M., Nielsen, J.P. and Bolancé, C., 2005, Kernel density estimation for heavy-tailed distributions using the Champernowne transformation, Statistics, 39, Cai, Z. and Wang, X., 2008, Nonparametric estimation of conditional VaR and expected shortfall, Journal of Econometrics, 147, Chen, S.X., 2008, Nonparametric Estimation of Expected Shortfall, Journal of Financial Econometrics, 6, Chen, S.X. and Tang, C.Y., 2005, Nonparametric inference of Value at Risk for dependent financial returns, Journal of Financial Econometrics, 3, Silverman, B.W., 1986, Density Estimation for Statistics and Data Analysis, Chapman & Hall/CRC Finance Series, London. Terrell, G.R., 1990, The maximal smoothing principle in density estimation, Journal of the American Statistical Association, 85, Alemany, Bolancé & Guillén (Riskcenter UB) ASTIN rd May / 35
35 Can we use kernel smoothing to estimate Value at Risk and Tail Value at Risk? Ramon Alemany, Catalina Bolancé and Montserrat Guillén Riskcenter - IREA Universitat de Barcelona ASTIN Colloquium The Hague 22nd-24th May 2013 Alemany, Bolancé & Guillén (Riskcenter UB) ASTIN rd May / 35
Nonparametric approach to analysing operational risk losses
Nonparametric approach to analysing operational risk losses Catalina Bolancé, Mercedes Ayuso and Montserrat Guillén Departament d Econometria, Estadística i Economia Espanyola. RFA-IREA University of Barcelona.
More information**BEGINNING OF EXAMINATION** A random sample of five observations from a population is:
**BEGINNING OF EXAMINATION** 1. You are given: (i) A random sample of five observations from a population is: 0.2 0.7 0.9 1.1 1.3 (ii) You use the Kolmogorov-Smirnov test for testing the null hypothesis,
More informationQuantitative modeling of operational risk losses when combining internal and external data sources
Quantitative modeling of operational risk losses when combining internal and external data sources Jens Perch Nielsen (Cass Business School, City University, United Kingdom) Montserrat Guillén (Riskcenter,
More informationFitting financial time series returns distributions: a mixture normality approach
Fitting financial time series returns distributions: a mixture normality approach Riccardo Bramante and Diego Zappa * Abstract Value at Risk has emerged as a useful tool to risk management. A relevant
More informationHomework Problems Stat 479
Chapter 10 91. * A random sample, X1, X2,, Xn, is drawn from a distribution with a mean of 2/3 and a variance of 1/18. ˆ = (X1 + X2 + + Xn)/(n-1) is the estimator of the distribution mean θ. Find MSE(
More informationReducing the mean squared error of quantile-based estimators by smoothing
Reducing the mean squared error of quantile-based estimators by smoothing Mia Hubert, Irène Gijbels and Dina Vanpaemel May 2, 202 Abstract Many univariate robust estimators are based on quantiles. As already
More informationPractice Exam 1. Loss Amount Number of Losses
Practice Exam 1 1. You are given the following data on loss sizes: An ogive is used as a model for loss sizes. Determine the fitted median. Loss Amount Number of Losses 0 1000 5 1000 5000 4 5000 10000
More informationHomework Problems Stat 479
Chapter 2 1. Model 1 is a uniform distribution from 0 to 100. Determine the table entries for a generalized uniform distribution covering the range from a to b where a < b. 2. Let X be a discrete random
More informationLecture 6: Non Normal Distributions
Lecture 6: Non Normal Distributions and their Uses in GARCH Modelling Prof. Massimo Guidolin 20192 Financial Econometrics Spring 2015 Overview Non-normalities in (standardized) residuals from asset return
More informationRandom Variables and Probability Distributions
Chapter 3 Random Variables and Probability Distributions Chapter Three Random Variables and Probability Distributions 3. Introduction An event is defined as the possible outcome of an experiment. In engineering
More informationProbability and Statistics
Kristel Van Steen, PhD 2 Montefiore Institute - Systems and Modeling GIGA - Bioinformatics ULg kristel.vansteen@ulg.ac.be CHAPTER 3: PARAMETRIC FAMILIES OF UNIVARIATE DISTRIBUTIONS 1 Why do we need distributions?
More informationSYLLABUS OF BASIC EDUCATION SPRING 2018 Construction and Evaluation of Actuarial Models Exam 4
The syllabus for this exam is defined in the form of learning objectives that set forth, usually in broad terms, what the candidate should be able to do in actual practice. Please check the Syllabus Updates
More informationA New Hybrid Estimation Method for the Generalized Pareto Distribution
A New Hybrid Estimation Method for the Generalized Pareto Distribution Chunlin Wang Department of Mathematics and Statistics University of Calgary May 18, 2011 A New Hybrid Estimation Method for the GPD
More informationMarket Risk Analysis Volume I
Market Risk Analysis Volume I Quantitative Methods in Finance Carol Alexander John Wiley & Sons, Ltd List of Figures List of Tables List of Examples Foreword Preface to Volume I xiii xvi xvii xix xxiii
More informationCase Study: Heavy-Tailed Distribution and Reinsurance Rate-making
Case Study: Heavy-Tailed Distribution and Reinsurance Rate-making May 30, 2016 The purpose of this case study is to give a brief introduction to a heavy-tailed distribution and its distinct behaviors in
More informationAnalysis of truncated data with application to the operational risk estimation
Analysis of truncated data with application to the operational risk estimation Petr Volf 1 Abstract. Researchers interested in the estimation of operational risk often face problems arising from the structure
More informationFrequency Distribution Models 1- Probability Density Function (PDF)
Models 1- Probability Density Function (PDF) What is a PDF model? A mathematical equation that describes the frequency curve or probability distribution of a data set. Why modeling? It represents and summarizes
More informationAsset Allocation Model with Tail Risk Parity
Proceedings of the Asia Pacific Industrial Engineering & Management Systems Conference 2017 Asset Allocation Model with Tail Risk Parity Hirotaka Kato Graduate School of Science and Technology Keio University,
More informationDistortion operator of uncertainty claim pricing using weibull distortion operator
ISSN: 2455-216X Impact Factor: RJIF 5.12 www.allnationaljournal.com Volume 4; Issue 3; September 2018; Page No. 25-30 Distortion operator of uncertainty claim pricing using weibull distortion operator
More informationREINSURANCE RATE-MAKING WITH PARAMETRIC AND NON-PARAMETRIC MODELS
REINSURANCE RATE-MAKING WITH PARAMETRIC AND NON-PARAMETRIC MODELS By Siqi Chen, Madeleine Min Jing Leong, Yuan Yuan University of Illinois at Urbana-Champaign 1. Introduction Reinsurance contract is an
More informationAn Improved Skewness Measure
An Improved Skewness Measure Richard A. Groeneveld Professor Emeritus, Department of Statistics Iowa State University ragroeneveld@valley.net Glen Meeden School of Statistics University of Minnesota Minneapolis,
More informationA Markov Chain Monte Carlo Approach to Estimate the Risks of Extremely Large Insurance Claims
International Journal of Business and Economics, 007, Vol. 6, No. 3, 5-36 A Markov Chain Monte Carlo Approach to Estimate the Risks of Extremely Large Insurance Claims Wan-Kai Pang * Department of Applied
More informationFinancial Risk Forecasting Chapter 9 Extreme Value Theory
Financial Risk Forecasting Chapter 9 Extreme Value Theory Jon Danielsson 2017 London School of Economics To accompany Financial Risk Forecasting www.financialriskforecasting.com Published by Wiley 2011
More informationHomework Problems Stat 479
Chapter 2 1. Model 1 in the table handed out in class is a uniform distribution from 0 to 100. Determine what the table entries would be for a generalized uniform distribution covering the range from a
More informationPoint Estimation. Some General Concepts of Point Estimation. Example. Estimator quality
Point Estimation Some General Concepts of Point Estimation Statistical inference = conclusions about parameters Parameters == population characteristics A point estimate of a parameter is a value (based
More informationIntroduction Models for claim numbers and claim sizes
Table of Preface page xiii 1 Introduction 1 1.1 The aim of this book 1 1.2 Notation and prerequisites 2 1.2.1 Probability 2 1.2.2 Statistics 9 1.2.3 Simulation 9 1.2.4 The statistical software package
More informationFinancial Time Series and Their Characteristics
Financial Time Series and Their Characteristics Egon Zakrajšek Division of Monetary Affairs Federal Reserve Board Summer School in Financial Mathematics Faculty of Mathematics & Physics University of Ljubljana
More informationA Comparison Between Skew-logistic and Skew-normal Distributions
MATEMATIKA, 2015, Volume 31, Number 1, 15 24 c UTM Centre for Industrial and Applied Mathematics A Comparison Between Skew-logistic and Skew-normal Distributions 1 Ramin Kazemi and 2 Monireh Noorizadeh
More informationCambridge University Press Risk Modelling in General Insurance: From Principles to Practice Roger J. Gray and Susan M.
adjustment coefficient, 272 and Cramér Lundberg approximation, 302 existence, 279 and Lundberg s inequality, 272 numerical methods for, 303 properties, 272 and reinsurance (case study), 348 statistical
More information1. You are given the following information about a stationary AR(2) model:
Fall 2003 Society of Actuaries **BEGINNING OF EXAMINATION** 1. You are given the following information about a stationary AR(2) model: (i) ρ 1 = 05. (ii) ρ 2 = 01. Determine φ 2. (A) 0.2 (B) 0.1 (C) 0.4
More informationSTRESS-STRENGTH RELIABILITY ESTIMATION
CHAPTER 5 STRESS-STRENGTH RELIABILITY ESTIMATION 5. Introduction There are appliances (every physical component possess an inherent strength) which survive due to their strength. These appliances receive
More information[D7] PROBABILITY DISTRIBUTION OF OUTSTANDING LIABILITY FROM INDIVIDUAL PAYMENTS DATA Contributed by T S Wright
Faculty and Institute of Actuaries Claims Reserving Manual v.2 (09/1997) Section D7 [D7] PROBABILITY DISTRIBUTION OF OUTSTANDING LIABILITY FROM INDIVIDUAL PAYMENTS DATA Contributed by T S Wright 1. Introduction
More informationChapter 4: Commonly Used Distributions. Statistics for Engineers and Scientists Fourth Edition William Navidi
Chapter 4: Commonly Used Distributions Statistics for Engineers and Scientists Fourth Edition William Navidi 2014 by Education. This is proprietary material solely for authorized instructor use. Not authorized
More informationModeling of Price. Ximing Wu Texas A&M University
Modeling of Price Ximing Wu Texas A&M University As revenue is given by price times yield, farmers income risk comes from risk in yield and output price. Their net profit also depends on input price, but
More informationMuch of what appears here comes from ideas presented in the book:
Chapter 11 Robust statistical methods Much of what appears here comes from ideas presented in the book: Huber, Peter J. (1981), Robust statistics, John Wiley & Sons (New York; Chichester). There are many
More informationTHE USE OF THE LOGNORMAL DISTRIBUTION IN ANALYZING INCOMES
International Days of tatistics and Economics Prague eptember -3 011 THE UE OF THE LOGNORMAL DITRIBUTION IN ANALYZING INCOME Jakub Nedvěd Abstract Object of this paper is to examine the possibility of
More informationESTIMATION OF MODIFIED MEASURE OF SKEWNESS. Elsayed Ali Habib *
Electronic Journal of Applied Statistical Analysis EJASA, Electron. J. App. Stat. Anal. (2011), Vol. 4, Issue 1, 56 70 e-issn 2070-5948, DOI 10.1285/i20705948v4n1p56 2008 Università del Salento http://siba-ese.unile.it/index.php/ejasa/index
More informationIntroduction to Algorithmic Trading Strategies Lecture 8
Introduction to Algorithmic Trading Strategies Lecture 8 Risk Management Haksun Li haksun.li@numericalmethod.com www.numericalmethod.com Outline Value at Risk (VaR) Extreme Value Theory (EVT) References
More informationPoint Estimation. Stat 4570/5570 Material from Devore s book (Ed 8), and Cengage
6 Point Estimation Stat 4570/5570 Material from Devore s book (Ed 8), and Cengage Point Estimation Statistical inference: directed toward conclusions about one or more parameters. We will use the generic
More informationFINITE SAMPLE DISTRIBUTIONS OF RISK-RETURN RATIOS
Available Online at ESci Journals Journal of Business and Finance ISSN: 305-185 (Online), 308-7714 (Print) http://www.escijournals.net/jbf FINITE SAMPLE DISTRIBUTIONS OF RISK-RETURN RATIOS Reza Habibi*
More informationCommonly Used Distributions
Chapter 4: Commonly Used Distributions 1 Introduction Statistical inference involves drawing a sample from a population and analyzing the sample data to learn about the population. We often have some knowledge
More informationAn Information Based Methodology for the Change Point Problem Under the Non-central Skew t Distribution with Applications.
An Information Based Methodology for the Change Point Problem Under the Non-central Skew t Distribution with Applications. Joint with Prof. W. Ning & Prof. A. K. Gupta. Department of Mathematics and Statistics
More informationQQ PLOT Yunsi Wang, Tyler Steele, Eva Zhang Spring 2016
QQ PLOT INTERPRETATION: Quantiles: QQ PLOT Yunsi Wang, Tyler Steele, Eva Zhang Spring 2016 The quantiles are values dividing a probability distribution into equal intervals, with every interval having
More informationA Skewed Truncated Cauchy Logistic. Distribution and its Moments
International Mathematical Forum, Vol. 11, 2016, no. 20, 975-988 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/imf.2016.6791 A Skewed Truncated Cauchy Logistic Distribution and its Moments Zahra
More informationA Robust Test for Normality
A Robust Test for Normality Liangjun Su Guanghua School of Management, Peking University Ye Chen Guanghua School of Management, Peking University Halbert White Department of Economics, UCSD March 11, 2006
More informationContinuous Distributions
Quantitative Methods 2013 Continuous Distributions 1 The most important probability distribution in statistics is the normal distribution. Carl Friedrich Gauss (1777 1855) Normal curve A normal distribution
More informationdiscussion Papers Some Flexible Parametric Models for Partially Adaptive Estimators of Econometric Models
discussion Papers Discussion Paper 2007-13 March 26, 2007 Some Flexible Parametric Models for Partially Adaptive Estimators of Econometric Models Christian B. Hansen Graduate School of Business at the
More informationStatistical Methodology. A note on a two-sample T test with one variance unknown
Statistical Methodology 8 (0) 58 534 Contents lists available at SciVerse ScienceDirect Statistical Methodology journal homepage: www.elsevier.com/locate/stamet A note on a two-sample T test with one variance
More informationEstimation of Large Insurance Losses: A Case Study
University of Nebraska - Lincoln DigitalCommons@University of Nebraska - Lincoln Journal of Actuarial Practice 1993-2006 Finance Department 2006 Estimation of Large Insurance Losses: A Case Study Tine
More informationProcess capability estimation for non normal quality characteristics: A comparison of Clements, Burr and Box Cox Methods
ANZIAM J. 49 (EMAC2007) pp.c642 C665, 2008 C642 Process capability estimation for non normal quality characteristics: A comparison of Clements, Burr and Box Cox Methods S. Ahmad 1 M. Abdollahian 2 P. Zeephongsekul
More informationFolded- and Log-Folded-t Distributions as Models for Insurance Loss Data
Folded- and Log-Folded-t Distributions as Models for Insurance Loss Data Vytaras Brazauskas University of Wisconsin-Milwaukee Andreas Kleefeld University of Wisconsin-Milwaukee Revised: September 009 (Submitted:
More informationExam STAM Practice Exam #1
!!!! Exam STAM Practice Exam #1 These practice exams should be used during the month prior to your exam. This practice exam contains 20 questions, of equal value, corresponding to about a 2 hour exam.
More informationINDIAN INSTITUTE OF SCIENCE STOCHASTIC HYDROLOGY. Lecture -5 Course Instructor : Prof. P. P. MUJUMDAR Department of Civil Engg., IISc.
INDIAN INSTITUTE OF SCIENCE STOCHASTIC HYDROLOGY Lecture -5 Course Instructor : Prof. P. P. MUJUMDAR Department of Civil Engg., IISc. Summary of the previous lecture Moments of a distribubon Measures of
More informationStatistical Modeling Techniques for Reserve Ranges: A Simulation Approach
Statistical Modeling Techniques for Reserve Ranges: A Simulation Approach by Chandu C. Patel, FCAS, MAAA KPMG Peat Marwick LLP Alfred Raws III, ACAS, FSA, MAAA KPMG Peat Marwick LLP STATISTICAL MODELING
More informationMortality Rates Estimation Using Whittaker-Henderson Graduation Technique
MATIMYÁS MATEMATIKA Journal of the Mathematical Society of the Philippines ISSN 0115-6926 Vol. 39 Special Issue (2016) pp. 7-16 Mortality Rates Estimation Using Whittaker-Henderson Graduation Technique
More informationKARACHI UNIVERSITY BUSINESS SCHOOL UNIVERSITY OF KARACHI BS (BBA) VI
88 P a g e B S ( B B A ) S y l l a b u s KARACHI UNIVERSITY BUSINESS SCHOOL UNIVERSITY OF KARACHI BS (BBA) VI Course Title : STATISTICS Course Number : BA(BS) 532 Credit Hours : 03 Course 1. Statistical
More informationStochastic model of flow duration curves for selected rivers in Bangladesh
Climate Variability and Change Hydrological Impacts (Proceedings of the Fifth FRIEND World Conference held at Havana, Cuba, November 2006), IAHS Publ. 308, 2006. 99 Stochastic model of flow duration curves
More information2.1 Random variable, density function, enumerative density function and distribution function
Risk Theory I Prof. Dr. Christian Hipp Chair for Science of Insurance, University of Karlsruhe (TH Karlsruhe) Contents 1 Introduction 1.1 Overview on the insurance industry 1.1.1 Insurance in Benin 1.1.2
More informationContinuous-Time Pension-Fund Modelling
. Continuous-Time Pension-Fund Modelling Andrew J.G. Cairns Department of Actuarial Mathematics and Statistics, Heriot-Watt University, Riccarton, Edinburgh, EH4 4AS, United Kingdom Abstract This paper
More informationChapter 2 Uncertainty Analysis and Sampling Techniques
Chapter 2 Uncertainty Analysis and Sampling Techniques The probabilistic or stochastic modeling (Fig. 2.) iterative loop in the stochastic optimization procedure (Fig..4 in Chap. ) involves:. Specifying
More informationTwo hours. To be supplied by the Examinations Office: Mathematical Formula Tables and Statistical Tables THE UNIVERSITY OF MANCHESTER
Two hours MATH20802 To be supplied by the Examinations Office: Mathematical Formula Tables and Statistical Tables THE UNIVERSITY OF MANCHESTER STATISTICAL METHODS Answer any FOUR of the SIX questions.
More informationEVA Tutorial #1 BLOCK MAXIMA APPROACH IN HYDROLOGIC/CLIMATE APPLICATIONS. Rick Katz
1 EVA Tutorial #1 BLOCK MAXIMA APPROACH IN HYDROLOGIC/CLIMATE APPLICATIONS Rick Katz Institute for Mathematics Applied to Geosciences National Center for Atmospheric Research Boulder, CO USA email: rwk@ucar.edu
More informationSOCIETY OF ACTUARIES EXAM STAM SHORT-TERM ACTUARIAL MATHEMATICS EXAM STAM SAMPLE QUESTIONS
SOCIETY OF ACTUARIES EXAM STAM SHORT-TERM ACTUARIAL MATHEMATICS EXAM STAM SAMPLE QUESTIONS Questions 1-307 have been taken from the previous set of Exam C sample questions. Questions no longer relevant
More informationMEASURING PORTFOLIO RISKS USING CONDITIONAL COPULA-AR-GARCH MODEL
MEASURING PORTFOLIO RISKS USING CONDITIONAL COPULA-AR-GARCH MODEL Isariya Suttakulpiboon MSc in Risk Management and Insurance Georgia State University, 30303 Atlanta, Georgia Email: suttakul.i@gmail.com,
More informationStatistical analysis and bootstrapping
Statistical analysis and bootstrapping p. 1/15 Statistical analysis and bootstrapping Michel Bierlaire michel.bierlaire@epfl.ch Transport and Mobility Laboratory Statistical analysis and bootstrapping
More informationContinuous random variables
Continuous random variables probability density function (f(x)) the probability distribution function of a continuous random variable (analogous to the probability mass function for a discrete random variable),
More informationIntroduction to Computational Finance and Financial Econometrics Descriptive Statistics
You can t see this text! Introduction to Computational Finance and Financial Econometrics Descriptive Statistics Eric Zivot Summer 2015 Eric Zivot (Copyright 2015) Descriptive Statistics 1 / 28 Outline
More informationDependence Modeling and Credit Risk
Dependence Modeling and Credit Risk Paola Mosconi Banca IMI Bocconi University, 20/04/2015 Paola Mosconi Lecture 6 1 / 53 Disclaimer The opinion expressed here are solely those of the author and do not
More informationWeek 2 Quantitative Analysis of Financial Markets Hypothesis Testing and Confidence Intervals
Week 2 Quantitative Analysis of Financial Markets Hypothesis Testing and Confidence Intervals Christopher Ting http://www.mysmu.edu/faculty/christophert/ Christopher Ting : christopherting@smu.edu.sg :
More informationSTAT 509: Statistics for Engineers Dr. Dewei Wang. Copyright 2014 John Wiley & Sons, Inc. All rights reserved.
STAT 509: Statistics for Engineers Dr. Dewei Wang Applied Statistics and Probability for Engineers Sixth Edition Douglas C. Montgomery George C. Runger 7 Point CHAPTER OUTLINE 7-1 Point Estimation 7-2
More informationA Skewed Truncated Cauchy Uniform Distribution and Its Moments
Modern Applied Science; Vol. 0, No. 7; 206 ISSN 93-844 E-ISSN 93-852 Published by Canadian Center of Science and Education A Skewed Truncated Cauchy Uniform Distribution and Its Moments Zahra Nazemi Ashani,
More informationModelling Environmental Extremes
19th TIES Conference, Kelowna, British Columbia 8th June 2008 Topics for the day 1. Classical models and threshold models 2. Dependence and non stationarity 3. R session: weather extremes 4. Multivariate
More informationSOCIETY OF ACTUARIES/CASUALTY ACTUARIAL SOCIETY EXAM C CONSTRUCTION AND EVALUATION OF ACTUARIAL MODELS EXAM C SAMPLE QUESTIONS
SOCIETY OF ACTUARIES/CASUALTY ACTUARIAL SOCIETY EXAM C CONSTRUCTION AND EVALUATION OF ACTUARIAL MODELS EXAM C SAMPLE QUESTIONS Copyright 2008 by the Society of Actuaries and the Casualty Actuarial Society
More informationModelling Environmental Extremes
19th TIES Conference, Kelowna, British Columbia 8th June 2008 Topics for the day 1. Classical models and threshold models 2. Dependence and non stationarity 3. R session: weather extremes 4. Multivariate
More informationTwo Hours. Mathematical formula books and statistical tables are to be provided THE UNIVERSITY OF MANCHESTER. 22 January :00 16:00
Two Hours MATH38191 Mathematical formula books and statistical tables are to be provided THE UNIVERSITY OF MANCHESTER STATISTICAL MODELLING IN FINANCE 22 January 2015 14:00 16:00 Answer ALL TWO questions
More informationStochastic Models. Statistics. Walt Pohl. February 28, Department of Business Administration
Stochastic Models Statistics Walt Pohl Universität Zürich Department of Business Administration February 28, 2013 The Value of Statistics Business people tend to underestimate the value of statistics.
More informationDistribution analysis of the losses due to credit risk
Distribution analysis of the losses due to credit risk Kamil Łyko 1 Abstract The main purpose of this article is credit risk analysis by analyzing the distribution of losses on retail loans portfolio.
More informationRobust Critical Values for the Jarque-bera Test for Normality
Robust Critical Values for the Jarque-bera Test for Normality PANAGIOTIS MANTALOS Jönköping International Business School Jönköping University JIBS Working Papers No. 00-8 ROBUST CRITICAL VALUES FOR THE
More informationEstimating the Parameters of Closed Skew-Normal Distribution Under LINEX Loss Function
Australian Journal of Basic Applied Sciences, 5(7): 92-98, 2011 ISSN 1991-8178 Estimating the Parameters of Closed Skew-Normal Distribution Under LINEX Loss Function 1 N. Abbasi, 1 N. Saffari, 2 M. Salehi
More informationModeling Credit Risk of Loan Portfolios in the Presence of Autocorrelation (Part 2)
Practitioner Seminar in Financial and Insurance Mathematics ETH Zürich Modeling Credit Risk of Loan Portfolios in the Presence of Autocorrelation (Part 2) Christoph Frei UBS and University of Alberta March
More informationStatistics 431 Spring 2007 P. Shaman. Preliminaries
Statistics 4 Spring 007 P. Shaman The Binomial Distribution Preliminaries A binomial experiment is defined by the following conditions: A sequence of n trials is conducted, with each trial having two possible
More informationSubject CS1 Actuarial Statistics 1 Core Principles. Syllabus. for the 2019 exams. 1 June 2018
` Subject CS1 Actuarial Statistics 1 Core Principles Syllabus for the 2019 exams 1 June 2018 Copyright in this Core Reading is the property of the Institute and Faculty of Actuaries who are the sole distributors.
More informationELEMENTS OF MONTE CARLO SIMULATION
APPENDIX B ELEMENTS OF MONTE CARLO SIMULATION B. GENERAL CONCEPT The basic idea of Monte Carlo simulation is to create a series of experimental samples using a random number sequence. According to the
More informationMeasuring Risk. Review of statistical concepts Probability distribution. Review of statistical concepts Probability distribution 2/1/2018
Measuring Risk Review of statistical concepts Probability distribution Discrete and continuous probability distributions. Discrete: Probability mass function assigns a probability to each possible out-come.
More informationChapter 8: Sampling distributions of estimators Sections
Chapter 8 continued Chapter 8: Sampling distributions of estimators Sections 8.1 Sampling distribution of a statistic 8.2 The Chi-square distributions 8.3 Joint Distribution of the sample mean and sample
More informationRisk Aggregation with Dependence Uncertainty
Risk Aggregation with Dependence Uncertainty Carole Bernard (Grenoble Ecole de Management) Hannover, Current challenges in Actuarial Mathematics November 2015 Carole Bernard Risk Aggregation with Dependence
More informationThe use of flexible quantile-based measures in risk assessment
Institut de Recerca en Economia Aplicada Regional i Pública Document de Treball 2011/07 pàg. 1 Research Institute of Applied Economics Working Paper 2011/07 pag.1 Institut de Recerca en Economia Aplicada
More informationدرس هفتم یادگیري ماشین. (Machine Learning) دانشگاه فردوسی مشهد دانشکده مهندسی رضا منصفی
یادگیري ماشین توزیع هاي نمونه و تخمین نقطه اي پارامترها Sampling Distributions and Point Estimation of Parameter (Machine Learning) دانشگاه فردوسی مشهد دانشکده مهندسی رضا منصفی درس هفتم 1 Outline Introduction
More informationAppendix A. Selecting and Using Probability Distributions. In this appendix
Appendix A Selecting and Using Probability Distributions In this appendix Understanding probability distributions Selecting a probability distribution Using basic distributions Using continuous distributions
More informationSYSM 6304 Risk and Decision Analysis Lecture 2: Fitting Distributions to Data
SYSM 6304 Risk and Decision Analysis Lecture 2: Fitting Distributions to Data M. Vidyasagar Cecil & Ida Green Chair The University of Texas at Dallas Email: M.Vidyasagar@utdallas.edu September 5, 2015
More informationComparing Downside Risk Measures for Heavy Tailed Distributions
Comparing Downside Risk Measures for Heavy Tailed Distributions Jón Daníelsson London School of Economics Mandira Sarma Bjørn N. Jorgensen Columbia Business School Indian Statistical Institute, Delhi EURANDOM,
More informationOperational Risk: Evidence, Estimates and Extreme Values from Austria
Operational Risk: Evidence, Estimates and Extreme Values from Austria Stefan Kerbl OeNB / ECB 3 rd EBA Policy Research Workshop, London 25 th November 2014 Motivation Operational Risk as the exotic risk
More informationComputational Statistics Handbook with MATLAB
«H Computer Science and Data Analysis Series Computational Statistics Handbook with MATLAB Second Edition Wendy L. Martinez The Office of Naval Research Arlington, Virginia, U.S.A. Angel R. Martinez Naval
More informationQuantitative Methods for Economics, Finance and Management (A86050 F86050)
Quantitative Methods for Economics, Finance and Management (A86050 F86050) Matteo Manera matteo.manera@unimib.it Marzio Galeotti marzio.galeotti@unimi.it 1 This material is taken and adapted from Guy Judge
More informationRisk Measure and Allocation Terminology
Notation Ris Measure and Allocation Terminology Gary G. Venter and John A. Major February 2009 Y is a random variable representing some financial metric for a company (say, insured losses) with cumulative
More informationModelling of Long-Term Risk
Modelling of Long-Term Risk Roger Kaufmann Swiss Life roger.kaufmann@swisslife.ch 15th International AFIR Colloquium 6-9 September 2005, Zurich c 2005 (R. Kaufmann, Swiss Life) Contents A. Basel II B.
More informationThe Two Sample T-test with One Variance Unknown
The Two Sample T-test with One Variance Unknown Arnab Maity Department of Statistics, Texas A&M University, College Station TX 77843-343, U.S.A. amaity@stat.tamu.edu Michael Sherman Department of Statistics,
More informationBasic notions of probability theory: continuous probability distributions. Piero Baraldi
Basic notions of probability theory: continuous probability distributions Piero Baraldi Probability distributions for reliability, safety and risk analysis: discrete probability distributions continuous
More informationModeling Co-movements and Tail Dependency in the International Stock Market via Copulae
Modeling Co-movements and Tail Dependency in the International Stock Market via Copulae Katja Ignatieva, Eckhard Platen Bachelier Finance Society World Congress 22-26 June 2010, Toronto K. Ignatieva, E.
More informationAssessing Regime Switching Equity Return Models
Assessing Regime Switching Equity Return Models R. Keith Freeland, ASA, Ph.D. Mary R. Hardy, FSA, FIA, CERA, Ph.D. Matthew Till Copyright 2009 by the Society of Actuaries. All rights reserved by the Society
More information