R/error_matrix.R

In juoe/spatialtoolbox: Collection of functions for doing spatial stuff.

#' 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)


}