Past exam of the mathematics course of the University of Cambridge / 2014 / iii / Paper 31 / 4

Write the year-$j$ count as a sum of its $m_j$ individual counts, with $m_j>0$. The specified conditional independence within that year givesDividing by $m_j$ and scaling the conditional variance by $m_j^{-2}$ therefore givesThis is the Bühlmann–Straub model: larger exposures reduce the process noise in a year's average.

Define the population parameters and total observed exposure byHere $v$ is the expected process variance and $a$ the variance of hypothetical means. Assume finite second moments. The law of total variance gives $\mathrm{Var}(X_j)=a+v/m_j$, and conditional independence between years gives $\mathrm{Cov}(X_j,X_k)=a$ for $j\ne k$: their only shared variation is the latent conditional mean.

To derive the Bühlmann–Straub credibility estimate, minimize mean squared error among affine estimates of $\mu (\Theta )$. For coefficients $b_j$ and $B={\sum }_j{b}_j$, the optimal intercept is $m_0(1-B)$. Centering at $m_0$ and conditioning on $\Theta$ shows that the resulting error isThe cross term vanishes because $X_j-\mu (\Theta )$ has conditional mean zero; the conditional noise cross terms vanish by the between-year independence assumption.

For fixed $B$, the Cauchy-Schwarz inequality giveswith equality when $b_j=B{m}_j/W$. Minimize the remaining quadratic $a(1-B{)}^2+(v/W)B^2$. Its minimizing value is the Bühlmann–Straub credibility factorThe experience term is exposure-weighted:Since the future conditional mean count is $m_{n+1}\mu (\Theta )$, multiplying the optimal estimate by the known future exposure givesIt is the best affine linear least-squares projection, not a claim that the exact conditional expectation given all data is always affine. If $a=0$, the conditional mean is a known constant almost surely and one takes $Z=0$; if $v=0<a$, the observations reveal it without process noise and $Z=1$. The ordinary case has $a,v>0$.

Use the gamma distribution shape-rate convention: $\Theta \sim \mathrm{Gamma}(\alpha ,\beta )$ with $\alpha ,\beta >0$. The conditional Poisson distribution has mean and variance both equal to $\Theta$, so the Bühlmann–Straub model parameters areFor $W={\sum }_j{m}_j$ and $y={\sum }_j{y}_j$, the Bühlmann–Straub credibility factor is $Z=W/(W+\beta )$. Thusand the required expected-count estimate isThe total exposure determines how much information the observed count carries; the number of years alone is not the appropriate denominator.

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.