res_scale {ecostatscale}

Resilience scaling function

Description

Extrapolate resilience observed at the scale of a single spatial or ecological scale b (e.g. a patch or species) to a larger scale, B (e.g. functional group or landscape).

Usage

res_scale(
  mvar_b,
  murho_b_abundance,
  mucov_b_abundance = NULL,
  msd_b,
  murho_b_disturbance,
  mucov_b_disturbance = NULL,
  b = 1,
  B,
  lambda
)

Arguments

mvar_b

Mean abundance variance observed at scale b

murho_b_abundance

Mean Pearson correlation coefficient of abundance values observed at scale b

mucov_b_abundance

Mean covariance of abundance values observed at scale b. Ignored unless murho_b_abundance is NULL Defaults to NULL.

msd_b

Mean disturbance standard deviation observed at scale b

murho_b_disturbance

Mean Pearson correlation coefficient of disturbances observed at scale b

mucov_b_disturbance

Mean covariance of disturbances observed at scale b. Ignored unless murho_b_abundance is NULL Defaults to NULL.

b

Size of observed scale. Defaults to 1.

B

Larger scale being extrapolated to (e.g. total number of species, or size of patch B relative to b)

lambda

Mean disturbance frequency.

Value

Extrapolated median resilience at scale of M species.

Examples

# extrapolate from scale of 1 species to 10 species
res_scale(mvar_b = 0.25, murho_b_abundance = -0.034, msd_b = sqrt(0.1),
           murho_b_disturbance = 0, B = 30, lambda=1)

# plot relationship for groups of 1 to 30 species
plot(1:30, res_scale(mvar_b = 0.25, murho_b_abundance = -0.034, msd_b = sqrt(0.1),
      murho_b_disturbance = 0, B = 1:30, lambda=1),
      xlab="ecological scale", ylab="resilience, r", type="b")