[R] Degrees of freedom for lm in logLik and AIC

Frank Harrell f.harrell at vanderbilt.edu
Mon Mar 28 05:36:23 CEST 2011


I have a question about the computation of the degrees of freedom in a linear
model:

x <- runif(20); y <- runif(20)
f <- lm(y ~ x)
logLik(f)
'log Lik.' -1.968056 (df=3)

The 3 is coming from f$rank + 1.  Shouldn't it be f$rank?  This affects
AIC(f).
Thanks
Frank

-----
Frank Harrell
Department of Biostatistics, Vanderbilt University
--
View this message in context: http://r.789695.n4.nabble.com/Degrees-of-freedom-for-lm-in-logLik-and-AIC-tp3410687p3410687.html
Sent from the R help mailing list archive at Nabble.com.



More information about the R-help mailing list