umx (version 1.2.7)

umxACEcov: umxACEcov

Description

Make a 2-group ACE Cholesky Twin model with covariates modeled (see Details below)

Usage

umxACEcov(name = "ACEcov", selDVs, selCovs, dzData, mzData, suffix = NULL, dzAr = 0.5, dzCr = 1, addStd = TRUE, addCI = TRUE, boundDiag = NULL, equateMeans = TRUE, bVector = FALSE, hint = c("none", "left_censored"), autoRun = getOption("umx_auto_run"))

Arguments

name

The name of the model (defaults to"ACE")

selDVs

The variables to include from the data (do not include suffixes)

selCovs

The covariates to include from the data (do not include suffixes)

dzData

The DZ dataframe

mzData

The MZ dataframe

suffix

suffix for twin 1 and twin 2, often "_T" Used to expand selDVs into full column names, i,e "dep" --> c("dep_T1", "dep_T2")

dzAr

The DZ genetic correlation (defaults to .5, vary to examine assortative mating)

dzCr

The DZ "C" correlation (defaults to 1: set to .25 to make an ADE model)

addStd

Whether to add the algebras to compute a std model (defaults to TRUE)

addCI

Whether to add intervals to compute CIs (defaults to TRUE)

boundDiag

= Whether to bound the diagonal of the a, c, and e matrices

equateMeans

Whether to equate the means across twins (defaults to TRUE)

bVector

Whether to compute row-wise likelihoods (defaults to FALSE)

hint

An analysis hint. Options include "none", (default) "left_censored". Default does nothing.

autoRun

Whether to run the model and return it, or just return it

Value

- mxModel of subclass mxModel.ACEcov

Details

umxACEcov supplements the umxACE Cholesky model with covariates.

References

- Neale, M. C., & Martin, N. G. (1989). The effects of age, sex, and genotype on self-report drunkenness following a challenge dose of alcohol. Behavior Genetics, 19(1), 63-78. doi:10.1007/BF01065884

See Also

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

Examples

require(umx)
data(twinData)
# add age 1 and age 2 columns
twinData$age1 = twinData$age2 = twinData$age
selDVs  = c("bmi")
selCovs = c("age")
selVars = umx_paste_names(c(selDVs, selCovs), textConstant = "", suffixes= 1:2)
# just top 80 so example runs in a couple of secs
mzData = subset(twinData, zyg == 1, selVars)[1:80, ]
dzData = subset(twinData, zyg == 3, selVars)[1:80, ]
# TODO update for new dataset variable zygosity
# mzData = subset(twinData, zygosity == "MZFF", selVars)[1:80, ]
# dzData = subset(twinData, zygosity == "DZFF", selVars)[1:80, ]
m1 = umxACEcov(selDVs = selDVs, selCovs = selCovs, dzData = dzData, mzData = mzData, 
	 suffix = "", autoRun = TRUE)
umxSummary(m1)
plot(m1)