# NOT RUN {
require(umx)
data(twinData)
twinData$wt1 = twinData$wt1/10 # help CSOLNP by putting wt on a similar scale to ht
twinData$wt2 = twinData$wt2/10 # help CSOLNP by putting wt on a similar scale to ht
selDVs = c("ht", "wt")
mzData <- subset(twinData, zygosity == "MZFF", umx_paste_names(selDVs, "", 1:2))
dzData <- subset(twinData, zygosity == "DZFF", umx_paste_names(selDVs, "", 1:2))
m1 = umxCP(selDVs = selDVs, dzData = dzData, mzData = mzData, suffix = "")
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 = TRUE,
extended = FALSE, showRg = TRUE, std = TRUE, CIs = TRUE);
umxSummaryCP(m1, ext = TRUE, file = "name")
umxSummaryCP(m1, file = "Figure 3", std = TRUE)
# }
Run the code above in your browser using DataLab