umxSummaryGxE: Summarize a GxE model

Description

Summarize a genetic moderation model, as returned by umxGxE(). Prints graphs of A, C, and E, standardized and raw.

Usage

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,
  ...
)

Value

optional mxModel()

Arguments

model: A fitted umxGxE() model to summarize
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

Details

Note: see also umxReduce() which knows how to reduce a GxE model.

References

https://github.com/tbates/umx, https://tbates.github.io

Examples

Run this code

if (FALSE) {
# 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]]),]
m1 = umxGxE(selDVs = "bmi", selDefs = "age", sep="", 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

The Learning Leader's Guide to AI Literacy