summary() needs ability to provide the df of the saturated model, as well as the saturated fit?

tbates's picture
Project:OpenMx
Component:Code
Category:feature request
Priority:normal
Assigned:Unassigned
Status:active
Description

Hi,
was just using this to compare a model to the saturated model:

> summary(twinSatSub1Fit, SaturatedLikelihood=1599.781)

It gets the χ² but not the df between the sub-model and the saturated model.

Do we need a parameter for summary to specifiy "sat model estimated parameters"?

ideally, one could simply do

> summary(subFit, satFit)

And summary would summarise the second model recursively then grab $Minus2LogLikelihood and $estimatedParameters

Cheers,
t