umx (version 1.7.5)

umxSummaryCP: umxSummaryCP

Description

Summarise a Commmon Pathway model, as returned by umxCP

Usage

umxSummaryCP(model, digits = 2, file = umx_set_auto_plot(silent = TRUE),
  returnStd = FALSE, extended = FALSE, showRg = TRUE, comparison = NULL,
  std = TRUE, CIs = FALSE, ...)

Arguments

model

A fitted umxCP model to summarize

digits

round to how many digits (default = 2)

file

The name of the dot file to write: NA = none; "name" = use the name of the model

returnStd

Whether to return the standardized form of the model (default = FALSE)

extended

how much to report (FALSE)

showRg

Whether to show the genetic correlations (FALSE) (Not implemented!)

comparison

Whether to run mxCompare on a comparison model (NULL)

std

Whether to show the standardized model (TRUE) (ignored: used extended = TRUE to get unstandardized)

CIs

Confidence intervals (F)

...

Optional additional parameters

Value

- optional mxModel

References

- http://www.github.com/tbates/umx, http://tbates.github.io

See Also

- umxCP(), plot(), umxSummary() work for IP, CP, GxE, SAT, and ACE models.

Other Twin Modeling Functions: plot.MxModel, umxACESexLim, umxACEcov, umxACE, umxCF_SexLim, umxCP, umxGxE_window, umxGxE, umxIP, umxPlotACEcov, umxPlotCP, umxPlotGxE, umxPlotIP, umxSummaryACEcov, umxSummaryACE, umxSummaryGxE, umxSummaryIP, umx

Examples

# NOT RUN {
require(umx)
data(twinData) 
zygList = c("MZFF", "MZMM", "DZFF", "DZMM", "DZOS")
twinData$ZYG = factor(twinData$zyg, levels = 1:5, labels = zygList)
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, ZYG == "MZFF", umx_paste_names(selDVs, "", 1:2))
dzData <- subset(twinData, ZYG == "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)
# }