Skip to content

Calculating AICc using conditional logistic regression

2 messages · K.Boughey at uea.ac.uk, Charles C. Berry

#
I have a case-control study that I'm analysing using the conditional
logistic regression function clogit from the survival package.

I would like to calculate the AICc of the models I fit using clogit.

I have a variety of scripts that can calculate AICc for models with a
logLik method, but clogit does not appear to use this method.

Is there a way I can calculate AICc from clogit in R?

Many thanks,
Katherine Boughey
#
On Tue, 12 Dec 2006, K.Boughey at uea.ac.uk wrote:

            
Review

 	?logLik
 	?coxph.object

and then roll your own logLik.coxph.

Here is one version:

logLik.coxph <-
     function(object,...) {
         y <-  -1 * (object$loglik[1] - object$loglik[2])
         class(y) <- "logLik"
         attr(y,'df')<- sum(!is.na(object$coef))
         y
}

Try it like this:
infert)
'log Lik.' 26.57712 (df=2)
[...]

Charles C. Berry                        (858) 534-2098
                                          Dept of Family/Preventive Medicine
E mailto:cberry at tajo.ucsd.edu	         UC San Diego
http://biostat.ucsd.edu/~cberry/         La Jolla, San Diego 92093-0717