Hi,
I don't know if I'm missing something basic here, but how does one apply constraints on means and variances across more than 2 groups (mxConstraint only accepts 2 arguments) in a saturated model, ie across sex and zygosity groups?
My data is setup such:
mzfData <- subset(twinsib.wide, ZygosityNumber ==1, selVars)
mzmData <- subset(twinsib.wide, ZygosityNumber ==2, selVars)
dzfData <- subset(twinsib.wide, ZygosityNumber ==3, selVars)
dzmData <- subset(twinsib.wide, ZygosityNumber ==4, selVars)
dzfmData <- subset(twinsib.wide, ZygosityNumber ==5, selVars)
dzmfData <- subset(twinsib.wide, ZygosityNumber ==6, selVars)
multiTwinSatModel <- mxModel("multiTwinSat",
mxModel("MZF",
...
),
mxModel("MZM",
...
),
mxModel("DZF",
...
),
etc...
),
I haven't nested individual models by zygosity or sex.
Thank you.
Paul