R/logLik_glm.R

Defines functions logLik_glm.fit

Documented in logLik_glm.fit

#' Calculating log-likelihood in generalized linear models
#' 
#' An utility function to extract log-likelihood based on AIC of glm.fit()
#' @param ... Arguments to be passed to \code{glm.fit()}.
#' @return Log-likelihood
logLik_glm.fit <- function(...) {
    fit <- stats::glm.fit(...)
    p <- fit$rank
    if (fit$family$family %in% c("gaussian", "Gamma", "inverse.gaussian"))
        p <- p + 1
    p - fit$aic / 2
}

Try the rNeighborQTL package in your browser

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

rNeighborQTL documentation built on April 9, 2021, 9:06 a.m.