Solution
= Solution
Use the shape-rate convention $R_k\sim\operatorname{Gamma}(\alpha_0,\beta_0)$. The factors involving $R_k$ in the <posterior density> are
$$
R_k^{I_k}e^{-s_kR_k}R_k^{\alpha_0-1}e^{-\beta_0R_k}
=R_k^{I_k+\alpha_0-1}e^{-(s_k+\beta_0)R_k}.
$$
By <Poisson-gamma conjugacy>,
$$
R_k\mid I_k,s_k\sim\operatorname{Gamma}(I_k+\alpha_0,s_k+\beta_0),
$$
so its <posterior mean> is $(I_k+\alpha_0)/(s_k+\beta_0)$.