Attachment | Size |
---|---|

AgeMod.txt | 9.78 KB |

Hi,

I adapted script from the Boulder workshop to test a univariate twin model with age moderating the variance components of observational data (one single global score per twin member). When I run the model, I get a non-zero status code 6, but I can get the model to run if I add in a main effect of age on the mean. However, this model gives me negative estimates, and I’m not sure if that’s reflective of something wrong with my code. Any help is appreciated!

Thanks!

Amanda

I don't think you should try to interpret the results of the first model that lacks the main effect of age on the mean (and which returned Status Red). One way to think about this kind of biometric moderation model is that the phenotype is being regressed onto latent A, latent C, latent E, an interaction between age and latent A, an interaction between age and latent C, and an interaction between age and latent E. Not including the main effect of age is wrong in the same way that it's wrong in ordinary multiple regression to include an interaction term without including the main effects of the two predictors.

While it's true you obtained some negative point estimates from the second model, that's not necessarily bad. Your parameters

`a11`

,`c11`

, and`e11`

represent what the loading (path coefficient) onto (respectively) A, C, and E will be when age equals zero. Your parameters`aM11`

,`cM11`

, and`eM11`

represent how much their corresponding loading changes for a 1-unit change in age. Note that the resulting loadings themselves are ambiguous with respect to sign. For instance, the A variance component is`(a+ D%`

, which in this case is equivalent to squaring a scalar, so it will always be non-negative.%aM) %% t(a+ D%*%aM)Thanks Rob, that's helpful!