Issues with confidence interval estimation for Cholesky Decomposition ACE model
Posted on
lyra9826
Joined: 10/20/2021
Attachment | Size |
---|---|
twoACE_Cholesky.R | 8.57 KB |
Forums
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.
output?
Log in or register to post comments
Invariance to sign
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.
Log in or register to post comments