umx (version 1.9.1)

umxReduceACE: Reduce an ACE model.

Description

This function can perform model reduction on umxACE models, testing dropping A and C, as well as an ADE or ACE model, displaying the results in a table, and returning the best model.

Usage

umxReduceACE(model, report = c("markdown", "inline", "html", "report"),
  baseFileName = "tmp", intervals = TRUE, ...)

Arguments

model

an ACE or ADE mxModel to reduce

report

How to report the results. "html" = open in browser

baseFileName

(optional) custom filename for html output (defaults to "tmp")

intervals

Recompute CIs (if any included) on the best model (default = TRUE)

...

Other parameters to control model summary

Value

Best fitting model

Details

It is designed for testing univariate models. You can offer up either the ACE or ADE base model.

Suggestions for more sophisticated automation welcomed!

References

- http://tbates.github.io

See Also

Other Core Modelling Functions: plot.MxModel, umxDiagnose, umxLatent, umxMatrix, umxPath, umxRAM, umxReduceGxE, umxRun, umxSuperModel, umx

Other Twin Modeling Functions: plot.MxModel, umxACESexLim, umxACE_cov_fixed, umxACEcov, umxACEv, umxACE, umxCF_SexLim, umxCP, umxGxE_window, umxGxE, umxIP, umxPlotACEcov, umxPlotCP, umxPlotGxE, umxPlotIP, umxSummaryACEcov, umxSummaryACEv, umxSummaryACE, umxSummaryCP, umxSummaryGxE, umxSummaryIP, umx_long2wide, umx_wide2long, umx, xmu_twin_check

Examples

# NOT RUN {
data(twinData)
mzData <- subset(twinData, zygosity == "MZFF")
dzData <- subset(twinData, zygosity == "DZFF")
m1 = umxACE(selDVs = "bmi", dzData = dzData, mzData = mzData, sep = "")
m2 = umxReduce(m1)
umxSummary(m2)
m1 = umxACE(selDVs = "bmi", dzData = dzData, mzData = mzData, sep = "", dzCr = .25)
m2 = umxReduce(m1)
# }