Copyright © 2007-2024 The OpenMx Project
Hi, The current trunk of OpenMx (999.0.0-3521) seems to be not estimating parameters correctly.
I've uploaded a gist that shows the issues here: It's a simple independence model with 3 variables in the built-in mtcars data set.
https://gist.github.com/tbates/d2261e8c5daf426e1c8e
You can see that the means are borderline OK, but the variances are way off.
# =================================== # = Means and variances in the data = # =================================== round(diag(cov(mtcars[,manifests])),2) mpg disp gear 36.32 15360.80 0.54 colMeans(mtcars[,manifests]) mpg disp gear 20.09 230.72 3.69 summary(m1)$parameters name matrix row col Estimate Std.Error lbound ubound 1 ind.S[1,1] S mpg mpg 6.955250e+05 NA 2 ind.S[2,2] S disp disp 1.089073e+08 NA 3 ind.S[3,3] S gear gear 1.879208e+04 NA 4 ind.M[1,1] M 1 mpg 2.001955e+01 145.28108 5 ind.M[1,2] M 1 disp 1.307563e+02 1833.42676 6 ind.M[1,3] M 1 gear 3.693204e+00 24.14331