JKI drift percentage data for field crops and pome/stone fruit

8 files changed, 128 insertions, 1 deletions

diff --git a/.gitignore b/.gitignore
index 4e7e1ab..614e3ac 100644
--- a/.gitignore
+++ b/.gitignore
@@ -4,3 +4,4 @@ pfm\.Rcheck
 pkg/ChangeLog
 test.R
 README.html
+pkg/inst/extdata/Tabelle\ der\ Abdrifteckwerte.xls
diff --git a/README.md b/README.md
index d0d08ab..1f70a86 100644
--- a/README.md
+++ b/README.md
@@ -92,3 +92,29 @@ PEC_soil(100, int = 0.25)
 ```
 ## [1] 0.1
 ```
+
+### Rautmann drift data
+
+Some of the drift percentage data published by the JKI are included. To
+see the data for one application:
+
+
+```r
+drift_data_JKI[1]
+```
+
+```
+## [[1]]
+##            crop
+## distance    Ackerbau Obstbau früh Obstbau spät
+##   1.000000      2.77           NA           NA
+##   3.000000        NA        29.20        15.73
+##   5.000000      0.57        19.89         8.41
+##   10.000000     0.29        11.81         3.60
+##   15.000000     0.20         5.55         1.81
+##   20.000000     0.15         2.77         1.09
+##   30.000000     0.10         1.04         0.54
+##   40.000000     0.07         0.52         0.32
+##   50.000000     0.06         0.30         0.22
+```
+
diff --git a/README.rmd b/README.rmd
index b8a5829..c4c514e 100644
--- a/README.rmd
+++ b/README.rmd
@@ -57,3 +57,14 @@ density of 1.5 kg/L, output in mg/kg:
 ```{r}
 PEC_soil(100, int = 0.25)
 ```
+
+### Rautmann drift data
+
+Some of the drift percentage data published by the JKI are included. To
+see the data for one application:
+
+```{r}
+drift_data_JKI[1]
+```
+
+
diff --git a/pkg/.Rbuildignore b/pkg/.Rbuildignore
index c56e58e..c48d8e3 100644
--- a/pkg/.Rbuildignore
+++ b/pkg/.Rbuildignore
@@ -1,2 +1,3 @@
 inst/testdata/SwashProjects/project_H_sw/MACRO/*/*.bin
 inst/testdata/Swash/
+inst/extdata/Tabelle\ der\ Abdrifteckwerte.xls
diff --git a/pkg/DESCRIPTION b/pkg/DESCRIPTION
index a6880f2..98e9142 100644
--- a/pkg/DESCRIPTION
+++ b/pkg/DESCRIPTION
@@ -2,7 +2,7 @@ Package: pfm
 Type: Package
 Title: Utilities for pesticide fate modelling
 Version: 0.1-3
-Date: 2015-05-07
+Date: 2015-06-11
 Authors@R: person("Johannes Ranke", email = "jranke@uni-bremen.de", 
                   role = c("aut", "cre", "cph"))
 Description: Utilities for simple PEC calculations and for dealing with data
diff --git a/pkg/R/drift_data_JKI.R b/pkg/R/drift_data_JKI.R
new file mode 100644
index 0000000..44a959b
--- /dev/null
+++ b/pkg/R/drift_data_JKI.R
@@ -0,0 +1,42 @@
+#' Deposition from spray drift expressed as percent of the applied dose as
+#' published by the JKI
+#' 
+#' Deposition from spray drift expressed as percent of the applied dose as
+#' published by the German Julius-Kühn Institute (JKI).
+#' 
+#' The data were extracted from the spreadsheet cited below using the R code
+#' given in the example section. The spreadsheet is not included in the package
+#' as its licence is not clear.
+#' 
+#' 
+#' @name drift_data_JKI
+#' @docType data
+#' @format A list currently containing matrices with spray drift percentage
+#' data for field crops (Ackerbau), and Pome/stone fruit, early and late
+#' (Obstbau früh, spät).
+#' @source JKI (2010) Spreadsheet 'Tabelle der Abdrifteckwerte.xls', retrieved
+#' from
+#' http://www.jki.bund.de/no_cache/de/startseite/institute/anwendungstechnik/abdrift-eckwerte.html
+#' on 2015-06-11
+#' @keywords datasets
+#' @examples
+#' 
+#' \dontrun{
+#'   # This is the code that was used to extract the data
+#'   library(readxl)
+#'   abdrift_path <- "../inst/extdata/Tabelle der Abdrifteckwerte.xls"
+#'   JKI_crops <- c("Ackerbau", "Obstbau früh", "Obstbau spät")
+#'   names(JKI_crops) <- c("Field crops", "Pome/stone fruit, early", "Pome/stone fruit, late")
+#'   drift_data_JKI <- list()
+#' 
+#'   for (n in 1:8) {
+#'     drift_data_raw <- read_excel(abdrift_path, sheet = n + 1, skip = 2)
+#'     drift_data <- as.matrix(drift_data_raw[1:9, 2:4]) 
+#'     dimnames(drift_data) <- list(distance = as.integer(drift_data_raw[1:9, 1]),
+#'                                             crop = JKI_crops)
+#'     drift_data_JKI[[n]] <- drift_data
+#'   }
+#'   save(drift_data_JKI, file = "../data/drift_data_JKI.RData")
+#' }
+#' 
+NULL
diff --git a/pkg/data/drift_data_JKI.RData b/pkg/data/drift_data_JKI.RData
new file mode 100644
index 0000000..061ca27
--- /dev/null
+++ b/pkg/data/drift_data_JKI.RData
diff --git a/pkg/man/drift_data_JKI.Rd b/pkg/man/drift_data_JKI.Rd
new file mode 100644
index 0000000..a2bbb77
--- /dev/null
+++ b/pkg/man/drift_data_JKI.Rd
@@ -0,0 +1,46 @@
+% Generated by roxygen2 (4.1.0.9001): do not edit by hand
+% Please edit documentation in R/drift_data_JKI.R
+\docType{data}
+\name{drift_data_JKI}
+\alias{drift_data_JKI}
+\title{Deposition from spray drift expressed as percent of the applied dose as
+published by the JKI}
+\format{A list currently containing matrices with spray drift percentage
+data for field crops (Ackerbau), and Pome/stone fruit, early and late
+(Obstbau früh, spät).}
+\source{
+JKI (2010) Spreadsheet 'Tabelle der Abdrifteckwerte.xls', retrieved
+from
+http://www.jki.bund.de/no_cache/de/startseite/institute/anwendungstechnik/abdrift-eckwerte.html
+on 2015-06-11
+}
+\description{
+Deposition from spray drift expressed as percent of the applied dose as
+published by the German Julius-Kühn Institute (JKI).
+}
+\details{
+The data were extracted from the spreadsheet cited below using the R code
+given in the example section. The spreadsheet is not included in the package
+as its licence is not clear.
+}
+\examples{
+\dontrun{
+  # This is the code that was used to extract the data
+  library(readxl)
+  abdrift_path <- "../inst/extdata/Tabelle der Abdrifteckwerte.xls"
+  JKI_crops <- c("Ackerbau", "Obstbau früh", "Obstbau spät")
+  names(JKI_crops) <- c("Field crops", "Pome/stone fruit, early", "Pome/stone fruit, late")
+  drift_data_JKI <- list()
+
+  for (n in 1:8) {
+    drift_data_raw <- read_excel(abdrift_path, sheet = n + 1, skip = 2)
+    drift_data <- as.matrix(drift_data_raw[1:9, 2:4])
+    dimnames(drift_data) <- list(distance = as.integer(drift_data_raw[1:9, 1]),
+                                            crop = JKI_crops)
+    drift_data_JKI[[n]] <- drift_data
+  }
+  save(drift_data_JKI, file = "../data/drift_data_JKI.RData")
+}
+}
+\keyword{datasets}
+