aboutsummaryrefslogtreecommitdiff
path: root/README.rmd
blob: f1c259920fda842963b0abe0570ff10af2238de5 (plain) (blame)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
# pfm
```{r, echo = FALSE}
options(device = function(file, width = 7, height = 7, ...) {
          cairo_pdf(tempfile(), width = width, height = height, ...)
})
```
The R package **pfm** provides some utilities for dealing with FOCUS pesticide fate modelling tools,
(currently only TOXSWA cwa and out files), made available under the GNU public license.
This means:

    This program is free software: you can redistribute it and/or modify it under
    the terms of the GNU General Public License as published by the Free Software
    Foundation, either version 3 of the License, or (at your option) any later
    version.

    This program is distributed in the hope that it will be useful, but WITHOUT
    ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
    FOR A PARTICULAR PURPOSE.  See the GNU General Public License for more
    details.

    You should have received a copy of the GNU General Public License along with
    this program. If not, see <http://www.gnu.org/licenses/>

## Installation

You can install the package from [github](http://github.com/jranke/pfm), e.g.
using the `devtools` package.  Using `quick = TRUE` skips docs,
multiple-architecture builds, demos, and vignettes, to make installation as
fast and painless as possible.

```{r, eval = FALSE}
library(devtools)
install_github("jranke/pfm", subdir = "pkg", quick = TRUE)
```

## Use

### Analyse TOXSWA output

Read in and analyse a cwa file:


```{r}
library(pfm, quietly = TRUE)
example_cwa <- read.TOXSWA_cwa("00003s_pa.cwa")
plot(example_cwa)
```

Get events above thresholds of 20 and 100 µg/L,
and do a moving window analysis for windows of 7 days
and 21 days, print the results:

```{r}
example_cwa$get_events(c(20, 100))
example_cwa$moving_windows(c(7, 21))
print(example_cwa)
```

This can also be done with out files, the function reads
out files from current TOXSWA versions as well as cwa files
from old TOXSWA versions.


### Calculate PEC soil

Simple PEC soil calculation for an application rate of 100 g/ha and
25% interception, assuming complete mixing into 5 cm and a soil bulk
density of 1.5 kg/L, output in mg/kg:


```{r}
PEC_soil(100, interception = 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]
```

### PEC surface water due to drift

Initial PEC values for an application of 100 g/ha in the vicinity of a 30 cm
deep water body are obtained using


```{r}
PEC_sw_drift(100, applications = 1)
```

Contact - Imprint