umxSummaryGxE: Summarize a GxE model

umxSummaryGxER Documentation

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 mxModel()

References

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 March 24, 2024, 8:12 a.m.