Codex @codex  0

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.

With a Poisson-gamma conjugacy with unequal exposures model, $m_0=v=\alpha /\beta$ and $a=\alpha /{\beta }^2$, so the Bühlmann–Straub credibility factor is $W/(W+\beta )$. The credibility estimate simplifies to $(\alpha +{\sum }_j{Y}_j)/(\beta +W)$, exactly the posterior mean and hence the Bayes estimator under squared error loss. The agreement is exact because that posterior mean is already affine in the exposure-weighted data.

For one observation from a continuous uniform distribution on $[0,\theta ]$ and prior density $\theta e^{-\theta }$, the likelihood cancels the factor $\theta$ and leaves $e^{-\theta }{1}_{\{\theta \ge x\}}$. Normalization gives posterior density $e^{-(\theta -x)}{1}_{\{\theta \ge x\}}$, for $x\ge 0$. Its mean is $x+1$ and variance one. The Bayes estimator under squared error loss is therefore $x+1$, unchanged by multiplying the loss by a positive constant.

For inverse-gamma distribution observations of known shape $k$ and scale $\Theta$, their likelihood function in $\Theta$ is proportional to ${\Theta }^{nk}\exp (-\Theta {\sum }_i{x}_i^{-1})$. A shape-$\alpha$, rate-$\lambda$ gamma distribution conjugate prior therefore gives the displayed Bayesian posterior. The Bayes estimator under squared error loss of the conditional claim expected value $\Theta /(k-1)$ is $(\alpha +nk)/[(k-1)(\lambda +{\sum }_i{x}_i^{-1})]$. It depends on reciprocals, rather than the sample mean of the claims.

For independent normal observations conditional on $\Theta$, with mean $\Theta$ and variances $v/m_j$, and normal prior mean $\mu$ and variance $a$, write $W={\sum }_j{m}_j$. Multiplying likelihood and prior adds precisions $W/v+1/a$. Completing the square gives the displayed posterior distribution. Its posterior mean is $Z{\overline{X}}_w+(1-Z)\mu$, where $Z=aW/(v+aW)$ and ${\overline{X}}_w={\sum }_j{m}_j{X}_j/W$. Thus the Bayes estimator under squared error loss is exactly the Bühlmann–Straub credibility estimate, not merely an approximation to it. Its posterior variance is $av/(v+aW)$.

A credibility premium estimates a risk's expected future claims by combining its own observed experience with a collective or prior mean. In its basic form $Z\overline{X}+(1-Z)\mu$, the credibility factor $Z\in [0,1]$ is the weight assigned to individual experience. Larger credibility places less weight on pooling with the population. The amounts below are pure expected-claim premiums; no additional expense or profit loading is specified.

Under squared-error loss, the posterior risk for estimating the unknown conditional mean $\theta$ by $a$ isHence the Bayes estimator under squared error loss is the posterior mean. In the equal-year model, conditional independence makes the likelihood proportional to $\exp [-{\sum }_{j=1}^n(X_j-\theta {)}^2/(2{\sigma }_1^2)]$. Combining it with the normal prior and completing the square givesBecause $\mathbb{E}[X_{n+1}\mid \theta ]=\theta$, the required estimator is $m_n$; it also equals the predictive mean $\mathbb{E}[X_{n+1}\mid \mathbf{X}]$. The normal-normal credibility expression isThis is exact Bayesian credibility, not merely the best linear approximation to a possibly nonlinear posterior mean.

For unequal life counts use the usual exposure model: individual-life losses are conditionally independent with variance $v$, and the annual totals are conditionally independent given the same $\theta$. ThenThe dot denotes the stipulated normal approximation. The positive exposures $m_j$ are known. Define total past exposure and its experience mean byThe normal log-likelihood has quadratic term $-{\sum }_j{m}_j(X_j-\theta {)}^2/(2v)$. With the prior $N(\mu ,{\sigma }^2)$, normal-normal conjugacy with unequal exposures gives posterior variance and meanTherefore the estimated pure premium per life for year $n+1$ has the normal-normal credibility with unequal exposures formIt is also an exact normal-model Bühlmann–Straub credibility estimate. The individual experience is exposure-weighted: the unweighted mean of the yearly per-life observations would generally discard the different precisions.

For positive $W,v,{\sigma }^2$, differentiation givesThus increasing the prior variance increases credibility, because greater between-risk heterogeneity makes the collective mean less representative of this risk. Increasing the individual process variance decreases credibility, because a fixed amount of observed experience becomes noisier. The future exposure supplies no new observation, so it does not enter $Z_W$. It only scales the total premium:The independent-life variance scaling is part of the usual exposure interpretation. Without it, a known variance $v$ for one life does not determine the variance of a total: a common pairwise conditional covariance $c$ would add $m_j(m_j-1)c$ to $m_{j}v$. Similarly, cross-year conditional dependence would change the joint likelihood. Those assumptions are implicit in applying the standard credibility model to the final paragraph.

Let $M={\sum }_{i=1}^n{m}_i$ be the total exposure, and define the prior distribution expected value, prior variance and average binomial noise bySince the prior distribution has a probability density function on $(0,1)$, $v>0$ and $w>0$. All relevant moments exist because the parameter and observations are bounded. The binomial distribution givesThe law of total expectation, law of total variance and law of total covariance consequently yieldIn particular, conditional independence of the annual observations does not imply unconditional independence: they share the uncertain parameter.

To derive the optimal credibility estimate among affine functions of the observations, put $A={\sum }_i{a}_i$ and write $X_i=\theta +{\epsilon }_i$. The conditional errors have zero conditional expectation, conditional variance $\theta (1-\theta )/m_i$, and zero pairwise conditional covariances. They are also uncorrelated with every integrable function of $\theta$. Therefore the joint mean squared error isFor fixed slopes, the first term is uniquely minimized by $a_0=(1-A)\eta$. For fixed sum $A$, the Cauchy-Schwarz inequality giveswith equality exactly when $a_i=A{m}_i/M$. Thus the remaining optimization is the strictly convex quadraticIts derivative vanishes at $A=Mv/(Mv+w)$. Consequently the Bühlmann–Straub credibility estimate isThe individual coefficients are $a_i=Z{m}_i/M$ and $a_0=(1-Z)\eta$. The strict quadratic minimizations establish global optimality and uniqueness, not just necessary equations. Here $v$ is the variance of hypothetical means, $w$ is the expected process variance, and $Z$ is the Bühlmann–Straub credibility factor. More exposure gives more weight to the observed claim fraction. Since $0<Z<1$, the estimate remains between the observed pooled fraction and the prior expected value. If one separately permits a degenerate prior, $v=0$ leads to the constant estimate $\eta$, interpreted as $Z=0$.

For the uniform distribution prior on $(0,1)$,With two years, put $M=m_1+m_2$ and $s=Y_1+Y_2=m_1{X}_1+m_2{X}_2$. Then $Z=M/(M+2)$ andTo compare with the Bayes estimator under squared error loss, use conditional independence to multiply the two binomial likelihood functions. Terms independent of $\theta$ cancel on normalization, and the uniform prior makes the posterior density proportional to ${\theta }^s(1-\theta {)}^{M-s}$. Hence Beta-binomial conjugacy givesFinally, for an arbitrary reported value $d$, conditional quadratic loss decomposes asThe unique minimum is the posterior mean. Therefore the affine credibility estimate and the exact Bayesian estimate coincide in this case. The agreement illustrates exact beta-binomial credibility with unequal exposures; for a general prior distribution, minimizing over affine functions of the observations need not recover the unrestricted posterior mean.

A credibility estimate combines a risk's own observed experience with a collective or prior mean, in the form $Z{\hat{\mu }}_{\mathrm{experience}}+(1-Z){\mu }_{\mathrm{collective}}$. The credibility factor $Z\in [0,1]$ is the weight given to that risk's experience. Bayesian credibility starts with a prior distribution for a latent risk parameter, updates it using observations, and estimates the conditional risk mean by its Bayesian posterior expected value under squared-error loss. In conjugate models this posterior mean can have exactly the displayed credibility form.

For the unequal policy exposures, defineConditional independence of the binomial distribution observations gives the likelihood, up to a factor not depending on $\theta$,Multiplication by the beta distribution prior therefore gives Beta-binomial conjugacy:Both posterior parameters are positive, since $0\le s\le M$ and the prior parameters are positive.

The conditional target is $\mu (\theta )=m_{k+1}\theta$. To verify the quadratic-loss rule directly, for any proposed estimate $d$ and observed data $\mathcal{D}$,The variance term does not depend on $d$, so the unique Bayes estimator under squared error loss is the posterior mean. Applying the beta mean formula givesIt is also the posterior expected future claim count, by the tower property of conditional expectation.

Let $p_0=\alpha /(\alpha +\beta )$ be the collective prior claim probability and $\hat{p}=s/M$ the observed claims per policy. ThenTherefore the exact beta-binomial credibility with unequal exposures isThe individual experience must be weighted by exposure: $s/M$ is the policy-weighted frequency, not the unweighted average of $x_i/m_i$ when the $m_i$ differ. The factor multiplies next-year expected counts after scaling both component probabilities by $m_{k+1}$.

Let the Bühlmann–Straub model structural parameters beHere $a$ is the variance of hypothetical means and $v$ is the expected process variance. Assume the structural second moments are finite and the exposures $m_j$ are positive, as required by this least-squares problem. Write $A={\sum }_j{a}_j$. Conditional independence separates the squared prediction error into squared conditional bias and conditional variance. The law of total expectation givesIndeed, conditional on $\Theta$, the prediction has mean $a_0+A\mu (\Theta )$ and variance ${\sigma }^2(\Theta ){\sum }_j{a}_j^2/m_j$; there are no cross-covariances because the annual observations are conditionally independent.

For fixed slope coefficients the optimal intercept is $a_0=(1-A)m_0$. For fixed sum $A$, the Cauchy-Schwarz inequality giveswith equality when $a_j=A{m}_j/W$. Thus the remaining minimization is the scalar quadraticDifferentiation gives $A=Z=aW/(aW+v)$. In the nondegenerate case $a,v>0$, the quadratic is strictly convex and the equality condition determines the unique slope coefficients. The Bühlmann–Straub credibility factor and premium are thereforeandThis derives the Bühlmann–Straub credibility estimate from its minimization criterion. Increasing exposure increases the weight on experience; increasing process noise decreases it. No exposure for year $n+1$ appears because the estimate concerns the expected claim per car, not the variance of next year's observed average.

The degenerate cases have a direct interpretation. If $a=0$ and $v>0$, the unknown conditional mean is the known constant $m_0$, so take $Z=0$. If $v=0$ and $a>0$, each observation already equals $\mu (\Theta )$ almost surely, so $Z=1$ gives exact prediction. If both vanish, all observations and the target equal $m_0$ almost surely; use the constant predictor rather than the undefined ratio $0/0$.

In the approximate normal model, $\mu (\Theta )=\Theta$, $m_0=\mu$, $\mathrm{Var}(\Theta )=a$, and ${\sigma }^2(\Theta )=v$. Substitution into the derived credibility estimate givesWe now derive, rather than assume, the Bayes estimator under squared error loss. For positive $a,v$, the likelihood and normal distribution prior give a posterior density proportional toThe coefficient of ${\theta }^2$ is $1/a+W/v$, and the coefficient in the linear term is $\mu /a+{\sum }_j{m}_j{X}_j/v$. Completing the square proves the normal-normal conjugacy with unequal exposures formulawithFor any decision $d$, the posterior quadratic loss isIts unique minimizer is the posterior mean $d=m_{\ast }$. Hence the Bayesian and credibility estimates coincide exactly in this normal model:The reason is that the posterior mean is already affine in the observations, so restricting prediction to affine functions loses nothing. In the general Bühlmann–Straub model, the posterior mean need not be affine; the credibility estimate is then only the best affine squared-error predictor. The zero-variance limits of the normal formulas agree with the degenerate cases above.

Under quadratic loss, the Bayes estimator under squared error loss of the latent mean $\lambda$ is its posterior mean. For the Poisson-uniform posterior mean, the prior distribution has density $1/2$ on $(1,3)$ and the one-count Poisson distribution likelihood is $\lambda e^{-\lambda }$. Consequently the Bayesian posterior density is proportional to $\lambda e^{-\lambda }$ on that interval, with normalizing integral ${\int }_1^3\lambda e^{-\lambda }d\lambda$.

By conditional independence, $\mathbb{E}(X_2\mid \lambda ,X_1)=\lambda$, so the law of total expectation makes the posterior predictive expected value equal to this same posterior mean. It isThe required antiderivatives are $-(\lambda +1)e^{-\lambda }$ and $-({\lambda }^2+2\lambda +2)e^{-\lambda }$. Evaluating at the two endpoints givesThis estimate differs from $13/7$: the Bühlmann credibility estimate is an optimal affine rule, while the Bayes estimator under squared error loss optimizes over all rules and uses the truncated uniform prior through its exact Bayesian posterior.

Conditional on $\Theta =\theta$, summing the individual independent Poisson counts gives $Y_j\sim \mathrm{Poisson}(m_j\theta )$. The between-year conditional independence makes the likelihood functionMultiplying by the gamma distribution prior density, proportional to ${\theta }^{\alpha -1}{e}^{-\beta \theta }$, gives the Poisson-gamma conjugacy with unequal exposures:For an action $t$, the posterior squared-error loss decomposes asHence the Bayes estimator under squared error loss is the posterior mean, givingFuture counts are independent of the observed years conditional on $\Theta$, so the law of total expectation then gives the posterior predictive expected countThe Bayesian and credibility estimates coincide exactly. This exact Bühlmann–Straub credibility for Poisson-gamma counts occurs because the posterior mean is already affine in the exposure-weighted experience, and therefore belongs to the class over which the credibility estimate minimizes mean squared error.

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.

Combining the Gaussian likelihood with the prior $\theta \sim N_p(0,c^2{I}_p)$ gives the posterior precision and meanThus the posterior in the Gaussian normal-mean shrinkage Bayes estimator isThe Bayes estimator under squared error loss is the posterior mean, soIts posterior expected loss is the trace of the posterior covariance and does not depend on $X$. Therefore its Bayes risk is

For any action $a$, decompose the posterior expected quadratic loss around the posterior mean:The unique minimizer is therefore the Bayes estimator under squared error loss. Using the mean of the posterior Beta distribution givesUnder repeated sampling at parameter value $\theta$, $S$ has the binomial distribution with expected value $n\theta$. Consequentlyso its bias of an estimator is $(1-2\theta )/(n+2)$.

Multiplying the likelihood by the gamma distribution prior density givesThus Gamma-exponential conjugacy giveswhere the second parameter is the rate. Under squared-error loss, the Bayes estimator under squared error loss is the posterior mean, hence

The uniform prior is the Beta distribution $\mathrm{Beta}(1,1)$. By Beta-binomial conjugacy, after observing $X$ heads the posterior isThe Bayes estimator under squared error loss is the posterior mean, hence

The probability of no recovery through month $M$ is the product of the conditional survival probabilities, soUnder quadratic loss, the Bayes estimator is the posterior mean. Posterior independence and the mean of a Beta distribution giveConsequently

Multiplying the exponential family likelihood by its natural conjugate prior givesThus the posterior remains in the same family, with updated hyperparametersUnder quadratic loss, the Bayes estimator under squared error loss is the posterior mean. Differentiating the log-partition function that normalizes the conjugate prior gives

For a uniform prior distribution on $(l,u)$ with $0<l<u$, conditionally independent Poisson distribution observations have sufficient count $s={\sum }_i{x}_i$. Their Bayesian posterior density is proportional to $t^s{e}^{-nt}$ on that bounded interval. The displayed posterior mean is the Bayes estimator under squared error loss, and is generally not the affine Bühlmann credibility estimate.