Attachment | Size |
---|---|
ACE_models.R [6] | 4.41 KB |
Hi OpenMx Community!
I just ran a single-trait ACE model, with one observed variable, and no definition variables. I used mxRefModels to get the Saturated and the Independence models. I wanted to see if I get similar values defining the Saturated model and the Independence model by hand.
My results differ basically in every aspect. Is my Saturated and Independence model definition wrong (see attached script), or is this a limitation of mxRefModels:
-Df is 6 vs 3: I guess mxRefModels includes A,C,E parameters, that is why it is 3 more. Should they be included in the Saturated model?
-My independence and saturated fit values is different, thus all fit statistics are different
#results of : summary(ACEFit, refModels = Sat_models) observed statistics: 198 estimated parameters: 4 degrees of freedom: 194 fit value ( units ): 736.6658 saturated fit value ( units ): 729.8184 number of observations: 99 chi-square: X2 ( df=6 ) = 6.847411, p = 0.3351895 Information Criteria: | df Penalty | Parameters Penalty | Sample-Size Adjusted AIC: 348.6658 744.6658 NA BIC: -154.7874 755.0463 742.4141 CFI: 0.9926078 TLI: 0.9975359 (also known as NNFI) RMSEA: 0.03777058 [95% CI (0, 0.1546974)] Prob(RMSEA <= 0.05): 0.4868861
#results of : summary(ACEFit, SaturatedLikelihood=SatSum$Minus2LogLikelihood , #SaturatedDoF=SatSum$degreesOfFreedom, #IndependenceLikelihood=IndepSum$Minus2LogLikelihood , #IndependenceDoF=IndepSum$degreesOfFreedom) observed statistics: 198 estimated parameters: 4 degrees of freedom: 194 fit value ( units ): 736.6658 saturated fit value ( units ): 731.998 number of observations: 99 chi-square: X2 ( df=3 ) = 4.667762, p = 0.1978056 Information Criteria: | df Penalty | Parameters Penalty | Sample-Size Adjusted AIC: 348.6658 744.6658 NA BIC: -154.7874 755.0463 742.4141 CFI: 0.9853227 TLI: 0.9902151 (also known as NNFI) RMSEA: 0.07493571 [95% CI (0, 0.2185226)] Prob(RMSEA <= 0.05): 0.2926299