Gets concentration(s) that protect specified percentage(s) of species.

ssd_hc(x, ...)

# S3 method for list
ssd_hc(x, percent = 5, hc = 5, ...)

# S3 method for fitdists
ssd_hc(
  x,
  percent = 5,
  hc = 5,
  ci = FALSE,
  level = 0.95,
  nboot = 1000,
  average = TRUE,
  delta = 7,
  min_pboot = 0.99,
  parametric = TRUE,
  control = NULL,
  ...
)

# S3 method for fitburrlioz
ssd_hc(
  x,
  percent = 5,
  ci = FALSE,
  level = 0.95,
  nboot = 1000,
  min_pboot = 0.99,
  parametric = FALSE,
  ...
)

Arguments

Value

A tibble of corresponding hazard concentrations.

Details

If ci = TRUE uses parameteric bootstrapping to get confidence intervals on the hazard concentrations(s).

Methods (by class)

  • list: Hazard Concentrations for Distributional Estimates

  • fitdists: Hazard Concentrations for fitdists Object

  • fitburrlioz: Hazard Concentrations for fitburrlioz Object '

Examples

fits <- ssd_fit_dists(ssddata::ccme_boron)
ssd_hc(fits)
#> # A tibble: 1 × 10
#>   dist    percent   est    se   lcl   ucl    wt method     nboot pboot
#>   <chr>     <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <chr>      <int> <dbl>
#> 1 average       5  1.24    NA    NA    NA     1 parametric     0    NA
ssd_hc(estimates(fits))
#> # A tibble: 6 × 9
#>   dist        percent   est    se   lcl   ucl    wt nboot pboot
#>   <chr>         <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <int> <dbl>
#> 1 gamma             5  1.07    NA    NA    NA     1     0    NA
#> 2 lgumbel           5  1.77    NA    NA    NA     1     0    NA
#> 3 llogis            5  1.56    NA    NA    NA     1     0    NA
#> 4 lnorm             5  1.68    NA    NA    NA     1     0    NA
#> 5 lnorm_lnorm       5  1.54    NA    NA    NA     1     0    NA
#> 6 weibull           5  1.09    NA    NA    NA     1     0    NA
ssd_hc(ssd_match_moments())
#> # A tibble: 6 × 9
#>   dist        percent   est    se   lcl   ucl    wt nboot pboot
#>   <chr>         <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <int> <dbl>
#> 1 gamma             5 0.481    NA    NA    NA     1     0    NA
#> 2 lgumbel           5 0.739    NA    NA    NA     1     0    NA
#> 3 llogis            5 0.562    NA    NA    NA     1     0    NA
#> 4 lnorm             5 0.558    NA    NA    NA     1     0    NA
#> 5 lnorm_lnorm       5 0.469    NA    NA    NA     1     0    NA
#> 6 weibull           5 0.501    NA    NA    NA     1     0    NA
fit <- ssd_fit_burrlioz(ssddata::ccme_boron)
ssd_hc(fit)
#> # A tibble: 1 × 10
#>   dist      percent   est    se   lcl   ucl    wt method         nboot pboot
#>   <chr>       <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <chr>          <int> <dbl>
#> 1 invpareto       5 0.387    NA    NA    NA     1 non-parametric     0    NA