Hi,
I've tried setting up an ACDE model as I have siblings (up to 3) data for my twins.
Sorry, basic question but I've managed to confuse myself with how to interpret best fit. I've embedded all submodels within the ACDE model. As per below, the AE model has the lowest AIC, so that is the best fit right? So the p value is lower than the ACE model because of the extra df compared to the main model?
> # Compare all models
> ACDENested <- list(multiCholACEFit,multiCholADEFit, multiCholAEFit, multiCholCEFit, multiCholEFit)
> tableFitStatistics(multiCholACDEFit,ACDENested)
Name ep -2LL df AIC diffLL diffdf p
Model 1 : multiCholACDE 18 10103.08 4076 1951.08 - - -
Model 2 : multiCholACE 15 10103.8 4079 1945.8 0.72 3 0.87
Model 3 : multiCholADE 15 10106.97 4079 1948.97 3.89 3 0.27
Model 4 : multiCholAE 12 10106.97 4082 1942.97 3.89 6 0.69
Model 5 : multiCholCE 12 10165.75 4082 2001.75 62.68 6 0
Model 6 : multiCholE 9 10459.53 4085 2289.53 356.45 9 0