vulnerability: Climatic vulnerability

Description Usage Arguments Details Value References See Also Examples

Description

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

Usage

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
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 σ and the exposure ε. If method = "arithmetic", they will be combined as

ν = (w_1σ + w_2ε) / (∑_i w_i).

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

ν = √(σ * ε).

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

1
2
3
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)

CENFA documentation built on Aug. 16, 2021, 9:06 a.m.