R/gmdl.dof.R

In plsRglm: Partial Least Squares Regression for Generalized Linear Models

#' @rdname aic.dof
#' @aliases gmdl.dof
#' @export gmdl.dof
gmdl.dof <- function (sigmahat, n, DoF, yhat) 
{
    SS <- sigmahat^2
    denominator <- DoF * SS
    FF <- (yhat)/(DoF * SS)
    FF[FF == 0] = Inf
    gmdl_temp <- (n/2) * log(SS) + (DoF/2) * log(FF) + (1/2) * log(n)
    return(gmdl_temp)
}