R/el2.lm.R

In markean/bayesEL: Bayesian Empicial Likelihood

el2.lm <- function(b, x, y, abstol = 1e-8, maxit = 100) {
  result <- el.lm(b, x, y, 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
}