Published on *OpenMx* (http://openmx.psyc.virginia.edu)

By *spanosal*

Created *10/04/2010 - 10:08*

Mon, 10/04/2010 - 10:08 — spanosal [1]

I am running a trivariate Cholesky decomposition using a script from the workshop. This script did not include CIs but I included it from another source. However, while the CI script allows me to run CIs for the genetic and environmental correlations, I can't figure out how to modify the script so I can run CIs for my standardized path estimates. The portion of my script that runs the CIs for correlations is below:

corComp_A <- mxAlgebra(solve(sqrt(ACE.I*ACE.A))%*%ACE.A%*% solve(sqrt(ACE.I*ACE.A)), name="corComp_A")

corComp_E <- mxAlgebra(solve(sqrt(ACE.I*ACE.E))%*%ACE.E%*% solve(sqrt(ACE.I*ACE.E)),name="corComp_E")

#Add GE correlations to the model, and ask the model to report CIs when it runs

CImodel <- mxModel(multAEFit, corComp_A, corComp_E, mxCI(c("corComp_A", "corComp_E")))

#Run the model

CImodelFit <- mxRun(CImodel, intervals=TRUE)

CIAEsumm <- summary(CImodelFit)CI

CIAEsumm

Any suggestions?

**Links:**

[1] http://openmx.psyc.virginia.edu/users/spanosal

[2] http://openmx.psyc.virginia.edu/thread/1055

[3] http://openmx.psyc.virginia.edu/thread/1006

[4] http://openmx.psyc.virginia.edu/forums/openmx-help/openmx-structural-equation-modeling