Our decision to use residualized scores, instead of explicitly modeling dummies for age in the twin models, was primarily related to feasibility (and comparability) of testing the various twin models. While there are limitations to using residualized scores (e.g. the impact of age cannot be simultaneously modeled), there are challenges to including age as a covariate in the some of the more complex models, due to power.
