Skip to contents

Gets hazard concentrations with confidence intervals that protect 1, 5, 10 and 20% of species using settings adopted by BC, Canada, Australia and New Zealand for official guidelines. This function can take several minutes to run with required 10,000 iterations.

Usage

ssd_hc_bcanz(x, nboot = 10000, delta = 10, min_pboot = 0.95)

Arguments

x

The object.

nboot

A count of the number of bootstrap samples to use to estimate the confidence limits. A value of 10,000 is recommended for official guidelines.

delta

A non-negative number specifying the maximum absolute AIC difference cutoff. Distributions with an absolute AIC difference greater than delta are excluded from the calculations.

min_pboot

A number between 0 and 1 of the minimum proportion of bootstrap samples that must successfully fit (return a likelihood) to report the confidence intervals.

Value

A tibble of corresponding hazard concentrations.

See also

ssd_hc().

Other BCANZ: ssd_fit_bcanz()

Examples

fits <- ssd_fit_bcanz(ssddata::ccme_boron)
ssd_hc_bcanz(fits, nboot = 100)
#> Warning: One or more pboot values less than 0.95 (decrease min_pboot with caution).
#> # A tibble: 4 × 11
#>   dist    proportion   est    se   lcl   ucl    wt method    nboot pboot samples
#>   <chr>        <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <chr>     <dbl> <dbl> <I<lis>
#> 1 average       0.01 0.267    NA    NA    NA     1 parametr…   100  0.85 <dbl>  
#> 2 average       0.05 1.26     NA    NA    NA     1 parametr…   100  0.85 <dbl>  
#> 3 average       0.1  2.38     NA    NA    NA     1 parametr…   100  0.85 <dbl>  
#> 4 average       0.2  4.81     NA    NA    NA     1 parametr…   100  0.85 <dbl>