R/utilities.R

In evwalz/uroc: Computes ROCM, UROC and CPA

# compute area under curve using trapezoidal rule

trap <- function(farate, hitrate) {

  diffs.far <- farate[-1] - farate[-length(farate)]
  means.vert <- (hitrate[-1] + hitrate[-length(hitrate)]) * 0.5

  return(sum(means.vert * diffs.far))
}