#' Calculate area-adjusted error matrix from model prediction
#'
#' @param map Predicted map
#' @param pred Vector with predicted values.
#' @param ref Vector with reference values.
#' @param labels Labels for map classes.
#' @param pixelsize Numeric. Size of pixels.
#' @param ci Numeric. Confidence interval to use.
#'
#' @return Area adjusted error matrix
#' @export
area_adjusted_error_matrix<- function(map, pred, ref, labels = NULL, pixelsize = .09, ci = .95) {
maparea <- as.numeric(table(raster::values(map)[which(!is.na(raster::values(map)))])) # in pixel counts
area <- maparea*pixelsize # in area unit
if(!is.null(labels)) {
names(area) <- labels
} else {
names(area) <- make.names(area)
}
### calculate error matrix
cm <- table(pred, ref)
dtwSat::twdtwAssess(cm, area = area)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.