= Solution
For a population-averaged rate trend, use the <marginal incidence trend with clustered observations> model
$$
\boxed{\log\mathbb E[Y_{it}]=\log P_{it}+\alpha+\beta t,\qquad t=0,\ldots,4.}
$$
Fit a <generalized estimating equation> with village as cluster, a log link and an appropriate working correlation; estimate uncertainty with a village-level <sandwich covariance matrix>. Then $e^\beta$ is the fitted annual marginal <rate ratio>, $100(e^\beta-1)$ is the annual percentage rate change, and $e^{4\beta}$ is the four-year endpoint ratio. Test $\beta=0$ using cluster-robust uncertainty, with a small-sample correction or village-resampling <bootstrap> given only 32 independent clusters. Also compare year as a factor with a linear trend, since a single slope can conceal nonlinear behaviour.
Report which population is averaged over. The directly pooled rate $\sum_iY_{it}/\sum_iP_{it}$ is exposure weighted; a mean of village rates gives each village equal weight. For standardized comparisons one can use the same village weights at every year, preventing changing population composition from masquerading as a rate trend.
A conditional random-slope coefficient cannot simply be substituted for the marginal slope: the <marginal mean of a Poisson random-slope model> contains time-dependent variance terms. \b[A direction, magnitude and confidence interval for the trend across all villages require the missing full dataset.] The increase in the two displayed villages does not determine the trend in the other thirty.
Back to article page
