3.6 Problems
Problem 3.6.1. Of \(250\) children examined in a school, \(3\) have a skin infection. Compute the Wald, Wilson and exact \(95\%\) intervals for the prevalence and comment. What is the exact upper limit if no child had been infected?
Show solution
Solution. \(\hat p = 0.012\). The Wald interval is \[ 0.012 \pm 1.96\sqrt {0.012 \times 0.988/250} = 0.012 \pm 0.0135 , \] that is \((-0.0015,\ 0.0255)\), which includes impossible negative values. The Wilson interval, from (3.1), is \((0.0041,\ 0.0347)\), and the Clopper–Pearson interval is \((0.0025,\ 0.0347)\). The two good intervals are asymmetric about \(\hat p\), as they should be near a boundary.
With \(x = 0\) the exact upper limit solves \((1 - p)^{250} = 0.025\), giving \(p_U = 1 - 0.025^{1/250} = 0.0146\), close to \(-\ln (0.025)/250 = 3.69/250\). The same argument for a one-sided \(95\%\) limit gives \(-\ln (0.05)/n \approx 3/n\), the “rule of three”.
Problem 3.6.2. In a cohort followed for a total of \(12{,}400\) person-years, no cases of a rare cancer occur. Give a \(95\%\) upper confidence limit for the incidence rate. How many person-years would be needed for the upper limit to fall below \(1\) per \(10{,}000\) person-years if again no cases were seen?
Show solution
Solution. With \(D = 0\), (3.2) gives \(\lambda _U = \chi ^2_{2,\,0.975}/(2T) = 7.378/(2 \times 12400) = 2.97\) per \(10{,}000\) person-years; equivalently \(P(D = 0) = e^{-\mu } = 0.025\) gives \(\mu = 3.689\) and \(\lambda _U = 3.689/T\). For \(3.689/T < 10^{-4}\) needs \(T > 36{,}890\) person-years.
Problem 3.6.3. People are subject to a constant hazard \(\lambda \) of a disease and a constant hazard \(\mu \) of death from other causes, independently. Show that the probability of developing the disease by time \(t\) is \[ \frac {\lambda }{\lambda + \mu }\bigl (1 - e^{-(\lambda + \mu )t}\bigr ) . \] Evaluate it for \(\lambda = 0.01\), \(\mu = 0.04\) per year and \(t = 20\), and compare with \(1 - e^{-\lambda t}\).
Show solution
Solution. A person is alive and disease-free at time \(u\) with probability \(e^{-(\lambda + \mu )u}\), since the two hazards add. The probability of developing the disease in \((u, u + du]\) is that times \(\lambda \,du\). Integrating, \[ \int _0^t \lambda e^{-(\lambda + \mu )u}\,du = \frac {\lambda }{\lambda + \mu }\bigl (1 - e^{-(\lambda + \mu )t}\bigr ) . \] With the given values: \(0.2(1 - e^{-1}) = 0.2 \times 0.632 = 0.126\), against \(1 - e^{-0.2} = 0.181\) without competing deaths. The second figure is the risk in a hypothetical population where nobody dies of anything else; the first is the proportion who actually develop the disease. Both are legitimate, but they answer different questions, and a risk reported in an elderly population must say which is meant.
Problem 3.6.4. Using Example 3.3.2, suppose instead that a researcher computed the risk as \(2/8 = 0.25\) over five years. Explain what is wrong. If the rate is constant, what five-year risk does the rate \(0.08\) per year imply?
Show solution
Solution. The denominator \(8\) treats all eight workers as observed for five years, but they contributed only \(25\) of the possible \(40\) person-years. Workers 3, 5 and 8 entered late, worker 4 moved away and worker 7 died, so their disease status at five years is unknown; treating them as disease-free underestimates the risk. Under a constant hazard the five-year risk is \(1 - e^{-0.08 \times 5} = 1 - e^{-0.4} = 0.33\).
Problem 3.6.5. A chronic disease has incidence \(2\) per \(1000\) person-years in non-smokers and \(4\) per \(1000\) in smokers, so the incidence rate ratio is \(2\). Among smokers the disease lasts (until death) \(4\) years on average, among non-smokers \(8\) years. In a steady state, find the prevalence in each group and the prevalence ratio. What would a cross-sectional study conclude?
Show solution
Solution. Prevalence odds are \(0.004 \times 4 = 0.016\) in smokers and \(0.002 \times 8 = 0.016\) in non-smokers, so both prevalences are \(0.016/1.016 = 1.57\%\) and the prevalence ratio is \(1\). Smoking doubles the incidence but halves the duration, and a cross-sectional study, which sees prevalence, would find no association. By Theorem 3.4.1 the prevalence odds ratio equals the incidence rate ratio times the ratio of mean durations, here \(2 \times \tfrac 12\).
Problem 3.6.6. In Example 3.5.2, suppose town B has lower death rates than town A in both age groups: \(1.5\) per \(1000\) under \(50\) and \(18\) per \(1000\) at \(50\) and over. Find town B’s crude rate and compare it with town A’s. Explain the result using Proposition 3.5.1.
Show solution
Solution. Deaths in B are \(40 \times 1.5 + 60 \times 18 = 60 + 1080 = 1140\), a crude rate of \(11.4\) per \(1000\), still twice town A’s \(5.6\). The crude rate weights the specific rates by the town’s own age distribution: \(0.4 \times 1.5 + 0.6 \times 18\) for B against \(0.8 \times 2 + 0.2 \times 20\) for A. Town B puts \(60\%\) of its weight on the high old-age rate, town A only \(20\%\), and that difference outweighs B’s lower rates in both groups. The crude comparison reverses the direction of the age-specific comparisons.
Questions on this section
Stuck on something here? Ask below and it stays attached to this topic.