This function computes overlap metrics between two rasters. It currently
implements Schoener's D and the inverse I of Hellinger's distance.
Usage
niche_overlap(x, y, method = c("Schoener", "Hellinger"))
Arguments
- x
a terra::SpatRaster with a single layer
- y
a terra::SpatRaster with a single layer
- method
a string (or vector of strings) taking values "Schoener" and
"Hellinger"
Value
a list of overlap metrics, with slots D and I (depending on method)
Details
Note that Hellinger's distance is
normalised by dividing by square root of 2 (which is the correct
asymptote for Hellinger's D),
rather than the incorrect 2 used originally in
Warren et al (2008), based on the Erratum for that paper.
References
Warren, D.L., Glor, R.E. & Turelli M. (2008) Environmental niche
equivalency versus conservativism: quantitative approaches to niche evolution.
Evolution 62: 2868-2883
