The Vasicek Distribution
|
|
- Joseph Lawson
- 6 years ago
- Views:
Transcription
1 The Vasicek Distribution Dirk Tasche Lloyds TSB Bank Corporate Markets Rating Systems Bristol / London, August 2008 The opinions expressed in this presentation are those of the author and do not necessarily reflect views of Lloyds TSB Bank. i
2 Outline What is special about the Vasicek distribution? (slide 1) Properties of the Vasicek distribution (slide 8) Modelling default rates (slide 17) Related topics (not covered) (slide 30) ii
3 What is special about the Vasicek distribution? (1) The Vasicek distribution is a two-parametric (0 < p < 1, 0 < ϱ < 1) continuous distribution on (0, 1) with density f p,ϱ (x) = 1 ϱ ϱ exp ( 1 2 { Φ 1 (x) 2 ( 1 ϱ Φ 1 (x) Φ 1 (p) ϱ ) 2 }). (1.1) Other two-parametric continuous distributions on (0, 1) include the beta-distribution with parameters α, β > 0 and density f α,β (x) = Γ(α+β) Γ(α) Γ(β) xα 1 (1 x) β 1, the Kumaraswamy distribution with parameters a, b > 0 and density f a,b (x) = a b x a 1 (1 x a ) b 1. 1
4 What is special about the Vasicek distribution? (2) Examples of densities 3.0 Mean 0.25, Stddev Density value Vasicek Beta Kumaraswamy % 25% 50% 75% 100% Loss percentage 2
5 What is special about the Vasicek distribution? (3) Examples of densities 5 Mean 0.25, Stddev Density value 3 2 Vasicek Beta Kumaraswamy 1 0 0% 25% 50% 75% 100% Loss percentage 3
6 What is special about the Vasicek distribution? (4) Observations: When fitted to the same mean and standard deviation, the Vasicek, beta, and Kumaraswamy densities do not differ much. Vasicek: Density looks sophisticated, moment matching of correlation parameter ϱ requires numerical root-finding. Kumaraswamy: Density looks simple, but moment matching requires numerical solution of two-dimensional optimisation problem. beta: Density looks relatively simple, moment matching is easy. Why not just working with the beta distribution? 4
7 What is special about the Vasicek distribution? (5) Kumaraswamy: Similar to beta but more efficient for simulation. beta: distribution Appears naturally as distribution of order statistics: X 1,..., X n i.i.d. uniformly distributed, X [k] = min I {1,...,n}, I =k max{x i : i I} (k-th smallest) X [k] β(k, n k + 1), i.e. P[X [k] x] = Γ(n) Γ(k) Γ(n k+1) x 0 uk 1 (1 u) n k du. (5.1) Dirichlet distribution is multi-variate extension of beta distribution. However, there is no economic motivation of the beta distribution. 5
8 What is special about the Vasicek distribution? (6) It s the economy, stupid Vasicek distribution can be interpreted in context of a trigger mechanism that is useful for modeling credit risk. Denote by A the change in value of an obligor s assets in a fixed time period (say 1 year). Denote by t the default threshold. The obligor s default is triggered if A t. Then P D = P[A t] is the unconditional probability of default. 6
9 What is special about the Vasicek distribution? (7) Oldrich Vasicek (1987) : Assume that A = ϱ S + 1 ϱ ξ where S is the systematic and ξ the idiosyncratic risk factor. S and ξ are independent and standard normal. Consider the conditional (on the systematic factor) probability of default P D(S) : P D(S) = P[A t S] = Φ ( t ϱ S ) 1 ϱ (7.1) Then P D(S) is a Vasicek-distributed random variable with p = Φ(t). Probability of Loss on Loan Portfolio, wp.html Φ denotes the standard normal distribution function: Φ(x) = e u2 /2 2 π du. 7
10 Properties of the Vasicek distribution (1) Important observation: If X has a Vasicek distribution with parameters p and ϱ, then by Slide 7 X can be represented as X = Φ ( t ϱ S 1 ϱ ), (8.1) where t = Φ 1 (p) and S is a standard normal random variable. Consequences: Distribution function: F p,ϱ (x) = P[X x] = Φ ( 1 ϱ Φ 1 ) (x) t ϱ Density: f p,ϱ (x) = x F p,ϱ(x) (see (1.1)) Quantiles: Fp,ϱ 1 (α) = Φ ( t+ ϱ Φ 1 (α) ) 1 ϱ 8
11 Properties of the Vasicek distribution (2) Moments Let n be a positive integer and ξ 1,..., ξ n i.i.d. standard normal. If X is Vasicek-distributed with parameters p = Φ(t) and ϱ, then E[X n ] = E [ Φ ( t ϱ S ) n ] 1 ϱ = E [ n i=1 P[ ϱ S + 1 ϱ ξ i t S] ] = E [ P[ ϱ S + 1 ϱ ξ 1 t,..., ϱ S + 1 ϱ ξ n t S] ] = P[ ϱ S + 1 ϱ ξ 1 t,..., ϱ S + 1 ϱ ξ n t] = P[Y 1 t,..., Y n t], (9.1) where (Y 1,..., Y n ) is a multi-variate normal vector with E[Y i ] = 0, var[y i ] = 1, and corr[y i, Y j ] = ϱ, i j. 9
12 Properties of the Vasicek distribution (3) Special cases: E[X] = P[Y t] = p var[x] = E[X 2 ] E[X] 2 = P[Y 1 t, Y 2 t] p 2 (10.1) Define Φ 2 (s, t; ϱ) = P[Y 1 s, Y 2 t], where (Y 1, Y 2 ) is a bivariate normal vector with E[Y i ] = 0, var[y i ] = 1, and corr[y 1, Y 2 ] = ϱ. Estimating p and ϱ from a sample x 1,..., x m (0, 1) by momentmatching: ˆp = 1 m m x i i=1 ˆϱ is unique solution of Φ 2 (t, t; ϱ) = 1 m m x 2 i i=1 10
13 Properties of the Vasicek distribution (4) Indirect moment matching Observation: If X is Vasicek-distributed with parameters p and ϱ, then Y = Φ 1 (X) is normally distributed with E[Y ] = µ = Φ 1 (p) 1 ϱ, var[y ] = σ 2 = ϱ 1 ϱ (11.1) Estimate µ and σ by ˆµ = 1 m m i=1 Φ 1 (x i ), ˆσ 2 = 1 m m i=1 Φ 1 (x i ) 2 ˆµ 2 (11.2) Solve then (11.1) for p and ϱ: p = Φ ( ) ˆµ 1+ˆσ 2, ϱ = ˆσ 2 1+ˆσ 2 (11.3) ( p, ϱ) is also the maximum likelihood estimator of (p, ϱ). 11
14 Properties of the Vasicek distribution (5) Quantile-based estimators Use again observation (11.1) on the transformation of the sample into a normally distributed sample. Choose probabilities 0 < α 1 < α 2 < 1 and calculate the empirical quantiles ˆq(α 1 ) < ˆq(α 2 ) of the transformed sample. Determine estimators ˆµ and ˆσ from the equations ˆq(α i ) = ˆµ + ˆσ Φ 1 (α i ), i = 1, 2 (12.1) Determine estimator p of p and ϱ of ϱ by solving (11.3) for p and ϱ. 12
15 Properties of the Vasicek distribution (6) Simulation study of estimation efficiency Simulate 1000 times a Vasicek sample of length 5, 10, 25, and 100, for parameter values p = 0.1 and ϱ = Compare results for direct moment matching (DMM) indirect moment matching (= maximum likelihood estimation, MLE) quantile-based estimation (QBE) Tabulate estimators (Est.), standard deviations (SD) of estimators, and standard errors (SE). 13
16 Properties of the Vasicek distribution (7) Results for estimation of p Sample size 5 10 Est. p SD SE Est. p SD SE DMM MLE QBE Sample size Est. p SD SE Est. p SD SE DMM MLE QBE
17 Properties of the Vasicek distribution (8) Results for estimation of ϱ Sample size 5 10 Est. ϱ SD SE Est. ϱ SD SE DMM MLE QBE Sample size Est. ϱ SD SE Est. ϱ SD SE DMM MLE QBE
18 Properties of the Vasicek distribution (9) Comments on the simulation study Estimation of p: Little difference between Direct Moment Matching (DMM) and Maximum Likelihood Estimation (MLE). Quantile-Based Estimation (QBE) is clearly less efficient (higher Standard Error). All estimators seem nearly unbiased (i.e. the mean of many estimates is close to the true value) for sample size 10 or greater. Estimation of ϱ: MLE is most efficient (as expected from statistical theory), DMM is second, QBE is worst. None of the estimators is unbiased. In average, they all underestimate the correlation parameter, even for relatively large sample size
19 Modelling default rates (1) Future default rates of a portfolio cannot be predicted with certainty. Based on appropriate probabilistic models, however, probabilities of interesting future events can be calculated (e.g. the probability to observe a default rate of more than 1%). Models should be able to explain clustering of defaults (i.e. very high observed rates in some years, very low rates in other years). Common approach: create default correlation by dependence on one or more systematic factors. 17
20 Modelling default rates (2) Binomial approach n independent borrowers Default event: D i = { Borrower i defaults } Default indicator: 1 Di = 1, if D i occurs 0, otherwise Identical default probability: P D = P[D i ] = P[1 Di = 1] Default rate: DR = 1 n n i=1 1 Di 18
21 Modelling default rates (3) Calculating the default rate distribution Exact binomial probability by incomplete beta integrals: ( ) k n = P D l (1 P D) n l l [ P DR k ] n Normal approximation: = l=0 1 P D 1 0 u k (1 u) n k 1 du u k (1 u) n k 1 du P [ DR k n ] Φ ( ) k n P D n P D (1 P D) (19.1) Poisson approximation: P [ DR k n ] e n P D k l=0 (n P D) l l! 19
22 Modelling default rates (4) Comments on binomial calculations On principle, there is no reason to use approximations for binomial distributions. Calculations are fast (even in MSExcel) when the incomplete beta representation is used. For PDs up to 10% the Poisson approximation is better than the normal approximation. The Poisson approximation is useful for algebraic calculations with default or loss rates (actuarial approach, CreditRisk + ). 20
23 Modelling default rates (5) General observations on the binomial approach E[DR] = P D, var[dr] = 1 n P D (1 P D) Hence lim var[dr] = 0 n Conclusion: If default events were independent, for large pools of borrowers with similar credit-worthiness the variation of observed default rates over time would be small. Contradiction to empirically observed credit cycles. Problem: Independence of default events P[Borrowers i and j default] = P[D i D j ] = P D 2 (21.1) 21
24 Modelling default rates (6) Vasicek approach: correlated binomial distribution Choose representation of default event D i such that defaults become dependent (cf. Slide 7): D i = { ϱ S + 1 ϱ ξ i t} (22.1) S is the systematic factor, ξ i is the idiosyncratic factor, and t = Φ 1 (P D). S and ξ i, i = 1,..., n, are independent and standard normal. Consequence: Defaults are no longer independent: P[Borrowers i and j default] = P[D i D j ] = Φ 2 (t, t; ϱ) > P D 2 (22.2) 22
25 Modelling default rates (7) Observations on the Vasicek approach E[DR] = P D, var[dr] = 1 n P D P D2 + n 1 n Φ 2(t, t; ϱ) Hence lim var[dr] = Φ n 2 (t, t; ϱ) P D 2 > 0 Conclusion: The Vasicek approach can depict cyclic variation of default rates even for very large portfolios. Expectation given systematic factor (from (7.1)) is Vasicek-distributed: E[DR S] = Φ ( t ϱ S ) 1 ϱ (23.1) Large portfolio behaviour (approximation by E[DR S]): lim n E[ (DR E[DR S]) 2] ( = lim 1 n n P D Φ2 (t, t; ϱ) ) = 0 (23.2) 23
26 Modelling default rates (8) Binomial and correlated binomial distributions Frequency PD = 0.1, rho = 0.12, n = Default rate Binomial Corr Binomial 24
27 Modelling default rates (9) Vasicek and correlated binomial distributions Frequency * n PD = 0.1, rho = 0.12, n = 100 Corr binomial Vasicek Default rate 25
28 Modelling default rates (10) Vasicek and correlated binomial distributions Frequency * n PD = 0.1, rho = 0.12, n = 250 Corr binomial Vasicek Default rate 26
29 Modelling default rates (11) Vasicek approach: Calculating the default rate distribution Exact probability by conditioning and integrating against density of systematic factor: P [ DR k n ] = E = [ P [ ϕ(s) DR k n k l=0 ]] S ( ) n Φ ( t ϱ s ) ( l l 1 ϱ 1 Φ ( t ϱ s ) ) n l 1 ϱ ds Alternative: by Monte-Carlo simulation (less accurate). Density of E[DR S]: by (1.1). ϕ denotes the standard normal density function: ϕ(s) = e s2 /2 2 π ds. (27.1) 27
30 Modelling default rates (12) Vasicek (2002) Fit the distribution of DR (correlated binomial) by a Vasicekdistributed variable X with parameters p = P D and ϱ such that Φ 2 (t, t; ϱ ) P D 2 = var[x] = var[dr] = 1 n P D P D2 + n 1 n Φ 2(t, t; ϱ) (28.1) Involves numerical root-finding and calculation of bivariate normal distribution. Kolmogorov-Smirnov distances for fitting correlated binomial distribution: (exact distribution by (27.1)) Approximation by Vasicek-distribution with P D, ϱ: Approximation by Vasicek-distribution with P D, ϱ : The distribution of loan portfolio value, RISK 15,
31 Modelling default rates (13) Vasicek fit and correlated binomial distributions Frequency * n PD = 0.1, rho = 0.12, rho* = 0.143, n = 100 Corr binomial Vasicek fit Vasicek Default rate 29
32 Related topics (not covered) Default rate modelling with negative binomial (conditional Poisson) distributions. Across-time default rate modelling. Multi-factor models. Portfolio-loss models (heterogeneity of exposures). Derivation of Basel II risk weight functions
Bayesian estimation of probabilities of default for low default portfolios
Bayesian estimation of probabilities of default for low default portfolios Dirk Tasche arxiv:1112.555v3 [q-fin.rm] 5 Apr 212 First version: December 23, 211 This version: April 5, 212 The estimation of
More informationDefinition 9.1 A point estimate is any function T (X 1,..., X n ) of a random sample. We often write an estimator of the parameter θ as ˆθ.
9 Point estimation 9.1 Rationale behind point estimation When sampling from a population described by a pdf f(x θ) or probability function P [X = x θ] knowledge of θ gives knowledge of the entire population.
More informationGRANULARITY ADJUSTMENT FOR DYNAMIC MULTIPLE FACTOR MODELS : SYSTEMATIC VS UNSYSTEMATIC RISKS
GRANULARITY ADJUSTMENT FOR DYNAMIC MULTIPLE FACTOR MODELS : SYSTEMATIC VS UNSYSTEMATIC RISKS Patrick GAGLIARDINI and Christian GOURIÉROUX INTRODUCTION Risk measures such as Value-at-Risk (VaR) Expected
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 informationMonte Carlo and Empirical Methods for Stochastic Inference (MASM11/FMSN50)
Monte Carlo and Empirical Methods for Stochastic Inference (MASM11/FMSN50) Magnus Wiktorsson Centre for Mathematical Sciences Lund University, Sweden Lecture 5 Sequential Monte Carlo methods I January
More informationUQ, STAT2201, 2017, Lectures 3 and 4 Unit 3 Probability Distributions.
UQ, STAT2201, 2017, Lectures 3 and 4 Unit 3 Probability Distributions. Random Variables 2 A random variable X is a numerical (integer, real, complex, vector etc.) summary of the outcome of the random experiment.
More informationPoint Estimators. STATISTICS Lecture no. 10. Department of Econometrics FEM UO Brno office 69a, tel
STATISTICS Lecture no. 10 Department of Econometrics FEM UO Brno office 69a, tel. 973 442029 email:jiri.neubauer@unob.cz 8. 12. 2009 Introduction Suppose that we manufacture lightbulbs and we want to state
More informationChapter 5. Statistical inference for Parametric Models
Chapter 5. Statistical inference for Parametric Models Outline Overview Parameter estimation Method of moments How good are method of moments estimates? Interval estimation Statistical Inference for Parametric
More informationsuch that P[L i where Y and the Z i ~ B(1, p), Negative binomial distribution 0.01 p = 0.3%, ρ = 10%
Irreconcilable differences As Basel has acknowledged, the leading credit portfolio models are equivalent in the case of a single systematic factor. With multiple factors, considerable differences emerge,
More informationLimit Theorems for the Empirical Distribution Function of Scaled Increments of Itô Semimartingales at high frequencies
Limit Theorems for the Empirical Distribution Function of Scaled Increments of Itô Semimartingales at high frequencies George Tauchen Duke University Viktor Todorov Northwestern University 2013 Motivation
More informationVersion A. Problem 1. Let X be the continuous random variable defined by the following pdf: 1 x/2 when 0 x 2, f(x) = 0 otherwise.
Math 224 Q Exam 3A Fall 217 Tues Dec 12 Version A Problem 1. Let X be the continuous random variable defined by the following pdf: { 1 x/2 when x 2, f(x) otherwise. (a) Compute the mean µ E[X]. E[X] x
More informationcontinuous rv Note for a legitimate pdf, we have f (x) 0 and f (x)dx = 1. For a continuous rv, P(X = c) = c f (x)dx = 0, hence
continuous rv Let X be a continuous rv. Then a probability distribution or probability density function (pdf) of X is a function f(x) such that for any two numbers a and b with a b, P(a X b) = b a f (x)dx.
More informationLecture 10: Point Estimation
Lecture 10: Point Estimation MSU-STT-351-Sum-17B (P. Vellaisamy: MSU-STT-351-Sum-17B) Probability & Statistics for Engineers 1 / 31 Basic Concepts of Point Estimation A point estimate of a parameter θ,
More informationGeneralized MLE per Martins and Stedinger
Generalized MLE per Martins and Stedinger Martins ES and Stedinger JR. (March 2000). Generalized maximum-likelihood generalized extreme-value quantile estimators for hydrologic data. Water Resources Research
More informationEstimating LGD Correlation
Estimating LGD Correlation Jiří Witzany University of Economics, Prague Abstract: The paper proposes a new method to estimate correlation of account level Basle II Loss Given Default (LGD). The correlation
More information6. Continous Distributions
6. Continous Distributions Chris Piech and Mehran Sahami May 17 So far, all random variables we have seen have been discrete. In all the cases we have seen in CS19 this meant that our RVs could only take
More informationBivariate Birnbaum-Saunders Distribution
Department of Mathematics & Statistics Indian Institute of Technology Kanpur January 2nd. 2013 Outline 1 Collaborators 2 3 Birnbaum-Saunders Distribution: Introduction & Properties 4 5 Outline 1 Collaborators
More 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 informationROM SIMULATION Exact Moment Simulation using Random Orthogonal Matrices
ROM SIMULATION Exact Moment Simulation using Random Orthogonal Matrices Bachelier Finance Society Meeting Toronto 2010 Henley Business School at Reading Contact Author : d.ledermann@icmacentre.ac.uk Alexander
More informationStatistical Analysis of Data from the Stock Markets. UiO-STK4510 Autumn 2015
Statistical Analysis of Data from the Stock Markets UiO-STK4510 Autumn 2015 Sampling Conventions We observe the price process S of some stock (or stock index) at times ft i g i=0,...,n, we denote it by
More informationHomework Assignments
Homework Assignments Week 1 (p. 57) #4.1, 4., 4.3 Week (pp 58 6) #4.5, 4.6, 4.8(a), 4.13, 4.0, 4.6(b), 4.8, 4.31, 4.34 Week 3 (pp 15 19) #1.9, 1.1, 1.13, 1.15, 1.18 (pp 9 31) #.,.6,.9 Week 4 (pp 36 37)
More informationModeling and Estimation of
Modeling and of Financial and Actuarial Mathematics Christian Doppler Laboratory for Portfolio Risk Management Vienna University of Technology PRisMa 2008 29.09.2008 Outline 1 2 3 4 5 Credit ratings describe
More informationEconomi Capital. Tiziano Bellini. Università di Bologna. November 29, 2013
Economi Capital Tiziano Bellini Università di Bologna November 29, 2013 Tiziano Bellini (Università di Bologna) Economi Capital November 29, 2013 1 / 16 Outline Framework Economic Capital Structural approach
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 informationPortfolio Models and ABS
Tutorial 4 Portfolio Models and ABS Loïc BRI François CREI Tutorial 4 Portfolio Models and ABS École ationale des Ponts et Chausées Département Ingénieurie Mathématique et Informatique Master II Loïc BRI
More informationPIVOTAL QUANTILE ESTIMATES IN VAR CALCULATIONS. Peter Schaller, Bank Austria Creditanstalt (BA-CA) Wien,
PIVOTAL QUANTILE ESTIMATES IN VAR CALCULATIONS Peter Schaller, Bank Austria Creditanstalt (BA-CA) Wien, peter@ca-risc.co.at c Peter Schaller, BA-CA, Strategic Riskmanagement 1 Contents Some aspects of
More informationCredit Portfolio Risk
Credit Portfolio Risk Tiziano Bellini Università di Bologna November 29, 2013 Tiziano Bellini (Università di Bologna) Credit Portfolio Risk November 29, 2013 1 / 47 Outline Framework Credit Portfolio Risk
More informationChapter 4: Asymptotic Properties of MLE (Part 3)
Chapter 4: Asymptotic Properties of MLE (Part 3) Daniel O. Scharfstein 09/30/13 1 / 1 Breakdown of Assumptions Non-Existence of the MLE Multiple Solutions to Maximization Problem Multiple Solutions to
More informationSlides for Risk Management
Slides for Risk Management Introduction to the modeling of assets Groll Seminar für Finanzökonometrie Prof. Mittnik, PhD Groll (Seminar für Finanzökonometrie) Slides for Risk Management Prof. Mittnik,
More informationMVE051/MSG Lecture 7
MVE051/MSG810 2017 Lecture 7 Petter Mostad Chalmers November 20, 2017 The purpose of collecting and analyzing data Purpose: To build and select models for parts of the real world (which can be used for
More informationLecture 23. STAT 225 Introduction to Probability Models April 4, Whitney Huang Purdue University. Normal approximation to Binomial
Lecture 23 STAT 225 Introduction to Probability Models April 4, 2014 approximation Whitney Huang Purdue University 23.1 Agenda 1 approximation 2 approximation 23.2 Characteristics of the random variable:
More informationLecture outline. Monte Carlo Methods for Uncertainty Quantification. Importance Sampling. Importance Sampling
Lecture outline Monte Carlo Methods for Uncertainty Quantification Mike Giles Mathematical Institute, University of Oxford KU Leuven Summer School on Uncertainty Quantification Lecture 2: Variance reduction
More informationCPSC 540: Machine Learning
CPSC 540: Machine Learning Monte Carlo Methods Mark Schmidt University of British Columbia Winter 2019 Last Time: Markov Chains We can use Markov chains for density estimation, d p(x) = p(x 1 ) p(x }{{}
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 informationThinking positively. Katja Pluto and Dirk Tasche. July Abstract
Thinking positively Katja Pluto and Dirk Tasche July 2005 Abstract How to come up with numerical PD estimates if there are no default observations? Katja Pluto and Dirk Tasche propose a statistically based
More informationChapter 6 Continuous Probability Distributions. Learning objectives
Chapter 6 Continuous s Slide 1 Learning objectives 1. Understand continuous probability distributions 2. Understand Uniform distribution 3. Understand Normal distribution 3.1. Understand Standard normal
More informationBayesian Approach to PD Calibration and Stress-testing in Low Default Portfolios
Journal of Applied Finance & Banking, vol. 7, no. 2, 2017, 83-98 ISSN: 1792-6580 (print version), 1792-6599 (online) Scienpress Ltd, 2017 Bayesian Approach to PD Calibration and Stress-testing in Low Default
More informationPreprint: Will be published in Perm Winter School Financial Econometrics and Empirical Market Microstructure, Springer
STRESS-TESTING MODEL FOR CORPORATE BORROWER PORTFOLIOS. Preprint: Will be published in Perm Winter School Financial Econometrics and Empirical Market Microstructure, Springer Seleznev Vladimir Denis Surzhko,
More informationPortfolio Credit Risk Models
Portfolio Credit Risk Models Paul Embrechts London School of Economics Department of Accounting and Finance AC 402 FINANCIAL RISK ANALYSIS Lent Term, 2003 c Paul Embrechts and Philipp Schönbucher, 2003
More informationCPSC 540: Machine Learning
CPSC 540: Machine Learning Monte Carlo Methods Mark Schmidt University of British Columbia Winter 2018 Last Time: Markov Chains We can use Markov chains for density estimation, p(x) = p(x 1 ) }{{} d p(x
More informationMAFS Computational Methods for Pricing Structured Products
MAFS550 - Computational Methods for Pricing Structured Products Solution to Homework Two Course instructor: Prof YK Kwok 1 Expand f(x 0 ) and f(x 0 x) at x 0 into Taylor series, where f(x 0 ) = f(x 0 )
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 informationEE641 Digital Image Processing II: Purdue University VISE - October 29,
EE64 Digital Image Processing II: Purdue University VISE - October 9, 004 The EM Algorithm. Suffient Statistics and Exponential Distributions Let p(y θ) be a family of density functions parameterized by
More informationP2.T6. Credit Risk Measurement & Management. Malz, Financial Risk Management: Models, History & Institutions
P2.T6. Credit Risk Measurement & Management Malz, Financial Risk Management: Models, History & Institutions Portfolio Credit Risk Bionic Turtle FRM Video Tutorials By David Harper, CFA FRM 1 Portfolio
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 informationEstimation Y 3. Confidence intervals I, Feb 11,
Estimation Example: Cholesterol levels of heart-attack patients Data: Observational study at a Pennsylvania medical center blood cholesterol levels patients treated for heart attacks measurements 2, 4,
More informationMath489/889 Stochastic Processes and Advanced Mathematical Finance Homework 5
Math489/889 Stochastic Processes and Advanced Mathematical Finance Homework 5 Steve Dunbar Due Fri, October 9, 7. Calculate the m.g.f. of the random variable with uniform distribution on [, ] and then
More informationSampling and sampling distribution
Sampling and sampling distribution September 12, 2017 STAT 101 Class 5 Slide 1 Outline of Topics 1 Sampling 2 Sampling distribution of a mean 3 Sampling distribution of a proportion STAT 101 Class 5 Slide
More informationReview of the Topics for Midterm I
Review of the Topics for Midterm I STA 100 Lecture 9 I. Introduction The objective of statistics is to make inferences about a population based on information contained in a sample. A population is the
More informationBonn Econ Discussion Papers
Bonn Econ Discussion Papers Discussion Paper 16/2001 Factor Models for Portofolio Credit Risk by Philipp J. Schönbucher December 2000 Bonn Graduate School of Economics Department of Economics University
More informationFinancial Risk Management
Financial Risk Management Professor: Thierry Roncalli Evry University Assistant: Enareta Kurtbegu Evry University Tutorial exercices #4 1 Correlation and copulas 1. The bivariate Gaussian copula is given
More informationAsset Correlation of German Corporate Obligors: Its Estimation, Its Drivers and Implications for Regulatory Capital
Asset Correlation of German Corporate Obligors: Its Estimation, Its Drivers and Implications for Regulatory Capital Klaus Düllmann and Harald Scheule March 2003 Abstract This paper addresses the gap between
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 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 informationChapter 7: Point Estimation and Sampling Distributions
Chapter 7: Point Estimation and Sampling Distributions Seungchul Baek Department of Statistics, University of South Carolina STAT 509: Statistics for Engineers 1 / 20 Motivation In chapter 3, we learned
More informationLecture 5: Fundamentals of Statistical Analysis and Distributions Derived from Normal Distributions
Lecture 5: Fundamentals of Statistical Analysis and Distributions Derived from Normal Distributions ELE 525: Random Processes in Information Systems Hisashi Kobayashi Department of Electrical Engineering
More informationDesirable properties for a good model of portfolio credit risk modelling
3.3 Default correlation binomial models Desirable properties for a good model of portfolio credit risk modelling Default dependence produce default correlations of a realistic magnitude. Estimation number
More informationParameter uncertainty for integrated risk capital calculations based on normally distributed subrisks
Parameter uncertainty for integrated risk capital calculations based on normally distributed subrisks Andreas Fröhlich and Annegret Weng March 7, 017 Abstract In this contribution we consider the overall
More informationSystematic Risk in Homogeneous Credit Portfolios
Systematic Risk in Homogeneous Credit Portfolios Christian Bluhm and Ludger Overbeck Systematic Risk in Credit Portfolios In credit portfolios (see [5] for an introduction) there are typically two types
More informationDecision theoretic estimation of the ratio of variances in a bivariate normal distribution 1
Decision theoretic estimation of the ratio of variances in a bivariate normal distribution 1 George Iliopoulos Department of Mathematics University of Patras 26500 Rio, Patras, Greece Abstract In this
More informationOption Pricing for Discrete Hedging and Non-Gaussian Processes
Option Pricing for Discrete Hedging and Non-Gaussian Processes Kellogg College University of Oxford A thesis submitted in partial fulfillment of the requirements for the MSc in Mathematical Finance November
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 informationLecture Notes 6. Assume F belongs to a family of distributions, (e.g. F is Normal), indexed by some parameter θ.
Sufficient Statistics Lecture Notes 6 Sufficiency Data reduction in terms of a particular statistic can be thought of as a partition of the sample space X. Definition T is sufficient for θ if the conditional
More informationVaR Estimation under Stochastic Volatility Models
VaR Estimation under Stochastic Volatility Models Chuan-Hsiang Han Dept. of Quantitative Finance Natl. Tsing-Hua University TMS Meeting, Chia-Yi (Joint work with Wei-Han Liu) December 5, 2009 Outline Risk
More informationSTAT/MATH 395 PROBABILITY II
STAT/MATH 395 PROBABILITY II Distribution of Random Samples & Limit Theorems Néhémy Lim University of Washington Winter 2017 Outline Distribution of i.i.d. Samples Convergence of random variables The Laws
More informationA GENERAL FORMULA FOR OPTION PRICES IN A STOCHASTIC VOLATILITY MODEL. Stephen Chin and Daniel Dufresne. Centre for Actuarial Studies
A GENERAL FORMULA FOR OPTION PRICES IN A STOCHASTIC VOLATILITY MODEL Stephen Chin and Daniel Dufresne Centre for Actuarial Studies University of Melbourne Paper: http://mercury.ecom.unimelb.edu.au/site/actwww/wps2009/no181.pdf
More informationModelling financial data with stochastic processes
Modelling financial data with stochastic processes Vlad Ardelean, Fabian Tinkl 01.08.2012 Chair of statistics and econometrics FAU Erlangen-Nuremberg Outline Introduction Stochastic processes Volatility
More informationStatistical Tables Compiled by Alan J. Terry
Statistical Tables Compiled by Alan J. Terry School of Science and Sport University of the West of Scotland Paisley, Scotland Contents Table 1: Cumulative binomial probabilities Page 1 Table 2: Cumulative
More informationChapter 5 Univariate time-series analysis. () Chapter 5 Univariate time-series analysis 1 / 29
Chapter 5 Univariate time-series analysis () Chapter 5 Univariate time-series analysis 1 / 29 Time-Series Time-series is a sequence fx 1, x 2,..., x T g or fx t g, t = 1,..., T, where t is an index denoting
More informationFirm Heterogeneity and Credit Risk Diversification
Firm Heterogeneity and Credit Risk Diversification Samuel G. Hanson* M. Hashem Pesaran Harvard Business School University of Cambridge and USC Til Schuermann* Federal Reserve Bank of New York and Wharton
More informationIntroduction to Statistics I
Introduction to Statistics I Keio University, Faculty of Economics Continuous random variables Simon Clinet (Keio University) Intro to Stats November 1, 2018 1 / 18 Definition (Continuous random variable)
More informationIntroduction to Probability Theory and Stochastic Processes for Finance Lecture Notes
Introduction to Probability Theory and Stochastic Processes for Finance Lecture Notes Fabio Trojani Department of Economics, University of St. Gallen, Switzerland Correspondence address: Fabio Trojani,
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 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 informationProbability Weighted Moments. Andrew Smith
Probability Weighted Moments Andrew Smith andrewdsmith8@deloitte.co.uk 28 November 2014 Introduction If I asked you to summarise a data set, or fit a distribution You d probably calculate the mean and
More 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 informationStochastic Volatility (SV) Models
1 Motivations Stochastic Volatility (SV) Models Jun Yu Some stylised facts about financial asset return distributions: 1. Distribution is leptokurtic 2. Volatility clustering 3. Volatility responds to
More informationExtend the ideas of Kan and Zhou paper on Optimal Portfolio Construction under parameter uncertainty
Extend the ideas of Kan and Zhou paper on Optimal Portfolio Construction under parameter uncertainty George Photiou Lincoln College University of Oxford A dissertation submitted in partial fulfilment for
More informationStrategies for Improving the Efficiency of Monte-Carlo Methods
Strategies for Improving the Efficiency of Monte-Carlo Methods Paul J. Atzberger General comments or corrections should be sent to: paulatz@cims.nyu.edu Introduction The Monte-Carlo method is a useful
More 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 informationExam M Fall 2005 PRELIMINARY ANSWER KEY
Exam M Fall 005 PRELIMINARY ANSWER KEY Question # Answer Question # Answer 1 C 1 E C B 3 C 3 E 4 D 4 E 5 C 5 C 6 B 6 E 7 A 7 E 8 D 8 D 9 B 9 A 10 A 30 D 11 A 31 A 1 A 3 A 13 D 33 B 14 C 34 C 15 A 35 A
More informationQualifying Exam Solutions: Theoretical Statistics
Qualifying Exam Solutions: Theoretical Statistics. (a) For the first sampling plan, the expectation of any statistic W (X, X,..., X n ) is a polynomial of θ of degree less than n +. Hence τ(θ) cannot have
More informationDependence Structure and Extreme Comovements in International Equity and Bond Markets
Dependence Structure and Extreme Comovements in International Equity and Bond Markets René Garcia Edhec Business School, Université de Montréal, CIRANO and CIREQ Georges Tsafack Suffolk University Measuring
More informationThe Bernoulli distribution
This work is licensed under a Creative Commons Attribution-NonCommercial-ShareAlike License. Your use of this material constitutes acceptance of that license and the conditions of use of materials on this
More information12 The Bootstrap and why it works
12 he Bootstrap and why it works For a review of many applications of bootstrap see Efron and ibshirani (1994). For the theory behind the bootstrap see the books by Hall (1992), van der Waart (2000), Lahiri
More 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 informationADVANCED OPERATIONAL RISK MODELLING IN BANKS AND INSURANCE COMPANIES
Small business banking and financing: a global perspective Cagliari, 25-26 May 2007 ADVANCED OPERATIONAL RISK MODELLING IN BANKS AND INSURANCE COMPANIES C. Angela, R. Bisignani, G. Masala, M. Micocci 1
More informationThe Binomial Lattice Model for Stocks: Introduction to Option Pricing
1/27 The Binomial Lattice Model for Stocks: Introduction to Option Pricing Professor Karl Sigman Columbia University Dept. IEOR New York City USA 2/27 Outline The Binomial Lattice Model (BLM) as a Model
More informationECON 214 Elements of Statistics for Economists 2016/2017
ECON 214 Elements of Statistics for Economists 2016/2017 Topic The Normal Distribution Lecturer: Dr. Bernardin Senadza, Dept. of Economics bsenadza@ug.edu.gh College of Education School of Continuing and
More informationmay be of interest. That is, the average difference between the estimator and the truth. Estimators with Bias(ˆθ) = 0 are called unbiased.
1 Evaluating estimators Suppose you observe data X 1,..., X n that are iid observations with distribution F θ indexed by some parameter θ. When trying to estimate θ, one may be interested in determining
More informationDiploma in Business Administration Part 2. Quantitative Methods. Examiner s Suggested Answers
Cumulative frequency Diploma in Business Administration Part Quantitative Methods Examiner s Suggested Answers Question 1 Cumulative Frequency Curve 1 9 8 7 6 5 4 3 1 5 1 15 5 3 35 4 45 Weeks 1 (b) x f
More informationMonte Carlo Simulations in the Teaching Process
Monte Carlo Simulations in the Teaching Process Blanka Šedivá Department of Mathematics, Faculty of Applied Sciences University of West Bohemia, Plzeň, Czech Republic CADGME 2018 Conference on Digital
More informationدرس هفتم یادگیري ماشین. (Machine Learning) دانشگاه فردوسی مشهد دانشکده مهندسی رضا منصفی
یادگیري ماشین توزیع هاي نمونه و تخمین نقطه اي پارامترها Sampling Distributions and Point Estimation of Parameter (Machine Learning) دانشگاه فردوسی مشهد دانشکده مهندسی رضا منصفی درس هفتم 1 Outline Introduction
More informationValuation of Forward Starting CDOs
Valuation of Forward Starting CDOs Ken Jackson Wanhe Zhang February 10, 2007 Abstract A forward starting CDO is a single tranche CDO with a specified premium starting at a specified future time. Pricing
More informationLecture 1 of 4-part series. Spring School on Risk Management, Insurance and Finance European University at St. Petersburg, Russia.
Principles and Lecture 1 of 4-part series Spring School on Risk, Insurance and Finance European University at St. Petersburg, Russia 2-4 April 2012 s University of Connecticut, USA page 1 s Outline 1 2
More informationWhat was in the last lecture?
What was in the last lecture? Normal distribution A continuous rv with bell-shaped density curve The pdf is given by f(x) = 1 2πσ e (x µ)2 2σ 2, < x < If X N(µ, σ 2 ), E(X) = µ and V (X) = σ 2 Standard
More informationFrequency of Price Adjustment and Pass-through
Frequency of Price Adjustment and Pass-through Gita Gopinath Harvard and NBER Oleg Itskhoki Harvard CEFIR/NES March 11, 2009 1 / 39 Motivation Micro-level studies document significant heterogeneity in
More informationPaper Series of Risk Management in Financial Institutions
- December, 007 Paper Series of Risk Management in Financial Institutions The Effect of the Choice of the Loss Severity Distribution and the Parameter Estimation Method on Operational Risk Measurement*
More informationDependence Modeling and Credit Risk
Dependence Modeling and Credit Risk Paola Mosconi Banca IMI Bocconi University, 20/04/2015 and 27/04/2015 Paola Mosconi Lecture 6 1 / 112 Disclaimer The opinion expressed here are solely those of the author
More informationProbability Theory and Simulation Methods. April 9th, Lecture 20: Special distributions
April 9th, 2018 Lecture 20: Special distributions Week 1 Chapter 1: Axioms of probability Week 2 Chapter 3: Conditional probability and independence Week 4 Chapters 4, 6: Random variables Week 9 Chapter
More information