umx (version 1.4.0)

umxPlotACEcov: umxPlotACEcov

Description

Make a graphical display of an ACE model

Usage

umxPlotACEcov(x = NA, file = "name", digits = 2, means = FALSE, std = TRUE, showMeans = NULL, ...)

Arguments

x

mxModel to plot (created by umxACE in order to inherit the MxModel.ACE class)

file

The name of the dot file to write: NA = none; "name" = use the name of the model

digits

How many decimals to include in path loadings (default is 2)

means

Whether to show means paths (default is FALSE)

std

Whether to standardize the model (default is TRUE)

showMeans

DEPRECATED use "means" instead

...

Additional (optional) parameters

Value

- optionally return the dot code

References

- http://tbates.github.io

See Also

Other Plotting functions: plot.MxModel, umxPlotACE, umxPlotCP, umxPlotGxE, umxPlotIP

Other Reporting functions: RMSEA.MxModel, RMSEA.summary.mxmodel, RMSEA, confint.MxModel, extractAIC.MxModel, loadings, logLik.MxModel, plot.MxModel, residuals.MxModel, umxCI_boot, umxCI, umxCompare, umxExpCov, umxExpMeans, umxFitIndices, umxPlotACE, umxPlotCP, umxPlotGxE, umxPlotIP, umxSummary.MxModel, umxSummaryACE, umx_drop_ok, umx_standardize_RAM

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

Examples

require(umx)
# BMI ?twinData from Australian twins. 
# Cohort 1 Zygosity 1 == MZ females 3 == DZ females
data(twinData)
# Pick the variables. We will use base names (i.e., "bmi") and set suffix.
selDVs  = c("bmi")
selCovs = c("age")
selVars = umx_paste_names(c(selDVs, selCovs), sep = "", suffixes= 1:2)
# just top few pairs so example runs quickly
mzData = subset(twinData, zyg == 1, selVars)[1:100, ]
dzData = subset(twinData, zyg == 3, selVars)[1:100, ]
# TODO update for new dataset variable zygosity
# mzData = subset(twinData, zygosity == "MZFF", selVars)[1:200, ]
# dzData = subset(twinData, zygosity == "DZFF", selVars)[1:200, ]
m1 = umxACEcov(selDVs = selDVs, selCovs = selCovs, dzData = dzData, mzData = mzData, 
	 suffix = "", autoRun = TRUE)
plot(m1)
plot(m1, std = FALSE) # don't standardize