plot.MxModelTwinMaker: Create and display a graphical path diagram for a path-based...

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

Create and display a graphical path diagram for a path-based twin model.

Description

Assumes the model has a group called "MZ" inside.

Usage

## S3 method for class 'MxModelTwinMaker'
plot(
  x = NA,
  std = FALSE,
  fixed = TRUE,
  means = TRUE,
  oneTwin = TRUE,
  sep = "_T",
  digits = 2,
  file = "name",
  labels = c("none", "labels", "both"),
  resid = c("circle", "line", "none"),
  strip_zero = FALSE,
  splines = TRUE,
  min = NULL,
  same = NULL,
  max = NULL,
  ...
)

Arguments

x	A umxTwinMaker() model from which to make a path diagram
std	Whether to standardize the model (default = FALSE)
fixed	Whether to show fixed paths (defaults to TRUE)
means	Whether to show means or not (default = TRUE)
oneTwin	(whether to plot a pair of twins, or just one (default = TRUE)
sep	The separator for twin variables ("_T")
digits	The number of decimal places to add to the path coefficients
file	The name of the dot file to write: NA = none; "name" = use the name of the model
labels	Whether to show labels on the paths. "none", "labels", or "both" (parameter + label).
resid	How to show residuals and variances default is "circle". Options are "line" & "none"
strip_zero	Whether to strip the leading "0" and decimal point from parameter estimates (default = FALSE)
splines	Whether to allow lines to curve: defaults to TRUE (nb: some models look better with FALSE)
min	optional list of objects to group at the top of the plot. Default (NULL) chooses automatically.
same	optional list of objects to group at the same rank in the plot. Default (NULL) chooses automatically.
max	optional list of objects to group at the bottom of the plot. Default (NULL) chooses automatically.
...	Optional parameters

Details

If you use umx_set_plot_format("graphviz"), they will open in a graphviz helper app (if installed). The commercial application “OmniGraffle” is great for editing these images. On unix and windows, plot() will create a pdf and open it in your default pdf reader.

See Also

• umx_set_plot_format(), plot.MxModel(), umxPlotACE(), umxPlotCP(), umxPlotIP(), umxPlotGxE()

Other Plotting functions: ggAddR(), plot.MxLISRELModel(), plot.MxModel(), umx, umxPlot(), umxPlotACE(), umxPlotACEcov(), umxPlotACEv(), umxPlotCP(), umxPlotDoC(), umxPlotFun(), umxPlotGxE(), umxPlotGxEbiv(), umxPlotIP(), umxPlotSexLim(), umxPlotSimplex()

Examples

## Not run: 
require(umx)
# 
# =====================
# = Make an ACE model =
# =====================
# 1. Clean data: Add separator and scale
data(twinData)
tmp = umx_make_twin_data_nice(data=twinData, sep="", zygosity="zygosity", numbering=1:2)
tmp = umx_scale_wide_twin_data(varsToScale= c("wt", "ht"), sep= "_T", data= tmp)
mzData = subset(tmp, zygosity %in%  c("MZFF", "MZMM"))
dzData = subset(tmp, zygosity %in%  c("DZFF", "DZMM"))

# 2. Define paths: You only need the paths for one person:
paths = c(
umxPath(v1m0 = c("a1", 'c1', "e1")),
umxPath(means = c("wt")),
umxPath(c("a1", 'c1', "e1"), to = "wt", values=.2)
)
m1 = umxTwinMaker("test", paths, mzData = mzData, dzData= dzData)
plot(m1, std= TRUE, means= FALSE)
plot(m1, means=FALSE, std=TRUE, strip=TRUE, splines="FALSE", max="intercept")

## End(Not run) # end dontrun

# =================
# = An ACEv model =
# =================
# Not complete

paths = c(
umxPath(v1m0 = c("A1", 'C1', "E1")),
umxPath(v1m0 = c("A2", 'C2', "E2")),
umxPath(v.m0 = c("l1", 'l2')),
umxPath(v.m. = c("wt", "ht")),
umxPath(c("A1", 'C1', "E1"), to = "l1", values= .2),
umxPath(c("A2", 'C2', "E2"), to = "l2", values= .2),
umxPath(c("l1", 'l2'), to = c("wt", "ht"), values= .2)
)

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