bootstrap_E: Bootstrap thermoregulation effectiveness

Description

Bootstrapping of the effectiveness of temperature regulation (E) from the original distributions of Te and Tb as described by Hertz et al. (1993). One can choose the number of resamples and has the option to calculate E as defined by Hertz et al. (1993) or Blouin-Demers & Weatherhead (2001). See calculate_E_hertz and calculate_E_blouin for more information about these two indices. The thermal quality of the habitat (de) and accuracy of temperature regulation (db) are calculated as part of this formula, so it is not necessary to run calculate_de and calculate_db before running this function.

Usage

bootstrap_E(te, tb, tset_low, tset_up, index, n)

Arguments

A vector containing operative temperatures.

A vector containing body temperature measurements.

tset_low

Lower boundary of a species or population set-point range that was determined through thermal preference trials in a temperature gradient. This may be a named double vector containing the lower boundary value, or simply the value itself.

tset_up

Upper boundary of the set-point range.

index

Either 'hertz' or 'blouin'.

The desired number of samples drawn with replacement.

Value

The mean E and its 95 percent confidence interval obtained through resampling with replacement n times.

References

Blouin-Demers, G., & Weatherhead, P. J. (2001). Thermal ecology of black rat snakes (Elaphe obsoleta) in a thermally challenging environment. Ecology, 82 (11), 3025-3043. Hertz, P. E., Huey, R. B., & Stevenson, R. D. (1993). Evaluating temperature regulation by field-active ectotherms: the fallacy of the inappropriate question. The American Naturalist, 142(5), 796-818.

Examples

Run this code

# NOT RUN {
te <- na.omit(bufbuf[,"te"])
tb <- na.omit(bufbuf[,"tb"])
E_bootstrapped <- bootstrap_E(te, tb,
                             19.35, 26.44,
                             'hertz',
                             1000)

# }