# NOT RUN {
require(umx)
data(twinData)
selDVs = c("ht1", "wt1", "ht2", "wt2")
mzData <- subset(twinData, zygosity == "MZFF")
dzData <- subset(twinData, zygosity == "DZFF")
m1 = umxIP(selDVs = selDVs, dzData = dzData, mzData = mzData)
m1 = umxRun(m1)
umxSummaryIP(m1)
plot(m1)
# }
# NOT RUN {
umxSummaryIP(m1, digits = 2, file = "Figure3", showRg = FALSE, CIs = TRUE);
# }
Run the code above in your browser using DataLab