@@ -190,96 +198,97 @@ distribution
The summary function returns a list with components, among others
version, RversionThe mkin and R versions used
-
date.fit, date.summaryThe dates where the fit and the summary were
- produced
-
diffsThe differential equations used in the model
-
use_of_ffWas maximum or minimum use made of formation fractions
-
bparOptimised and backtransformed
- parameters
-
dataThe data (see Description above).
-
startThe starting values and bounds, if applicable, for optimised
- parameters.
-
fixedThe values of fixed parameters.
-
errmin The chi2 error levels for
- each observed variable.
-
bparms.odeAll backtransformed ODE
- parameters, for use as starting parameters for related models.
-
errparmsError model parameters.
-
ffThe estimated formation fractions derived from the fitted
- model.
-
distimesThe DT50 and DT90 values for each observed variable.
-
SFORBIf applicable, eigenvalues of SFORB components of the model.
- The print method is called for its side effect, i.e. printing the summary.
+
date.fit, date.summaryThe dates where the fit and the summary were
+produced
+
diffsThe differential equations used in the model
+
use_of_ffWas maximum or minimum use made of formation fractions
+
bparOptimised and backtransformed
+parameters
+
dataThe data (see Description above).
+
startThe starting values and bounds, if applicable, for optimised
+parameters.
+
fixedThe values of fixed parameters.
+
errmin The chi2 error levels for
+each observed variable.
+
bparms.odeAll backtransformed ODE
+parameters, for use as starting parameters for related models.
+
errparmsError model parameters.
+
ffThe estimated formation fractions derived from the fitted
+model.
+
distimesThe DT50 and DT90 values for each observed variable.
+
SFORBIf applicable, eigenvalues of SFORB components of the model.
+The print method is called for its side effect, i.e. printing the summary.
References
FOCUS (2006) “Guidance Document on Estimating Persistence
- and Degradation Kinetics from Environmental Fate Studies on Pesticides in
- EU Registration” Report of the FOCUS Work Group on Degradation Kinetics,
- EC Document Reference Sanco/10058/2005 version 2.0, 434 pp,
- http://esdac.jrc.ec.europa.eu/projects/degradation-kinetics
+and Degradation Kinetics from Environmental Fate Studies on Pesticides in
+EU Registration” Report of the FOCUS Work Group on Degradation Kinetics,
+EC Document Reference Sanco/10058/2005 version 2.0, 434 pp,
+
http://esdac.jrc.ec.europa.eu/projects/degradation-kinetics
Examples
#> mkin version used for fitting: 0.9.49.9
-#> R version used for fitting: 3.6.3
-#> Date of fit: Mon Mar 30 21:43:00 2020
-#> Date of summary: Mon Mar 30 21:43:00 2020
+
summary(
mkinfit(
mkinmod(
parent = mkinsub(
"SFO")),
FOCUS_2006_A,
quiet = TRUE))
#> mkin version used for fitting: 0.9.50.1
+#> R version used for fitting: 4.0.0
+#> Date of fit: Tue May 12 08:07:03 2020
+#> Date of summary: Tue May 12 08:07:03 2020
#>
#> Equations:
-#> d_parent/dt = - k_parent_sink * parent
+#> d_parent/dt = - k_parent * parent
#>
#> Model predictions using solution type analytical
#>
-#> Fitted using 131 model solutions performed in 0.306 s
+#> Fitted using 131 model solutions performed in 0.026 s
#>
#> Error model: Constant variance
#>
#> Error model algorithm: OLS
#>
#> Starting values for parameters to be optimised:
-#> value type
-#> parent_0 101.24 state
-#> k_parent_sink 0.10 deparm
+#> value type
+#> parent_0 101.24 state
+#> k_parent 0.10 deparm
#>
#> Starting values for the transformed parameters actually optimised:
-#> value lower upper
-#> parent_0 101.240000 -Inf Inf
-#> log_k_parent_sink -2.302585 -Inf Inf
+#> value lower upper
+#> parent_0 101.240000 -Inf Inf
+#> log_k_parent -2.302585 -Inf Inf
#>
#> Fixed parameter values:
#> None
#>
+#> Results:
+#>
+#> AIC BIC logLik
+#> 55.28197 55.5203 -24.64099
+#>
#> Optimised, transformed parameters with symmetric confidence intervals:
-#> Estimate Std. Error Lower Upper
-#> parent_0 109.200 3.70400 99.630 118.700
-#> log_k_parent_sink -3.291 0.09176 -3.527 -3.055
-#> sigma 5.266 1.31600 1.882 8.649
+#> Estimate Std. Error Lower Upper
+#> parent_0 109.200 3.70400 99.630 118.700
+#> log_k_parent -3.291 0.09176 -3.527 -3.055
+#> sigma 5.266 1.31600 1.882 8.649
#>
#> Parameter correlation:
-#> parent_0 log_k_parent_sink sigma
-#> parent_0 1.000e+00 5.428e-01 1.642e-07
-#> log_k_parent_sink 5.428e-01 1.000e+00 2.507e-07
-#> sigma 1.642e-07 2.507e-07 1.000e+00
+#> parent_0 log_k_parent sigma
+#> parent_0 1.000e+00 5.428e-01 1.642e-07
+#> log_k_parent 5.428e-01 1.000e+00 2.507e-07
+#> sigma 1.642e-07 2.507e-07 1.000e+00
#>
#> Backtransformed parameters:
#> Confidence intervals for internally transformed parameters are asymmetric.
#> t-test (unrealistically) based on the assumption of normal distribution
#> for estimators of untransformed parameters.
-#> Estimate t value Pr(>t) Lower Upper
-#> parent_0 109.20000 29.47 4.218e-07 99.6300 118.70000
-#> k_parent_sink 0.03722 10.90 5.650e-05 0.0294 0.04712
-#> sigma 5.26600 4.00 5.162e-03 1.8820 8.64900
+#> Estimate t value Pr(>t) Lower Upper
+#> parent_0 109.20000 29.47 4.218e-07 99.6300 118.70000
+#> k_parent 0.03722 10.90 5.650e-05 0.0294 0.04712
+#> sigma 5.26600 4.00 5.162e-03 1.8820 8.64900
#>
#> FOCUS Chi2 error levels in percent:
#> err.min n.optim df
#> All data 8.385 2 6
#> parent 8.385 2 6
#>
-#> Resulting formation fractions:
-#> ff
-#> parent_sink 1
-#>
#> Estimated disappearance times:
#> DT50 DT90
#> parent 18.62 61.87
@@ -296,17 +305,10 @@ distribution
#> 118 parent 0.39 1.351 -0.9613
-
@@ -317,7 +319,7 @@ distribution