Closed tbates closed 3 years ago
closing point 1: a table of parameters
require(umx); data(twinData)
twinData$age1 = twinData$age2 = twinData$age
mzData = subset(twinData, zygosity == "MZFF")
dzData = subset(twinData, zygosity == "DZFF")
m1 = umxGxE(selDVs = "bmi", selDefs = "age", dzData = dzData, mzData = mzData, sep= "", try="yes")
selDVs: 'bmi1' and 'bmi2'
1 dz rows dropped due to missing def var for Twin 1 or Twin 2 or both
Running G_by_E with 9 parameters
name | Estimate | Std.Error |
---|---|---|
a_r1c1 | 0.713 | 0.042 |
c_r1c1 | 0.000 | 0.300 |
e_r1c1 | 0.292 | 0.030 |
am_r1c1 | 0.002 | 0.001 |
cm_r1c1 | 0.000 | 0.007 |
em_r1c1 | 0.006 | 0.001 |
mean | 20.590 | 0.134 |
lin11 | 0.038 | 0.007 |
quad11 | 0.000 | 0.000 |
version 4.2.0 adds CI ribbons
Thanks to #139 the user can now also customise the output:
tmp[[2]]+ylim(c(0,1.5))
Yesssss. Great stuff!
custom moderator range
umxGxE