This function calculates several metrics of niche overlap based on predictions for two species (or for the same species but different models) at the same sites.
nicheOverlapMetrics(
x1,
x2,
method = c("meanDiff", "meanAbsDiff", "rmsd", "d", "i", "esp", "cor", "rankCor"),
w = rep(1, length(x1)),
na.rm = FALSE,
...
)List object with one element per value specified by the argument in method.
Numeric. Vector of predictions from a model.
Numeric. Vector of predictions from another model.
Character vector, indicates type of metric to calculate:
meanDiff: Average difference
meanAbsDiff: Average of absolute values of difference
rmsd: Root-mean square deviation
d: Schoener's D
i: Warren's I
esp: Godsoe's ESP
cor: Pearson correlation between x1 and x2 (will apply logitAdj() first unless logit=FALSE).
rankCor: Spearman rank correlation.
Numeric vector. Weights of predictions in x1 and x2.
Logical. If TTRUE then remove elements in x1 and 2 that are NA in either x1 or x2.
Other arguments (not used).
Warren, D.L., Glor, R.E., and Turelli, M. 2008. Environmental niche equivalency versus conservatism: Quantitative approaches to niche evolution. Evolution 62:2868-2883. tools:::Rd_expr_doi("10.1111/j.1558-5646.2008.00482.x")
Warren, D.L., Glor, R.E., and Turelli, M. 2008. Erratum. Evolution 62:2868-2883. tools:::Rd_expr_doi("10.1111/j.1558-5646.2010.01204.x")
Godsoe, W. 2014. Inferring the similarity of species distributions using Species' Distribution Models. Ecography 37:130-136. tools:::Rd_expr_doi("10.1111/j.1600-0587.2013.00403.x")
compareResponse
x1 <- seq(0, 1, length.out=100)
x2 <- x1^2
nicheOverlapMetrics(x1, x2)
Run the code above in your browser using DataLab