Previous measures of overlap between suitability rasters that were implemented in ENMTools were based on measures of similarity between probability distributions. That's fine as far as it goes, but my feeling is that it's more useful as a measure of species' potential spatial distributions than as a measure of similarity of the underlying model. Here's a quick example:
# Perfect positive correlation
sp1 = seq(0.1, 1.0, 0.001)
sp2 = seq(0.1, 1.0, 0.001)
Okay, that's all well and good - perfectly positively correlated responses get a 1 from all metrics. Now what if they're perfectly negatively correlated?
# Perfect negative correlation
sp1 = seq(0.1, 1.0, 0.001)
sp2 = seq(1.0, 0.1, 0.001)
What's going on here? Spearman rank correlation tells us that they are indeed negatively correlated, but D and I both have somewhat high values! The reason is that the values of the two functions are fairly similar across a fairly broad range of environments, even though the functions themselves are literally as different as they could possibly be. Thinking about what this means in terms of species occurrence is quite informative; if the threshold for suitability for a species to occur is low (e.g., .0005 in this cartoon example), they might co-occur across a fairly broad range of environments; both species would find env values from 250 to 750 suitable and might therefore overlap across about 2/3 of their respective ranges. That's despite them having completely opposite responses to that environmental gradient, strange though that may seem.
So do you want to measure the potential for your species to occupy the same environments, or do you want to measure the similarity in their estimated responses to those environments? That's entirely down to what question you're trying to answer!
Okay, one more reason I kinda like correlations:
sp1 = abs(rnorm(1000))
sp2 = abs(rnorm(1000))
Again we've got fairly high overlaps between species using D and I, but Spearman rank correlation is really close to zero. That's exactly what we'd expect if there's no signal at all. Of course the fact that species distributions and the environment are both spatially autocorrelated means that we'll likely have higher-than-zero (at least in absolute value) correlations even if there is no causal relationship between the environment and species distributions, but at least it's nice to know that we do have a clear expected value when chaos reigns.
Code for this is here:
Interesting, thanks! I came upon this today as I was thinking about some similar ideas. I wrote up a few notes, FYI: https://plantarum.ca/2021/12/02/schoenersd/ReplyDelete