require(umx)
data(twinData)
labList = c("MZFF", "MZMM", "DZFF", "DZMM", "DZOS")
twinData$ZYG = factor(twinData$zyg, levels = 1:5, labels = labList)
selDVs = c("ht", "wt") # will be expanded into "ht1", "wt1", "ht2", "wt2"
mzData <- subset(twinData, ZYG == "MZFF")
dzData <- subset(twinData, ZYG == "DZFF")
m1 = umxCP(selDVs = selDVs, dzData = dzData, mzData = mzData, suffix = "")
m1 = mxRun(m1)
umxSummaryCP(m1, file = NA) # suppress plot creation with file
umxSummary(m1, file = NA) # generic summary is the same
stdFit = umxSummaryCP(m1, digits = 2, file = NA, returnStd = FALSE,
extended = FALSE, showRg = TRUE, showStd = TRUE, CIs = TRUE);
examples which will create graphical output
umxSummaryCP(fit);
umxSummaryCP(fit, file = "Figure 3", showStd = TRUE)Run the code above in your browser using DataLab