Displaying 1 result from an estimated 1 matches for "lrmaic".
Did you mean:
limbic
2005 Aug 12
0
HowTo derive a correct likelihood-ratio chi-squared statistic from lrm() with a rsc() ?
...dAIC <- function(fit){
logl <- oos.loglik(fit) # derive -2*logL of the model
k <- 2
edf <- fit$stats[4]
# 'edf' is the equivalent degrees of freedom (i.e., the number of
parameters for usul parametric models) of 'fit'.
lrmAIC <- logl + k*edf
return(lrmAIC)
}
Best regards,
Jan
_______________________________________________________________________
Ir. Jan Verbesselt
Research Associate
Group of Geomatics Engineering
Department Biosystems ~ M³-BIORES
Vital Decosterstraat 102, 3000 Leuven...