Example datasets from the NAFTA SOP published 2015

Data taken from US EPA (2015), p. 19 and 23.

NAFTA_SOP_2015

Format

2 datasets with observations on the following variables.

name

a factor containing the name of the observed variable

time

a numeric vector containing time points

value

a numeric vector containing concentrations

Source

NAFTA (2011) Guidance for evaluating and calculating degradation kinetics in environmental media. NAFTA Technical Working Group on Pesticides https://www.epa.gov/pesticide-science-and-assessing-pesticide-risks/guidance-evaluating-and-calculating-degradation accessed 2019-02-22

US EPA (2015) Standard Operating Procedure for Using the NAFTA Guidance to Calculate Representative Half-life Values and Characterizing Pesticide Degradation https://www.epa.gov/pesticide-science-and-assessing-pesticide-risks/standard-operating-procedure-using-nafta-guidance

Examples

  nafta_evaluation <- nafta(NAFTA_SOP_Appendix_D, cores = 1)
#> Warning: Could not calculate correlation; no covariance matrix
#> Warning: Could not calculate correlation; no covariance matrix
#> Warning: Could not calculate correlation; no covariance matrix
#> The SFO model is rejected as S_SFO is equal or higher than the critical value S_c
#> The representative half-life of the IORE model is longer than the one corresponding
#> to the terminal degradation rate found with the DFOP model.
#> The representative half-life obtained from the DFOP model may be used
  print(nafta_evaluation)
#> Sums of squares:
#>       SFO      IORE      DFOP 
#> 1378.6832  615.7730  517.8836 
#> 
#> Critical sum of squares for checking the SFO model:
#> [1] 717.4598
#> 
#> Parameters:
#> $SFO
#>               Estimate Pr(>t) Lower Upper
#> parent_0       83.7558     NA    NA    NA
#> k_parent_sink   0.0017     NA    NA    NA
#> sigma           8.7518     NA    NA    NA
#> 
#> $IORE
#>                     Estimate Pr(>t) Lower Upper
#> parent_0            9.69e+01     NA    NA    NA
#> k__iore_parent_sink 8.40e-14     NA    NA    NA
#> N_parent            6.68e+00     NA    NA    NA
#> sigma               5.85e+00     NA    NA    NA
#> 
#> $DFOP
#>          Estimate Pr(>t) Lower Upper
#> parent_0 9.76e+01     NA    NA    NA
#> k1       4.24e-02     NA    NA    NA
#> k2       8.24e-04     NA    NA    NA
#> g        2.88e-01     NA    NA    NA
#> sigma    5.36e+00     NA    NA    NA
#> 
#> 
#> DTx values:
#>      DT50    DT90 DT50_rep
#> SFO   407    1350      407
#> IORE  541 5190000  1560000
#> DFOP  429    2380      841
#> 
#> Representative half-life:
#> [1] 841.41
  plot(nafta_evaluation)