Nothing
resampleSummary <- function(obs, resampled, index = NULL, keepData = TRUE)
{
numPred <- apply(resampled, 2, function(u) sum(!is.na(u)))
# for everything but LOO, we should get multiple predictions per resample
if(all(numPred >= 2))
{
# calculate performance metrics for each resample
performanceStats <- apply(
resampled,
2,
postResample,
obs = obs)
#summarize resample dists
out <- c(
apply(performanceStats, 1, mean, na.rm = TRUE),
apply(performanceStats, 1, sd, na.rm = TRUE))
# optionally returen the data in "vertical" format
# to conserve space remove the many missing values
if(keepData)
{
# stack has issues when there are a lot of missing values,
# so we'll use lapply to stack the columns of resampled
if(is.factor(obs))
{
outResample <- data.frame(
obs = rep(obs, dim(resampled)[2]),
pred = factor(unlist(lapply(resampled, as.character)),
levels = levels(obs)),
group = paste(
"Resample",
rep(
1:dim(resampled)[2],
each = dim(resampled)[1], sep = "")))
} else {
outResample <- data.frame(
obs = rep(obs, dim(resampled)[2]),
pred = unlist(lapply(resampled, I)),
group = paste(
"Resample",
rep(
1:dim(resampled)[2],
each = dim(resampled)[1], sep = "")))
}
} else outResample <- NULL
} else {
pred <- apply(resampled, 2, function(u) u[!is.na(u)])
if(is.factor(obs)) pred <- factor(as.character(pred), levels = levels(obs))
tmp <- postResample(pred, obs)
tmp2 <- tmp * 0
out <- c(
tmp,
tmp * 0)
outResample <- data.frame(
obs = obs,
pred = pred,
group = "Resample1")
}
if(keepData) outResample <- outResample[!is.na(outResample$pred),]
list(metrics = out, data = outResample)
}
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.