ioptimum | R Documentation |
Estimates the sampling effort at which the improvement in precision (MultSE) per additional sampling unit becomes sub-optimal or redundant, based on predefined cut-off thresholds.
ioptimum(xx, multi.site = TRUE, c1 = 10, c2 = 5, c3 = 2.5)
xx |
A data frame generated by |
multi.site |
Logical. Indicates whether multiple sites were simulated. |
c1 |
First cut threshold. Default is 10% improvement over the highest MultSE. |
c2 |
Second cut threshold. Default is 5% improvement over the highest MultSE. |
c3 |
Third cut threshold. Default is 2.5% improvement over the highest MultSE. |
Sampling efforts between the minimum (e.g. 2 samples) and c1
represent the necessary effort to achieve acceptable precision.
Efforts between c1
and c2
reflect sub-optimal gains, and those between c2
and c3
are considered optimal.
Beyond c3
, any additional effort results in marginal improvements in MultSE and may be considered redundant.
This classification helps support cost-benefit decisions in ecological survey design (see Underwood, 1990).
If c3
is not reached within the simulated range, the maximum available effort is returned with a warning.
A vector or matrix indicating the sampling sizes corresponding to each cut-off point.
The cut-off thresholds are arbitrary and should be adjusted based on the ecological question and resource availability.
In some cases, c3
may not be reached within the range of simulated sampling efforts.
Underwood, A. J. (1990). Experiments in ecology and management: Their logics, functions and interpretations. Australian Journal of Ecology, 15, 365–389.
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: \Sexpr[results=rd]{tools:::Rd_expr_doi("10.1111/ecog.05284")}
plot_ssp
, summary_ssp
, sampsd
## Single site example
data(micromollusk)
par.mic <- assempar(data = micromollusk, type = "P/A", Sest.method = "average")
sim.mic <- simdata(par.mic, cases = 3, N = 20, 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)
opt.mic <- ioptimum(xx = summ.mic, multi.site = FALSE)
## Multiple sites example
data(sponges)
par.spo <- assempar(data = sponges, type = "counts", Sest.method = "average")
sim.spo <- simdata(par.spo, cases = 3, N = 10, sites = 3)
sam.spo <- sampsd(dat.sim = sim.spo,
Par = par.spo,
transformation = "square root",
method = "bray",
n = 10,
m = 3,
k = 3)
summ.spo <- summary_ssp(results = sam.spo, multi.site = TRUE)
opt.spo <- ioptimum(xx = summ.spo, multi.site = TRUE)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.