-
Notifications
You must be signed in to change notification settings - Fork 18
/
Copy pathwqg.R
68 lines (66 loc) · 2.47 KB
/
wqg.R
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
# Copyright 2015 Province of British Columbia
#
# Licensed under the Apache License, Version 2.0 (the "License");
# you may not use this file except in compliance with the License.
# You may obtain a copy of the License at
#
# https://www.apache.org/licenses/LICENSE-2.0
#
# Unless required by applicable law or agreed to in writing, software
# distributed under the License is distributed on an "AS IS" BASIS,
# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
# See the License for the specific language governing permissions and
# limitations under the License.
#' Water Quality Guideline for British Columbia
#'
#' Calculates the 5% Hazard Concentration for British Columbia
#' after rescaling the data
#' based on the log-logistic, log-normal and gamma distributions
#' using the parametric bootstrap and AICc model averaging.
#'
#' Returns a tibble the model averaged 5% hazard concentration with
#' standard errors, 95% lower and upper confidence limits
#' and the number of bootstrap samples as well as the proportion of bootstrap
#' samples that successfully returned a likelihood
#' (convergence of the bootstrap sample is not required).
#'
#' @inheritParams params
#'
#' @return A tibble of the 5% hazard concentration with 95% confidence intervals.
#' @seealso [`ssd_fit_dists()`] and [`ssd_hc()`]
#' @export
#'
#' @examples
#' \dontrun{
#' ssd_wqg_bc(ssdtools::boron_data)
#' }
ssd_wqg_bc <- function(data, left = "Conc") {
dists <- ssd_dists("bc")
fits <- ssd_fit_dists(data, left = left, rescale = TRUE)
ssd_hc(fits, ci = TRUE, delta = 7, nboot = 10000)
}
#' Water Quality Guideline for Burrlioz
#'
#' Calculates the 5% Hazard Concentration
#' based on either the BurrIII3, Inverse Pareto or Log-Gumbel distribution.
#'
#' Returns a tibble the model averaged 5% hazard concentration with
#' standard errors, 95% lower and upper confidence limits
#' and the number of bootstrap samples as well as the proportion of bootstrap
#' samples that successfully returned a likelihood
#' (convergence of the bootstrap sample is not required).
#'
#' @inheritParams params
#'
#' @return A tibble of the 5% hazard concentration with 95% confidence intervals.
#' @seealso [`ssd_fit_dists()`] and [`ssd_hc()`]
#' @export
#'
#' @examples
#' \dontrun{
#' ssd_wqg_burrlioz(ssdtools::boron_data)
#' }
ssd_wqg_burrlioz <- function(data, left = "Conc") {
fit <- ssd_fit_burrlioz(data, left = left)
ssd_hc(fit, ci = TRUE, average = FALSE, nboot = 10000)
}