bcgov/ssdtools

Use reciprocal of current burrIII distribution scale parameter so that increases with higher concentrations

Closed this issue · 1 comments

library(ssdtools)
#> Please replace the following in your scripts:
#> - `ssdtools::boron_data` with `ssddata::ccme_boron`
#> - `ssdtools::ccme_data` with `ssddata::ccme_data`

data <- ssddata::ccme_chloride

ssd_fit_dists(data, dists = "burrIII3")
#> Distribution 'burrIII3'
#>   scale 0.000707885
#>   shape1 0.585703
#>   shape2 1.7658
#> 
#> Parameters estimated from 28 rows of data.

data$Conc <- data$Conc / 1000

ssd_fit_dists(data, dists = "burrIII3", rescale = TRUE)
#> Distribution 'burrIII3'
#>   scale 4.83061
#>   shape1 0.585696
#>   shape2 1.7658
#> 
#> Parameters estimated from 28 rows of rescaled (6.824) data.

Created on 2021-10-29 by the reprex package (v2.0.1)

library(ssdtools)
#> Please replace the following in your scripts:
#> - `ssdtools::boron_data` with `ssddata::ccme_boron`
#> - `ssdtools::ccme_data` with `ssddata::ccme_data`

data <- ssddata::ccme_chloride

ssd_fit_dists(data, dists = "burrIII3")
#> Distribution 'burrIII3'
#>   scale 1412.85
#>   shape1 0.585714
#>   shape2 1.76589
#> 
#> Parameters estimated from 28 rows of data.


data$Conc <- data$Conc / 1000

ssd_fit_dists(data, dists = "burrIII3", rescale = TRUE)
#> Distribution 'burrIII3'
#>   scale 0.207014
#>   shape1 0.585694
#>   shape2 1.7658
#> 
#> Parameters estimated from 28 rows of rescaled (6.824) data.

Created on 2021-10-31 by the reprex package (v2.0.1)