Past exam of the mathematics course of the University of Cambridge / 2006 / iii / Paper 44 / 1

Use a longitudinal study with repeated measurements of the same individuals. Recruit people across several starting ages and measure the same ear, with a standardized anatomical definition and instrument, at baseline and subsequent visits over several years. Repeat measurements at each visit to estimate measurement error; where practical, conceal previous measurements from the observer. Record sex, body size, birth cohort and observer, and document reasons for loss to follow-up.

A simple random-intercept linear mixed model for ear length $Y_{ij}$ at elapsed follow-up time $t_{ij}$ isThe individual intercept ${\alpha }_i$ absorbs persistent differences between people. Alternatively, differences $Y_{ij}-Y_{i0}$ remove that intercept directly. Estimate the within-person regression coefficient $\beta$ and its confidence interval, accounting for the correlation of repeated measurements. Compare it with the cross-sectional slope rather than treating the two slopes as automatically identical.

This longitudinal study directly checks whether an individual's ears tend to lengthen as time passes. Several entry-age cohorts help assess whether the change varies with starting age. Calendar-period changes and selective dropout can still affect interpretation, so repeated observation alone does not eliminate every possible source of confounding or selection bias.

The association is strongly statistically significant under the reported regression assumptions, but age alone has only moderate predictive ability. The reported confidence interval for the age regression coefficient excludes zero by a wide margin. If it is the usual two-sided 95% interval from a normal linear model with an intercept and one predictor, there are $206-2=204$ degrees of freedom. Taking $t_{0.975,204}\simeq 1.97$ givesThus the two-sided p-value against zero slope is far below $0.001$. This conclusion concerns a statistical association in the sampled population; it does not establish an individual growth rate. Clustering by practitioner, nonconstant variance or nonindependent observations would require an appropriate standard error rather than blind reliance on this calculation.

The scatter plot has substantial vertical spread at each age. A narrow slope confidence interval describes uncertainty about the average trend, whereas a prediction interval for a new person's ear length must also include the residual variance. The reported interval even permits an approximate numerical assessment: write $S_{xx}={\sum }_i(x_i-\bar{x}{)}^2$ and let $s_e^2=\mathrm{RSS}/(n-2)$. In simple linear regression, $t^2={\hat{\beta }}^2{S}_{xx}/s_e^2$, while the fitted and residual sums of squares are ${\hat{\beta }}^2{S}_{xx}$ and $(n-2)s_e^2$. Hence the coefficient of determination isUnder these assumptions, age explains only about 27% of the sample variation; roughly 73% remains unexplained. The calculation is approximate because the printed slope and interval are rounded. For a new individual aged $x_0$, the usual prediction interval is based onThe leading one represents individual residual variation. The provided summaries do not determine $s_e$ and $S_{xx}$ separately, so they do not supply a numerical prediction interval.

Record sex, height or another measure of body size, ancestry, family characteristics and practitioner or measurement method. Use these covariates in a multiple linear regression, allowing scientifically justified nonlinear age effects or interactions, and assess prediction error on held-out observations. A guessed overall sex proportion cannot replace individual sex measurements or identify sex-adjusted effects.

There is a source inconsistency: the PDF's numerical ear-length summary is ten times the scale of its plotted measurements and fitted line. Indeed, the fitted line at the reported mean age gives about $67.7$ mm, whereas the printed mean is about ten times larger. Missing decimal points in the summary are a plausible explanation. The significance and approximate $R^2$ above use the mutually consistent slope, interval and graph; they do not treat the inconsistent summary as valid data.

A cross-sectional age association need not equal within-person growth. The following alternatives each suggest a concrete extension of the study.

These explanations are hypotheses to be checked, not established facts about the participants. Neither a small p-value nor an approximately linear scatter plot excludes them.