umx (version 1.4.0)

umxSummaryGxE: umxSummaryGxE

Description

Summarise a Moderation model, as returned by umxGxE

Usage

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("1", "2", "html"), ...)

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

default = F

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

"1" = regular, "2" = add descriptive sentences; "html" = open a browser and copyable tables

...

Optional additional parameters

Value

- optional mxModel

References

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

See Also

- umxGxE(), plot(), umxSummary() work for IP, CP, GxE, and ACE models.

Other Twin Modeling Functions: plot.MxModel, umxACESexLim, umxACEcov, umxACE, umxCF_SexLim, umxCP, umxGxE_window, umxGxE, umxIP, umxPlotACEcov, umxPlotCP, umxPlotGxE, umxPlotIP, umxSummaryACEcov, umxSummaryACE, umxSummaryCP, umxSummaryIP, umx_make_TwinData, umx

Examples

# 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) 
labList = c("MZFF", "MZMM", "DZFF", "DZMM", "DZOS")
twinData$ZYG = factor(twinData$zyg, levels = 1:5, labels = labList)
twinData$age1 = twinData$age2 = twinData$age
selDVs  = c("bmi1", "bmi2")
selDefs = c("age1", "age2")
selVars = c(selDVs, selDefs)
mzData  = subset(twinData, ZYG == "MZFF", selVars)
dzData  = subset(twinData, ZYG == "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 = selDVs, selDefs = selDefs, dzData = dzData, mzData = mzData)
m1 = umxRun(m1)
# Plot Moderation
umxSummaryGxE(m1)
umxSummaryGxE(m1, location = "topright")
umxSummaryGxE(m1, separateGraphs = FALSE)