diff options
Diffstat (limited to 'R/llhist.R')
| -rw-r--r-- | R/llhist.R | 20 |
1 files changed, 9 insertions, 11 deletions
@@ -1,8 +1,7 @@ #' Plot the distribution of log likelihoods from multistart objects #' -#' Produces a histogram of log-likelihoods, and an overlayed kernel density -#' estimate. In addition, the likelihood of the original fit is shown as -#' a red vertical line. +#' Produces a histogram of log-likelihoods. In addition, the likelihood of the +#' original fit is shown as a red vertical line. #' #' @param object The [multistart] object #' @param breaks Passed to [hist] @@ -10,9 +9,10 @@ #' @param main Title of the plot #' @param \dots Passed to [hist] #' @seealso [multistart] -#' @importFrom KernSmooth bkde #' @export -llhist <- function(object, breaks = "Sturges", lpos = "topleft", main = "", ...) { +llhist <- function(object, breaks = "Sturges", lpos = "topleft", main = "", + ...) +{ oldpar <- par(no.readonly = TRUE) on.exit(par(oldpar, no.readonly = TRUE)) @@ -26,19 +26,17 @@ llhist <- function(object, breaks = "Sturges", lpos = "topleft", main = "", ...) } ll <- stats::na.omit(sapply(object, llfunc)) - kde <- KernSmooth::bkde(ll) par(las = 1) h <- hist(ll, freq = TRUE, - xlim = range(kde$x), xlab = "", main = main, ylab = "Frequency of log likelihoods", breaks = breaks, ...) freq_factor <- h$counts[1] / h$density[1] - lines(kde$x, freq_factor * kde$y) + abline(v = logLik(attr(object, "orig")), col = 2) + legend(lpos, inset = c(0.05, 0.05), bty = "n", - lty = 1, col = c(2, 1), - legend = c("original log likelihood", - "kernel density estimate")) + lty = 1, col = c(2), + legend = "original fit") } |