G*E interaction with a moderator

Posted on
No user picture. Qiuzhi Xie Joined: 12/27/2018
I am currently doing G*E interaction with SES as a moderator. The plot (attached here) suggests that SES likely moderates with the A factor. To analyze whether the moderation effect is significant on A factor, I drop the moderator on the Path A only and compare the model fit (diffLL). I find that the p value is insignificant. I want to ask whether this is the correct way to see the significance of moderator effect.

Also, is it appropriate to see the unstandardized variance components by SES moderation, rather than the standardized variance components?

Replied on Wed, 03/13/2019 - 09:46
Picture of user. tbates Joined: 07/31/2009

The plot can be misleading: if the moderator on the Path A doesn't improve fit significantly, then that tells us the plot apparently sloping around reflected the parameter being estimated at a non-zero value, but including other values including zero in its CI.

Not sure what you mean by is it appropriate to see the unstandardized variance components by SES moderation. it's often very informative.

umxGxE does this model and the plots for you.

Also umxReduce() will reduce the model in ways that typically make sense for a table in a journal.

Replied on Wed, 03/13/2019 - 10:14
No user picture. Qiuzhi Xie Joined: 12/27/2018

In reply to by tbates

Many thanks for the reply! In another example that the slope in the plot is not so obvious (aMLL value is smaller) but exactly linear, the Δ-2LL is significant. Therefore, I wonder whether linear moderator effect is be more easily recognized and quadratic moderator effect is more likely be recognized as insignificant by comparing model fit?
Replied on Wed, 03/13/2019 - 13:25
Picture of user. tbates Joined: 07/31/2009

The moderation is a linear effect in path coefficient terms, and the effect on variance is squared, hence the common curvature.
(PS: umxGxE also includes quadratic effects on means).

If you want to explore non-linear moderator effects, I would use losem [1], implemented in umx as umxGxE_window()

If one had a particular model in mind (e.g. Eric Turkheimer and the late Irv Gottesman proposed a "bent chicken wire" model in which h^2 increased until a "good enough environment" was reached, after which h^2 was no longer coupled to environment. One could estimate that inflection point in the model, changing the algebra to include this new form of moderated effect. Likewise for other models of GxE/CxE and other shapes of function (growth, inverted-U etc.)

[1] Briley, D. A., Harden, K. P., Bates, T. C., & Tucker-Drob, E. M. (2015). Nonparametric Estimates of Gene x Environment Interaction Using Local Structural Equation Modeling. Behavior Genetics, 45(5), 581-596. doi:10.1007/s10519-015-9732-8
[pubmed](http://www.ncbi.nlm.nih.gov/pubmed/26318287)

[pdf](https://labs.la.utexas.edu/harden/files/2013/05/Briley-et-al-2015-LOSEM-for-GxE.pdf)