R/aic_logregr.R

Defines functions AIC.logregr.madlib.grps logLik.logregr.madlib.grps extractAIC.logregr.madlib.grps logLik.logregr.madlib extractAIC.logregr.madlib

Documented in AIC.logregr.madlib.grps extractAIC.logregr.madlib extractAIC.logregr.madlib.grps logLik.logregr.madlib logLik.logregr.madlib.grps

extractAIC.glm.madlib <- extractAIC.logregr.madlib <- function(fit, scale=0, k=2, ...)
{
    edf <- length(fit$coef)
    c(edf, -2 * fit$log_likelihood + k * edf)
}

## ----------------------------------------------------------------------

logLik.glm.madlib <- logLik.logregr.madlib <- function(object, ...)
{
    ll <- object$log_likelihood
    attr(ll, "df") <- length(object$coef)
    class(ll) <- "logLik"
    ll
}

## ----------------------------------------------------------------------

extractAIC.glm.madlib.grps <- extractAIC.logregr.madlib.grps <- function(fit, scale=0, k=2, ...)
    lapply(fit, extractAIC, scale=scale, k=k, ...)

## ----------------------------------------------------------------------

logLik.glm.madlib.grps <- logLik.logregr.madlib.grps <- function(object, ...)
    lapply(object, logLik, ...)

## ----------------------------------------------------------------------

AIC.glm.madlib.grps <- AIC.logregr.madlib.grps <- function(object, ..., k=2)
    sapply(object, AIC, k=k, ...)

Try the PivotalR package in your browser

Any scripts or data that you put into this service are public.

PivotalR documentation built on March 13, 2021, 1:06 a.m.