aboutsummaryrefslogtreecommitdiff
path: root/R
diff options
context:
space:
mode:
authorJohannes Ranke <jranke@uni-bremen.de>2020-06-15 09:05:27 +0200
committerJohannes Ranke <jranke@uni-bremen.de>2020-06-15 09:05:27 +0200
commit068b168d78c9550a77514e0c39923f8d40b48b03 (patch)
treebea0c6801e687d4e0f326ddbdbf5a0a27df8f355 /R
parente6f9e9ca89e35e610d9895b979f1351a47451db0 (diff)
Fix a bug, update for current vdiffr
Vignette FOCUS_L failed as I had introduced a bug in the handling of warnings. Current vdiffr only runs visual tests if R < 4.1.0, skipping r-devel for now, see https://github.com/r-lib/vdiffr/commit/630a29d013361fd63fea242f531e2db6aef37919
Diffstat (limited to 'R')
-rw-r--r--R/mkinfit.R2
-rw-r--r--R/summary_DFOP_FOCUS_C.txt82
2 files changed, 1 insertions, 83 deletions
diff --git a/R/mkinfit.R b/R/mkinfit.R
index ec2d3412..154c2a18 100644
--- a/R/mkinfit.R
+++ b/R/mkinfit.R
@@ -853,7 +853,7 @@ mkinfit <- function(mkinmod, observed,
if (fit$convergence != 0) {
convergence_warning = paste0("Optimisation did not converge:\n", fit$message)
- summary_warnings <- c(warnings, convergence_warning)
+ summary_warnings <- c(summary_warnings, convergence_warning)
warning(convergence_warning)
} else {
if(!quiet) message("Optimisation successfully terminated.\n")
diff --git a/R/summary_DFOP_FOCUS_C.txt b/R/summary_DFOP_FOCUS_C.txt
deleted file mode 100644
index ab64a588..00000000
--- a/R/summary_DFOP_FOCUS_C.txt
+++ /dev/null
@@ -1,82 +0,0 @@
-mkin version used for fitting: Dummy 0.0 for testing
-R version used for fitting: Dummy R version for testing
-Date of fit: Dummy date for testing
-Date of summary: Dummy date for testing
-
-Equations:
-d_parent/dt = - ((k1 * g * exp(-k1 * time) + k2 * (1 - g) * exp(-k2 *
- time)) / (g * exp(-k1 * time) + (1 - g) * exp(-k2 * time)))
- * parent
-
-Model predictions using solution type analytical
-
-Fitted using test 0 model solutions performed in test time 0 s
-
-Error model: Constant variance
-
-Error model algorithm: OLS
-
-Starting values for parameters to be optimised:
- value type
-parent_0 85.10 state
-k1 0.10 deparm
-k2 0.01 deparm
-g 0.50 deparm
-
-Starting values for the transformed parameters actually optimised:
- value lower upper
-parent_0 85.100000 -Inf Inf
-log_k1 -2.302585 -Inf Inf
-log_k2 -4.605170 -Inf Inf
-g_ilr 0.000000 -Inf Inf
-
-Fixed parameter values:
-None
-
-Results:
-
- AIC BIC logLik
- 29.02372 30.00984 -9.511861
-
-Optimised, transformed parameters with symmetric confidence intervals:
- Estimate Std. Error Lower Upper
-parent_0 85.0000 0.66620 83.1500 86.8500
-log_k1 -0.7775 0.03380 -0.8713 -0.6836
-log_k2 -4.0260 0.13100 -4.3890 -3.6620
-g_ilr 1.2490 0.05811 1.0870 1.4100
-sigma 0.6962 0.16410 0.2406 1.1520
-
-Parameter correlation:
-[1] "Correlation matrix is platform dependent, not tested"
-
-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 85.00000 127.600 1.131e-08 83.15000 86.85000
-k1 0.45960 29.580 3.887e-06 0.41840 0.50480
-k2 0.01785 7.636 7.901e-04 0.01241 0.02568
-g 0.85390 83.310 6.221e-08 0.82310 0.88020
-sigma 0.69620 4.243 6.618e-03 0.24060 1.15200
-
-FOCUS Chi2 error levels in percent:
- err.min n.optim df
-All data 2.661 4 5
-parent 2.661 4 5
-
-Estimated disappearance times:
- DT50 DT90 DT50_k1 DT50_k2
-parent 1.887 21.25 1.508 38.83
-
-Data:
- time variable observed predicted residual
- 0 parent 85.1 85.003 0.09726
- 1 parent 57.9 58.039 -0.13912
- 3 parent 29.9 30.054 -0.15351
- 7 parent 14.6 13.866 0.73388
- 14 parent 9.7 9.787 -0.08657
- 28 parent 6.6 7.532 -0.93205
- 63 parent 4.0 4.033 -0.03269
- 91 parent 3.9 2.447 1.45348
- 119 parent 0.6 1.484 -0.88424

Contact - Imprint