R/logL.R

In npmr: Nuclear Penalized Multinomial Regression

Defines functions logL

Documented in logL

logL <-
function(B, b, X, Y) {
    eta = matrix(1, nrow(X), 1) %*% t(b) + X %*% B
    sum(Y*(eta - log(rowSums(exp(eta)))))
}

Try the npmr package in your browser