Nothing
#' @export
hist.evmOpt <-
function(x, xlab, ylab, main, ...){
# Want parameters as a matrix with one row for passing
# through to family$rng etc.
co <- t(x$coefficients)
th <- x$threshold
if (!is.finite(th))
th <- min(x$data$y)
UpperEndPoint <- x$family$endpoint(co, x)
dat <- x$data$y
dfun <- x$family$density
h <- hist(dat, plot = FALSE)
xx <- seq(th, min(UpperEndPoint, max(h$breaks)), length = 100)
y <- dfun(xx + .Machine$double.eps, co, x)
breaks <- seq(from=min(dat), to=max(dat), len=nclass.Sturges(dat) + 1)
res <- list(dat=dat, dens=cbind(x=xx, y=y), breaks=breaks)
oldClass(res) <- "hist.evmOpt"
res
}
#' @export
plot.hist.evmOpt <- function(x, xlab=NULL, ylab=NULL, main=NULL, ...){
if (missing(xlab) || is.null(xlab)) xlab <- "Data"
if (missing(ylab) || is.null(ylab)) ylab <- ""
if (missing(main) || is.null(main)) main <- "Histogram and density"
hist(x$dat, prob = TRUE, ylim = c(0, max(x$dens[, 2])),
xlab=xlab, ylab=ylab, main=main, breaks = x$breaks, ...)
lines(x$dens[, 1], x$dens[, 2], col = 4)
rug(x$dat)
invisible()
}
#' @export
print.hist.evmOpt <- function(x, xlab=NULL, ylab=NULL, main=NULL, ...){
plot.hist.evmOpt(x, xlab=xlab, ylab=ylab, main=main, ...)
invisible(x)
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.