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 recommended 10,000 iterations.

Usage

ssd_hc_bcanz(x, nboot = 10000, 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.

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

Examples

fits <- ssd_fit_bcanz(ssddata::ccme_boron)
ssd_hc_bcanz(fits, nboot = 100)
#> # 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 0.449 0.0654  1.73     1 paramet…   100     1 <dbl>  
#> 2 average       0.05 1.26  0.878 0.538   3.65     1 paramet…   100     1 <dbl>  
#> 3 average       0.1  2.38  1.28  1.25    5.74     1 paramet…   100     1 <dbl>  
#> 4 average       0.2  4.81  1.99  2.50    9.77     1 paramet…   100     1 <dbl>