Learn R Programming

SSP (version 1.1.0)

summary_ssp: Summary of MultSE for Each Sampling Effort in Simulated Data Sets

Description

Computes the average MultSE (pseudo-multivariate standard error) for each sampling effort across simulated datasets, and estimates associated variation and rate of change.

Usage

summary_ssp(results, multi.site)

Value

A data frame summarizing MultSE for each sampling effort, including the mean, quantiles, relativized values, and estimated derivative.

Arguments

results

A matrix generated by sampsd containing MultSE values for each simulation and sampling configuration.

multi.site

Logical. Indicates whether multiple sites were simulated.

Details

For each sampling effort in each simulated data set, the average MultSE is computed (Anderson & Santana-Garcon, 2015). The function then calculates the overall mean and associated lower and upper quantiles of these averages. To evaluate how precision improves with effort, the average MultSE values are relativized to the maximum (typically at the lowest effort), and a numerical forward finite difference derivative is calculated to approximate the rate of change.

This output is used to support the identification of optimal and redundant sampling efforts based on precision gain.

References

Anderson, M. J., & Santana-Garcon, J. (2015). Measures of precision for dissimilarity-based multivariate analysis of ecological communities. Ecology Letters, 18(1), 66–73.

Guerra-Castro, E.J., Cajas, J.C., Simões, N., Cruz-Motta, J.J., & Mascaró, M. (2021). SSP: an R package to estimate sampling effort in studies of ecological communities. Ecography 44(4), 561-573. doi: tools:::Rd_expr_doi("10.1111/ecog.05284")

See Also

sampsd, ioptimum

Examples

Run this code
## Single site example
data(micromollusk)
par.mic <- assempar(data = micromollusk, type = "P/A", Sest.method = "average")
sim.mic <- simdata(par.mic, cases = 2, N = 10, sites = 1)
sam.mic <- sampsd(dat.sim = sim.mic, Par = par.mic, transformation = "P/A",
                  method = "jaccard", n = 10, m = 1, k = 3)
summ.mic <- summary_ssp(results = sam.mic, multi.site = FALSE)

## See the full multi-site workflow in the vignette

Run the code above in your browser using DataLab