The AICs do not seem right to me either. Unless I am missing something, it appears that the formula:
AIC= -2x logLik -2k is being applied, rather than: AIC= -2x logLik +2k Meaning models with fewer degrees of freedom are being penalised. So in your example I make the degrees of freedom 9.61 +2*2= 13.61 for Integrated logLik 23.78 +2*13.7 for Penalised logLik. Someone please let me know if I am mistaken here! -- View this message in context: http://r.789695.n4.nabble.com/coxme-AIC-score-and-p-value-mismatch-tp2333980p2550348.html Sent from the R help mailing list archive at Nabble.com. ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.