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