I am fitting logistic regression models, by defining my own link function, and would like to get AICc values. Using the glm command gives a value for AIC, but I haven't been able to get R to convert that to AICc. Is there a code that has already been written for this? Right now I am just putting the AIC values into an excel spreadsheet and calculating AICc, likelihood, and AIC weights that way, but it would be much more efficient if I could do this in R. Thanks!
Katie Stumpf wrote:> > I am fitting logistic regression models, by defining my own link > function, and would like to get AICc values. Using the glm command > gives a value for AIC, but I haven't been able to get R to convert > that to AICc. Is there a code that has already been written for > this? Right now I am just putting the AIC values into an excel > spreadsheet and calculating AICc, likelihood, and AIC weights that > way, but it would be much more efficient if I could do this in R. > >RSiteSearch("AICc", restrict="functions") -- View this message in context: http://www.nabble.com/AICc-tp23303142p23304120.html Sent from the R help mailing list archive at Nabble.com.
Katie- I think you can calculate it easily in R using the definitions: AIC<- -2*LL + 2*K AICc<-AIC + 2*K*(K+1)/(n-K-1) where LL is the logLik, K is the the number of parameters in the fitted model, and n is the number of observation, so if you know AIC then it is easy to calculate AICc Ron Katie Stumpf wrote:> I am fitting logistic regression models, by defining my own link > function, and would like to get AICc values. Using the glm command > gives a value for AIC, but I haven't been able to get R to convert > that to AICc. Is there a code that has already been written for > this? Right now I am just putting the AIC values into an excel > spreadsheet and calculating AICc, likelihood, and AIC weights that > way, but it would be much more efficient if I could do this in R. > > Thanks! > > ______________________________________________ > R-help at 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. >-- R. R. Burns Oceanside, CA