Tutorials › CAS Exam 8 - Advanced Ratemaking › Choosing Prior Parameters From Actuarial Information

CAS Exam 8 · Bayesian Credibility Foundations · Tutorial 433 of 989

Choosing Prior Parameters From Actuarial Information

Translate actuarial assumptions about a claim rate or probability into prior-distribution parameters, then check that the parameters reproduce the intended mean and dispersion.

Intermediate 11 min read

What You'll Learn

  • Match Gamma shape and rate parameters to an assumed claim-rate mean and variance.
  • Use a coefficient of variation to express and calibrate Gamma prior dispersion.
  • Match Beta parameters to an assumed claim probability mean and variance.
  • Check whether proposed Beta moments are mathematically feasible.
  • Detect shape-rate versus shape-scale errors by verifying the implied moments.

Why Prior Parameters Need Calibration

Actuarial information about heterogeneous risks is often expressed in practical terms: an average claim frequency, a range of plausible risk levels, a standard deviation, or a coefficient of variation. A Bayesian model, however, requires parameters for a probability distribution. Choosing those parameters is a moment-matching task: select a distributional family, translate the supplied assumptions into its parameters, and verify that the resulting distribution has the intended moments.

This tutorial considers two common choices. A Gamma prior can describe a positive, uncertain claim rate, while a Beta prior can describe an uncertain probability between zero and one. The focus is calibration of the prior distribution, not updating it with observed claims. The likelihood and posterior calculations are separate steps.

Definition: Prior calibration by moment matching means choosing distribution parameters so that the distribution's theoretical mean and variance equal actuarial assumptions for the uncertain quantity. The assumed moments must be in the range the selected distribution can attain.

Gamma Priors for Claim Rates

Let \(\lambda\) be an unknown claim rate, such as expected claims per policy-year. Use a Gamma distribution with shape parameter \(\alpha>0\) and rate parameter \(\beta>0\), denoted \(\lambda\sim\operatorname{Gamma}(\alpha,\beta)\). In this convention, the density is proportional to \(\lambda^{\alpha-1}e^{-\beta\lambda}\), and the rate parameter \(\beta\), not its reciprocal, appears in the exponential.

$$ E[\lambda]=\frac{\alpha}{\beta}, \qquad \operatorname{Var}(\lambda)=\frac{\alpha}{\beta^2}. $$

Suppose actuarial information supplies a positive prior mean \(m\) and variance \(s^2\). Equating the Gamma moments to these assumptions gives

$$ m=\frac{\alpha}{\beta}, \qquad s^2=\frac{\alpha}{\beta^2} \quad\Longrightarrow\quad \alpha=\frac{m^2}{s^2}, \qquad \beta=\frac{m}{s^2}. $$

The equations have a unique positive solution when \(m>0\) and \(s^2>0\). A standard deviation can be used after squaring it to obtain the variance. If the assumptions instead give a coefficient of variation \(c=s/m\), then \(s^2=c^2m^2\), so the same formulas reduce to \(\alpha=1/c^2\) and \(\beta=1/(c^2m)\).

Worked Example: Calibrating a Gamma Prior From a Rate and Standard Deviation

A hypothetical portfolio has an uncertain claim frequency \(\lambda\), measured in claims per policy-year. Actuarial assumptions specify \(E[\lambda]=0.08\) and a prior standard deviation of \(0.04\) claims per policy-year. Find the Gamma shape and rate parameters. Assume no interest discounting or mortality component; this is a frequency model, and the prior family is Gamma.

Step 1: Convert the standard deviation to a variance.

$$ m=0.08,\qquad s=0.04,\qquad s^2=(0.04)^2=0.0016. $$

Step 2: Match the mean and variance. For the shape-rate convention,

$$ \alpha=\frac{m^2}{s^2} =\frac{(0.08)^2}{0.0016} =\frac{0.0064}{0.0016} =4, \qquad \beta=\frac{m}{s^2} =\frac{0.08}{0.0016} =50. $$

Thus \(\lambda\sim\operatorname{Gamma}(4,50)\), with the second parameter interpreted as a rate. The verification is \(E[\lambda]=4/50=0.08\) and \(\operatorname{Var}(\lambda)=4/50^2=4/2500=0.0016\), whose square root is \(0.04\). Both assumed moments are reproduced.

Method and trap avoided. This is direct moment matching. Calling \(0.02=s^2/m\) the rate would confuse the scale \(s^2/m\) with the rate \(m/s^2\); the implied mean would then be \(4/0.02=200\), not \(0.08\).

Why the Gamma Formulas Work

The moment-matching equations are obtained by solving the distribution's moment equations, rather than by treating shape and rate as informal labels. For a Gamma variable in the stated shape-rate convention, its mean is \(\alpha/\beta\) and its variance is \(\alpha/\beta^2\). The mean equation gives \(\alpha=m\beta\). Substitute this into the variance equation:

$$ s^2=\frac{m\beta}{\beta^2}=\frac{m}{\beta} \quad\Longrightarrow\quad \beta=\frac{m}{s^2}. $$

Then \(\alpha=m\beta=m^2/s^2\). This also explains the role of dispersion: at a fixed mean, a larger variance produces a smaller shape and a smaller rate. Both parameters change, but their ratio continues to preserve the assumed mean.

Units provide a useful independent check. If \(\lambda\) is measured in claims per policy-year, \(\beta\) has units of policy-years per claim-rate unit, so that \(\alpha/\beta\) has the units of \(\lambda\). Shape \(\alpha\) is dimensionless. If exposure is rescaled—for example, from policy-years to thousands of policy-years—the numerical rate parameter must be consistent with the new numerical unit for \(\lambda\).

Worked Example: Using a Coefficient of Variation for a Gamma Prior

For a hypothetical commercial vehicle segment, the uncertain expected frequency is \(\lambda\) claims per 1,000 vehicle-years. The prior mean is 30 claims per 1,000 vehicle-years and the coefficient of variation is \(1/3\). Find the shape-rate parameters and check the moments. Assume no interest discounting or mortality component; the modeled quantity is a claim rate, and the prior family is Gamma.

Step 1: Recover the standard deviation and variance. The coefficient of variation is standard deviation divided by mean, so

$$ s=\frac{1}{3}(30)=10, \qquad s^2=10^2=100. $$

Step 2: Calculate the parameters.

$$ \alpha=\frac{30^2}{100}=\frac{900}{100}=9, \qquad \beta=\frac{30}{100}=0.3. $$

The calibrated prior is \(\lambda\sim\operatorname{Gamma}(9,0.3)\), with \(\lambda\) measured numerically in claims per 1,000 vehicle-years. The check gives \(E[\lambda]=9/0.3=30\) and \(\operatorname{Var}(\lambda)=9/(0.3)^2=9/0.09=100\), so its standard deviation is 10 and its coefficient of variation is \(10/30=1/3\).

Method and trap avoided. The coefficient of variation is relative dispersion, not variance. Squaring \(1/3\) gives the squared coefficient of variation \(1/9\), and \(s^2=(1/9)(30^2)=100\). Using \(1/3\) directly as a variance would not reproduce the intended dispersion.

Beta Priors for Claim Probabilities

Let \(p\) be an unknown claim probability, such as the probability that a hypothetical insured has at least one claim during a defined period. A Beta prior has parameters \(a>0\) and \(b>0\), written \(p\sim\operatorname{Beta}(a,b)\). Its support is \(0<p<1\), and its moments are

$$ E[p]=\frac{a}{a+b}, \qquad \operatorname{Var}(p)=\frac{ab}{(a+b)^2(a+b+1)}. $$

If the assumed mean is \(m\) and variance is \(s^2\), define the concentration \(k=a+b\). Since \(a=mk\) and \(b=(1-m)k\), substitution into the variance gives

$$ s^2=\frac{m(1-m)}{k+1} \quad\Longrightarrow\quad k=\frac{m(1-m)}{s^2}-1, \qquad a=mk,\quad b=(1-m)k. $$

For a non-degenerate Beta distribution, \(0<m<1\) and \(0<s^2<m(1-m)\). The upper bound follows because \(k\) must be positive; equality would require \(k=0\), which is not an allowed Beta parameter sum. If the assumed mean is at an endpoint, or the variance is zero, a point mass rather than an ordinary Beta distribution is needed. If the variance violates the bound, the assumptions cannot be represented by a Beta prior.

Worked Example: Calibrating a Beta Prior From Probability Assumptions

A hypothetical underwriting study treats the probability \(p\) that a small contractor has a claim during one year as uncertain. The assumed prior mean is \(0.10\), and the prior standard deviation is \(0.06\). Find the Beta parameters and verify both moments. Assume no interest discounting or mortality component; \(p\) is a one-year event probability, and the prior family is Beta.

Step 1: Find the variance and check feasibility.

$$ m=0.10,\qquad s^2=(0.06)^2=0.0036,\qquad m(1-m)=0.10(0.90)=0.09. $$

Since \(0<0.0036<0.09\), these moments are feasible for a non-degenerate Beta distribution.

Step 2: Calculate the concentration and parameters.

$$ k=\frac{0.09}{0.0036}-1=25-1=24, \qquad a=0.10(24)=2.4, \qquad b=0.90(24)=21.6. $$

Thus \(p\sim\operatorname{Beta}(2.4,21.6)\). The mean check is \(2.4/(2.4+21.6)=2.4/24=0.10\). The variance check is

$$ \operatorname{Var}(p) =\frac{(0.10)(0.90)}{24+1} =\frac{0.09}{25} =0.0036, $$

so the standard deviation is \(\sqrt{0.0036}=0.06\). Both assumptions are matched.

Method and trap avoided. The calculation first uses the variance, not the standard deviation, and then solves for \(a+b\). Setting \(a+b\) equal to \(m(1-m)/s^2\) without subtracting one would give 25 rather than 24 and would produce the wrong variance.

Feasibility Checks and Common Exam Traps

Parameter calculations are not complete until the resulting distribution is checked against the original actuarial assumptions. For Gamma calibration, verify that the mean and variance are positive and use a single parameter convention consistently. For Beta calibration, check the feasibility bound before calculating \(a\) and \(b\). These checks catch impossible inputs as well as arithmetic and interpretation errors.

  • Confusing rate with scale. In the Gamma shape-rate convention used here, \(\beta=m/s^2\). The scale is \(s^2/m\), its reciprocal. Always substitute the proposed parameters into the stated moment equations.
  • Using a standard deviation where a variance is required. Both Gamma and Beta formulas use \(s^2\). Square the supplied standard deviation first.
  • Ignoring probability bounds. For a Beta prior, calculate \(m(1-m)\) and check that the assumed variance is strictly smaller. A mathematically infeasible pair of moments cannot be fixed by rounding parameters.
  • Changing units without changing interpretation. A numerical claim rate per policy-year is not the same numerical quantity as claims per 1,000 policy-years. State the exposure unit alongside the rate and parameter values.
  • Stopping after solving for parameters. Substitute them back into the distribution's mean and variance formulas. This final check is part of the calibration, not an optional presentation detail.

For example, a proposed Beta prior with mean \(0.08\) and standard deviation \(0.30\) is impossible: its variance is \(0.09\), whereas \(m(1-m)=0.08(0.92)=0.0736\), and a proper Beta distribution requires variance strictly below \(0.0736\). The problem is not a computational mistake; the requested moments are incompatible with the chosen family.

Key takeaway: Match the assumed mean and variance to the moments of the chosen prior, respect the distribution's parameter convention and feasibility conditions, and verify the resulting moments before using the prior in a model. This calibration skill is aligned with the CAS Exam 8 (Advanced Ratemaking) syllabus.

Check Your Understanding

For each problem, identify the distribution and parameter convention before calculating. Verify feasibility and units where applicable.

  1. A hypothetical liability segment has claim rate mean \(0.05\) claims per policy-year and prior standard deviation \(0.025\) claims per policy-year. Find the Gamma shape and rate. (Answer: shape \(4\), rate \(80\).)
  2. A hypothetical property portfolio's claim rate per 1,000 insured-years has mean 18 and coefficient of variation \(0.5\). Find the Gamma shape and rate. (Answer: shape \(4\), rate \(2/9\).)
  3. A hypothetical driver group's one-year claim probability has mean \(0.40\) and standard deviation \(0.12\). Find the Beta parameters. (Answer: \(a=94/15\approx6.2667\), \(b=47/5=9.4\).)
  4. A proposed Beta prior for a hypothetical coverage has mean \(0.08\) and standard deviation \(0.30\). Decide whether a proper Beta distribution can match these moments. (Answer: No; the variance \(0.09\) exceeds \(m(1-m)=0.0736\).)
  5. A hypothetical portfolio has claim rate mean \(0.06\) claims per policy-year and standard deviation \(0.03\) claims per policy-year. Find the Gamma shape-rate parameters and verify the mean. (Answer: shape \(4\), rate \(200/3\approx66.6667\); \(4/(200/3)=0.06\).)
Clinical Biostats