aboutsummaryrefslogtreecommitdiff
path: root/R
diff options
context:
space:
mode:
Diffstat (limited to 'R')
-rw-r--r--R/llhist.R20
1 files changed, 9 insertions, 11 deletions
diff --git a/R/llhist.R b/R/llhist.R
index 9ddf5b10..22e3aa08 100644
--- a/R/llhist.R
+++ b/R/llhist.R
@@ -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")
}

Contact - Imprint