require(umx)
data(twinData)
labList = c("MZFF", "MZMM", "DZFF", "DZMM", "DZOS")
twinData$ZYG = factor(twinData$zyg, levels = 1:5, labels = labList)
selDVs = c("bmi1", "bmi2")
mzData <- subset(twinData, ZYG == "MZFF", selDVs)
dzData <- subset(twinData, ZYG == "DZFF", selDVs)
m1 = umxACE(selDVs = selDVs, dzData = dzData, mzData = mzData)
m1 = umxRun(m1)
umxSummaryACE(m1)
## Not run:
# umxSummaryACE(m1, file = NA);
# umxSummaryACE(m1, file = "name", std = TRUE)
# stdFit = umxSummaryACE(m1, returnStd = TRUE);
# ## End(Not run)
Run the code above in your browser using DataLab