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.IACE.A))%%ACE.A%% solve(sqrt(ACE.IACE.A)), name="corComp_A")
corComp_E <- mxAlgebra(solve(sqrt(ACE.IACE.E))%%ACE.E%% solve(sqrt(ACE.IACE.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?