> ### To get around that, I computed 2*log(likelihood ratio) manually: > > lglk0 <- logLik(fit0) > lglk.ID1. <- logLik(Fit.ID1.) > pchisq(as.numeric(chisq.ID.), 1, lower=FALSE) > > [1] 0.008545848
(I think you're missing a line in there) But isn't this rather perilous unless you are confident that the two models are using exactly the same formulation of the likelihood? (ie. that they are truly nested) Hadley ______________________________________________ [email protected] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
