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?