umxSummaryGxE: Summarize a GxE model

In tbates/umx: Structural Equation Modeling and Twin Modeling in R

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

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.

Value

• optional OpenMx::mxModel()

References

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

See Also

• umxGxE(), umxReduce(), plot(), [umxSummary)] all work for IP, CP, GxE, and ACE models.

[umxSummary)]: R:umxSummary)

Other Summary functions: umxSummary.MxModel(), umxSummaryACEcov(), umxSummaryCP(), umxSummaryIP(), umxSummaryMRDoC()

Examples

## 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]]),]
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)

tbates/umx documentation built on Dec. 14, 2024, 11:28 a.m.