R/Lhat_eta.r

In logcondens: Estimate a Log-Concave Probability Density from Iid Observations

Defines functions Lhat_eta

Documented in Lhat_eta

Lhat_eta <- function(x, w, eta) {
  n <- length(x)
  dx <- diff(x)
  phi <- etaphi(x, eta)
  ll <- sum(w * phi) - sum(dx * J00(phi[1:(n - 1)], phi[2:n]))
  return(list(ll = ll))
}

Try the logcondens package in your browser