You are here

multivariate (4) ACE model confidence intervals for standardized A,C,E

3 posts / 0 new
Last post
jwee2's picture
Offline
Joined: 07/01/2014 - 12:57
multivariate (4) ACE model confidence intervals for standardized A,C,E

Greetings, fairly new student of OpenMx here. I've successfully run a univariate twin model with confidence intervals thanks to mxCI, but in the quadvariate script below (obtained from http://ibg.colorado.edu/cdrom2014/bartels/Multivariate/Trivariate.R), I don't know where to put the mxCI command to request CIs for A/V, C/V, and E/V. The output is in the form of a matrix too, so I am unsure how to proceed. Any assistance would be greatly appreciated.

# Matrices declared to store a, c, and e Path Coefficients
pathA values=svPa, labels=labLower("a",nv), lbound=lbPa, name="a" )
pathC values=svPa, labels=labLower("c",nv), lbound=lbPa, name="c" )
pathE values=svPa, labels=labLower("e",nv), lbound=lbPa, name="e" )

# Matrices generated to hold A, C, and E computed Variance Components
covA covC covE

# Algebra to compute total variances and standard deviations (diagonal only)
covP matI invSD

# Algebra for expected Mean and Variance/Covariance Matrices in MZ & DZ twins
meanG values=svMe, labels=labFull("me",1,nv), name="expMean" )
covMZ covDZ

# Data objects for Multiple Groups
dataMZ dataDZ

# Objective objects for Multiple Groups
objMZ objDZ

# Combine Groups
pars modelMZ modelDZ minus2ll obj CholAceModel

# ------------------------------------------------------------------------------
# RUN GENETIC MODEL

# Run Cholesky Decomposition ACE model
CholAceFit

jwee2's picture
Offline
Joined: 07/01/2014 - 12:57
An update to my situation, I

An update to my situation, I discovered that the summary of the multivariate fit path model through summary(CholeAceFit) provides 'free parameter' estimates for all paths, including their standard error. My advisor says that presenting un-standardized estimates in a path model with confidence intervals calculated using those standard errors provided is fine. So although the original question was not solved, my problems overall have been remedied! :)

mhunter's picture
Offline
Joined: 07/31/2009 - 15:26
Try this

The following code should get you confidence intervals on the StandardizedA algebra in the submodel MZ.

avMZ <- mxAlgebra(A/V, name="StandardizedA")
avMZ_CI <- mxCI("StandardizedA")
modelMZ <- mxModel( pars, covMZ, dataMZ, objMZ, avMZ, avMZ_CI, name="MZ" )

Similar expressions would be possible for other algebras. Confidence intervals work on "named entities" within a model. So you can get a confidence interval on any algebra, matrix, or free parameter by creating an mxCI object that refers to the entity by name and putting the mxCI in the same model as the entity.

Log in or register to post comments