UEA Actuarial Methods CMP-5001B 2015/2016 Actuarial methods β Calculations Lecture 1 Introduction These six lectures run alongside the sessions at Aviva. They develop in more detail some of the calculations that may be used in General Insurance. They also develop applications for some of the statistical techniques you studied last semester. The material in these six lectures will be examined in a two hour exam in May/June. There should be no conflict between the material taught here and that taught in other courses. If you feel there is a difference in emphasis please talk to me. Aviva are approaching the subject from the viewpoint of practical experience whereas I am sticking more rigorously to the approach in the institute/faculty of actuaries syllabus and course (CT6). This first lecture will look at some statistical distributions and their applications in general insurance. In the second lecture we look at risk models. The third lecture will take a simple look at ruin theory which will be developed more fully in the final year. In the fourth lecture we will develop some techniques based on Bayesian statistics. The fifth lecture will continue looking at Bayesian methods and specifically credibility theory. The sixth lecture will talk about run off triangles. This session is where there will be a degree of overlap with the Aviva course and I will be relying on the fact that you have been taught the chain ladder method in their session eight, and spend more time on developing other methods. Background to General Insurance Features of Insurance The policyholder must have an interest in the risk being insured. This is a primary difference between insurance and a bet. So you can take out a policy to cover your own home but not to cover the home of your accident prone friend. The risk must be relatively easy to calculate in financial terms. The relevant calculation basis will be set out in the contract eg if your TV explodes you may get a new replacement of a similar type or you may get only the money to replace with a similar model of a similar age, depending on the wording of your contract. Ideally all risks insured will be independent of one another. This will not always be possible, an insurer will not refuse to insure a property because they already insure the neighbour. Events that are being insured should be relatively rare. You cannot insure against the grass in your garden growing. Insurance groups large numbers of similar (but still ideally independent) risks this should reduce the variance and enable you to estimate the potential cost more easily. If you are insuring static caravans 1 UEA Actuarial Methods CMP-5001B 2015/2016 against fire your ideal portfolio would include a large number of caravans spread over a large number of caravan parks rather than the same number of caravans spread over 10 caravan parks in one county. The insurer needs to set a limit on the liability they can incur under each contract. This may be by writing a limit in the contract or by the use of reinsurance which we will talk about later on today. The limit could be the replacement cost of an appliance. The insurer will try to prevent moral hazards. So the terms of a motor insurance policy will discourage you from driving into a wall to write the car off, for example they will not replace your 10 year old car with a brand new one. Clearly it will not be possible to meet all the above criteria in every insurance contract. Contracts will normally be for a fixed time period, very often this will be a year but for example travel insurance may only cover one trip, although even that will usually cover you for some risks from the time you book the trip. Usually the claims are not for a fixed amount as losses can vary, for example if your home is burgled, payment will be based on what has actually been stolen not on a set amount for any burglary. A claim occurring does not usually bring the policy to an end, in fact you may claim more than once in a year, you may have a stone hitting the car windscreen and claim for new glass and then be involved in an accident a few weeks later. Although the policy itself is usually for a fixed term claims may occur at any time during the term or even after this term. This is particularly the case for say industrial injuries, it may be many years before damage due to inhalation of asbestos comes to light for example. Estimating these future claims is a large part of the work we will be looking at in this course. Reinsurance Reinsurance is a type of insurance that the insurer takes out to protect itself against excessive payouts to its policyholders. There are special insurance companies known as reinsurance companies who specialise in this business. We will in this course be especially interested in two types of reinsurance known as excess of loss reinsurance and proportional reinsurance. It is important you understand how they work because it is possible they may crop up in questions across the range of calculation techniques we develop. Excess of loss reinsurance The insurer and reinsurer agree what is known as a retention limit. The insurer pays the full claim value if the value is less than the retention limit. If the claim amount is more than the retention limit then the insurer pays the retention limit and the reinsurer the excess over this limit. So for example the insurer has a retention limit of £40,000 on motor insurance claims, the majority of these will be less than £40,000 so the insurer pays the full amount but if somebody crashes a brand new Ferrari and writes it off then the insurer will pay out £40,000 and the reinsurer the value say £150,000 less the retention so in this case they will pay £110,000. 2 UEA Actuarial Methods CMP-5001B 2015/2016 Proportional reinsurance The other form of reinsurance we will consider is proportional where it is agreed that the insurer and reinsurer will each pay the same proportion of every claim. So if the proportion is 80% insurer 20% reinsurer then for a repair costing £1500 the insurer will pay £1200 and the reinsurer £300. For the £150,000 Ferrari the insurer will pay £120,000 and the reinsurer £30,000. Policy excess An excess on a policy is a way insurers discourage small claims. The insured has to pay the first £X of any claim. The excess will vary according to the type of policy but typically it will be at least £50. This means it is not worth claiming for very minor repairs. If you as the insured want to reduce the cost of your insurance one way of doing this is to agree to a larger excess but you will only save if you do not make a claim. Loss Distributions In general insurance there are a number of situations where we need to models to enable us to make estimates and this is often done by applying statistical distributions. We will look at some of the distributions which may be useful. Some of these you will have already met in your statistics modules, others are more specific to Actuarial work you will find key information for 19 distributions (7 discrete, 12 continuous) in the tables (yellow pages 6-15). We will not look at all of them in this course and some you will already be familiar with anyway. Whilst it is useful to know they are there in the constraints of an exam you will not want to spend your time leafing through these. Probably the most fundamental need of any insurer is estimating the total amount that will be paid out in claims during say a year or a month. It might be possible to model this directly. Alternatively we could model separately the number of claims and the amount of individual claims. Immediately we see why we need both discrete and continuous distributions. The number of claims is clearly discrete whereas the amount may be regarded as continuous. We could assume we have claims from a known distribution but often we will need to estimate the parameters of the distribution from the data we have. It is usually assumed that claims distributions are positively skewed and long tailed. Defining our distributions Discrete distributions have a probability function π(π₯) which gives the probability that π₯ takes each possible discrete value where the values may be finite or infinite. If we look at tossing a coin 10 times and observing the number of heads then x can only take integer values between 0 and 10. Whereas if we look at the number of claims on a very large insurance portfolio we may regard this as infinite especially for general insurance when usually a policyholder can make any number of claims on a policy. ββ ββ π(π₯) = 1 Continuous distributions have a probability density function (PDF) π(π₯) this can be regarded as analogous to the probability function. We tend to regard distributions of claim amounts as continuous even if in practice an amount can only be defined down to pence. Again it is possible to have distributions which are finite or infinite if there is a maximum total claim per annum then this is finite if there is no maximum β then it may be regarded as infinite. β«ββ π(π₯) = 1 3 UEA Actuarial Methods CMP-5001B 2015/2016 Both discrete and continuous distributions have what the tables describe as the distribution function πΉ(π₯). This is perhaps more usefully (and commonly) known as the cumulative distribution function. It gives the probability that a function has a value of π or less. Because we know the total probability will be 1 we can also find the probability that it is greater than π as 1 β πΉ(π). The formula for the cumulative π distribution function may be complex it is βπ₯ββ π(π₯) or β«ββ π(π₯) in these cases values may be tabulated. You will be familiar with the use of normal distribution tables. If we want the probability a value lies π between a and b we want to find πΉ(π) β πΉ(π) or βππ π(π₯) or β«π π(π₯). To complete the definition of the distribution we need to specify a parameter or parameters. One of the major areas we will be looking at in this subject is estimating these parameters. The tables also include probability generating functions we will not use these in this subject: β Once we have a distribution of claims then we can use it to calculate the average claim β«0 π₯ π(π₯) or to find π the probability that a claim is between a & b β«π π(π₯). This latter calculation could also be done using the cumulative distribution function πΉ(π) β πΉ(π). Moments We are very familiar with the first two moments of a distribution, the mean and the variance. The tables give formulae for these in terms of the parameters. We may need to calculate these from first principles in which case we use: π πΈ(π₯) = β«π π₯π(π₯)ππ₯ or πΈ(π₯) = βππ π₯π(π₯) as appropriate And πππ(π₯) = πΈ(π₯ 2 ) β (πΈ(π₯))2 For a number of distributions they also give us a moment generating function. A moment generating function is ππ (π‘) = πΈ(π π‘π ) = 1 + π‘π1 + moment about the origin. π‘ 2 π2 2! + π‘ 3 π3 3! + β― where ππ is the π π‘β Moment generating functions do not always exist but if they do then to get the π π‘β moment about the origin we differentiate π times with respect to π‘ and set π‘ = 0. Useful distributions We start by looking at some useful continuous distributions we will look at some more distributions including discrete ones next time. In general we will use these to look at claim size. 4 UEA Actuarial Methods CMP-5001B 2015/2016 Exponential distribution π(π₯) = ππ βππ₯ πΉ(π₯) = 1 β π βππ₯ 1 1 The mean of the exponential is π and variance π2 exponential lambda = 0.5 0.5 0.4 0.3 exponential lambda = 0.5 0.2 0.1 0 0.10.91.72.53.34.14.95.76.57.38.18.99.7 The shape of the exponential is appropriate for larger claims but in reality people do not tend to claim very small amounts. In fact many policies contain a policy excess clause as mentioned earlier meaning that they do not accept claims less than £50 or £100 say. In this case they will deduct the excess from any claim they do pay (for this reason it is sometimes known as a deductible). There may also be reinsurance which means the insurance has an upper limit on the amount they pay. The exponential distribution may not be the distribution of choice for insurers but it is the distribution of choice when asking questions as it is relatively easy to calculate with. Many questions in this course will require you to find the mean (and sometimes variance) of what is known as a curtate exponential distribution, in other words just a portion of the distribution. To do this you will need to use integration by parts which you should be familiar with: β« π’ππ£ = π’π£ β π£ππ’ Let us consider the situation where an insurer wants to find the mean claim size. They assume that claims are exponentially distributed and if they settled every claim in full the mean claim size would be £200. In fact there excess of loss insurance meaning that the reinsurer will pay the value of the claim over and above £500. Find the expected value of the insurers portion of the claim. 500 In this case we can find the expected value as β«0 β π₯π(π₯) + β«500 500π(π₯) Which for the exponential distribution becomes: 5 UEA Actuarial Methods CMP-5001B 2015/2016 500 β« π₯ππ β βππ₯ + β« 500ππ βππ₯ 0 500 1 1 Now we know the mean claim size is £200 and πΈ(π) = π so 200 = π Giving π = 0.005 500 β« β 0.005π₯π β0.005π₯ + β« 500 β 0.005π β0.005π₯ 0 500 The second integral is straightforward and we will leave it for now except for bringing constant outside. The first integral can be evaluated using integration by parts: Let π’ = π₯ and ππ£ = 0.005π β0.005π₯ so π£ = β« ππ£ = β« 0.005π β0.005π₯ = β π β0.005π₯ And ππ’ = 1 So our formula becomes 500 β 500 [βπ₯π β0.005π₯ ] + β« π β0.005π₯ + 2.5 β« π β0.005π₯ 0 0 500 We can then proceed with our new integral and the second integral from our original equation 500 1 500 β1 β0.005π₯ β + [β π β0.005π₯ ] + 2.5 [ π ] 0 0.005 0 0.005 500 1 1 β1 ββ 1 [β500π β0.005β500 + 0] + [β π β0.005β500 + π 0 ] + 2.5 [ π + π β0.005β500 ] 0.005 0.005 0.005 0.005 1 1 1 [β500π β2.5 ] + [β π β2.5 + π β2.5 ] ] + 2.5 [ 0.005 0.005 0.005 1 1 β500π β2.5 β π β2.5 + + 500π β2.5 0.005 0.005 First and last term cancel out 1 1 β π β2.5 = 200 β 16.417 = 183.58 0.005 0.005 [βπ₯π β0.005π₯ ] You will normally end up with terms cancelling out as above. If you are given a policy excess E then the lower limit on your first integral will be E and the amount payable will be π₯ β πΈ. If you want a variance you will need to calculate πΈ(π₯ 2 ) which will mean you need to integrate by part twice. 6 UEA Actuarial Methods CMP-5001B 2015/2016 Normal Distribution The normal distribution is not in reality a particularly useful distribution for claim amounts. However as it is relatively easy to use it is often used in questions to show understanding. π(π₯) = 1 πβ2π π 1 π₯βπ 2 β ( ) 2 π Where ππππ = π and ππππππππ = π 2 πΉ(π₯) is tabulated for a standard normal distribution where π = 0 and π 2 = 1 You will be pleased to know you will not be expected to integrate this to find the mean of a curtate distribution Normal mu=0 sigma =1 0.5 0.4 0.3 Normal 0.2 0.1 0 -3 -2.6-2.2-1.8-1.4 -1 -0.6-0.2 0.2 0.6 1 1.4 1.8 2.2 2.6 3 The type of question you may get asked is what retention limit should be set that if we want 95% of claims to be less than the retention limit. 7 UEA Actuarial Methods CMP-5001B 2015/2016 Lognormal distribution The lognormal distribution is defined as X has a lognormal distribution if Log(X) is normally distributed. This is potentially a realistic claim distribution and can be used in a similar way to the normal distribution. π(π₯) = 1 1 β1(ππ(π₯)βπ) π 2 π π β2π π₯ 2 you would not be expected to remember this. In this case the parameters do not represent the mean and variance which can be calculated as 1 2 2 2 Mean = π π+2π variance = π 2π+π (π π β 1) lognormal mu=0.5 sigma=0.5 1 0.8 0.6 lognormal 0.4 0.2 0.1 0.7 1.3 1.9 2.5 3.1 3.7 4.3 4.9 5.5 6.1 6.7 7.3 7.9 8.5 9.1 9.7 0 8 UEA Actuarial Methods CMP-5001B 2015/2016 Gamma distribution The Gamma distribution is defined as ππΌ π(π₯) = π₯ πΌβ1 π βππ₯ Ξ(πΌ) Reminder from your statistics: Ξ(1) = 1, Ξ(πΌ) = πΌ β 1! if πΌ is integer F(x) needs to be tabulated. πΌ πΌ Mean =π Variance = π2 Gamma alpha=1 gamma =0.5 0.5 0.4 0.3 Gamma alpha=1 gamma =0.5 0.2 0.1 0 0.10.71.31.92.53.13.74.34.95.56.16.77.37.98.59.19.7 There are few questions in this particular part of the course using the Gamma distribution but we will use it later on. 9 UEA Actuarial Methods CMP-5001B 2015/2016 Uniform Distribution Whilst the uniform distribution is not really a realistic claims distribution it can be useful to demonstrate calculations because of its simplicity. It has two parameters π and π where these are the lower and upper limits respectively. 1 π(π₯) = πβπ π₯βπ πΉ(π₯) = πβπ 1 1 Mean = 2 (π + π) Variance = 12 (π β π)2 Uniform a=0 b=100 0.012 0.01 0.008 Uniform a=0 b=100 0.006 0.004 0.002 0 8 16 24 32 40 48 56 64 72 80 88 96 0 10 UEA Actuarial Methods CMP-5001B 2015/2016 Pareto distribution The Pareto distribution can be described with varying numbers of parameters. The tables show both two and three parameters we are only concerned with the two parameter version in this course. The course actually specifies the two parameter version differently from the book. I set out below the course version i.e. the one you need to know. πΌππΌ π(π₯) = (π+π₯)πΌ+1 π πΌ ) π+π₯ πΉ(π₯) = 1 β ( The tail of the Pareto distribution is heavier than the tail of the exponential this can be useful to an insurer who wishes to be cautious as it errs on the side of there being a higher claim payments. The mean of the Pareto is π and the πΌβ1 πΌπ2 (πΌβ2) variance is (πΌβ1)2 Pareto alpha= 0.3 & lambda = 0.6 0.45 0.4 0.35 0.3 0.25 0.2 0.15 0.1 0.05 0 0.1 0.7 1.3 1.9 2.5 3.1 3.7 4.3 4.9 5.5 6.1 6.7 7.3 7.9 8.5 9.1 9.7 Pareto alpha= 0.3 & lambda = 0.6 11 UEA Actuarial Methods CMP-5001B 2015/2016 Weibull distribution The Weibull distribution gives a tail somewhere between the exponential and the Pareto if the parameter πΎ < 1. It is a very flexible distribution. If πΎ > 1 then the tail is lighter than the exponential The distribution is defined as πΎ π(π₯) = ππΎπ₯ πΎβ1 π βππ₯ πΎ πΉ(π₯) = 1 β π βππ₯ The moments of the Weibull distribution are π 1 πΈ(π π ) = Ξ (1 + πΎ) π for r= 1,2,3..... ππΎ Weibull c=0.3 gamma=0.6 0.45 0.4 0.35 0.3 0.25 0.2 0.15 0.1 0.05 0 Weibull c=0.3 gamma=0.6 0.1 0.7 1.3 1.9 2.5 3.1 3.7 4.3 4.9 5.5 6.1 6.7 7.3 7.9 8.5 9.1 9.7 12 UEA Actuarial Methods CMP-5001B 2015/2016 Fitting a distribution Whilst questions may say claims are distributed according to a known distribution with known parameters in reality we may assume a distribution and are unlikely to know the actual parameters. We will therefore need to use our data to estimate the parameters. We will look at three main methods used for fitting a distribution, the method of moments, maximum likelihood and method of percentiles. Different methods will be more suitable for fitting different distributions. Method of moments In this intuitively straightforward method of estimation we can calculate the parameters of any distribution by setting sample moments about 0, in other words we do not use the variance as our second moment but πΈ(π 2 ) to distribution moments. If the distribution has several parameters then you will need to do this for several moments for which you could use the moment generating function and then solve simultaneous equations. This method can be used very simply for the exponential distribution although it may not be a particularly good estimator for a one parameter distribution. It can also be used for many of the other distributions as an illustration we will look at the lognormal distribution as this is more difficult than some. Example We believe claims are log normally distributed but do not know the parameters of the distribution. In the last year the mean claim size has been 200 and the variance of claims has been 2500. Use the method of moments to fit a lognormal distribution. So we have 1 2 200 = π π+2π and 2 2 2500 = π 2π+π (π π β 1) 2 1 2 2 Now we note that π 2π+π = (π π+2π ) 2 So 2500 = 2002(π π β 1) 2 2500 = 40000π π β 40000 42500 2 = ππ 40000 Taking natural logs 0.06062 = π 2 π = 0.24622 Substituting into first equation 1 200 = π π+20.06062 Taking natural logs 5.29832 = π + 0.03031 π = 5.26801 13 UEA Actuarial Methods CMP-5001B 2015/2016 Method of Percentiles This method is similar to the method of moments, but rather than equating moments you equate percentiles for example you could equate the quartiles. For certain distributions where method of moments and maximum likelihood are difficult this can prove useful. If there is only one parameter then you would equate the median. If there are two parameters then you would normally use the upper and lower quartile, but the method is not uniquely defined so other percentiles may be used. We will look at an example for the normal distribution. Example We have observed a number of claims on a class of business where we believe claim values are normally distributed. 25% of the claims are below £300 and 25% are above £500. Estimate the parameters of the distribution using the method of percentiles. Using the standard normal distribution tables on pages 160-161 of the tables we find the value of π₯ for which 75% of observations are to the left of π₯. This falls somewhere between 0.67 and 0.68. We could then interpolate to find a closer value. However very helpfully if you turn over to page 162 they have actually done this for us for common percentages and we can look up 25% in the tail as 0.6745. (It is up to you which of these methods you use. The latter is quicker but as you may not have this table available in your statistics exam you may prefer to use the former method for consistency) 300βπ 500βπ So we have π = β0.6745 and π = 0.6745 300 β π = β0.6745π and 500 β π = 0.6745π Adding these two equations gives 800 β 2π = 0 So π = 400 as we might expect Therefore 500 β 400 = 0.6745π And hence π = 148 14 UEA Actuarial Methods CMP-5001B 2015/2016 Method of maximum likelihood Put simply if we want to estimate π for each observed value we calculate the probability of observing π₯π given that the distribution has parameter π. Technically this is equivalent to the likelihood of π given π₯π . If we have a sample size of n we multiply these n probabilities together and find π to give this a maximum probability. In formulae we calculate a likelihood function π πΏ(π) = β π(π₯π | π) 1 It can be more practical to deal with addition rather than multiplication so we can take logs π ππππΏ(π) = β ππππ(π₯π | π) 1 Maximising either of these will give the maximum likelihood estimator for π. To maximise we need to differentiate the expression and set the differential equal to zero. It is important to check that we have found a maximum and not a minimum. As an example we will look at estimating the parameter of an exponential distribution. Example We believe claim values are exponentially distributed we want to estimate the parameter π of the distribution. This year we have had 200 claims for a total of £100,000. The Likelihood function will be πΏ = ππ βππ₯1 β ππ βππ₯2 β ππ βππ₯3 β β¦ ππ βππ₯200 200 πΏ = β ππ βππ₯π 1 200 200 ππππΏ = β ππππ β β ππ₯π 1 1 200 ππππΏ = 200ππππ β π β β π₯π 1 ππππΏ = 200ππππ β π β 100000 Differentiating we get 200 β 100000 π Equating this to zero gives 200 = 100000 π Or π = 0.002 15 UEA Actuarial Methods CMP-5001B 2015/2016 1 This gives a mean of 0.002 = 500 which is what we would expect. Consider now the case where there is an excess of loss reinsurance treaty with the reinsurer paying claims above £1000. In this case we have 180 claims below £1000 totalling £80,000 and 20 claims where the insurer pays £1000 (and reinsurer pays balance) The probability that a claim is greater than £1000 is π βπ1000 (from 1 β πΉ(π₯)) So likelihood function becomes πΏ = ππ βππ₯1 β ππ βππ₯2 β ππ βππ₯3 β β¦ ππ βππ₯180 * π βπ1000 β β¦ π βπ1000 180 πΏ = β ππ βππ₯π β π β20π1000 1 180 180 ππππΏ = β ππππ β β ππ₯π β 20000π 1 1 180 ππππΏ = 180ππππ β π β β π₯π β 20000π 1 ππππΏ = 180ππππ β π β 80000 β 20000π Differentiating we get 180 β 100000 π Equating this to zero gives 180 = 100000 π Or π = 0.0018 16 UEA Actuarial Methods CMP-5001B 2015/2016 Effects of reinsurance Without reinsurance we can say expected level of claim is β β«0 π₯ π(π₯) If we have excess of loss reinsurance with retention R then expected level of claim will be π β β«0 π₯ π(π₯) +R β«π π(π₯) If we also allow for k% inflation in claim sizes but leave the retention unchanged we have expected claims of π /(1+π) β (1 + π)π₯ π(π₯) +R β«π /(1+π) π(π₯) β«0 In the case of excess of loss reinsurance method of moments cannot be employed so we need to use method of percentiles or maximum likelihood. If we have proportional reinsurance where the insurer is liable for L% the expected level of claim will be β πΏ β« π₯ π(π₯) 100 0 Effects of policy excess This works very similarly to excess of loss reinsurance expected claims if excess is E will be πΈ β 0 β«0 π(π₯) + β«πΈ (π₯ β πΈ)π(π₯) Exercises 1. April 2005 question 5 2. April 2007 question 3 (This question is not difficult in fact it is similar to 1. but the language is a bit different) 3. September 2009 question 5 4. September 2010 question 7 5. An insurer believes claims on a particular class of business are normally distributed with mean £600 and variance £80,000. They have found that it costs them £100 to administer each claim and on this basis they decide to impose a policy excess of £100 what percentage of claims will now not be paid. If you want something a bit more challenging try these too 6. September 2006 question 8 7. September 2008 question 11 17
© Copyright 2026 Paperzz