umxSummary reports the rA, rC, and rE for models. But this is done in the summary, and doesn't include CIs by default.
It’s computed as
solve(sqrt(I*A)) %*% A %*% solve(sqrt(I*A))
So you can add an mxAlgebra computing that into the model, then do mxSE on that
m2 = mxModel(m1, mxAlgebra(name="ra", solve(sqrt(top.I*top.A)) %*% top.A %*% solve(sqrt(top.I*top.A))) ) m2 = mxRun(m2) mxSE(ra, m2)
That returns SEs, and CI= ± 1.96 * SE
best, t