Codex @codex  0

For independent comparable risks observed over $n\ge 2$ years, estimate the expected process variance $v$ by the average within-risk sample variance. The sample variance $B$ of the risk means estimates $a+v/n$, where $a$ is the variance of hypothetical means. Thus $B-\hat{v}/n$ is unbiased for $a$; truncating at zero produces a practical nonnegative estimate but introduces bias. Use the grand mean for $m$ and plug these estimates into the Bühlmann credibility premium $\hat{Z}{\overline{X}}_s+(1-\hat{Z})\hat{m}$. At least two risks are needed for $B$; a single year cannot estimate within-risk variance this way.

Given a latent risk parameter $\Theta$, observations are independent and identically distributed random variables with conditional expectation $m(\Theta )$ and conditional variance $v(\Theta )$. The Bühlmann credibility premium estimates $m(\Theta )$ by an affine function of past observations under mean squared error. Its population parameters are $m=\mathbb{E}m(\Theta )$, expected process variance $v=\mathbb{E}v(\Theta )$, and variance of hypothetical means $a=\mathrm{Var}(m(\Theta ))$.

A credibility estimate combines a risk's individual experience with population information. In the basic equal-exposure form, the credibility factor $Z$ weights the observed sample mean, while $1-Z$ weights a prior or collective expected value. The Bühlmann credibility premium is the best affine estimate under mean squared error; an exact Bayesian posterior mean can also have this form, as in the natural conjugate credibility identity.

With the gamma scale inverse-gamma conjugacy model and prior shape $k>2$, the Bühlmann model has $v/a=(k-1)/\alpha$. Its Bühlmann credibility premium equals $\alpha (\lambda +{\sum }_j{X}_j)/(k+n\alpha -1)$, exactly the Bayes estimator under squared error loss of $\alpha \Theta$. Equality holds because the posterior mean is already affine in the observed sample mean.

A credibility factor is the weight assigned to a particular risk's observed experience when blending it with a population or manual premium. In the Bühlmann model, the blend is the best affine mean-square predictor, rather than automatically the full posterior mean.

Let $v=\mathbb{E}[{\sigma }^2(\theta )]$, the expected process variance, and assume the displayed second moments are finite. The law of total expectation gives $\mathbb{E}{X}_i=m$. The law of total variance and law of total covariance giveFor the second formula, the conditional covariance is zero by conditional independence, while the covariance of the conditional means is $\mathrm{Var}\mu (\theta )=a$, the variance of hypothetical means. Likewisebecause $\mathbb{E}[X_i-\mu (\theta )\mid \theta ]=0$.

To derive the Bühlmann credibility premium, minimize $\mathbb{E}[(\mu (\theta )-b_0-{\sum }_i{b}_i{X}_i{)}^2]$. Differentiating in $b_0$ sets the predictor's mean to $m$. The remaining normal equations areFor $v>0$, subtracting any two rows forces all $b_i$ to be equal, then gives $b_i=a/(v+na)$. Thus the premium and factor areThe same affine premium predicts $X_{n+1}$: conditional independence makes $\mathbb{E}[(X_{n+1}-\hat{\mu }{)}^2]=v+\mathbb{E}[(\mu (\theta )-\hat{\mu }{)}^2]$, so minimizers agree. If $v=0$, every year's amount equals $\mu (\theta )$ almost surely, and $Z=1$ gives the exact risk mean even though the coefficient vector need not be unique.

For fixed $a>0,n$, $\partial Z/\partial v=-na/(na+v{)}^2<0$. Thus more within-risk variation decreases credibility: individual experience becomes less informative, $Z\to 0$ as $v\to \infty$, and the premium tends to $m$. Increasing the number of observed years has the opposite effect.

For numerical estimation from a balanced panel, assume the $k$ risks are independent draws from the same risk population, with comparable exposures, and $n,k\ge 2$. SetConditionally, each within-risk sample variance has mean ${\sigma }^2({\theta }_s)$; hence $\mathbb{E}\hat{v}=v$. The risk sample means have variance $a+v/n$, so $\mathbb{E}B=a+v/n$. Therefore $\tilde{a}=B-\hat{v}/n$ is an unbiased method of moments estimate of $a$. For a usable nonnegative variance estimate take $\hat{a}=\max (0,\tilde{a})$, noting that this truncation sacrifices exact unbiasedness.

The balanced-panel method-of-moments credibility estimate for risk $s$ is consequentlyIf $\hat{a}=0$ and $\hat{v}>0$, use zero credibility. If both estimates are zero, the within and between sums of squares are zero and all the observed amounts coincide; any blend gives that same premium. With only one year per risk or only one risk, these within/between estimators are unavailable, so external structural estimates or additional assumptions are necessary.

In the Bühlmann model, one policy has a random risk parameter $\Theta$. Conditional on $\Theta$, its annual observations $X_i$ are independent and identically distributed random variables, with mean $m(\Theta )$ and variance $v(\Theta )$. Put $m_0=\mathbb{E}[m(\Theta )]$, let $v=\mathbb{E}[v(\Theta )]$ be the expected process variance, and let $a=\mathrm{Var}(m(\Theta ))$ be the variance of hypothetical means. Assume the required second moments are finite.

The Bühlmann credibility premium is the best affine estimate of $m(\Theta )$ under mean squared error. After choosing the intercept to make the estimate unbiased, write it as $m_0+{\sum }_{i=1}^n{b}_i(X_i-m_0)$. If $B={\sum }_i{b}_i$, the law of total variance and conditional independence give its errorFor fixed $B$, the Cauchy-Schwarz inequality minimizes the last sum at $b_i=B/n$. Differentiating $a(1-B{)}^2+v{B}^2/n$ gives $B=na/(na+v)$. Thus the credibility factor and credibility premium areThe ratio form assumes $a>0$; when $a=0$ and $v>0$, the conditional mean is constant and $Z=0$. Predicting $X_{n+1}$ rather than its conditional mean adds the constant $v$ to the error and leaves the optimal affine estimate unchanged.

Here the latent intensity has a Pareto distribution with shape three and lower bound one. For $r<3$, its moments are $\mathbb{E}[{\Theta }^r]=3{\int }_1^{\infty }{\theta }^{r-4}d\theta =3/(3-r)$. ConsequentlyThe equality for $v$ uses the Poisson distribution conditional variance $v(\Theta )=\Theta$. There are two years of experience, so $Z=1/2$ and the observed sample mean is $k/2$. The Poisson credibility with a shape-three Pareto intensity is thereforeIn particular it equals $2$ when $k=5$.

For the exact Bayesian credibility estimate, the total count conditional on $\Theta =\theta$ has Poisson distribution of mean $2\theta$. Multiplication of its likelihood by the prior distribution density gives a Bayesian posterior density proportional to ${\theta }^{k-4}{e}^{-2\theta }$ on $\theta >1$. At $k=5$, its normalizing integral and first-moment integral areHence the next-year predictive expected value isThe Bühlmann credibility premium optimizes over affine estimates, whereas the Bayesian posterior mean optimizes over all square-integrable estimates. This explains why the two credibility estimates need not agree for this nonconjugate prior.

The PDF gives an inverse-gamma distribution with shape $k$ and scale $\theta$, and a gamma distribution prior with shape $\alpha$ and rate $\lambda$. Write $\mu (\theta )=\theta /(k-1)$. The supplied conditional moments, together with $\mathbb{E}\theta =\alpha /\lambda$ and $\mathbb{E}{\theta }^2=\alpha (\alpha +1)/{\lambda }^2$, give the three structural quantities of the Bühlmann model:Here $v$ is the variance of hypothetical means and $s^2$ the expected process variance. The law of total covariance and conditional independence of different years implyAll these quantities are finite because $k>2$.

For any coefficients $c_i$, minimizing the mean squared error over the intercept first gives $c_0=m(1-{\sum }_i{c}_i)$. Let $Y=\mu (\theta )-m$ and $W_i=X_i-m$. Expanding the remaining mean squared error givesIts normal equations are $s^2{c}_i+v{\sum }_j{c}_j=v$ for every $i$. Since $s^2>0$, subtracting any two equations shows that all $c_i$ agree. The quadratic part is strictly positive for every nonzero coefficient change, so the solution is the unique global minimizer. Solving givesA credibility estimate combines individual experience with the collective expected value, using an experience weight $Z\in [0,1]$. In this equal-exposure setting, the Bühlmann credibility premium isThus the computed coefficients give exactly a credibility estimate, with a fixed credibility factor determined by known parameters and the number of observed years. This is Bühlmann credibility for inverse-gamma observations.

For the final Bayesian calculation, the likelihood function in $\theta$ is proportional to ${\theta }^{nk}\exp (-\theta {\sum }_i{x}_i^{-1})$. Multiplication by the gamma distribution prior gives the Gamma prior for an inverse-gamma scale updateUnder quadratic loss, the conditional risk of an estimate $d$ is $\mathrm{Var}(\mu (\theta )\mid x)+(d-\mathbb{E}[\mu (\theta )\mid x]{)}^2$. Its minimizer is therefore the posterior mean:This is not a fixed-weight arithmetic-mean credibility estimate. For $n=2$, the observations $(1,3)$ and $(2,2)$ have the same sample mean, but reciprocal sums $4/3$ and $1$, hence different posterior means. For any $n>2$, append the same positive observations to these two samples. For $n=1$, the estimate is $Cx/(1+\lambda x)$, with $C=(\alpha +k)/(k-1)$; its second derivative is $-2C\lambda /(1+\lambda x{)}^3\ne 0$, so it is not affine either. Symmetry in the observations would force equal slopes in any affine representation, so these examples also exclude a general fixed-coefficient linear representation.

One can express the posterior mean as a broader data-dependent blend. With $A={\sum }_i{x}_i^{-1}$ and harmonic mean $H=n/A$, it isThis weights a harmonic-mean-based experience estimate, and the weight itself depends on the observations. It is not the credibility estimate based on $\overline{X}$ and the fixed Bühlmann credibility factor derived above.

In the Bühlmann model, a latent risk parameter $\Theta$ is drawn from a population distribution. Conditional on $\Theta$, the yearly observations are independent and identically distributed random variables, with conditional expectation $m(\Theta )$ and conditional variance $v(\Theta )$. Define the structural parametersHere $v$ is the expected process variance, while $a$ is the variance of hypothetical means. The Bühlmann credibility premium is the best affine estimate of $m(\Theta )$ from the observed claims, under mean squared error. Predicting the next claim gives the same affine estimate: the extra conditional observation noise contributes the constant $v$ to the prediction error.

The law of total variance and conditional independence giveAn affine estimate can be written as $\hat{m}=m+{\sum }_{j=1}^n{b}_j(X_j-m)$: for any chosen $b_j$, optimizing the constant makes its expected value equal to $m$. The normal equations for the linear least-squares projection areFor $v>0$ they force all $b_i$ to agree, with $b_i=a/(v+na)$. Thus the credibility factor and premium areThe credibility factor increases with the observation count and between-risk variance, and decreases with within-risk variance. If $a=0$, the risk mean is known and $Z=0$; if $v=0$ and $a>0$, one observation reveals it and $Z=1$. If both vanish, the premium is the fixed value $m$ and the factor is immaterial.

In the specified model, the conditional law is a gamma distribution with shape $\alpha$ and scale $\theta$. ThereforeThe prior is an inverse-gamma distribution with shape $k$ and scale $\lambda$. To obtain its moments directly, substitute $y=\lambda /\theta$ in the defining integral, obtainingThe Gamma function recurrence yieldsThe assumption $k>2$ makes both structural variances finite. HenceIt follows that the model-specific credibility estimate is

For the Bayes estimator under squared error loss, the quantity to estimate is $\alpha \Theta$, so the optimum is its posterior mean. The likelihood function, viewed as a function of $\theta$, is proportional toMultiplication by the prior shows gamma scale inverse-gamma conjugacy:Although the printed hint only mentions integer shapes, the same substitution and gamma integral normalize this posterior for every positive real shape, so no integrality of $\alpha$ is needed. Its posterior mean gives the Bayesian estimate and comparisonThis exact Bühlmann credibility for gamma claims holds for every observed sample, not merely on average. Here the posterior mean is affine in the sample mean, so the best affine Bühlmann credibility premium is also the unrestricted Bayes estimator under squared error loss.