R/el2.mean.R

In markean/bayesEL: Bayesian Empicial Likelihood

el2.mean <- function(theta, x, abstol = 1e-8, maxit = 50) {
  result <- el.mean(theta, x, abstol, maxit)
  result$logLR <- ifelse(all.equal(sum(result$w), 1) == T,
                         result$logLR,
                         -.Machine$double.xmax)
  result$logL <- ifelse(all.equal(sum(result$w), 1) == T,
                        result$logL,
                        -.Machine$double.xmax)
  result
}