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.