← Back to Funcspace Scripts

density_threshold(z, p)

Find the density value at a given cumulative probability threshold from a KDE grid.

functional diversitystatistics
Args:z — KDE density matrixp — target probability (e.g. 0.95)
density_threshold <- function(z, p) {
  dx <- sort(as.vector(z), decreasing = TRUE)
  cs <- cumsum(dx) / sum(dx)
  dx[which(cs >= p)[1]]
}