Nothing
predPlot <- function(model = NULL, obs = NULL, pred = NULL, thresh = "preval", main = "Classified predicted values", legend.pos = "n", pch = 1, cex = 0.5, col = c("black", "grey"), na.rm = TRUE, rm.dup = FALSE, interval = 0.01, quant = 0) {
# version 1.7 (26 Oct 2022)
obspred <- inputMunch(model, obs, pred, na.rm = na.rm, rm.dup = rm.dup)
if (!is.null(obs) || !is.null(model)) obs <- obspred[ , "obs"]
pred <- obspred[ , "pred"]
pred0 <- pred[obs == 0]
pred1 <- pred[obs == 1]
thresh_exists <- !is.na(thresh) && !is.null(thresh)
if (thresh_exists) {
if (!(is.numeric(thresh) || thresh %in% modEvAmethods("getThreshold")))
stop("'thresh' must be either NA, NULL, a numeric value, or one of the options obtained with modEvAmethods('getThreshold')")
if (thresh %in% modEvAmethods("getThreshold"))
# if (thresh == "meanPred") obs <- 0 # so it works even if 'obs' not provided (it's not used)
thresh <- getThreshold(obs = obs, pred = pred, threshMethod = thresh, interval = interval, quant = quant, na.rm = na.rm)
}
# opar <- par(no.readonly = TRUE)
# on.exit(par(opar))
par_mar <- par()$mar
on.exit(par(mar = par_mar))
par(mar = c(5, 5.2, 3, 1))
plot(x = c(0, 1), y = c(-0.5, 1.5), xlab = "Predicted value", type = "n", ylab = "", yaxt = "n", main = main)
axis(side = 2, at = c(0, 1), tick = FALSE, labels = c("Observed\nabsences", "Observed\npresences"), las = 1)
if (thresh_exists) {
abline(v = thresh, lty = 2)
points(x = pred0, y = sapply(rep(0, length(pred0)), jitter, 10), pch = pch, cex = cex, col = ifelse(pred0 < thresh, col[2], col[1]))
points(x = pred1, y = sapply(rep(1, length(pred1)), jitter, 10), pch = pch, cex = cex, col = ifelse(pred1 < thresh, col[2], col[1]))
if (!is.na(legend.pos) && legend.pos != "n") legend(legend.pos, legend = c("Predicted presence", "Predicted absence"), pch = pch, col = col)
} else { # end if thresh_exists
points(x = pred0, y = sapply(rep(0, length(pred0)), jitter, 10), pch = pch, cex = cex, col = col[1])
points(x = pred1, y = sapply(rep(1, length(pred1)), jitter, 10), pch = pch, cex = cex, col = col[1])
if (!is.na(legend.pos) && legend.pos != "n") legend(legend.pos, legend = "Prediction", pch = pch, col = col[1])
} # end if thresh NA or null
}
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.