It is recommened to specify the arguments rate, water_depth and water_width using units::units from the units package.

Since pfm version 0.6.5, the function is vectorised with respect to rates, +applications, water depth, crop groups and distances

@@ -208,6 +214,43 @@ if the specified rate does not have #> Units: [µg/L] #> 1 m 5 m 10 m 20 m #> 0.60169999 0.18937304 0.10402698 0.05517095 + +# The function is vectorised with respect to rates, applications, water depth, +# crop groups and distances +PEC_sw_drift( + rate = rep(100, 6), + applications = c(1, 2, rep(1, 4)), + water_depth = c(30, 30, 30, 60, 30, 30), + crop_group_JKI = c(rep("Ackerbau", 4), rep("Obstbau frueh", 2)), + distances = c(rep(5, 4), 10, 5)) +#> Units: [µg/L] +#> 5 m 5 m 5 m 5 m 10 m 5 m +#> 0.1900000 0.1566667 0.1900000 0.0950000 3.9366667 6.6300000 + +# Try the same with the Rautmann formula +PEC_sw_drift( + rate = rep(100, 6), + applications = c(1, 2, rep(1, 4)), + water_depth = c(30, 30, 30, 60, 30, 30), + drift_data = "RF", + crop_group_RF = c(rep("arable", 4), rep("fruit, early", 2)), + distances = c(rep(5, 4), 10, 5)) +#> Units: [µg/L] +#> 5 m 5 m 5 m 5 m 10 m 5 m +#> 0.19064473 0.15991216 0.19064473 0.09532236 3.93566026 6.62814740 + +# And with the FOCUS variant +PEC_sw_drift( + rate = rep(100, 6), + applications = c(1, 2, rep(1, 4)), + water_depth = c(30, 30, 30, 60, 30, 30), + drift_data = "RF", + formula = "FOCUS", + crop_group_RF = c(rep("arable", 4), rep("fruit, early", 2)), + distances = c(rep(5, 4), 10, 5)) +#> Units: [µg/L] +#> 5 m 5 m 5 m 5 m 10 m 5 m +#> 0.1741454 0.1456444 0.1741454 0.0870727 3.7957683 6.1809560

See also