13.3 A heterogeneous example
Example 13.3.1. For the British doctors’ data (Example 10.6.1), where the stratum rate ratios fall from \(5.74\) to \(0.90\) with age:
| Method | Rate ratio | \(95\%\) interval |
| Crude | 1.72 | (1.39, 2.12) |
| Direct, combined standard | 1.42 | (1.15, 1.75) |
| Direct, smokers as standard | 1.42 | (1.15, 1.75) |
| Direct, non-smokers as standard | 1.43 | (1.16, 1.76) |
| SMR, non-smokers’ rates as standard | 1.42 | (1.31, 1.53) |
| Mantel–Haenszel | 1.42 | (1.15, 1.76) |
| Inverse-variance | 1.36 | (1.10, 1.68) |
| Poisson maximum likelihood | 1.43 | (1.15, 1.76) |
Here too most methods agree, because their weights happen to be similar: all put about \(85\%\) of the weight on the three oldest strata, where most of the deaths are. The inverse-variance estimate is lower because it averages on the log scale. A weighted geometric mean is never larger than the corresponding arithmetic mean, and the gap is large when the values are spread out: the ratio \(5.74\) in the youngest stratum pulls an arithmetic average up far more than a geometric one. The agreement of the other methods should not be mistaken for evidence that \(1.42\) is a meaningful summary. The deviance of \(12.1\) on \(4\) degrees of freedom shows that the common ratio model fails, and the proper report is the age-specific ratios or the trend model of Example 12.4.3.
Questions on this section
Stuck on something here? Ask below and it stays attached to this topic.