Re: [R] Getting AIC from lrm in Design package

Dear David, Thank you for the reference to Frank Harrell's excellent text. I will read up to correct my statistical deficiencies offline. Thank you. On Sun, Oct 25, 2009 at 1:24 PM, David Winsemius wrote: > > On Oct 25, 2009, at 12:55 PM, Kyle Werner wrote: > >> David, >> >> Thank you for your

Re: [R] Getting AIC from lrm in Design package

On Oct 25, 2009, at 12:55 PM, Kyle Werner wrote: David, Thank you for your reply. I am not using glm, but instead lrm. Does not matter. "lrm" is giving you the same output as would glm with a logistic link. I am consulting the documentation to try to parse out what the output "Model L.R

Re: [R] Getting AIC from lrm in Design package

David, Thank you for your reply. I am not using glm, but instead lrm. I am consulting the documentation to try to parse out what the output "Model L.R." actually means: http://lib.stat.cmu.edu/S/Harrell/help/Design/html/lrm.fit.html ("model likelihood ratio chi-square") >From my read of the docum

Re: [R] Getting AIC from lrm in Design package

On Oct 25, 2009, at 9:24 AM, Kyle Werner wrote: I am trying to obtain the AICc after performing logistic regression using the Design package. For simplicity, I'll talk about the AIC. I tried building a model with lrm, and then calculating the AIC as follows: likelihood.ratio <- unname(lrm(succ

[R] Getting AIC from lrm in Design package

I am trying to obtain the AICc after performing logistic regression using the Design package. For simplicity, I'll talk about the AIC. I tried building a model with lrm, and then calculating the AIC as follows: likelihood.ratio <- unname(lrm(succeeded~var1+var2,data=scenario,x=T,y=T)$stats["Model