Thanks. It is another milestone!
Unfortunately, it broke my metaSEM package. I am trying to fix it. I am testing a random-effects meta-analysis model. This model is equivalent to the one in the sources /openmx/models/passing/RAM-FIML-defvars.R There are only two parameters---one mean and one variance. The models are specified in OpenMx and in metaSEM that generates the code for OpenMx. The R syntax, outputs, and R objects are test2.0.R, “test2.0 output.txt”, and test2.0.txt, respectively.
As shown in “test2.0 output.txt”, the -2log likelihood and the parameter estimates are the same for them. However, there is no standard errors reported in the version created by metaSEM.
When I check the output slots, the calculated Hessian are different:
Specified in OpenMx:
test.Tau[1,1] 371.31703 -16.31094
test.M[1,1] -16.31094 181.77602
Specified in metaSEM:
Intercept1 -207769350869 144970978723
Tau2_1_1 144970978723 -404613761388
I have no idea on how to proceed. Any suggestions are highly appreciated! Thanks in advance.
|test2.0 output.txt ||9.62 KB|
|test2.0.R ||1.37 KB|
|test2.0.txt ||4.56 KB|