From 35eaf4f7bac89b249ff488e4e3021d4226ed9e23 Mon Sep 17 00:00:00 2001 From: Johannes Ranke Date: Wed, 15 Sep 2021 11:22:36 +0200 Subject: Correction of the dimethenamid dataset The data given in previous mkin versions for "Borstel 2" (Staudenmaier 2013) are in fact a subset of the "Borstel 1" data, with the "Borstel 2" data simply being further analysed for isomer composition by Staudenmaier (2013). While correcting this, the URL of the RAR published by EFSA was updated, as the EFSA Register of Questions was replaced by the Open EFSA portal. This also affects the recent publication in "Environments" [1], where these two datasets were pooled. However, the influence on the final results can be expected to be minimal. Out of the 16 observations for the Borstel soil, five were erroneously duplicated which just changes the weight of the different sampling points and modifies the variance. The conclusions of the publication are qualitative and are not considered to be affected. [1] https://doi.org/10.3390/environments8080071 --- inst/dataset_generation/dimethenamid_2018.R | 24 ++++++++---------------- 1 file changed, 8 insertions(+), 16 deletions(-) (limited to 'inst/dataset_generation') diff --git a/inst/dataset_generation/dimethenamid_2018.R b/inst/dataset_generation/dimethenamid_2018.R index 34a9dad1..52eef995 100644 --- a/inst/dataset_generation/dimethenamid_2018.R +++ b/inst/dataset_generation/dimethenamid_2018.R @@ -15,15 +15,7 @@ dimethenamid_2018 <- mkindsg$new( ) ) ), - mkinds$new("Borstel 1", # p. 22 - mkin_wide_to_long( - data.frame( - t = c(0, 28, 58, 89, 119), - DMTAP = c(100.5, 51.4, 26.8, 15.7, 7.9) - ) - ) - ), - mkinds$new("Borstel 2", # p. 19/20 + mkinds$new("Borstel", # p. 19/20 mkin_wide_to_long( data.frame( t = rep(c(0, 2, 7, 14, 28, 58, 89, 119), each = 2), @@ -134,20 +126,20 @@ dimethenamid_2018 <- mkindsg$new( # but the FOCUS generic guidance suggests them to be equivalent to the USDA # classes with the same name meta = data.frame( - study = c("Unsworth 2014", "Staudenmaier 2013", "Staudenmaier 2009", + study = c("Unsworth 2014", "Staudenmaier 2009", rep("Wendt 1997", 2), "König 1996", rep("König 1995", 2)), - usda_soil_type = c("Sandy loam", rep("Sand", 2), + usda_soil_type = c("Sandy loam", "Sand", rep("Clay loam", 2), "Sandy clay loam", "Loamy sand", "Sandy loam"), - study_moisture_ref_type = c("pF2", rep("pF1", 2), + study_moisture_ref_type = c("pF2", "pF1", rep("pF2.5", 2), rep("pF1", 3)), - rel_moisture = c(1, rep(0.5, 2), + rel_moisture = c(1, 0.5, rep(0.75, 2), rep(0.4, 3)), - study_ref_moisture = c(NA, rep(23, 2), rep(33.37, 2), rep(NA, 3)), - temperature = c(20, rep(20, 2), + study_ref_moisture = c(NA, 23, rep(33.37, 2), rep(NA, 3)), + temperature = c(20, 20, rep(23, 2), rep(20, 3)) ) ) -f_time_norm_focus(dimethenamid_2018) +f_time_norm_focus(dimethenamid_2018) # This modifies the R6 object in place save(dimethenamid_2018, file = "data/dimethenamid_2018.rda", version = 2) -- cgit v1.2.1