Dear umx Experts,
I ran an ADE Cholesky model with 4 variables. When trying to obtain CIs for the one-headed paths (using the command umxConfint), I got weird CIs and the note "alpha level not reached" near some of the CIs. After consultation in the behavioral genetics forum, I understood that this the result of an optimization problem.
AdminRobK suggested that instead of looking at the CIs, I will try to get robust standard errors, then use mxSE() to get standard errors for the elements and form a CI from those standard errors
My question is- how can I obtain robust standard errors for a model that was built in umx (as opposed to openMx)?
Note: I tried to use the imxRobustSE() command, but I got the following error:
Error: Unknown expected covariance name 'top.expCovMZ' detected in the expectation function of model 'MZ'
My guess is that I need to do some sort of adaptation to the model before running this command. Any suggestion will be very appreciated.
Thank you very much,
Lior