p-values from summary and mxCompare do not match

Posted on Fri, 11/18/2011 - 04:07

JWiley Joined: 03/25/2011

Attachment

Size

pvalExample.R

1.82 KB

Forums

OpenMx General Help

Specifically, I fit a model as well as a saturated model. When calling summary, I pass the -2LL and DF from the saturated model, with mxCompare, I compare my model with the saturated. Here are the results, everything else is as expected. Am I missing something?
from summary():
chi-square: 5.35
p: 1

from mxCompare():
minus2LL df AIC diffLL diffdf p
1 5318 742 3834 NA NA NA
2 5323 744 3835 5.35 2 0.0691

A reproducible example is provided in the attached script. I am using r57586 on Windows x64 with OpenMx_999.0.0-1661.

Replied on Thu, 12/01/2011 - 13:42

mspiegel Joined: 07/31/2009

mxCompare() invokes the

mxCompare() invokes the summary() function on the 'base' argument, and then invokes the summary() function on the 'comparison' argument, but in either case it does not pass any of the optional arguments (SaturatedLikelihood, SaturatedDoF, IndependenceLikelihood, IndependenceDoF) to the summary() function.

Whether or not mxCompare() should be passing optional arguments to the summary() function is a question for the other people on the development team.

Log in or register to post comments

Replied on Thu, 12/01/2011 - 13:55

tbates Joined: 07/31/2009

support mxCompare(numbs, numStats, SaturatedLikelihood, indep)

Makes sense to me that numbs, numStats, SaturatedLikelihood, and indep should be passed on to the omxLapply calls to summary if they are provided to mxCompare()