this function s returning the compliment of the hp value ie. It’s giving the fraction affected rather than protected.
ssd_hp() should be a decreasing function of increasing concentration.
Eg:
> m<-ssd_fit_bcanz(aims_aluminium_marine)
> hc<-ssd_hc(m,c(1,5,10,20))$est # HCx estimates
Nominal fraction protected at HC1,HC5,HC10, and HC20
> round(ssd_hp(m,hc)$est,2)
[1] 1.61 5.13 10.00 20.09 # aside: why aren’t HC1 and HC5 fractions closer to nominal 1% and 5% ??
fraction protected at 2 x HCx
> round(ssd_hp(m,hc*2)$est,2)
[1] 2.90 9.54 16.98 30.69
Simplest solution is to replace “protected” to “affected” in Description above.
joethorley
commented
2 years ago
this function s returning the compliment of the hp value ie. It’s giving the fraction affected rather than protected.