Past exam of the mathematics course of the University of Cambridge / 2013 / iii / Paper 30 / 4 / c / Solution

Allow each subject to have a separate baseline and slope through a correlated random-intercept and random-slope model:The subject random effects are independent of all measurement errors, and subjects are independent. The within-person covariance between times $s,t$ is ${\tau }_0^2+(s+t)\rho {\tau }_0{\tau }_1+st{\tau }_1^2$, with an additional ${\sigma }^2$ when the same reading is used twice.

Prefer the model with a random slope. The maximum likelihood estimation refits improve twice the log-likelihood by $227.9604$ while adding two covariance statistical parameters. Both the Akaike information criterion ($5165.779$ to $4941.819$) and Bayesian information criterion ($5183.367$ to $4968.200$) strongly favour it. The printed nominal likelihood-ratio test is overwhelming. The usual reference chi-squared distribution with two statistical degrees of freedom is not an exact regular calibration, since zero slope variance is a boundary and the intercept–slope correlation is then unidentified; a design-specific parametric bootstrap could calibrate it. This qualification does not undermine the substantial descriptive improvement shown by both information criteria.

The preferred fit estimates population mean baseline strength ${\hat{\beta }}_0=61.09286$ kg and weekly gain ${\hat{\beta }}_1=2.83143$ kg/week. The between-person baseline standard deviation is ${\hat{\tau }}_0=15.785582$ kg, and the between-person slope standard deviation is ${\hat{\tau }}_1=2.738666$ kg/week. Their estimated correlation $\hat{\rho }=0.117$ is weakly positive, giving random-effect covariance about $5.058$ kg${}^2$/week; its uncertainty is not supplied. The measurement-error standard deviation is $\hat{\sigma }=9.923280$ kg. The printed $-0.038$ instead describes the correlation between the estimated fixed effects, not between the subject random effects. These statistical parameter estimates come from restricted maximum likelihood; the model comparison uses the separate maximum likelihood estimation fits.

The population mean fitted trajectory is $61.09286+2.83143t$ kg. HenceUsing the reported slope standard error, an approximate 95% confidence interval for the population mean gain is $28.3143\pm 1.96(2.984464)=(22.46,34.16)$ kg; a Student t confidence interval with 499 statistical degrees of freedom is almost identical.

Successful individuals can improve far more than the average. Their latent ten-week gains have fitted normal distributionA clear estimate for an upper-performing group uses a between-person slope quantile: the 95th percentile is $28.3143+1.645(27.38666)\simeq 73.4$ kg, and the 97.5th percentile is about $82.0$ kg. The upper 5% of fitted underlying gains begin around 73 kg. These are person-to-person performance quantiles, not confidence intervals for the population mean. Identifying the best observed trainee would require that person's data or fitted subject random effects; the aggregate output cannot identify a literal maximum. If performance means the observed difference of endpoint readings, add the measurement-error variance $2{\sigma }^2$ to the latent gain variance.