Actuarial Statistics · Credibility theory: Bayesian and empirical Bayes credibility premiums
Bayesian Credibility Premium: Formula, Derivation and Poisson-Gamma Example
Updated 11 October 2026 · Fact-checked
The Bayesian credibility premium is the predictive mean of the next claim, E[X_{n+1} | x1,...,xn], found using the posterior distribution of the risk parameter. For conjugate pairs such as Poisson-gamma it equals Z × sample mean + (1 − Z) × prior mean, where Z = n ÷ (n + k).
Understand Bayesian Credibility Premium
An insurer has a risk with an unknown risk parameter θ. Claims X1, X2, ... are independent given θ, but θ itself is random. This is why past claims of one risk tell you something about its future claims.
Before seeing data you hold a prior distribution for θ. After seeing claims x1,...,xn you update it to the posterior distribution using Bayes' theorem: posterior ∝ likelihood × prior. The Bayesian premium is the best estimate of the next claim under squared error loss. That best estimate is the predictive mean: E[X_{n+1} | x] = E[ E[X | θ] | x ] = E[μ(θ) | x], where μ(θ) = E[X | θ]. In words, you take the mean claim for a given θ and average it over the posterior of θ.
This premium is sometimes a credibility estimate: a weighted average Z × x̄ + (1 − Z) × E[μ(θ)], where x̄ is the sample mean and E[μ(θ)] is the prior mean claim. Z is the credibility factor, and it lies between 0 and 1. This form is exact for certain likelihood and conjugate prior pairs. The main ones are Poisson-gamma, exponential-gamma, normal-normal (known variance) and binomial-beta. For these, the posterior mean is linear in the data.
For other pairs the Bayesian premium is not linear in x̄. Then the credibility form is only an approximation. The Bühlmann premium is the best linear approximation to the Bayesian premium. If the Bayesian premium is exactly linear, it equals the Bühlmann premium.
Z rises towards 1 as n grows. More data means more weight on the risk's own experience and less on the prior. A more spread-out prior (more uncertainty about θ) also raises Z.
Key rules to remember
- Bayesian premium (predictive mean)
- E[X_{n+1} | x] = ∫ μ(θ) f(θ | x) dθ
- Use a sum for a discrete θ. μ(θ) = E[X | θ]. X_{n+1} is independent of the past given θ.
- Posterior density
- f(θ | x) ∝ f(θ) × Π f(xᵢ | θ)
- Drop every factor not involving θ, then recognise the distribution from its kernel.
- Credibility form
- P = Z × x̄ + (1 − Z) × E[μ(θ)]
- Valid exactly only when the posterior mean of μ(θ) is linear in the data, as in conjugate cases above.
- Poisson-gamma model
- X | λ ~ Poisson(λ), λ ~ Gamma(α, β) (mean α ÷ β). Posterior: Gamma(α + Σxᵢ, β + n)
- Here β is the rate parameter. Check which parametrisation the question uses.
- Poisson-gamma premium and Z
- P = (α + Σxᵢ) ÷ (β + n) = Z x̄ + (1 − Z)(α ÷ β), with Z = n ÷ (n + β)
- So k = β. Z is found by matching the coefficient of x̄.
- Normal-normal (known variance)
- X | θ ~ N(θ, σ²), θ ~ N(μ, s²). Z = n ÷ (n + σ² ÷ s²)
- Posterior mean = Z x̄ + (1 − Z) μ.
- Binomial-beta (m trials each)
- X | q ~ Bin(m, q), q ~ Beta(a, b). Posterior mean of q = (a + Σxᵢ) ÷ (a + b + nm)
- Premium per trial is this value, so Z = nm ÷ (a + b + nm).
How to solve Bayesian Credibility Premium questions
Use this method for any Bayesian credibility question, whether it asks for the premium, Z or a check of credibility form.
- 1Write the model: the distribution of X given θ, and the prior of θ. Note the parametrisation.
- 2Write μ(θ) = E[X | θ], the quantity you want to predict.
- 3Find the posterior: multiply likelihood by prior, keep only θ terms, and identify the distribution and its parameters.
- 4Compute the posterior mean of μ(θ). This is the Bayesian premium.
- 5Rewrite it as Z × x̄ + (1 − Z) × prior mean. Read Z from the coefficient of x̄.
- 6Check that Z is between 0 and 1 and that the prior mean equals E[μ(θ)] from the original prior.
- 7Substitute the data and state the premium with units. If asked, comment on how Z changes with n.
Quickest way: Conjugate shortcut
When to use it: Use the add-the-sum, add-the-count shortcut only for the Poisson-gamma pair. For exponential-gamma, binomial-beta and normal-normal, use the update rule given in step 2. Then ask for a premium or Z.
- For Poisson-gamma (rate β), the posterior is Gamma(α + Σxᵢ, β + n): add the data sum to α and the count to β.
- For the other pairs use their own rules. Exponential-gamma (X | θ ~ Exp(θ) with mean 1 ÷ θ, θ ~ Gamma(α, β) rate): posterior Gamma(α + n, β + Σxᵢ), and the premium is (β + Σxᵢ) ÷ (α + n − 1), which needs α + n > 1. Binomial-beta: posterior Beta(a + Σxᵢ, b + nm − Σxᵢ). Normal-normal: do not add counts. The posterior mean is the precision-weighted mean Z x̄ + (1 − Z) μ with Z = n ÷ (n + σ² ÷ s²).
- Write the posterior mean as a ratio and divide top and bottom to isolate x̄.
- Read Z = n ÷ (n + k). For Poisson-gamma, k = β (rate). For exponential-gamma, k = α − 1. For binomial-beta, Z = nm ÷ (a + b + nm).
- Sanity check: with n = 0 you must get the prior mean, and as n grows the premium moves to x̄.
Common mistakes in Bayesian Credibility Premium
Using the posterior mean of θ instead of the posterior mean of μ(θ).
In Poisson-gamma μ(λ) = λ, so the two agree, and students assume that is always true.
Fix: Always write μ(θ) = E[X | θ] first. For exponential with mean 1 ÷ θ you need E[1 ÷ θ | x], not E[θ | x].
Mixing up the gamma rate and scale parameters.
Texts and questions use different conventions.
Fix: Check the stated mean. If the prior mean is α ÷ β, β is a rate and k = β. If it is αθ, convert.
Taking Z = n ÷ (n + α).
Students remember a shape parameter and not the correct one.
Fix: Derive it: (α + nx̄) ÷ (β + n) gives Z = n ÷ (n + β).
Claiming every Bayesian premium is a credibility premium.
The textbook examples are all conjugate and linear.
Fix: State that the form is exact only when the posterior mean is linear in the data. Otherwise it is only approximated by the Bühlmann premium.
Using the sum of claims where the mean is needed in Z × x̄.
The posterior formula contains Σxᵢ, so the sum gets plugged into the credibility form.
Fix: Use x̄ = Σxᵢ ÷ n in the credibility form. Check Z is below 1.
Forgetting the prior mean, or using the posterior mean, in the (1 − Z) term.
Confusion between prior and posterior.
Fix: The complement weight always multiplies the prior mean of μ(θ), here α ÷ β.
Worked examples
Example 1
Claim numbers per year for a risk are Poisson with mean λ. The prior for λ is gamma with α = 3 and rate β = 6 (prior mean 0.5). Over 4 years the claims were 1, 0, 2, 1. Find the Bayesian premium for year 5 and the credibility factor Z.
Show the solution
- Σxᵢ = 1 + 0 + 2 + 1 = 4, n = 4, x̄ = 1.
- Posterior is gamma with α' = 3 + 4 = 7 and β' = 6 + 4 = 10.
- μ(λ) = λ, so the premium is the posterior mean = 7 ÷ 10 = 0.7.
- Z = n ÷ (n + β) = 4 ÷ 10 = 0.4.
- Check: 0.4 × 1 + 0.6 × 0.5 = 0.4 + 0.3 = 0.7.
Answer: Premium = 0.7 claims; Z = 0.4.
Example 2
Annual claim counts for a risk are Poisson with mean λ, where λ has a gamma prior with mean 2 and variance 0.5 (shape α, rate β). How many years of data are needed for Z to be at least 0.8?
Show the solution
- For gamma with rate β: mean = α ÷ β = 2 and variance = α ÷ β² = 0.5.
- Divide: (α ÷ β²) ÷ (α ÷ β) = 1 ÷ β = 0.5 ÷ 2 = 0.25, so β = 4. Then α = 8.
- Z = n ÷ (n + 4) ≥ 0.8.
- n ≥ 0.8n + 3.2, so 0.2n ≥ 3.2, n ≥ 16.
Answer: At least 16 years of data are needed (Z = 16 ÷ 20 = 0.8 exactly at n = 16).
Exam tips
- Write the model, μ(θ) and the posterior in separate lines. Marks are given for each stage even if a number slips.
- If asked to show the premium is a credibility estimate, finish by naming Z and the prior mean explicitly.
- In written derivations, show the kernel step: drop constants and recognise the gamma, normal or beta form.
- Always verify with n = 0 and with the credibility form to catch parametrisation errors.
- Expect MCQs asking for Z or the premium from given prior parameters. Check whether the prior is given by rate, scale, or mean and variance.
Practice questions from Credibility theory: Bayesian and empirical Bayes credibility premiums
- In a Bühlmann model, the expected process variance is 300 and the variance of hypothetical means is 100. How many years of experience n are …
- Claim counts for a policyholder in one year are Poisson with mean λ. The prior for λ is Gamma with shape α = 3 and rate β = 6. The policyhol…
- Empirical Bühlmann estimation uses data on 3 risks, each with 4 years of claims. Risk means are 10, 14 and 12. Within-risk sample variances …
- Given θ, claim sizes are Exponential with mean 1/θ. The prior for θ is Gamma(3, 600). Four claims total 1,000 rupees (in thousands, use the …
- Four risks each have three years of data. The risk means are 20, 24, 16 and 28. The sums of squared deviations of the observations from thei…
Bayesian Credibility Premium: frequently asked questions
Why is the Bayesian premium the predictive mean?
Under squared error loss, the best estimate of a future value is its conditional expectation given the data. That is the mean of the predictive distribution of X_{n+1}. It equals the posterior expectation of μ(θ).
When is the Bayesian estimate exactly a credibility estimate?
When the posterior mean of μ(θ) is a linear function of the sample mean. This happens for conjugate pairs from the exponential family such as Poisson-gamma, normal-normal and binomial-beta. For other models it is not exact.
How do I calculate Z for the Poisson-gamma model?
Write the premium as (α + Σxᵢ) ÷ (β + n) and divide through by n to match Z x̄ + (1 − Z) α ÷ β. This gives Z = n ÷ (n + β), where β is the gamma rate.
How is this different from the Bühlmann model?
The Bayesian premium uses a full prior and the exact posterior. The Bühlmann premium only uses the prior mean, the expected process variance and the variance of the hypothetical means. The two coincide in the exact-credibility cases.