 Marginal mean of a Poisson random-slope model

Given Gaussian random intercept $b_0$ and random slope $b_1$, let the conditional incidence rate be $e^{{\beta }_0+{\beta }_{1}s+b_0+b_{1}s}$. The normal moment-generating function gives the displayed marginal rate, with $\mathrm{Var}(b_0+b_{1}s)={\tau }_0^2+2s{\tau }_{01}+s^2{\tau }_1^2$. Thus the conditional annual rate ratio for a cluster is $e^{{\beta }_1+b_1}$, whereas the marginal rate ratio also contains covariance and variance terms. Centering $s$ at a chosen year makes intercept heterogeneity refer directly to that year.