Issues with confidence interval estimation for Cholesky Decomposition ACE model

Posted on
No user picture. lyra9826 Joined: 10/20/2021
I am estimating a Cholesky decomposition ACE model and deriving the rA, rC, and rE estimates. The confidence intervals are not estimating for many of these parameters. There appear to be varying diagnoses when I use verbose=T (i.e., active box constraint, alpha level not reached infeasible non-linear constraint). I have tried choosing a different optimizer, re-scaling the phenotypes, using mxTryHard, and changing start values. The univariate variance components as well as their correlations are significant. I would appreciate any guidance regarding how to obtain these confidence intervals.

Please note we will also look at the correlated ACEs with direct variance but this is in response to a reviewer request for the Cholesky.

Replied on Wed, 03/01/2023 - 11:52
Picture of user. AdminRobK Joined: 01/24/2014

Could you provide your verbose `summary()` output?
Replied on Fri, 03/03/2023 - 10:31
Picture of user. AdminNeale Joined: 03/01/2013

Hi

You may be running into issues with a path-specified ACE model where an estimate of -.5 for a results in the same fit as +.5. In these circumstances, the lower CI may end up at negative the upper CI. To avoid this you can bound a c and e parameters to be positive. It is also worth noting that Cholesky is not a recommended multivariate model, because it has implicit constraints (to keep things positive definite) and incorrect Type I error rates when, e.g., dropping the C matrix. See Verhulst et al 2019 for an alternative specification that does not have these problems.