= Solution
The command fits a Gaussian <generalized additive mixed model>
$$
\boxed{Y_{ij}=\alpha+s(t_{ij})+b_jt_{ij}+\varepsilon_{ij},\qquad b_j\overset{\rm iid}{\sim}N(0,\tau^2),\quad \varepsilon_{ij}\overset{\rm iid}{\sim}N(0,\sigma^2).}
$$
The population growth curve $s$ is a centered <penalized regression spline> with a cubic regression basis. The incubator term remains a <random slope> without a <random intercept>, because `~0+hours` removes the intercept. A population linear component is allowed within the smooth's unpenalized part; curvature is penalized rather than being required.
A `gamm` fit returns a list with `gam` and `lme` components. Plot the population smooth and pointwise uncertainty bands with
``
plot(fly.model2$gam, se=TRUE, residuals=TRUE)
``
Inspect whether the fitted curve departs materially from an <affine function>, particularly in regions supported by observations. The <effective degrees of freedom> near one suggest that the penalization has selected an approximately linear centered smooth. Bands and partial residuals help distinguish convincing curvature from noisy deviations, although pointwise bands are not a simultaneous test of linearity. A formal comparison would need to account for both smoothing selection and the incubator dependence. \b[The smooth plot tests the visual plausibility of a common straight growth curve while retaining incubator-specific slope variation.]
Back to article page
