R/c_.R

c_ <- function (d, beta, y, X) {

  beta_1 <- log(mean(y * exp(-X %*% cbind(beta))))

  mean(log(y + d * exp(beta_1 + X %*% cbind(beta))) - beta_1 - X %*% cbind(beta))
}

Try the IOLS package in your browser

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

IOLS documentation built on April 8, 2023, 1:15 a.m.