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,
file = getOption("umx_auto_plot"),
returnStd = NULL,
std = NULL,
reduce = FALSE,
CIs = NULL,
report = c("markdown", "html"),
show = NULL,
...
)
A fitted umxGxE()
model to summarize
round to how many digits (default = 2)
label for the x-axis of plot
default = "topleft"
If TRUE, both std and raw plots in one figure (default FALSE)
The name of the dot file to write: NA = none; "name" = use the name of the model
Whether to return the standardized form of the model (default = FALSE)
Whether to show the standardized model (not implemented! TRUE)
Whether run and tabulate a complete model reduction...(Defaults to FALSE)
Confidence intervals (FALSE)
"markdown" or "html" = open a browser for copyable tables
not doing anything yet (required for all summary functions)
Optional additional parameters
optional mxModel()
Note: see also umxReduce()
which knows how to reduce a GxE model.
umxGxE()
, umxReduce()
, plot()
, [umxSummary)] all work for IP, CP, GxE, and ACE models.
[umxSummary)]: R:umxSummary)
Other Twin Reporting Functions:
umxPlotCP()
,
umxPlotDoC()
,
umxReduceACE()
,
umxReduceGxE()
,
umxReduce()
,
umxSummarizeTwinData()
,
umxSummaryACEcov()
,
umxSummaryACEv()
,
umxSummaryACE()
,
umxSummaryCP()
,
umxSummaryDoC()
,
umxSummaryGxEbiv()
,
umxSummaryIP()
,
umxSummarySexLim()
,
umxSummarySimplex()
,
umx
# NOT RUN {
# 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 = selDVs, selDefs = selDefs, dzData = dzData, mzData = mzData)
# Plot Moderation
umxSummaryGxE(m1)
umxSummaryGxE(m1, location = "topright")
umxSummaryGxE(m1, separateGraphs = FALSE)
# }
Run the code above in your browser using DataLab