Experience Rating in General Insurance by Credibility Estimation
|
|
- Brian Underwood
- 5 years ago
- Views:
Transcription
1 Experience Rating in General Insurance by Credibility Estimation Xian Zhou Department of Applied Finance and Actuarial Studies Macquarie University, Sydney, Australia Abstract This work presents a new credibility estimation of the probability distributions of risks under Bayes settings in a nonparametric framework without any assumption on the prior distribution. The method is then applied to general insurance premium pricing, a procedure commonly known as experience rating, which utilizes the insured s past claim experience to calculate the premium under a given premium principle (also referred to as a risk measure). By estimating the entire probability distributions of losses, our method provides a unified nonparametric framework to experience rating under arbitrary premium principles. This encompasses the advantages of the well-known approaches introduced by Bühlmann and Ferguson, while overcomes their drawbacks. This is a joint work with Prof. Xiaoqiang Cai at the Chinese University of Hong Kong, Dr. Limin Wen at Jiangxi Normal University, and Prof Xianyi Wu at East China Normal University. 1
2 1 Introduction Pricing or measuring risks is one of the central tasks of financial enterprises and regulators in risk management. Numerous risk measures have been proposed for this purpose, such as value-at-risk, conditional value-at-risk, coherent risk measures, and distortion risk measures. In compensating insured losses in property, business, employment, health, etc., pricing risks is carried out by premium principle, which is considered as the translation of risk measures into the general insurance context. In this talk, the terms risk measure and premium principle are alternately used to indicate a rule H(X) that assigns a premium to a risk X in terms of a functional of its distribution function. In practice, however, to determine an appropriate premium principle is only the first step in the process of risk pricing, because the involved distributions are not known so that H(X) can only be estimated by using the insured s claim experience with certain effective statistical methods. These procedures form the so-called experience rating. We are concerned with solutions to this problem under the Bayesian setting: the distribution of a risk X is identified by an unknown and unobservable parameter (vector) θ, via S(x, θ) = Pr(X > x θ). 2
3 The unobservable parameter θ is modeled by a random vector. Its distribution π(θ) is referred to as the prior distribution in statistics or structural function in actuarial science. The majority of the literature has mainly focussed on fully specified S(x, θ) given θ and known (or assumed) π(θ). Thus the problem can be solved with the standard parametric Bayesian methodology through posterior update. Although the applications of the parametric Bayesian methodology have been extensively investigated for general insurance, the reality is that π(θ) is rarely known, making these applications impractical. To deal with such situations, one approach is to allow some unknown parameters in π(θ), but maintain some assumed mathematical forms of S(x, θ) and π(θ). Then the risk pricing can be carried out by the empirical Bayes analysis. In many insurance practices, however, the mathematical forms of S(x, θ) and π(θ) are difficult to specify or unavailable. Under such circumstances, the analysis can only be done on a distribution-free or nonparametric basis. The best-known solution in the actuarial community is probably the Bühlmann s approach to net premium principle (known as the credibility theory). 3
4 Under the net premium principle, the credibility premium of a future claim X n+1 can be expressed as a weighted average between the empirical individual mean (indicating the information delivered by historical data of the risk itself) and the collective mean (indicating the aggregated information obtained from all possible insureds): Ĥ n+1 = z X + (1 z)µ 0, (1.1) where X = n 1 n X i is the empirical mean of the claims data {X 1, X 2,..., X n }, µ 0 = E[X] is known as the collective mean of the risk, and z is the credibility factor. This (linear) weighted average has been so well received today that it is almost regarded as a synonym of the credibility theory. However, it is not easy to directly transplant Bühlmann s method to other premium principles. Consequently most contributions in this area have been devoted to the net premium principle. Another important solution stream to this problem, originated from Ferguson (1973) 1, is the Bayesian nonparametrics. Consider S(x, θ) as a Dirichlet process with time x. Its distribution is determined by a finite Borel measure α( ) on the real line R. After losses X 1,..., X n are observed, the posterior distribution of θ is also a Dirichlet process but with α being updated to α+ n δ X i, where δ Xi is a probability measure degenerated at X i. 1 Ferguson, T. A Bayesian analysis of some non-parametric problems. Annals of Statistics, 1,
5 Consequently, S(x, θ) can be estimated by the posterior mean Ŝ(x, θ) = α(r) n α([x, )) + α(r) + n where S n (x) = n 1 n I(X i > x). α(r) + n S n(x), (1.2) Then the experience rating can be carried out by inserting the estimated distribution of the loss X into any premium principle. While the approach of Ferguson (1973) provides a unified solution for all premium principles, it requires the assumption of Dirichlet prior. This assumption is mainly for mathematical convenience, but hardly justifiable in practice. Without the Dirichlet assumption, the estimator in (1.2) would lose its justification, and the corresponding risk premium H(X n+1 ) would not be credible when the true prior deviates far from the Dirichlet. In this talk, we will introduce a new distribution-free approach for experience rating under arbitrary premium principles without any assumption on the prior distribution. It has two direct advantages: (i) truly distribution-free settings as in the Bühlman s credibility theory, and (ii) generating experience rating for arbitrary premium principles as achieved by Ferguson s Bayesian nonparametric method. 5
6 The main idea is to first derive an estimate of S(x, θ) by minimizing the L 2 -distance, and then embed its estimate, say Ŝ(x, θ), into the premium principle functional to obtain an estimate of the corresponding premium, i.e., H(X, θ) = H( Ŝ(x, θ)). Specifically, we consider the following two models: Model I: The first two moments of S(x, θ) with respect to the prior can be specified. In this case, the data are from X 1, X 2,..., which are (conditionally) i.i.d. copies of X given θ. This model brings insight into and motivates the estimation methods of the more practically meaningful model II next. Model II: This is a general model with totally unknown S(x, θ). Suppose that we have a portfolio of risks X i, i = 1, 2,..., K, where each X i is identified with a parameter (vector) θ i and contributes a series of data X i1, X i2,...x ini,..., which are i.i.d. copies of X i, given θ i. The purpose is to estimate the distributions and then assign proper premiums to the future claims X i,ni +1, i = 1, 2,..., K under a given premium principle. In order to overcome the restriction of Dirichlet priors, we develop a new approach in a quite different way from Ferguson (1973). It turns out, however, when the prior distribution is Dirichlet, our estimation happens to coincide with that of Ferguson (1973). This shows that our results actually generalize those of Ferguson (1973). 6
7 2 The Linear Bayes Estimation of Distributions We consider Model I in this section. Write S 0 (x) = E θ [S(x, θ)], τ 2 0 (x) = Var θ [S(x, θ)] (2.1) and σ0 2 (x) = E [ ( θ Var I(X>x) θ )], (2.2) where E θ and Var θ indicate that the expectations are computed with respect to the distribution of θ. We also denote µ(θ) = E[X θ], σ 2 = E θ [Var(X θ)], τ 2 = Var θ (E(X θ)). The performance measure to be optimized As S(x, θ) = E [I(X > x) θ], a pointwise credibility estimator of S(x, θ) for a fixed x following Bühlmann s classical credibility theory can be given by S(x, θ) = Z(x)S n (x) + (1 Z(x)) S 0 (x), (2.3) where S n (x) = n 1 n I(X i > x) is the empirical estimate of S(x, θ) based on the data X 1, X 2,..., X n, and Z(x) = nτ 2 0(x) σ 2 0 (x) + nτ 2 0 (x) is the credibility factor. By (2.1) (2.2), σ 2 0 (x) = S 0(x) E θ [ S 2 (x, θ) ], τ0 2 (x) = E [ θ S 2 (x, θ) ] S0 2 (x). (2.4) 7
8 The naive estimate in (2.3) can be more precisely expressed as S (x, θ) = n(e θ[s 2 (x, θ)] S 2 0 (x))s n(x) + (S 0 (x) E θ [S 2 (x, θ)])s 0 (x) n(e θ [S 2 (x, θ)] S 2 0(x)) + (S 0 (x) E θ [S 2 (x, θ)]) = n+1 (n i + 1)E θ [S 2 (x, θ)] (n i)s0 2(x) S 0(x)E θ [S 2 (x, θ)] I (n 1)E θ [S 2 (x, θ)] ns0(x) 2 [X(i 1),X + S 0 (x) (i) )(x), (2.5) where = X (0) < X (1) X (2) X (n) < X (n+1) = are the order statistics of the sample. This naive credibility estimate (2.3) (or (2.5)) would require: checking monotonicity of S(x, θ) in x, making an intuitive interpretation for S(x, θ), and computing a premium under S(x, θ). These tasks are however not feasible due to the dependence of the credibility factor Z(x) on x. A natural remedy is to consider Z(x) = Z to be independent of x. More specifically, we seek a function of the form n β 0 (x) + λ i I(X i > x), (2.6) which minimizes the square of the L 2 -distance: min E β 0 (x),λ 1,...,λ n + ( S(x,θ) β 0 (x) 2 n λ i I(X i > x)) dx, (2.7) where β 0 (x) is a nonincreasing function independent of the claims history and λ 1,..., λ n are real-valued decision variables. 8
9 Estimation of distribution functions The following theorem gives an estimator of S(x, θ) by solving the optimization problem (2.7). Theorem 2.1 The credibility estimator of S(x, θ), which minimizes (2.7), is given by Ŝ(x, θ) = ZS n (x) + (1 Z) S 0 (x), (2.8) where with τ 2 0 = τ 2 0 (x)dx Z = nτ 2 0 nτ σ2 0 and σ2 0 = σ 2 0 (x)dx. (2.9) Moreover, the integrated mean-squared loss for the optimal estimator given by (2.8) is + E ) 2 τ (Ŝ(x, 2 θ) S(x, θ) dx = 0 σ2 0 σ0 2 + nτ 0 2. (2.10) Note that Theorem 2.1 does not require to specify the prior distribution. If the prior is not Dirichlet, then Ferguson s method cannot be applied, but the estimator defined in (2.8) and (2.9) remains valid. In this sense, our new estimation provides a more general solution for the estimation of the distributions than that of Ferguson (1973). 9
10 Properties of estimator If S(x, θ) is a Dirichlet process with a finite Borel measure α( ) on R, then with density S(x, θ) Beta(α(x), α(r) α(x)) Γ(α(R)) Γ(α(x))Γ(α(R) α(x)) tα(x) 1 (1 t) α(r) α(x) 1, where α(x) = α((x, )). It is easy to show that in such a case, the credibility factor Z in (2.9) becomes Z = nτ 0 2 nτ0 2 + = α(r) σ2 0 α(r) + n. Thus the estimate in (2.8) reduces to Ferguson s estimate (1.2). The estimator in (2.8) is obviously a survival distribution function. By (2.9), 0 Z 1, Z 1 as n and Z 0 as n 0, (2.11) which allows the traditional credibility interpretation: the more data, the more credible the empirical distribution S n (x). The estimator (2.8) is uniformly strongly consistent in the sense that lim sup S(x, θ) Ŝ(x, θ) = 0 n x almost surely (a.s.). 10
11 3 The Empirical Bayes Estimation We now consider Model II, in which both the conditional distribution S(x, θ) and the structure function π(θ) are completely unknown, including the structure parameters S 0 (x), τ 2 0 and σ2 0. We propose an empirical Bayes method to estimate S 0 (x), τ 2 0 and σ 2 0 based on the claim experiences over a number of risks in the same portfolio. Let X 1,..., X K denote K risks under observation. The distribution of each X i is characterized by a risk parameter θ i and contributes a sequence of claim experiences X i = (X i1,...x ini ) over n i time periods, subject to the following usual assumptions. Assumption 3.1 (i) Conditional on θ i, the random variables X ij (j = 1, 2,..., n i ) are i.i.d. with common unknown distribution and moments: ( ) S(x, θ i ) = Pr (X ij > x θ i ), σ 2 (x, θ i ) = Var I ( X ij >x) θ i, i = 1, 2,..., K, j = 1, 2,..., n i. (ii) The random variables θ 1,..., θ K are i.i.d. with a common prior distribution π(θ). Write S 0 (x) = E θ [S(x, θ i )], τ 2 0 (x) = Var θ [S(x, θ i )] and σ 2 0(x) = E θ [ σ 2 (x, θ i ) ]. 11
12 Homogeneous estimation of the distribution Write S i (x) = 1 n i n i j=1 I (X ij > x), i = 1, 2,..., K. (3.1) To obtain the credibility estimator of S(x, θ i ), consider the class of distribution functions { K n s L = a st I (X st > x),a st R, s=1 t=1 K s=1 n s t=1 a st = 1 }, (3.2) Solve the minimization problem { + } min E [g S(x, θ i )] 2 dx. (3.3) g L The solution is stated in the next theorem. Theorem 3.1 The homogeneous credibility estimator of S(x, θ i ) as the solution to (3.3) is S(x, θ i ) = Z i S i (x) + (1 Z i ) S(x),, i = 1,..., K, (3.4) where Z i = n iτ 2 0 σ n iτ 2 0 and S(x) = 1 K Z i K Z i S i (x). (3.5) We can also derive an inhomogeneous estimator of the form K n s β 0 (x) + a st I (X st > x). s=1 t=1 The result, however, is the same as that presented by (2.8) (2.9). This is a typical phenomenon in credibility theory. 12
13 Estimation of structure parameters We now proceed to the estimation of the structure parameters τ 2 0 and σ 2 0, which needs the estimation of S 0 (x). At this moment, however, S(x) is not a suitable estimate of S 0 (x) because Z i s involve τ 2 0 and σ2 0 (cf. (3.5)). Moreover, if we use an estimator of the form S(x) = K w is i (x) by minimizing (S(x) S 0 (x)) 2 dx, it leads to S(x) = S(x). Therefore, we simply select S(x) = 1 N K n i S i (x), where N = K n i, (3.6) which is unbiased and strongly consistent in the sense that S(x) converges to S 0 (x) almost surely as K, under the condition n 2 K <. (3.7) (N K) 2 K=1 In view of the definitions of σ 2 0(x) and τ 2 0(x), the parameters σ 2 0 and τ 2 0 and can be estimated using SSE(x) = K SSA(x) = n i j=1 [I (X ij > x) S i (x)] 2 K 2 n i [S i (x) S(x)]. 13
14 The estimators can be formally defined as and τ 2 0 = N N 2 K n2 i σ 2 0 = 1 N K [ SSA(x)dx K 1 N K By inserting the estimates of τ 2 0 and σ2 0 SSE(x)dx (3.8) ] SSE (x) dx. (3.9) into (3.4) and (3.5), we can get empirical Bayes estimators for the distribution S(x, θ i ) as ( S(x, Ẑi) θ i ) = ẐiS i (x) + 1 Ŝ0 (x), i = 1,..., K, (3.10) where and Ẑ i = n iˆτ 0 2 ˆσ n iˆτ 0 2, i = 1,..., K, (3.11) Ŝ 0 (x) = 1 K Ẑi K Ẑ i S i (x). (3.12) The properties of the estimators in (3.8) (3.9) are provided in the following theorem. Theorem 3.2 σ 2 0 and τ 2 0 have the following properties. (1) σ 2 0 and τ 2 0 are unbiased estimators of σ2 0 and τ 2 0, respectively. (2) Under condition (3.7), σ2 0 σ 2 0 and τ 2 0 τ 2 0 almost surely as K. 14
15 Asymptotic optimality For simplicity, we discuss the balanced case with n i n, thus Z i and Ẑi are all the same. Consequently, Ŝ 0 (x) = 1 K K S i (x). An additional condition is needed for the asymptotic optimality needs the following extra: there is a constant δ > 0 such that + ( σ 2 0 (x) n ) δ/(2+δ) + τ 0 2 (x) dx <. (3.13) Write and Q(x, θ i ) = Q(x, θ i ) = + + ( S(x, θi ) S(x, θ i )) 2 dx ( S(x, θi ) S(x, θ i )) 2 dx. The asymptotic optimality is stated in the theorem below. Theorem 3.3 Under conditions (3.7) and (3.13), the S(x, θ i ) in (3.10) (3.12) are asymptotically optimal in the sense that [ lim max E Q(x, θi ) Q(x, θ i )] = 0, (3.14) K 1 i K where S(x, θ i ) is the linear estimator defined by (2.8) and (2.9). 15
16 4 Applications to Experience Ratemaking Having estimated the distributions, the application to the experience ratemaking is apparent. For a given premium principle H, write H (X θ) for the individual premium of a risk X with distribution S(x, θ), referred also to as the risk premium of X. By replacing S(x, θ) with Ŝ(x, θ) or S(x, θ i ), an estimator H(X θ) of the individual premium H(X θ) is obtained. We consider the experience premium by inserting Ŝ(x, θ) for S(x, θ) and compare it with existing methods. The version using S(x, θ i ) can be discussed in a similar way. The consistency of H (X θ) is presented in the theorem below, which follows directly from Theorems 2.1 and 3.2. Theorem 4.1 If the premium principle H is continuous with respect to the L -norm of distribution functions, then H(X θ) is strongly consistent for H(X θ). According to Theorem 4.1, the strong consistency of H(X θ) for H(X θ) is guaranteed as long as the premium principle H satisfies the condition of the theorem. 16
17 This is in contrast to existing results in the literature, where either the credibility estimator is not generally (strongly) consistent, or the consistency needs to be proved separately in each case. There are a number of well known premium principles for which strong consistency of the experience ratemaking can be obtained directly or by checking the condition of Theorem 4.1, including the net premium H(X θ) = E[X θ], variance premium H(X θ) = E[X θ] + αvarh(x θ), modified variance premium standard deviation premium Esscher premiums H(X θ) = E[X θ] + α Var(X θ), E[X θ] H(X θ) = E[X θ] + α Var(X θ), H(X θ) = E[XehX θ], E[e hx θ] conditional tail expectation premium H(X θ) = E[X X > α, θ], exponential premium H(X θ) = α 1 log [ E ( e αx θ )], and Kamp s premium H(X θ) = E[X(1 eαx ) θ]. E[(1 e αx ) θ] 17
18 The following are two exceptions where the proof of the consistency require additional work and conditions. 1. Dutch s premium principle: H(X θ) = E[X θ] + ηe [ (X αe[x θ]) + θ ], where α 1 and 0 < η 1. Observe that [ ] 1 n E Sn (X αe[x θ])+ = (X i α µ(θ)) +. n The estimate of H(X θ) is then given by H(X θ) = µ(θ) + η { Z n n ( X i α µ(θ) ) } ( + (1 Z) x α µ(θ) ) df 0 (x), + + where µ(θ) = ZX + (1 Z)µ Distortion premium principle: H(X θ) = The premium estimate is H(X θ) = 0 0 g (S(x, θ)) dx g (ZS n (x) + (1 Z)S 0 (x))dx 0 0 where the credibility factor is given by (2.9). g(1 S (x, θ))dx. g (ZF n (x) + (1 Z)F 0 (x))dx, The consistency in this case can be proved under certain additional regularity condition on g, such as the Lipschitz condition of g on [0, 1] if X has bounded support, or of g on [0, 1] if the support is unbounded. 18
19 5 Performance Evaluation We evaluate the performance of the proposed approach by comparing it with some well-known representative methods in the literature. Evaluation under the net premium principle Under the net premium principle, the optimal premium is given by the Bühlmann s credibility premium µ(θ) c = Z c X n + (1 Z c ) µ 0, where the superscript c stands for classical. The premium using our method is µ(θ) = H(X θ) = ZX + (1 Z) µ 0. (5.1) The expected quadratic losses of µ(θ) approaches that of µ(θ) c : lim n E[ µ(θ) µ(θ)] 2 E[ µ(θ) c µ(θ)] 2 = 1. (5.2) This shows that µ(θ) is asymptotically equivalent to the optimum. Assume S(x, θ) = e θx and θ Gamma(α, β). Then the expected squared losses of µ(θ), µ(θ) c and µ 0 (the collective premium) can be expressed by β 2 V cu (α), β 2 V c (α) and β 2 V col respectively. The curves of V cu (α), V c (α) and V col (α) as functions of α for n = 10 are plotted in Figure 1 below, which shows that µ(θ) is close to the optimal µ(θ) c. 19
20 Figure 1. The curves of V cu (α), V c (α) and V col (α) as functions of α. 12 V col V cu V c α 20
21 Evaluation under the variance premium principle Consider the Poisson-exponential model, where S(x, θ) is Poisson with mean θ and π(θ) is exponential with rate λ. Bühlmann suggested a credibility premium formula that involves the forth moments of X i given θ. In comparison, we only need the first two moments of the risk distribution. Let V & denote the expected squared loss with & = B, cu, Bul, col corresponding to Bayes, proposed, Bühlmann, collective premiums respectively under the variance premium principle. Figures 2 below shows the numerical results of V & for λ = 0.2 to 1.0 with sample size n = 100. Figure 2. The curves of V B,V cu and V Bul as functions of λ with n = V B V cu V Bul λ 21
22 Evaluation under the Esscher premium principle Consider the Poisson-Gamma model with S(x, θ) P oisson(θ) and θ Gamma(α, β) with α ranging from 2 to 8 and β = 6. Let V & denote the expected squared loss with & = B, cu, col, P under the Esscher premium principle, where P denotes the method introduced in Pan et.al (2008) 2. Figure 3 below shows the numerical results of V & for different values of α and sample size 100. Figure 3. Curves of V & for different values of α and sample size V B V cu 3 V P V col α 2 Pan, M., Wang, R. and Wu, X. On the consistency of credibility premiums regarding Esscher principle. Insurance: Mathematics and Economics, 42, ,
The Credibility Estimators with Dependence Over Risks
Applied Mathematical Sciences, Vol. 8, 2014, no. 161, 8045-8050 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2014.410803 The Credibility Estimators with Dependence Over Risks Qiang Zhang
More informationBayesian Point Process Modeling for Extreme Value Analysis, with an Application to Systemic Risk Assessment in Correlated Financial Markets
Bayesian Point Process Modeling for Extreme Value Analysis, with an Application to Systemic Risk Assessment in Correlated Financial Markets Athanasios Kottas Department of Applied Mathematics and Statistics,
More informationA few basics of credibility theory
A few basics of credibility theory Greg Taylor Director, Taylor Fry Consulting Actuaries Professorial Associate, University of Melbourne Adjunct Professor, University of New South Wales General credibility
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 informationCourse 4 Solutions November 2001 Exams
Course 4 Solutions November 001 Exams November, 001 Society of Actuaries Question #1 From the Yule-Walker equations: ρ φ + ρφ 1 1 1. 1 1+ ρ ρφ φ Substituting the given quantities yields: 0.53 φ + 0.53φ
More informationNonlife Actuarial Models. Chapter 4 Risk Measures
Nonlife Actuarial Models Chapter 4 Risk Measures Learning Objectives 1. Risk measures based on premium principles 2. Risk measures based on capital requirements 3. Value-at-Risk and conditional tail expectation
More informationParameter Estimation
Parameter Estimation Chapters 13-15 Stat 477 - Loss Models Chapters 13-15 (Stat 477) Parameter Estimation Brian Hartman - BYU 1 / 23 Methods for parameter estimation Methods for parameter estimation Methods
More informationREGRESSION TREE CREDIBILITY MODEL
LIQUN DIAO AND CHENGGUO WENG Department of Statistics and Actuarial Science, University of Waterloo Advances in Predictive Analytics Conference, Waterloo, Ontario Dec 1, 2017 Overview Statistical }{{ Method
More informationNon-Life Insurance: Mathematics and Statistics
ETH Zürich, D-MATH HS 08 Prof. Dr. Mario V. Wüthrich Coordinator Andrea Gabrielli Non-Life Insurance: Mathematics and Statistics Solution sheet 9 Solution 9. Utility Indifference Price (a) Suppose that
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 informationPractice Exam 1. (A) (B) (C) (D) (E) You are given the following data on loss sizes:
Practice Exam 1 1. Losses for an insurance coverage have the following cumulative distribution function: F(0) = 0 F(1,000) = 0.2 F(5,000) = 0.4 F(10,000) = 0.9 F(100,000) = 1 with linear interpolation
More informationMarkov Chain Monte Carlo (MCMC)
Markov Chain Monte Carlo (MCMC Dependent Sampling Suppose we wish to sample from a density π, and we can evaluate π as a function but have no means to directly generate a sample. Rejection sampling can
More informationFoundations of Nonparametric Bayesian Methods
1 / 27 Foundations of Nonparametric Bayesian Methods Part II: Models on the Simplex Peter Orbanz http://mlg.eng.cam.ac.uk/porbanz/npb-tutorial.html 2 / 27 Tutorial Overview Part I: Basics Part II: Models
More informationQualifying Exam in Probability and Statistics. https://www.soa.org/files/edu/edu-exam-p-sample-quest.pdf
Part 1: Sample Problems for the Elementary Section of Qualifying Exam in Probability and Statistics https://www.soa.org/files/edu/edu-exam-p-sample-quest.pdf Part 2: Sample Problems for the Advanced Section
More informationApplications of axiomatic capital allocation and generalized weighted allocation
6 2010 11 ( ) Journal of East China Normal University (Natural Science) No. 6 Nov. 2010 Article ID: 1000-5641(2010)06-0146-10 Applications of axiomatic capital allocation and generalized weighted allocation
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 informationStandard Error of Technical Cost Incorporating Parameter Uncertainty
Standard Error of Technical Cost Incorporating Parameter Uncertainty Christopher Morton Insurance Australia Group This presentation has been prepared for the Actuaries Institute 2012 General Insurance
More informationEcon 2148, spring 2019 Statistical decision theory
Econ 2148, spring 2019 Statistical decision theory Maximilian Kasy Department of Economics, Harvard University 1 / 53 Takeaways for this part of class 1. A general framework to think about what makes a
More informationNon-Life Insurance: Mathematics and Statistics
Exercise sheet 1 Exercise 1.1 Discrete Distribution Suppose the random variable N follows a geometric distribution with parameter p œ (0, 1), i.e. ; (1 p) P[N = k] = k 1 p if k œ N \{0}, 0 else. (a) Show
More informationBayesian estimation of the discrepancy with misspecified parametric models
Bayesian estimation of the discrepancy with misspecified parametric models Pierpaolo De Blasi University of Torino & Collegio Carlo Alberto Bayesian Nonparametrics workshop ICERM, 17-21 September 2012
More informationMATH c UNIVERSITY OF LEEDS Examination for the Module MATH2715 (January 2015) STATISTICAL METHODS. Time allowed: 2 hours
MATH2750 This question paper consists of 8 printed pages, each of which is identified by the reference MATH275. All calculators must carry an approval sticker issued by the School of Mathematics. c UNIVERSITY
More informationStatistics: Learning models from data
DS-GA 1002 Lecture notes 5 October 19, 2015 Statistics: Learning models from data Learning models from data that are assumed to be generated probabilistically from a certain unknown distribution is a crucial
More informationQualifying Exam in Probability and Statistics. https://www.soa.org/files/edu/edu-exam-p-sample-quest.pdf
Part : Sample Problems for the Elementary Section of Qualifying Exam in Probability and Statistics https://www.soa.org/files/edu/edu-exam-p-sample-quest.pdf Part 2: Sample Problems for the Advanced Section
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 informationSTA 732: Inference. Notes 10. Parameter Estimation from a Decision Theoretic Angle. Other resources
STA 732: Inference Notes 10. Parameter Estimation from a Decision Theoretic Angle Other resources 1 Statistical rules, loss and risk We saw that a major focus of classical statistics is comparing various
More informationTail Conditional Expectations for. Extended Exponential Dispersion Models
Tail Conditional Expectations for Extended Exponential Dispersion Models A Thesis Presented to the Faculty of the Department of Mathematical Sciences Middle Tennessee State University In Partial Fulfillment
More informationEfficient rare-event simulation for sums of dependent random varia
Efficient rare-event simulation for sums of dependent random variables Leonardo Rojas-Nandayapa joint work with José Blanchet February 13, 2012 MCQMC UNSW, Sydney, Australia Contents Introduction 1 Introduction
More informationMixture distributions in Exams MLC/3L and C/4
Making sense of... Mixture distributions in Exams MLC/3L and C/4 James W. Daniel Jim Daniel s Actuarial Seminars www.actuarialseminars.com February 1, 2012 c Copyright 2012 by James W. Daniel; reproduction
More informationClaims Reserving under Solvency II
Claims Reserving under Solvency II Mario V. Wüthrich RiskLab, ETH Zurich Swiss Finance Institute Professor joint work with Michael Merz (University of Hamburg) April 21, 2016 European Congress of Actuaries,
More informationRMSC 2001 Introduction to Risk Management
RMSC 2001 Introduction to Risk Management Tutorial 4 (2011/12) 1 February 20, 2012 Outline: 1. Failure Time 2. Loss Frequency 3. Loss Severity 4. Aggregate Claim ====================================================
More informationLecture 8: Information Theory and Statistics
Lecture 8: Information Theory and Statistics Part II: Hypothesis Testing and I-Hsiang Wang Department of Electrical Engineering National Taiwan University ihwang@ntu.edu.tw December 23, 2015 1 / 50 I-Hsiang
More informationExpectation Propagation Algorithm
Expectation Propagation Algorithm 1 Shuang Wang School of Electrical and Computer Engineering University of Oklahoma, Tulsa, OK, 74135 Email: {shuangwang}@ou.edu This note contains three parts. First,
More informationSolutions of the Financial Risk Management Examination
Solutions of the Financial Risk Management Examination Thierry Roncalli January 9 th 03 Remark The first five questions are corrected in TR-GDR and in the document of exercise solutions, which is available
More informationEcon 2140, spring 2018, Part IIa Statistical Decision Theory
Econ 2140, spring 2018, Part IIa Maximilian Kasy Department of Economics, Harvard University 1 / 35 Examples of decision problems Decide whether or not the hypothesis of no racial discrimination in job
More informationCOS513 LECTURE 8 STATISTICAL CONCEPTS
COS513 LECTURE 8 STATISTICAL CONCEPTS NIKOLAI SLAVOV AND ANKUR PARIKH 1. MAKING MEANINGFUL STATEMENTS FROM JOINT PROBABILITY DISTRIBUTIONS. A graphical model (GM) represents a family of probability distributions
More informationExtreme Value Analysis and Spatial Extremes
Extreme Value Analysis and Department of Statistics Purdue University 11/07/2013 Outline Motivation 1 Motivation 2 Extreme Value Theorem and 3 Bayesian Hierarchical Models Copula Models Max-stable Models
More informationIntroduction to Self-normalized Limit Theory
Introduction to Self-normalized Limit Theory Qi-Man Shao The Chinese University of Hong Kong E-mail: qmshao@cuhk.edu.hk Outline What is the self-normalization? Why? Classical limit theorems Self-normalized
More informationReinsurance and ruin problem: asymptotics in the case of heavy-tailed claims
Reinsurance and ruin problem: asymptotics in the case of heavy-tailed claims Serguei Foss Heriot-Watt University, Edinburgh Karlovasi, Samos 3 June, 2010 (Joint work with Tomasz Rolski and Stan Zachary
More informationOn Backtesting Risk Measurement Models
On Backtesting Risk Measurement Models Hideatsu Tsukahara Department of Economics, Seijo University e-mail address: tsukahar@seijo.ac.jp 1 Introduction In general, the purpose of backtesting is twofold:
More informationLecture 7 Introduction to Statistical Decision Theory
Lecture 7 Introduction to Statistical Decision Theory I-Hsiang Wang Department of Electrical Engineering National Taiwan University ihwang@ntu.edu.tw December 20, 2016 1 / 55 I-Hsiang Wang IT Lecture 7
More informationAsymptotic distribution of the sample average value-at-risk
Asymptotic distribution of the sample average value-at-risk Stoyan V. Stoyanov Svetlozar T. Rachev September 3, 7 Abstract In this paper, we prove a result for the asymptotic distribution of the sample
More informationLosses Given Default in the Presence of Extreme Risks
Losses Given Default in the Presence of Extreme Risks Qihe Tang [a] and Zhongyi Yuan [b] [a] Department of Statistics and Actuarial Science University of Iowa [b] Smeal College of Business Pennsylvania
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 informationRMSC 2001 Introduction to Risk Management
RMSC 2001 Introduction to Risk Management Tutorial 4 (2011/12) 1 February 20, 2012 Outline: 1. Failure Time 2. Loss Frequency 3. Loss Severity 4. Aggregate Claim ====================================================
More informationUniversity Of Calgary Department of Mathematics and Statistics
University Of Calgary Department of Mathematics and Statistics Hawkes Seminar May 23, 2018 1 / 46 Some Problems in Insurance and Reinsurance Mohamed Badaoui Department of Electrical Engineering National
More informationExam C Solutions Spring 2005
Exam C Solutions Spring 005 Question # The CDF is F( x) = 4 ( + x) Observation (x) F(x) compare to: Maximum difference 0. 0.58 0, 0. 0.58 0.7 0.880 0., 0.4 0.680 0.9 0.93 0.4, 0.6 0.53. 0.949 0.6, 0.8
More informationEXPERIENCE-BASED LONGEVITY ASSESSMENT
EXPERIENCE-BASED LONGEVITY ASSESSMENT ERMANNO PITACCO Università di Trieste ermanno.pitacco@econ.units.it p. 1/55 Agenda Introduction Stochastic modeling: the process risk Uncertainty in future mortality
More informationCreating New Distributions
Creating New Distributions Section 5.2 Stat 477 - Loss Models Section 5.2 (Stat 477) Creating New Distributions Brian Hartman - BYU 1 / 18 Generating new distributions Some methods to generate new distributions
More informationChapter 4 HOMEWORK ASSIGNMENTS. 4.1 Homework #1
Chapter 4 HOMEWORK ASSIGNMENTS These homeworks may be modified as the semester progresses. It is your responsibility to keep up to date with the correctly assigned homeworks. There may be some errors in
More informationChapter 3: Unbiased Estimation Lecture 22: UMVUE and the method of using a sufficient and complete statistic
Chapter 3: Unbiased Estimation Lecture 22: UMVUE and the method of using a sufficient and complete statistic Unbiased estimation Unbiased or asymptotically unbiased estimation plays an important role in
More informationCalculation of Bayes Premium for Conditional Elliptical Risks
1 Calculation of Bayes Premium for Conditional Elliptical Risks Alfred Kume 1 and Enkelejd Hashorva University of Kent & University of Lausanne February 1, 13 Abstract: In this paper we discuss the calculation
More informationEco517 Fall 2004 C. Sims MIDTERM EXAM
Eco517 Fall 2004 C. Sims MIDTERM EXAM Answer all four questions. Each is worth 23 points. Do not devote disproportionate time to any one question unless you have answered all the others. (1) We are considering
More informationEmpirical Bayes Quantile-Prediction aka E-B Prediction under Check-loss;
BFF4, May 2, 2017 Empirical Bayes Quantile-Prediction aka E-B Prediction under Check-loss; Lawrence D. Brown Wharton School, Univ. of Pennsylvania Joint work with Gourab Mukherjee and Paat Rusmevichientong
More informationConjugate Predictive Distributions and Generalized Entropies
Conjugate Predictive Distributions and Generalized Entropies Eduardo Gutiérrez-Peña Department of Probability and Statistics IIMAS-UNAM, Mexico Padova, Italy. 21-23 March, 2013 Menu 1 Antipasto/Appetizer
More informationDecision principles derived from risk measures
Decision principles derived from risk measures Marc Goovaerts Marc.Goovaerts@econ.kuleuven.ac.be Katholieke Universiteit Leuven Decision principles derived from risk measures - Marc Goovaerts p. 1/17 Back
More informationA new approach for stochastic ordering of risks
A new approach for stochastic ordering of risks Liang Hong, PhD, FSA Department of Mathematics Robert Morris University Presented at 2014 Actuarial Research Conference UC Santa Barbara July 16, 2014 Liang
More informationTime Series and Dynamic Models
Time Series and Dynamic Models Section 1 Intro to Bayesian Inference Carlos M. Carvalho The University of Texas at Austin 1 Outline 1 1. Foundations of Bayesian Statistics 2. Bayesian Estimation 3. The
More informationQualifying Exam in Probability and Statistics.
Part 1: Sample Problems for the Elementary Section of Qualifying Exam in Probability and Statistics https://www.soa.org/files/edu/edu-exam-p-sample-quest.pdf Part 2: Sample Problems for the Advanced Section
More informationA Very Brief Summary of Statistical Inference, and Examples
A Very Brief Summary of Statistical Inference, and Examples Trinity Term 2008 Prof. Gesine Reinert 1 Data x = x 1, x 2,..., x n, realisations of random variables X 1, X 2,..., X n with distribution (model)
More informationGeneral approach to the optimal portfolio risk management
General approach to the optimal portfolio risk management Zinoviy Landsman, University of Haifa Joint talk with Udi Makov and Tomer Shushi This presentation has been prepared for the Actuaries Institute
More informationMarket Risk. MFM Practitioner Module: Quantitiative Risk Management. John Dodson. February 8, Market Risk. John Dodson.
MFM Practitioner Module: Quantitiative Risk Management February 8, 2017 This week s material ties together our discussion going back to the beginning of the fall term about risk measures based on the (single-period)
More informationNonparametric estimation of extreme risks from heavy-tailed distributions
Nonparametric estimation of extreme risks from heavy-tailed distributions Laurent GARDES joint work with Jonathan EL METHNI & Stéphane GIRARD December 2013 1 Introduction to risk measures 2 Introduction
More informationEmpirical Risk Minimization is an incomplete inductive principle Thomas P. Minka
Empirical Risk Minimization is an incomplete inductive principle Thomas P. Minka February 20, 2001 Abstract Empirical Risk Minimization (ERM) only utilizes the loss function defined for the task and is
More informationLecture 4: September Reminder: convergence of sequences
36-705: Intermediate Statistics Fall 2017 Lecturer: Siva Balakrishnan Lecture 4: September 6 In this lecture we discuss the convergence of random variables. At a high-level, our first few lectures focused
More informationBayesian Nonparametric Regression for Diabetes Deaths
Bayesian Nonparametric Regression for Diabetes Deaths Brian M. Hartman PhD Student, 2010 Texas A&M University College Station, TX, USA David B. Dahl Assistant Professor Texas A&M University College Station,
More informationIntroduction: exponential family, conjugacy, and sufficiency (9/2/13)
STA56: Probabilistic machine learning Introduction: exponential family, conjugacy, and sufficiency 9/2/3 Lecturer: Barbara Engelhardt Scribes: Melissa Dalis, Abhinandan Nath, Abhishek Dubey, Xin Zhou Review
More informationBayesian Methods for Machine Learning
Bayesian Methods for Machine Learning CS 584: Big Data Analytics Material adapted from Radford Neal s tutorial (http://ftp.cs.utoronto.ca/pub/radford/bayes-tut.pdf), Zoubin Ghahramni (http://hunch.net/~coms-4771/zoubin_ghahramani_bayesian_learning.pdf),
More informationST5215: Advanced Statistical Theory
Department of Statistics & Applied Probability Wednesday, October 19, 2011 Lecture 17: UMVUE and the first method of derivation Estimable parameters Let ϑ be a parameter in the family P. If there exists
More informationRare event simulation for the ruin problem with investments via importance sampling and duality
Rare event simulation for the ruin problem with investments via importance sampling and duality Jerey Collamore University of Copenhagen Joint work with Anand Vidyashankar (GMU) and Guoqing Diao (GMU).
More informationEcon 2148, fall 2017 Gaussian process priors, reproducing kernel Hilbert spaces, and Splines
Econ 2148, fall 2017 Gaussian process priors, reproducing kernel Hilbert spaces, and Splines Maximilian Kasy Department of Economics, Harvard University 1 / 37 Agenda 6 equivalent representations of the
More informationStatistical Theory MT 2007 Problems 4: Solution sketches
Statistical Theory MT 007 Problems 4: Solution sketches 1. Consider a 1-parameter exponential family model with density f(x θ) = f(x)g(θ)exp{cφ(θ)h(x)}, x X. Suppose that the prior distribution has the
More informationUncertain Inference and Artificial Intelligence
March 3, 2011 1 Prepared for a Purdue Machine Learning Seminar Acknowledgement Prof. A. P. Dempster for intensive collaborations on the Dempster-Shafer theory. Jianchun Zhang, Ryan Martin, Duncan Ermini
More informationStatistical Theory MT 2006 Problems 4: Solution sketches
Statistical Theory MT 006 Problems 4: Solution sketches 1. Suppose that X has a Poisson distribution with unknown mean θ. Determine the conjugate prior, and associate posterior distribution, for θ. Determine
More informationExplicit Bounds for the Distribution Function of the Sum of Dependent Normally Distributed Random Variables
Explicit Bounds for the Distribution Function of the Sum of Dependent Normally Distributed Random Variables Walter Schneider July 26, 20 Abstract In this paper an analytic expression is given for the bounds
More informationLecture notes on statistical decision theory Econ 2110, fall 2013
Lecture notes on statistical decision theory Econ 2110, fall 2013 Maximilian Kasy March 10, 2014 These lecture notes are roughly based on Robert, C. (2007). The Bayesian choice: from decision-theoretic
More informationResearch Article Optimal Portfolio Estimation for Dependent Financial Returns with Generalized Empirical Likelihood
Advances in Decision Sciences Volume 2012, Article ID 973173, 8 pages doi:10.1155/2012/973173 Research Article Optimal Portfolio Estimation for Dependent Financial Returns with Generalized Empirical Likelihood
More informationActuarial models. Proof. We know that. which is. Furthermore, S X (x) = Edward Furman Risk theory / 72
Proof. We know that which is S X Λ (x λ) = exp S X Λ (x λ) = exp Furthermore, S X (x) = { x } h X Λ (t λ)dt, 0 { x } λ a(t)dt = exp { λa(x)}. 0 Edward Furman Risk theory 4280 21 / 72 Proof. We know that
More informationLecture 2. We now introduce some fundamental tools in martingale theory, which are useful in controlling the fluctuation of martingales.
Lecture 2 1 Martingales We now introduce some fundamental tools in martingale theory, which are useful in controlling the fluctuation of martingales. 1.1 Doob s inequality We have the following maximal
More informationA Dynamic Contagion Process with Applications to Finance & Insurance
A Dynamic Contagion Process with Applications to Finance & Insurance Angelos Dassios Department of Statistics London School of Economics Angelos Dassios, Hongbiao Zhao (LSE) A Dynamic Contagion Process
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 informationGeneral Bayesian Inference I
General Bayesian Inference I Outline: Basic concepts, One-parameter models, Noninformative priors. Reading: Chapters 10 and 11 in Kay-I. (Occasional) Simplified Notation. When there is no potential for
More information9 Bayesian inference. 9.1 Subjective probability
9 Bayesian inference 1702-1761 9.1 Subjective probability This is probability regarded as degree of belief. A subjective probability of an event A is assessed as p if you are prepared to stake pm to win
More informationFactor-Adjusted Robust Multiple Test. Jianqing Fan (Princeton University)
Factor-Adjusted Robust Multiple Test Jianqing Fan Princeton University with Koushiki Bose, Qiang Sun, Wenxin Zhou August 11, 2017 Outline 1 Introduction 2 A principle of robustification 3 Adaptive Huber
More informationCompetitive Equilibria in a Comonotone Market
Competitive Equilibria in a Comonotone Market 1/51 Competitive Equilibria in a Comonotone Market Ruodu Wang http://sas.uwaterloo.ca/ wang Department of Statistics and Actuarial Science University of Waterloo
More informationSOCIETY OF ACTUARIES EXAM STAM SHORT-TERM ACTUARIAL MODELS EXAM STAM SAMPLE SOLUTIONS
SOCIETY OF ACTUARIES EXAM STAM SHORT-TERM ACTUARIAL MODELS EXAM STAM SAMPLE SOLUTIONS Questions -37 have been taken from the previous set of Exam C sample questions. Questions no longer relevant to the
More informationProbability Transforms with Elliptical Generators
Probability Transforms with Elliptical Generators Emiliano A. Valdez, PhD, FSA, FIAA School of Actuarial Studies Faculty of Commerce and Economics University of New South Wales Sydney, AUSTRALIA Zurich,
More informationA Bayesian Analysis of Some Nonparametric Problems
A Bayesian Analysis of Some Nonparametric Problems Thomas S Ferguson April 8, 2003 Introduction Bayesian approach remained rather unsuccessful in treating nonparametric problems. This is primarily due
More informationIntroduction to Machine Learning. Lecture 2
Introduction to Machine Learning Lecturer: Eran Halperin Lecture 2 Fall Semester Scribe: Yishay Mansour Some of the material was not presented in class (and is marked with a side line) and is given for
More informationPrequential Analysis
Prequential Analysis Philip Dawid University of Cambridge NIPS 2008 Tutorial Forecasting 2 Context and purpose...................................................... 3 One-step Forecasts.......................................................
More information2 Statistical Estimation: Basic Concepts
Technion Israel Institute of Technology, Department of Electrical Engineering Estimation and Identification in Dynamical Systems (048825) Lecture Notes, Fall 2009, Prof. N. Shimkin 2 Statistical Estimation:
More informationLECTURE 5 NOTES. n t. t Γ(a)Γ(b) pt+a 1 (1 p) n t+b 1. The marginal density of t is. Γ(t + a)γ(n t + b) Γ(n + a + b)
LECTURE 5 NOTES 1. Bayesian point estimators. In the conventional (frequentist) approach to statistical inference, the parameter θ Θ is considered a fixed quantity. In the Bayesian approach, it is considered
More informationCoherent Risk Measures. Acceptance Sets. L = {X G : X(ω) < 0, ω Ω}.
So far in this course we have used several different mathematical expressions to quantify risk, without a deeper discussion of their properties. Coherent Risk Measures Lecture 11, Optimisation in Finance
More informationTail Dependence of Multivariate Pareto Distributions
!#"%$ & ' ") * +!-,#. /10 243537698:6 ;=@?A BCDBFEHGIBJEHKLB MONQP RS?UTV=XW>YZ=eda gihjlknmcoqprj stmfovuxw yy z {} ~ ƒ }ˆŠ ~Œ~Ž f ˆ ` š œžÿ~ ~Ÿ œ } ƒ œ ˆŠ~ œ
More informationLecture 2: Basic Concepts of Statistical Decision Theory
EE378A Statistical Signal Processing Lecture 2-03/31/2016 Lecture 2: Basic Concepts of Statistical Decision Theory Lecturer: Jiantao Jiao, Tsachy Weissman Scribe: John Miller and Aran Nayebi In this lecture
More informationStatistics. Statistics
The main aims of statistics 1 1 Choosing a model 2 Estimating its parameter(s) 1 point estimates 2 interval estimates 3 Testing hypotheses Distributions used in statistics: χ 2 n-distribution 2 Let X 1,
More informationProbability and Measure
Part II Year 2018 2017 2016 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2018 84 Paper 4, Section II 26J Let (X, A) be a measurable space. Let T : X X be a measurable map, and µ a probability
More informationCourse 1 Solutions November 2001 Exams
Course Solutions November Exams . A For i =,, let R = event that a red ball is drawn form urn i i B = event that a blue ball is drawn from urn i. i Then if x is the number of blue balls in urn, ( R R)
More informationFlexible Estimation of Treatment Effect Parameters
Flexible Estimation of Treatment Effect Parameters Thomas MaCurdy a and Xiaohong Chen b and Han Hong c Introduction Many empirical studies of program evaluations are complicated by the presence of both
More informationModern Methods of Statistical Learning sf2935 Auxiliary material: Exponential Family of Distributions Timo Koski. Second Quarter 2016
Auxiliary material: Exponential Family of Distributions Timo Koski Second Quarter 2016 Exponential Families The family of distributions with densities (w.r.t. to a σ-finite measure µ) on X defined by R(θ)
More information1 of 7 7/16/2009 6:12 AM Virtual Laboratories > 7. Point Estimation > 1 2 3 4 5 6 1. Estimators The Basic Statistical Model As usual, our starting point is a random experiment with an underlying sample
More information