glm.fit() gave me the same AIC's regardless of TRUE or FALSE intercept option.
> myX <- as.matrix(1:10) > myY <- 3+5*myX > foo <- glm.fit(x=myX, y=myY, family = gaussian(link = "identity"), > intercept=TRUE) > foo$aic [1] 38.94657 > foo <- glm.fit(x=myX, y=myY, family = gaussian(link = "identity"), > intercept=FALSE) > foo$aic [1] 38.94657 > AIC(lm(myY~0+myX, data=data.frame(myY,myX))) [1] 38.94657 > AIC(lm(myY~1+myX, data=data.frame(myY,myX))) [1] -650.9808 -Luke ______________________________________________ R-devel@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-devel
