Summarize a GxE model
Summarize a genetic moderation model, as returned by umxGxE(). Prints graphs of A, C, and E, standardized and raw.
umxSummaryGxE(
model = NULL,
digits = 2,
xlab = NA,
location = "topleft",
separateGraphs = FALSE,
gg = TRUE,
file = getOption("umx_auto_plot"),
returnStd = NULL,
std = NULL,
reduce = FALSE,
CIs = NULL,
report = c("markdown", "html"),
show = NULL,
...
)model |
A fitted |
digits |
round to how many digits (default = 2) |
xlab |
label for the x-axis of plot |
location |
default = "topleft" |
separateGraphs |
If TRUE, both std and raw plots in one figure (default FALSE) |
gg |
Whether to use ggplot to create the graphs (default TRUE) |
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) |
std |
Whether to show the standardized model (not implemented! TRUE) |
reduce |
Whether run and tabulate a complete model reduction...(Defaults to FALSE) |
CIs |
Confidence intervals (FALSE) |
report |
"markdown" or "html" = open a browser for copyable tables |
show |
not doing anything yet (required for all summary functions) |
... |
Optional additional parameters |
Note: see also umxReduce() which knows how to reduce a GxE model.
optional mxModel()
umxGxE(), umxReduce(), plot(), [umxSummary)] all work for IP, CP, GxE, and ACE models.
[umxSummary)]: R:umxSummary)
Other Summary functions:
umxSummary.MxModel(),
umxSummaryACEcov(),
umxSummaryCP(),
umxSummaryIP()
# The total sample has been subdivided into a young cohort,
# aged 18-30 years, and an older cohort aged 31 and above.
# Cohort 1 Zygosity is coded as follows 1 == MZ females 2 == MZ males
# 3 == DZ females 4 == DZ males 5 == DZ opposite sex pairs
require(umx)
data(twinData)
twinData$age1 = twinData$age2 = twinData$age
selDVs = c("bmi1", "bmi2")
selDefs = c("age1", "age2")
selVars = c(selDVs, selDefs)
mzData = subset(twinData, zygosity == "MZFF", selVars)
dzData = subset(twinData, zygosity == "DZMM", selVars)
# Exclude cases with missing Def
mzData = mzData[!is.na(mzData[selDefs[1]]) & !is.na(mzData[selDefs[2]]),]
dzData = dzData[!is.na(dzData[selDefs[1]]) & !is.na(dzData[selDefs[2]]),]
## Not run:
m1 = umxGxE(selDVs = "bmi", selDefs = "age", sep="", dzData = dzData, mzData = mzData)
# Plot Moderation
umxSummaryGxE(m1)
umxSummaryGxE(m1, location = "topright")
umxSummaryGxE(m1, separateGraphs = FALSE)
## End(Not run)Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.