Attachment | Size |
---|---|

test2.0 output.txt [6] | 9.62 KB |

test2.0.R [7] | 1.37 KB |

test2.0.txt [8] | 4.56 KB |

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:

$calculatedHessian

test.Tau[11] test.M[1,1]

test.Tau[1,1] 371.31703 -16.31094

test.M[1,1] -16.31094 181.77602

Specified in metaSEM:

$calculatedHessian

Intercept1 Tau2_1_1

Intercept1 -207769350869 144970978723

Tau2_1_1 144970978723 -404613761388

I have no idea on how to proceed. Any suggestions are highly appreciated! Thanks in advance.