I am very new with ACE twin modeling as well as R, therefore I need some help.
I am trying to run a univariate ACE model using the attached script (UnivACE.R) in OpenMx. It would be very grateful to know whether the script does not have any mistakes. And one more thing, what modifications are needed in the script to get the p-values of the estimated h2,c2 and e2 in ACE model? Furthermore, when using the same dataset to run ACE model, Mplus reported different CIs for the path coefficients compare to that in OpenMX results. Which one is correct, and why?
Another problem is how to set the start value, since I found this sometimes has a great influence on the final results.
Thanks a lot for your help!