vulnerability: Climatic vulnerability

In rinnan/CENFA: Climate and Ecological Niche Factor Analysis

Climatic vulnerability

Description

Calculates the climatic vulnerability of a species using a cnfa and departure object.

Usage

vulnerability(
  cnfa,
  dep,
  method = "geometric",
  w,
  parallel = FALSE,
  n = 1,
  filename = "",
  ...
)

## S4 method for signature 'cnfa,departure'
vulnerability(
  cnfa,
  dep,
  method = "geometric",
  w,
  parallel = FALSE,
  n = 1,
  filename = "",
  ...
)

Arguments

cnfa	Object of class cnfa
dep	Object of class departure
method	character. What type of mean should be used to combine sensitivity and exposure. Choices are "arithmetic" and "geometric"
w	numeric. Optional vector of length two specifying the relative weights of sensitivity and exposure. See Details
parallel	logical. If TRUE then multiple cores are utilized
n	numeric. Number of cores to use for calculation
filename	character. Output filename (optional)
...	Additional arguments for file writing as for writeRaster

Details

The values of the vulnerability raster are calculated by combining the sensitivity \sigma and the exposure \epsilon. If method = "arithmetic", they will be combined as

\nu = (w_1\sigma + w_2\epsilon) / (\sum_i w_i).

If method = "geometric", they will be combined as

\nu = \sqrt(\sigma * \epsilon).

Value

Returns an S4 object of class vulnerability with the following slots:

call

Original function call

vf

Vulnerability factor. Vector of length p that describes the amount of vulnerability in each climate variable

vulnerability

Magnitude of the vulnerability factor

ras

RasterLayer of climate vulnerability

weights

Raster layer of weights used for departure calculation

References

Rinnan, D. Scott and Lawler, Joshua. Climate-niche factor analysis: a spatial approach to quantifying species vulnerability to climate change. Ecography (2019): <doi:10.1111/ecog.03937>.

See Also

departure

Examples

mod1 <- cnfa(x = climdat.hist, s.dat = ABPR, field = "CODE")
dep <- departure(x = climdat.hist, y = climdat.fut, s.dat = ABPR)
vuln <- vulnerability(cnfa = mod1, dep = dep)

rinnan/CENFA documentation built on July 19, 2023, 12:58 p.m.