R/loglik.R

Defines functions loglik

loglik <- function(X, y, beta, family) {
    link = as.vector(X %*% beta)
    n = length(y) 
    if (family == "gaussian") 
      return(n * log(mean((y - link)^2)))
    if (family == "poisson") 
        return(-2 * sum(exp(link) + 2 * y * link))
    if (family == "binomial") 
        return(2 * sum(log(1 + exp(link)) - y * link))
       
}

Try the RAMP package in your browser

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

RAMP documentation built on Jan. 16, 2020, 5:02 p.m.