Tutorials › CAS Exam 8 - Advanced Ratemaking › Bühlmann Credibility for Claim Frequencies

CAS Exam 8 · Bühlmann Credibility Factor · Tutorial 362 of 989

Bühlmann Credibility for Claim Frequencies

Apply Bühlmann credibility to claim frequencies, estimate its variance components from a balanced dataset, and calculate a credibility-weighted frequency in consistent exposure units.

Intermediate 10 min read

What You'll Learn

  • Define claim frequency as claims per stated, consistent exposure unit.
  • Identify when the standard Bühlmann factor applies to frequency observations.
  • Estimate EPV and VHM from a balanced set of repeated frequency observations.
  • Calculate \(K\), the credibility factor \(Z\), and a weighted frequency estimate.
  • Convert between claim counts and frequencies without changing exposure units.
  • Recognize when unequal exposure makes the simple \(n/(n+K)\) formula inapplicable.

Why Apply Credibility to Claim Frequency?

A risk’s claim frequency can vary from one observation period to another, even when its underlying tendency to generate claims is stable. A short run of unusually high or low frequency should therefore not automatically determine the estimate for future periods. Bühlmann credibility balances the risk’s observed frequency against the collective average: more consistent, relevant experience receives more weight, while the collective mean supplies a complement.

Here, frequency means the number of claims divided by a stated exposure amount, such as claims per policy-year or per 1,000 vehicle-years. The exposure basis must be specified and used consistently. The standard Bühlmann factor \(Z=n/(n+K)\) applies when the \(n\) observations are comparable frequency measurements with the same exposure basis and model variance structure. It is not a formula for weighting raw claim counts from periods with different exposure.

Definition: Let \(X_{it}\) be the claim frequency for risk \(i\) in observation period \(t\), measured on a common exposure basis. In the Bühlmann model, each risk has an unobserved risk parameter \(\Theta_i\). Conditional on \(\Theta_i\), the observations have mean \(m(\Theta_i)=E[X_{it}\mid\Theta_i]\), variance \(s^2(\Theta_i)=\operatorname{Var}(X_{it}\mid\Theta_i)\), and are independent with the same conditional distribution. The EPV is \(E[s^2(\Theta)]\), the expected process variance. The VHM is \(\operatorname{Var}(m(\Theta))\), the variance of hypothetical means across risks. The collective mean is \(\mu=E[m(\Theta)]=E[X]\).

For a risk with \(n\) observations, let \(\bar X\) be its arithmetic mean frequency. When EPV and VHM are positive, define \(K=\mathrm{EPV}/\mathrm{VHM}\). The Bühlmann credibility-weighted estimate of that risk’s expected frequency is \(Z\bar X+(1-Z)\mu\), where \(Z=n/(n+K)\). The estimate is a frequency, not a prediction that the next period will contain a fractional claim count.

$$ K=\frac{\mathrm{EPV}}{\mathrm{VHM}},\qquad Z=\frac{n}{n+K},\qquad \widehat{m}=Z\bar X+(1-Z)\mu $$

The units matter. If \(X\) is claims per 1,000 policy-years, then \(\mu\), \(\bar X\), EPV, and VHM all refer to observations on that frequency basis. \(K\) is measured in observation-period units, consistent with the definition of \(n\). A calculated \(\widehat m\) can be converted to an expected claim count for a specified future exposure by multiplying by that exposure in the matching units.

Why the Credibility Weight Has This Form

The formula is the best linear estimate of the hypothetical mean under squared-error loss. Consider estimates of the form \(a\bar X+(1-a)\mu\), where \(a\) is a weight between zero and one. The risk-specific target is \(m(\Theta)\). Conditional independence gives \(\operatorname{Var}(\bar X\mid\Theta)=s^2(\Theta)/n\). Averaging over risks gives \(\operatorname{Var}(\bar X-m(\Theta))=\mathrm{EPV}/n\); the difference between the collective mean and the risk-specific mean has variance VHM. The cross-product has expectation zero because the conditional sampling error has mean zero given \(\Theta\).

Thus the mean squared error of the linear estimate is

$$ E\left[\left\{a\bar X+(1-a)\mu-m(\Theta)\right\}^2\right] =a^2\frac{\mathrm{EPV}}{n}+(1-a)^2\mathrm{VHM}. $$

Differentiate with respect to \(a\) and set the derivative equal to zero:

$$ 2a\frac{\mathrm{EPV}}{n}-2(1-a)\mathrm{VHM}=0. $$

Rearranging gives \(a(\mathrm{EPV}/n+\mathrm{VHM})=\mathrm{VHM}\), so

$$ a=\frac{\mathrm{VHM}}{\mathrm{VHM}+\mathrm{EPV}/n} =\frac{n}{n+\mathrm{EPV}/\mathrm{VHM}} =\frac{n}{n+K}. $$

The second derivative of the mean squared error is \(2(\mathrm{EPV}/n+\mathrm{VHM})\), which is positive when the variance components are positive. This establishes that the weight is the minimizing linear weight. The result does not say that the weighted estimate is the exact posterior mean for every possible distribution of risk parameters; it is the Bühlmann linear credibility estimate.

Worked Calculations for Frequency Experience

The following examples use hypothetical insurance portfolios. For the standard factor, each observation is one comparable period on a fixed exposure basis. When data are used to estimate variance components, the balanced setup has the same number of periods for each risk.

Worked Example: Weighting Eight Years of Claim Frequency

A hypothetical insurer observes one risk for eight years, with one policy-year of exposure in each year. The annual claim counts are \(0,1,2,1,0,1,1,2\), so each count is also the frequency per policy-year. Assume the collective mean is \(0.90\) claims per policy-year, EPV is \(0.24\), and VHM is \(0.06\). Estimate the risk’s expected frequency. Name the method and identify the common error to avoid.

Method and units. Apply the Bühlmann frequency estimate. Each annual observation is on the same one-policy-year basis, so \(n=8\). First calculate the risk’s sample mean:

$$ \bar X=\frac{0+1+2+1+0+1+1+2}{8} =\frac{8}{8}=1.00\text{ claim per policy-year}. $$

Calculate \(K\) and \(Z\).

$$ K=\frac{0.24}{0.06}=4,\qquad Z=\frac{8}{8+4}=\frac{8}{12}=\frac{2}{3}. $$

Credibility-weight the frequency. The complement weight is \(1-Z=1/3\), so

$$ \widehat m=\frac{2}{3}(1.00)+\frac{1}{3}(0.90) =0.6667+0.3000=0.9667 $$

The estimate is \(0.9667\) claims per policy-year. As a check, it lies between the observed mean of \(1.00\) and collective mean of \(0.90\), closer to the observed mean because \(Z=2/3\). The common error is to use \(K=0.24/0.06\) incorrectly as the credibility weight; \(K=4\) is the volume parameter, and it must be used in \(n/(n+K)\).

Worked Example: Estimating EPV and VHM from a Balanced Frequency Dataset

A hypothetical portfolio has four risks, each observed for three equal-exposure periods. The frequencies per exposure unit are shown below. Estimate EPV and VHM using the balanced method-of-moments formulas, then calculate \(K\), \(Z\), and the credibility-weighted frequency for Risk B. Assume the frequency observations satisfy the Bühlmann model.

RiskPeriod 1Period 2Period 3Risk mean
A010
B112
C212
D323

Step 1: Calculate the risk means and collective mean. The four risk totals are \(1,4,5,8\), each divided by three periods. Thus the risk means are \(1/3,4/3,5/3,8/3\). The total of all 12 observations is \(18\), so the grand mean is \(18/12=1.50\).

Step 2: Estimate EPV. For each risk, the sum of squared deviations from its own mean is \(2/3\). For example, Risk A has squared deviations \((0-1/3)^2+(1-1/3)^2+(0-1/3)^2=1/9+4/9+1/9=2/3\). The other three risks have the same sum. The pooled within-risk estimate divides the total within-risk sum of squares by \(J(n-1)=4(3-1)=8\):

$$ \widehat{\mathrm{EPV}} =\frac{4(2/3)}{8} =\frac{1}{3}. $$

Step 3: Estimate VHM. The squared deviations of the four risk means from \(1.50\) are \(49/36,1/36,1/36,49/36\), which sum to \(100/36=25/9\). Their sample variance is \((25/9)/(4-1)=25/27\). For three observations per risk, subtract the EPV contribution to the variance of a risk mean:

$$ \widehat{\mathrm{VHM}} =\frac{25}{27}-\frac{1/3}{3} =\frac{25}{27}-\frac{1}{9} =\frac{22}{27}. $$

Step 4: Calculate the estimate for Risk B. Its observed mean is \(4/3\), and the collective mean is \(3/2\). Therefore,

$$ \widehat K=\frac{1/3}{22/27}=\frac{9}{22}, \qquad \widehat Z=\frac{3}{3+9/22}=\frac{66}{75}=0.88, $$ $$ \widehat m_B=0.88\left(\frac{4}{3}\right)+0.12\left(\frac{3}{2}\right) =\frac{88}{75}+\frac{9}{50} =\frac{203}{150} \approx 1.3533. $$

Risk B’s estimated frequency is \(1.3533\) per exposure unit. It lies between Risk B’s mean \(1.3333\) and the collective mean \(1.50\). The common error is to calculate the between-risk variance and call it VHM without correction. A risk’s observed mean still contains process variation, here \(\widehat{\mathrm{EPV}}/3=1/9\), which must be removed.

Balanced-data result: With \(J\) risks and \(n\) equal-exposure observations per risk, let \(\bar X_i\) be risk \(i\)'s mean, \(\bar X\) the grand mean, and \(s_i^2\) the sample variance within risk \(i\). The usual method-of-moments estimates are \(\widehat{\mathrm{EPV}}=\sum_i(n-1)s_i^2/[J(n-1)]\) and \(\widehat{\mathrm{VHM}}=s_{\bar X}^2-\widehat{\mathrm{EPV}}/n\), where \(s_{\bar X}^2=\sum_i(\bar X_i-\bar X)^2/(J-1)\). These formulas rely on the balanced, common-model setup; they are not a general recipe for unequal exposure or changing frequency definitions.

Worked Example: Converting Counts to a Standardized Frequency

A hypothetical commercial program records a risk’s claim counts over four periods, each with 2,000 policy-years of exposure: \(3,1,4,2\). The frequency basis is claims per 1,000 policy-years. Assume the collective frequency is \(1.10\) per 1,000 policy-years, EPV is \(0.16\), and VHM is \(0.04\), both defined for one 2,000-policy-year observation after conversion to the stated frequency basis. Estimate the risk’s expected frequency and its expected claims over a future 2,000-policy-year period.

Method and exposure conversion. Each observation contains two units of 1,000 policy-years, so divide each count by 2 to obtain its frequency per 1,000 policy-years. The four frequencies are \(1.5,0.5,2.0,1.0\). Their mean is

$$ \bar X=\frac{1.5+0.5+2.0+1.0}{4}=\frac{5.0}{4}=1.25. $$

Calculate credibility. There are four comparable observation periods, so \(n=4\). Then

$$ K=\frac{0.16}{0.04}=4,\qquad Z=\frac{4}{4+4}=0.50. $$

Estimate frequency and count.

$$ \widehat m=0.50(1.25)+0.50(1.10) =0.625+0.550=1.175 $$

The estimated frequency is \(1.175\) claims per 1,000 policy-years. A future 2,000-policy-year period contains two units of 1,000 policy-years, so its estimated claim count is \(1.175(2)=2.35\). This is an expected count, not a claim-count guarantee. The common mistake is to treat the raw counts \(3,1,4,2\) as frequencies per 1,000 policy-years; doing so doubles the observed rate and mixes units with the collective mean.

Conditions, Interpretation, and Exam Traps

The frequency application uses the same credibility structure as other Bühlmann applications, but the observation must be defined carefully. For instance, counts recorded over unequal exposure periods do not have the same variance merely because each is one calendar year. In a conditional Poisson model, a count over exposure \(e\) has conditional mean and variance proportional to \(e\); dividing by \(e\) changes the variance with exposure. Consequently, applying \(n/(n+K)\) to unadjusted counts, or to exposure-normalized rates with materially different exposure without accounting for their differing variances, is not justified by the equal-observation formula. Use comparable exposure units or an appropriate unequal-exposure method.

Other assumptions deserve attention: the risk’s underlying frequency tendency is treated as stable over the observations; periods are conditionally independent; and the collective mean and variance components describe the same portfolio basis as the risk being estimated. A trend, coverage change, or shift in claim definition can make historical frequency observations non-comparable. The credibility calculation does not correct such changes automatically.

The weighted estimate is a convex combination when \(0\leq Z\leq1\), so it lies between the observed mean and collective mean. As \(n\) grows with \(K\) fixed, \(Z\) increases and the risk’s own experience has greater influence. If VHM is zero, there is no modeled variation in hypothetical means across risks, and the ratio defining \(K\) is not usable in its ordinary form; the frequency data provide no between-risk differentiation under that model. If variance-component estimates are unstable or outside plausible bounds, investigate their estimation and modeling assumptions rather than treating the resulting factor as self-validating.

  • Do not confuse frequency with count. A count is tied to its exposure; frequency is count divided by stated exposure.
  • Do not confuse \(K\) with \(Z\). \(K\) is EPV divided by VHM; \(Z=n/(n+K)\) is the weight on the observed mean.
  • Do not omit the VHM correction. The variation among observed risk means includes process variation in those means.
  • Do not mix units. The observed mean, collective mean, and variance components must correspond to the same defined observation basis.

The practical purpose is to stabilize a risk’s estimated frequency without discarding its information. In ratemaking, frequency can contribute to a broader estimate of expected claim cost, but this calculation alone does not supply claim severity, expenses, or a final premium. The frequency estimate is one model component, and its interpretation depends on the exposure definition and assumptions stated above. This tutorial is aligned with the CAS Exam 8 (Advanced Ratemaking) syllabus.

Check Your Understanding

For each problem, use the stated frequency basis and assumptions. Answers are given for checking.

  1. A risk has six equal-exposure frequency observations with mean \(1.4\). The collective mean is \(1.2\), EPV is \(0.30\), and VHM is \(0.10\). Calculate \(K\), \(Z\), and the credibility-weighted frequency. (Answer: \(K=3\), \(Z=2/3\), estimate \(1.3333\) per exposure unit.)
  2. A hypothetical risk has four equal periods, each with 500 policy-years of exposure, and counts \(2,0,1,3\). Express frequency per 500 policy-years. If EPV is \(0.50\), VHM is \(0.25\), and the collective mean is \(1.20\) per 500 policy-years, calculate the weighted estimate. (Answer: observed mean \(1.50\), \(K=2\), \(Z=2/3\), estimate \(1.40\) claims per 500 policy-years.)
  3. A risk’s five comparable frequency observations are \(0.8,1.2,1.0,0.6,1.4\). Given collective mean \(0.90\), EPV \(0.09\), and VHM \(0.03\), calculate the weighted expected frequency. (Answer: observed mean \(1.00\), \(K=3\), \(Z=5/8=0.625\), estimate \(0.9625\) per exposure unit.)
  4. A program’s three equal periods each have 6,000 policy-years of exposure, with counts \(6,9,3\). Frequency is reported per 1,000 policy-years. Given collective mean \(1.20\), EPV \(0.12\), and VHM \(0.04\) for these standardized period rates, calculate the frequency estimate and expected count over 6,000 future policy-years. (Answer: period rates \(1,1.5,0.5\), mean \(1.00\), \(K=3\), \(Z=1/2\), estimate \(1.10\) per 1,000 policy-years, expected count \(6.60\).)
Clinical Biostats