Mathematical models for insurance payments – part 3 – other modifications

Posted on Updated on

This post is a continuation of the discussion on models of insurance payments initiated in two previous posts. Part 1 focuses on the models of insurance payments in which the insurance policy imposes a policy limit. Part 2 continues the discussion by introducing models in which the insurance policy imposes an ordinary deductible. In each of these two previous posts, the insurance coverage has only one coverage modification. A more interesting and more realistic scenario would be insurance coverage that contains a combination of several coverage modifications. This post is to examine the effects on the insurance payments as a result of having some or all of these coverage modifications – policy limit, ordinary deductible, franchise deductible and inflation. Additional topics: expected payment per loss versus expected payment per payment and loss elimination ratio.

\text{ }

Ordinary Deductible and Policy Limit

The previous two posts discuss the expectations E[X \wedge u] and E[(X-d)_+]. The first is the expected insurance payment when the coverage has a policy limit u. The second is the expected insurance payment when the coverage has an ordinary deductible d. They are the expected values of the following two variables.

    \displaystyle  X \wedge u=\left\{ \begin{array}{ll}                     \displaystyle  X &\ X \le u \\           \text{ } & \text{ } \\           \displaystyle  u &\ X > u           \end{array} \right.

    \displaystyle  (X-d)_+=\left\{ \begin{array}{ll}                     \displaystyle  0 &\ X \le d \\           \text{ } & \text{ } \\           \displaystyle  X-d &\ X > d           \end{array} \right.

It is easy to verify that X=(X-d)_+ + X \wedge d. Buying a coverage with an ordinary deductible d and another coverage with policy limit d equals full coverage. Thus we have the following relation.

(1) \ \ \ \ E[(X-d)_+]=E[X]-E[X \wedge d]

The limited expectation E[X \wedge d] has expressions in closed form in some cases or has expressions in terms of familiar functions (e.g. gamma function) in other cases. Thus the expectation E[(X-d)_+] can be computed by knowing the original expected value E[X] and the limited expectation E[X \wedge d].

Inflation

Suppose that losses (or claims) in the next period are expected to increase uniformly by 100r%. For example, r=0.10 means 10%. What would be the effect of inflation on the expectations E[X \wedge u] and E[(X-d)_+]?

First, the effect on the limited expectation E[X \wedge u]. As usual, X is the random loss and u is the policy limit. With inflation rate r, the loss variable for the next period would be (1+r)X. One approach is to derive the distribution for the inflated loss variable and use the new distribution to calculate E[(1+r)X \wedge u]. Another approach is to express it in terms of the limited expectation of the pre-inflated loss X. The following is the expectation E[(1+r)X \wedge u], assuming that there is no change in the policy limit.

\displaystyle (2) \ \ \ \ E[(1+r) X \wedge u]=(1+r) \ E \biggl[X \wedge \frac{u}{1+r} \biggr]

Relation (2) relates the limited expectation of the inflated variable to the limited expectation of the pre-inflated loss X. It says that the limited expectation of the inflated variable (1+r) X is obtained by inflating the limited expectation of the pre-inflated loss but at a smaller policy limit \frac{u}{1+r}.

The following is the expectation E[(1+r) (X-d)_+]=E[(1+r) X]-E[(1+r) X \wedge d].

\displaystyle \begin{aligned} (3) \ \ \ \ E[((1+r) X-d)_+]&=E[(1+r) X]-E[(1+r) X \wedge d] \\&=(1+r) E[X]-(1+r) \ E \biggl[X \wedge \frac{d}{1+r} \biggr] \\&=(1+r) \biggl( E[X]- \ E \biggl[X \wedge \frac{d}{1+r} \biggr] \biggr) \\&=(1+r) E \biggl[ \biggl(X-\frac{d}{1+r} \biggr)_+ \biggr] \end{aligned}

Similarly, (3) expresses the expected payment on the inflated loss in terms of the expected payment of the pre-inflated loss. It says that when the loss is inflated, the expected payment per loss is obtained by inflating the expected payment on the pre-inflated loss but at a smaller deductible.

Insurance Payment Per Loss versus Per Payment

The previous post (Part 2) shows how to evaluate the average amount paid to the insured when the coverage has an ordinary deductible. The average payment discussed in Part 2 is the average per loss (over all losses). As a simple illustration, let’s say the amounts of losses in a given period for an insured are 7, 4, 33 and 17 subject to an ordinary deductible of 5. Then the insurance payments are: 2, 0, 28 and 12. The average payment per loss would be (2+0+28+12)/4 = 10.5. If we only count the losses that require a payment, the average is (2+28+12)/3 = 14. Thus the average payment per payment is greater than the average payment per loss since only the losses exceeding the deductible are counted in the average payment per payment. In the calculation discussed here, the average payment per payment is obtained by dividing the average payment per loss by the probability that the loss exceeds the deductible. Note that 10.5/0.75 = 14.

Suppose that the random variable X is the size of the loss (if a loss occurs). Under an insurance policy with an ordinary deductible d, the first d dollars of a loss is responsible by the insured and the amount of the loss in excess of d is paid by the insurer. Under such an arrangement, a certain number of losses are not paid by the insurer, precisely those losses that are less than or equal to d. If we only count the losses that are paid by the insurer, the payment amount is the conditional random variable X-d \lvert X>d. The expected value of this conditional random variable is denoted by e_X(d) or e(d) if the loss X is understood.

Given that P(X>d)>0, the variable X-d \lvert X>d is called the excess loss variable. Its expected value e_X(d) is called the mean excess loss function. Other names are mean residual life and complete expectation of life (when the context is that of a mortality study).

For the discussion in this post and other posts in the same series, we use Y_P to denote X-d \lvert X>d. The P stands for payment so that its expected value would be average insurance payment per payment, i.e. the expected amount paid given that the loss exceeds the deductible. When the random variable X is the age at death, e_X(x) would be the expected remaining time until death given that the life has survived to age x.

The expected value e_X(d) is thus the expected payment per payment (or expected cost per payment) under an ordinary deductible. In contrast, the expected value E[(X-d)_+] is the expected payment per loss (or expected cost per loss), which is discussed in the previous post. The two expected values are related. The calculation of one will give the other. The following compares the two calculation. Let f_X(x) and F_X(x) be the PDF and CDF of X, respectively.

\displaystyle (4) \ \ \ \ E[(X-d)_+]=\int_d^\infty (x-d) \ f_X(x) \ dx

\text{ }

\displaystyle (5) \ \ \ \ e_X(d)=\int_d^\infty (x-d) \ \frac{f_X(x)}{S_X(d)} \ dx=\frac{E[(X-d)_+]}{S_X(d)}

Note that e_X(d) is calculated using a conditional density function. As a result, e_X(d) can be obtained by dividing the expected payment per loss divided by the probability that there is a payment. Thus the expected payment per payment is the expected payment per loss divided by the probability that there is a payment. This is described in the following relation.

\displaystyle (6) \ \ \ \ e_X(d)=\frac{E[(X-d)_+]}{S_X(d)}=\frac{E[X]-E[X \wedge d]}{S_X(d)}

\text{ }

\text{ }

Distributions of Insurance Payment Variables

Relation (1) and Relation (4) are used for the calculation of expected insurance payment when the coverage has an ordinary deductible. For deriving other information about insurance payment in the presence of an ordinary deductible, it is helpful to know the distributions of the insurance payment (per loss and per payment).

Let Y_L=(X-d)_+ (payment per loss) and let Y_P=X-d \lvert X>d (payment per payment). We now discuss the distribution for Y_P. First, the following gives the PDF, CDF and the survival function of Y_L.

Payment Per Loss
PDF \displaystyle  f_{Y_L}(y)=\left\{ \begin{array}{ll}                     \displaystyle  F_X(d) &\ \ \ \ y=0 \\           \text{ } & \text{ } \\           \displaystyle  f_X(y+d) &\ \ \ \ y > 0           \end{array} \right.
CDF \displaystyle  F_{Y_L}(y)=\left\{ \begin{array}{ll}                     \displaystyle  0 &\ \ \ \ y<0 \\           \text{ } & \text{ } \\           \displaystyle  F_X(y+d) &\ \ \ \ y \ge 0                      \end{array} \right.
Survival Function \displaystyle  S_{Y_L}(y)=\left\{ \begin{array}{ll}                     \displaystyle  1 &\ \ \ \  y<0 \\           \text{ } & \text{ } \\           \displaystyle  S_X(y+d) &\ \ \ \ y \ge 0                      \end{array} \right.

Note that the above functions have a point mass at y=0 to account for the losses that are not paid. For Y_P=X-d \lvert X>d, there is no point mass at y=0. We only need to consider y>0. Normalizing the function f_X(y+d) would give the PDF of Y_P=X-d \lvert X>d. Thus the following gives the PDF, the survival function and the CDF of Y_P=X-d \lvert X>d.

Payment Per Payment
PDF \displaystyle f_{Y_P}(y)=\frac{f_X(y+d)}{S_X(d)} \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ y > 0
CDF \displaystyle F_{Y_P}(y)=\frac{F_X(y+d)-F_X(d)}{S_X(d)} \ \ \ \ \ \ \ \ \ \ \ \ \ y > 0
Survival Function \displaystyle S_{Y_P}(y)=\frac{S_X(y+d)}{S_X(d)} \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ y > 0

The following calculates the two averages using the respective PDFs, as a result deriving the same relationship between the two expected values.

\displaystyle (7) \ \ \ \ E(Y_L)=E[(X-d)_+]=\int_0^\infty y \ f_X(y+d) \ dy

\text{ }

\displaystyle (8) \ \ \ \ E(Y_P)=e_X(d)=\int_0^\infty y \ \frac{f_X(y+d)}{S_X(d)} \ dy

Relation (1) and Relation (7) are identical. The former is calculated using the distribution of the original loss X and the latter is calculated using the distribution of the payment per loss variable. Similarly, compare Relation (5) and Relation (8). The former is computed from the distribution of the unmodified loss X and the latter is computed using the distribution of the payment per payment variable.

Note that Relation (7) and Relation (8) also lead to Relation (6), which relates the expected payment E[(X-d)_+] and the expected payment e_X(d).

With the PDFs and CDFs for the payment per loss and the payment per payment developed, other distributional quantities can be derived, e.g. hazard rate, variance, skewness and kurtosis.

Loss Elimination Ratio

When the insurance coverage has an ordinary deductible, a natural question is (from the insurance company’s perspective), what is the impact of the deductible on the payment that is made to the insured? More specifically, on average what is the reduction of the payment? The loss elimination ratio is the proportion of the expected loss that is not paid to the insured by the insurer. For example, if the expected loss before application of the deductible is 45 (of which 36 is expected to be paid by the insurer), the loss elimination ratio is 9/45 = 0.20 (20%). In this example, the insurer has reduced its obligation by 20%. More formally, the loss elimination ratio (LER) is defined as:

    \displaystyle LER=\frac{E(X)-E[(X-d)_+]}{E(X)}=\frac{E[X \wedge d]}{E(X)}

LER is the ratio of the expected reduction in payment as a result of imposing the ordinary deductible to the expected payment without the deductible. Though it is possible to define LER as the ratio of the reduction in expected payment as a result of a coverage modification to the expected payment without the modification for a modification other than an ordinary deductible, we do not attempt to further generalize this concept.

Policy with a Limit and a Deductible

In part 1, expected insurance payment under a policy limit is developed. In part 2, expected insurance payment under a policy with an ordinary deductible is developed. We now combine both provisions in the same insurance policy. First, let’s define two terms. A policy limit is the maximum amount that will be paid by a policy. For example, if the policy limit is 10,000, the policy will paid at most 10,000 per loss. If the actual loss is 15,000, then the policy will paid 10,000 and the insured will have to be responsible for the remaining 5,000. On the other hand, maximum covered loss is the level above which no loss will be paid. For example, Suppose that the policy covers up to 10,000 per loss subject to a deductible of 1,000. If the actual loss is 20,000, then the maximum covered loss is 10,000 with the policy limit being 9,000. This is because the policy only covers the first 10,000 with the first 1,000 paid by the insured.

Suppose that an insurance policy has an ordinary deductible d, a maximum covered loss u with d<u and no other coverage modifications. Any loss below d is not paid by the insurer. For any loss exceeding d, the insurer pays the loss amount in excess of d up to the maximum covered loss u, with the policy limit being u-d. The following describes the payment rule more explicitly.

    \displaystyle  (X \wedge u)-(X \wedge d)=\left\{ \begin{array}{ll}                     \displaystyle  0 &\ \ \ \ X \le d \\           \text{ } & \text{ } \\           \displaystyle  X-d &\ \ \ \ d < X \le u \\           \text{ } & \text{ } \\           \displaystyle  u-d &\ \ \ \ X > u           \end{array} \right.

Under such a policy, the maximum amount paid (policy limit) is u-d, which is the maximum covered loss minus the deductible. The policy limit is reached when X > u. When d < X \le u, the payment X \wedge u is X and the payment X \wedge d is d. Then the expected payment per loss under such a policy is:

\displaystyle  (9) \ \ \ \ E[X \wedge u]-E[X \wedge d]

There is not special notation for the expected payment. In words, it is the limited expected value at u (the maximum covered loss) minus the limited expected value at the ordinary deductible d. The higher moments can be derived by evaluating E(Y^k) using the PDF of the loss X where Y=(X \wedge u)-(X \wedge d).

The payment Y=(X \wedge u)-(X \wedge d) is on a per loss basis. It is also possible to consider payment or payment by removing the point mass at zero. The expected payment per payment is obtained by dividing (9) by the probability of a positive payment.

\text{ }

Franchise Deductible

An alternative to the ordinary deductible is the franchise deductible. It works like an ordinary deductible except that when the loss exceeds the deductible, the policy pays the loss in full. The following gives the payment rule.

    \displaystyle  Y=\left\{ \begin{array}{ll}                     \displaystyle  0 &\ \ \ \ X \le d \\           \text{ } & \text{ } \\           \displaystyle  X &\ \ \ \ X > d           \end{array} \right.

Note that when the loss exceeds the deductible (X > d), the policy with a franchise deductible pays more than a policy with the same ordinary deductible. By how much? By the amount d. Thus the expected payment per loss under a franchise deductible is E[(X-d)_+]+d \ S_X(d). The addition of d \cdot S_X(d) reflects the additional benefit when X > d. Instead of deriving calculation specific to franchise deductible, we can derive the payments under franchise deductible by adding the additional benefit appropriately.

\displaystyle  (10) \ \ \ \ E[(X-d)_+]+d \ S_X(d)

\text{ }

Benefit Tables

The previous two posts and this post discuss coverage with two types of deductible (ordinary and franchise) as well as coverage that may include a limit. In order to keep the expected payments straight, the following table organizes the different combinations of coverage options discussed up to this point.

A = Ordinary Deductible Only
B = Franchise Deductible Only
C = Ordinary Deductible + Limit
D = Franchise Deductible + Limit

Expected Cost Per Loss Expected Cost Per Payment
A \displaystyle E[(X-d)_+]=E(X)-E(X \wedge d) \displaystyle \frac{E(X)-E(X \wedge d)}{S_X(d)}=e_X(d)
B \displaystyle E(X)-E(X \wedge d)+d \ S_X(d) \displaystyle \frac{E(X)-E(X \wedge d)}{S_X(d)}+d=e_X(d)+d
C \displaystyle E(X \wedge u)-E(X \wedge d ) \displaystyle \frac{E(X \wedge u)-E(X \wedge d)}{S_X(d)}
D \displaystyle E(X \wedge u)-E(X \wedge d )+d \ S_X(d) \displaystyle \frac{E(X \wedge u)-E(X \wedge d)}{S_X(d)}+d

Row A shows the expected payment in a coverage with an ordinary deductible for both per loss and per payment. The expected per loss payment in Row A is the subject of the previous post and the expected payment per payment is discussed earlier in this post. Row B is for expected payments in a coverage with a franchise deductible. Note that the coverage with a franchise deductible pays more than a coverage with an identical ordinary deductible. Thus Row B is Row A plus the added benefit, which is the deductible d.

Row C is for the coverage of an ordinary deductible with a policy limit. The per loss expected payment is \displaystyle E(X \wedge u)-E(X \wedge d ). The per payment expected value is a conditional one, conditional on the loss greater than the deductible. Thus the per payment expected value is the per loss expected value divided by S_X(d).

Note that Row D is Row C plus the amount of d, the additional benefit as a result of having a franchise deductible instead of an ordinary deductible.

Rather than memorizing these formulas, focus on the general structure of the table. For example, understand the payments involving ordinary deductible (Row A and Row C). Then the payments for franchise deductible are obtained by adding an appropriate additional benefit.

The following table shows the expected payments under the influence of claim cost inflation. The maximum covered loss u and the deductible d are identical to the above benefit table. The losses are subject to a (100r)% inflation in the next exposure period. Note that u^* and d^* are the modified maximum covered loss and deductible that will make the formulas work.

E = Ordinary Deductible Only
F = Franchise Deductible Only
G = Ordinary Deductible + Limit
H = Franchise Deductible + Limit

\displaystyle u^*=\frac{u}{1+r}

\displaystyle d^*=\frac{d}{1+r}

Expected Cost Per Loss (with Inflation) Expected Cost Per Payment (with Inflation)
E \displaystyle (1+r) \biggl \{ E(X)-E [X \wedge d^* ] \bigg \} \displaystyle (1+r) \biggl \{ \frac{\displaystyle  E(X)-E [X \wedge d^* ]}{\displaystyle S_{X} (d^*)} \biggr \}
F \displaystyle (1+r) \biggl \{ E(X)-E [X \wedge d^* ]+d^* \ S_{X}(d^* ) \bigg \} \displaystyle (1+r) \biggl \{ \frac{\displaystyle  E(X)-E [X \wedge d^* ]}{\displaystyle S_{X} (d^* )}+d^* \biggr \}
G \displaystyle (1+r) \biggl \{ E [X \wedge u^* ]-E [X \wedge d^* ] \bigg \} \displaystyle (1+r) \biggl \{ \frac{\displaystyle  E [X \wedge u^* ]-E [X \wedge d^*]}{\displaystyle S_{X} (d^* )} \biggr \}
H \displaystyle (1+r) \biggl \{ E [X \wedge u^* ]-E [X \wedge d^* ]+d^* \ S_{X}   (d^* ) \bigg \} \displaystyle (1+r) \biggl \{ \frac{\displaystyle  E[X \wedge u^* ]-E [X \wedge d^* ]}{\displaystyle S_{X} (d^* )}+d^* \biggr \}

There is really no need to mathematically derive the formulas for the second insurance benefit table (Rows E through H). Recall the effect of inflation on the expected payments E[X \wedge u] and E[ (X-d)_+] (see Relation (2) and Relation (3)). Then apply the inflation effect on the first insurance benefit table. Note that the second table is obtained by inflating the first table by 1+r but with smaller upper limit u^* and deductible d^*. Also note the relation between Row E and Row F (same relation between Row A and Row B) and the relation between Row G and Row H (same relation between Row C and Row D).

\text{ }

Calculation

One approach of modeling insurance payments is to assume that the unmodified insurance losses are from a catalog of parametric distributions. For certain parametric distributions, the limited expectations E[X \wedge d] have convenient forms. Examples are: exponential distribution, Pareto distribution and lognormal distribution. Other distributions do not close form for E[X \wedge d] but the limited expectation can be expressed as a function that can be numerically calculated. One example is the gamma function.

The table in this link is an inventory of distributions that may be useful in modeling insurance losses. Included in the table are the limited expectations for various distributions. The following table shows the limited expectations for exponential, Pareto and lognormal.

Limited Expectation
Exponential E[X \wedge x]=\theta (1-e^{-x/\theta})
Pareto \displaystyle E[X \wedge x]=\frac{\theta}{\alpha-1} \biggl[1-\biggl(\frac{\theta}{x+\theta} \biggr)^{\alpha-1} \biggr], \alpha \ne 1
Lognormal \displaystyle E[(X \wedge x)^k]=\exp(k \mu+k^2 \sigma^2/2) \Phi \biggl(\frac{\log(x)-\mu-k \sigma^2}{\sigma} \biggr)+x^k \biggl[1-\Phi \biggl(\frac{\log(x)-\mu}{\sigma} \biggr) \biggr]

See the table in this link for the limited expectations of the distributed not shown in the above table. Once E[X \wedge d] is computed or estimated, the expected payment for various types of coverage provisions can be derived.

Examples

The first example demonstrates the four categories of expected payments in the table in the preceding section.

Example 1
Suppose that the loss distribution is described by the PDF \displaystyle f_X(x)=\frac{1}{5000} (100-x) with support 0<x<100. Walk through all the expected payments discussed in the above table (Rows A through D) assuming d=12 and u=60.

First, the basic calculation.

    \displaystyle E(X)=\int_0^{100} x \ \frac{1}{5000} (100-x) \ dx=\frac{100}{3}=33.3333

    \displaystyle F_X(x)=\frac{1}{5000} (100x- \frac{1}{2} x^2); \ \ \ 0<x<100

    F_X(12)=0.2256

    S_X(12)=1-0.2256=0.7744

    F_X(60)=0.84

    S_X(60)=1-0.84=0.16

    \displaystyle E(X \wedge 12)=\int_0^{12} x \ \frac{1}{5000} (100-x) \ dx +12 \cdot S_X(12)=\frac{6636}{625}

    \displaystyle E(X \wedge 60)=\int_0^{60} x \ \frac{1}{5000} (100-x) \ dx +12 \cdot S_X(60)=\frac{12864}{625}

Four categories of expected payments are derived and are shown in the following table.

A = Ordinary Deductible Only
B = Franchise Deductible Only
C = Ordinary Deductible + Limit
D = Franchise Deductible + Limit

Expected Cost Per Loss Expected Cost Per Payment
A \displaystyle \frac{100}{3}-\frac{6636}{625}=\frac{42592}{1875}=22.7157 \displaystyle e_X(d)=\frac{42592}{1875} \frac{1}{0.7744}=\frac{10648}{363}=29.3333
B \displaystyle \frac{42592}{1875}+12 \cdot 0.7744=\frac{60016}{1875}=32.0085 \displaystyle \frac{10648}{363}+12=\frac{15004}{363}=41.3333
C \displaystyle \frac{12864}{625}-\frac{6636}{625}=\frac{6228}{625}=9.9648 \displaystyle \frac{6228}{625} \frac{1}{0.7744}=\frac{6228}{484}=12.8678
D \displaystyle \frac{6228}{625}+12 \cdot 0.7744=\frac{12036}{625}=19.2576 \displaystyle \frac{6228}{484}+12=\frac{12036}{484}=24.8678

\text{ }

Example 2
Suppose that a coverage has an ordinary deductible and suppose that the CDF of the insurance payment per loss is given by:

    \displaystyle  G(y)=\left\{ \begin{array}{ll}                     \displaystyle  0 &\ \ \ \ y<0 \\           \text{ } & \text{ } \\           \displaystyle  1-e^{-\frac{y+20}{100}} &\ \ \ \ y \ge 0                      \end{array} \right.

Determine the expected value and the variance of the insurance cost per loss.

Note that there is a jump in the CDF at y=0. Thus y=0 is a point mass. The following is the PDF of the insurance payment.

    \displaystyle  g(y)=\left\{ \begin{array}{ll}                     \displaystyle  1-e^{-1/5} &\ \ \ \ y=0 \\           \text{ } & \text{ } \\           \displaystyle  \frac{1}{100} e^{-\frac{y+20}{100}}=e^{-1/5} \ \frac{1}{100} e^{-y/100} &\ \ \ \ y > 0                      \end{array} \right.

It is clear that g(y)=f(y+20) where f(y) is the PDF of the exponential distribution with mean 100. Thus g(y) is the PDF of the insurance payment per loss when the coverage has an ordinary deductible of 20 where the loss has exponential distribution with mean 100. We can use g(y) to calculate the mean and variance of Y. Using g(y), the mean and variance are:

    \displaystyle \begin{aligned} E(Y)&=\int_0^\infty y \ e^{-1/5} \ \frac{1}{100} e^{-y/100} \ dy \\&=e^{-1/5} \ \int_0^\infty y \ \frac{1}{100} e^{-y/100} \ dy \\&=e^{-1/5} \ (100)=81.873 \end{aligned}

    \displaystyle \begin{aligned} E(Y^2)&=\int_0^\infty y^2 \ e^{-1/5} \ \frac{1}{100} e^{-y/100} \ dy \\&=e^{-1/5} \ \int_0^\infty y^2 \ \frac{1}{100} e^{-y/100} \ dy \\&=e^{-1/5} \ 2(100)^2 \end{aligned}

    \displaystyle Var(Y)=e^{-1/5} \ 2(100)^2-\biggl( e^{-1/5} \ (100) \biggr)^2=9671.414601

Note that the calculation for E(Y) and E(Y^2) is done by multiplying e^{-1/5} by mean and second moment of the exponential distribution with mean 100.

\text{ }

Example 3
You are given the following:

  • Losses follow a Pareto distribution with parameters \alpha=3 and \theta=500.
  • The coverage has an ordinary deductible of 100.

Determine the PDF and CDF of the payment to the insured per payment. Determine the expected cost per payment.

Note that the setting of this example is identical to Example 3 in this previous post. Since this only concerns the insurance when the loss exceeds the deductible, the insurance payment is the conditional distribution Y=X-100 \lvert X>100 where X is the Pareto loss distribution.

The following gives the PDF and CDF and other information of the loss X.

    \displaystyle f_X(x)=\frac{3 \ 500^3}{(x+500)^4} \ \ \ \ \ \ \ \ \ \ \ \ \ \ x> 0

    \displaystyle F_X(x)=1-\biggl(\frac{500}{x+500} \biggr)^3 \ \ \ \ \ \ x> 0

    \displaystyle S_X(x)=\biggl(\frac{500}{x+500} \biggr)^3 \ \ \ \ \ \ \ \ \ \ \ \ x> 0

    \displaystyle S_X(100)=\biggl(\frac{500}{100+500} \biggr)^3=\biggl(\frac{500}{600} \biggr)^3=\frac{125}{216}

The following shows the PDF and CDF for the payment per payment variable Y.

    \displaystyle f_Y(y)=\frac{f_X(y+100)}{S_X(100)}=\frac{3 \ 500^3}{(y+100+500)^4} \ \biggl(\frac{600}{500} \biggr)^3=\frac{3 \cdot 600^3}{(y+600)^4} \ \ \ \ \ \ \ \ \ y> 0

    \displaystyle S_Y(y)=\frac{S_X(y+100)}{S_X(100)}=\biggl(\frac{500}{y+100+500} \biggr)^3 \ \biggl(\frac{600}{500} \biggr)^3=\biggl(\frac{600}{y+600} \biggr)^3 \ \ \ \ \ \ y> 0

    \displaystyle F_Y(y)=1-S_Y(y)=1-\biggl(\frac{600}{y+600} \biggr)^3 \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ y> 0

Note that these PDF and CDF are for a Pareto distribution with \alpha=3 and \theta=600. Thus the expected payment per payment and the variance are:

    \displaystyle E(Y)=\frac{600}{3-1}=300

    \displaystyle E(Y^2)=\frac{600^2 \cdot 2}{(3-1) (3-2)}=600^2

    Var(Y)=600^2 - 300^2=270000

Example 4
Suppose that losses X follow a Pareto distribution with shape parameter \alpha and scale parameter \theta. Suppose that an insurance coverage pays each loss subject to an ordinary deductible d. Derive a formula for the expected payment per payment e_X(d).

Let Y=X-d \lvert X>d be the payment per payment. In Example 3, we see that the CDF and PDF of Y are also Pareto but with shape parameter \alpha (same as for X) and scale parameter \theta+d (the original scale parameter plus the deductible). Thus the following gives the expected payment per payment.

    \displaystyle e_X(d)=\frac{\theta+d}{\alpha-1}=\frac{\theta}{\alpha-1}+\frac{1}{\alpha-1} d

Note that e_X(d) in this case is an increasing linear function of the deductible d. The higher the deductible, the larger the expected payment. This is a clear sign that the Pareto distribution is a heavy tailed distribution.

Example 5
Suppose that losses follow a lognormal distribution with parameters \mu=5 and \sigma=0.6. An insurance coverage has an ordinary deductible of 100. Compute the expected payment per loss. Compute the expected payment per loss if the deductible is a franchise deductible.

The following calculates E(X) and E[X \wedge 100].

    E(X)=e^{5+0.6^2/2}=e^{5.18}=177.682811

    \displaystyle \begin{aligned} E[X \wedge 100]&=e^{5+0.6^2/2} \ \Phi \biggl(\frac{\log(100)-5-0.6^2}{0.6} \biggr)+100 \ \biggl[1-\Phi \bigg(\frac{\log(100)-5}{0.6} \biggr) \biggr] \\&=e^{5.18} \ \Phi (-1.26)+100 \ [1-\Phi (-0.66 )]\\&=e^{5.18} \ 0.1038+100(0.7454)=e^{5.18} \ 0.1038+74.54=92.98347578 \end{aligned}

As a result, the expected payment per loss is E(X)-E[X \wedge 100]=84.69933522=84.70. If the deductible of 100 is a franchise deductible, the expected payment per loss is obtained by adding an additional insurance payment, which is d \ S(d). The following is the added payment.

    \displaystyle \begin{aligned} \text{added payment}&=100 \ \biggl[1-\Phi \bigg(\frac{\log(100)-5}{0.6} \biggr) \biggr] \\&=100 \ \biggl[1-\Phi (-0.66 ) \biggr] \\&=100 (0.7454) \\&=74.54  \end{aligned}

Thus the expected payment per loss with a franchise deductible of 250 is 84.70 + 74.54 = 159.24.

Practice Problems

The following practice problem sets are to reinforce the concepts discussed here.

Practice Problem Set 7
Practice Problem Set 8
Practice Problem Set 9

\text{ }

\text{ }

\text{ }

models of insurance payment

ordinary deductible, franchise deductible

maximum covered loss

Daniel Ma Math

Daniel Ma Mathematics

Dan Ma Actuarial

\copyright 2017 – Dan Ma

Advertisements

4 thoughts on “Mathematical models for insurance payments – part 3 – other modifications

    […] problem set is to reinforce the 3-part discussion on insurance payment models (Part 1, Part 2 and Part 3). The practice problems in this post are basic problems on calculating average insurance payment […]

    […] problem set is to reinforce the 3-part discussion on insurance payment models (Part 1, Part 2 and Part 3). The practice problems in this post are additional practice problems in addition to Practice […]

    […] problem set is to reinforce the 3-part discussion on insurance payment models (Part 1, Part 2 and Part 3). The practice problems in this post are additional practice problems in addition to Practice […]

    […] this blog post in a companion blog for more information on mean excess loss function and limited expectation […]

Leave a Reply

Fill in your details below or click an icon to log in:

WordPress.com Logo

You are commenting using your WordPress.com account. Log Out /  Change )

Google+ photo

You are commenting using your Google+ account. Log Out /  Change )

Twitter picture

You are commenting using your Twitter account. Log Out /  Change )

Facebook photo

You are commenting using your Facebook account. Log Out /  Change )

Connecting to %s