From 8a3475c59f3d91ce5ce7d980d6de09360617e7fe Mon Sep 17 00:00:00 2001 From: Johannes Ranke Date: Tue, 7 May 2019 08:12:27 +0200 Subject: After the OLS step, use OLS parameter estimates - Fix the respective error in the code - Static documentation rebuilt by pkgdown --- docs/reference/test_data_from_UBA_2014.html | 54 ++++++++++------------------- 1 file changed, 18 insertions(+), 36 deletions(-) (limited to 'docs/reference/test_data_from_UBA_2014.html') diff --git a/docs/reference/test_data_from_UBA_2014.html b/docs/reference/test_data_from_UBA_2014.html index bc988340..11cb6420 100644 --- a/docs/reference/test_data_from_UBA_2014.html +++ b/docs/reference/test_data_from_UBA_2014.html @@ -159,20 +159,13 @@ # model). m_ws <- mkinmod(parent_w = mkinsub("SFO", "parent_s"), parent_s = mkinsub("SFO", "parent_w"))
#> Successfully compiled differential equation model from auto-generated C code.
f_river <- mkinfit(m_ws, test_data_from_UBA_2014[[1]]$data, quiet = TRUE)
#> Warning: Observations with value of zero were removed from the data
plot_sep(f_river)
- summary(f_river)$bpar
#> Estimate se_notrans t value Pr(>t) -#> parent_w_0 9.598567e+01 2.12351789 4.520126e+01 9.476190e-18 -#> k_parent_w_sink 3.603743e-01 0.03149282 1.144306e+01 4.126593e-09 -#> k_parent_w_parent_s 6.031371e-02 0.01603582 3.761186e+00 9.436275e-04 -#> k_parent_s_sink 5.108964e-11 0.09482736 5.387647e-10 5.000000e-01 -#> k_parent_s_parent_w 7.419672e-02 0.10737376 6.910135e-01 2.500560e-01 -#> sigma 2.982879e+00 0.50545649 5.901356e+00 1.454535e-05 -#> Lower Upper -#> parent_w_0 91.48420503 100.4871438 -#> k_parent_w_sink 0.30668904 0.4234571 -#> k_parent_w_parent_s 0.03423904 0.1062455 -#> k_parent_s_sink 0.00000000 Inf -#> k_parent_s_parent_w 0.02289956 0.2404043 -#> sigma 2.00184022 3.9639169
mkinerrmin(f_river)
#> err.min n.optim df + summary(f_river)$bpar
#> Warning: Could not calculate correlation; no covariance matrix
#> Estimate se_notrans t value Pr(>t) Lower Upper +#> parent_w_0 9.598567e+01 NA NA NA NA NA +#> k_parent_w_sink 3.603743e-01 NA NA NA NA NA +#> k_parent_w_parent_s 6.031371e-02 NA NA NA NA NA +#> k_parent_s_sink 5.108964e-11 NA NA NA NA NA +#> k_parent_s_parent_w 7.419672e-02 NA NA NA NA NA +#> sigma 2.982879e+00 NA NA NA NA NA
mkinerrmin(f_river)
#> err.min n.optim df #> All data 0.09246946 5 6 #> parent_w 0.06377096 3 3 #> parent_s 0.20882325 2 3
@@ -183,28 +176,17 @@ M2 = mkinsub("SFO", "M3"), M3 = mkinsub("SFO"), use_of_ff = "max")
#> Successfully compiled differential equation model from auto-generated C code.
- f_soil <- mkinfit(m_soil, test_data_from_UBA_2014[[3]]$data, quiet = TRUE)
#> Warning: Observations with value of zero were removed from the data
plot_sep(f_soil, lpos = c("topright", "topright", "topright", "bottomright"))
summary(f_soil)$bpar
#> Estimate se_notrans t value Pr(>t) Lower -#> parent_0 76.55425584 0.859186619 89.1008474 1.113866e-26 74.755959756 -#> k_parent 0.12081956 0.004601922 26.2541544 1.077373e-16 0.111561582 -#> k_M1 0.84258649 0.806231419 1.0450926 1.545475e-01 0.113839803 -#> k_M2 0.04210878 0.017083049 2.4649453 1.170195e-02 0.018013807 -#> k_M3 0.01122919 0.007245890 1.5497322 6.885127e-02 0.002909463 -#> f_parent_to_M1 0.32240199 0.240803555 1.3388589 9.820820e-02 NA -#> f_parent_to_M2 0.16099854 0.033691991 4.7785403 6.531224e-05 NA -#> f_M1_to_M3 0.27921501 0.269443514 1.0362655 1.565440e-01 0.022992937 -#> f_M2_to_M3 0.55641331 0.595125445 0.9349513 1.807725e-01 0.008003317 -#> sigma 1.14005399 0.149696423 7.6157731 1.727024e-07 0.826735778 -#> Upper -#> parent_0 78.35255192 -#> k_parent 0.13084582 -#> k_M1 6.23641265 -#> k_M2 0.09843279 -#> k_M3 0.04333950 -#> f_parent_to_M1 NA -#> f_parent_to_M2 NA -#> f_M1_to_M3 0.86443083 -#> f_M2_to_M3 0.99489847 -#> sigma 1.45337221
mkinerrmin(f_soil)
#> err.min n.optim df + f_soil <- mkinfit(m_soil, test_data_from_UBA_2014[[3]]$data, quiet = TRUE)
#> Warning: Observations with value of zero were removed from the data
plot_sep(f_soil, lpos = c("topright", "topright", "topright", "bottomright"))
summary(f_soil)$bpar
#> Warning: Could not calculate correlation; no covariance matrix
#> Estimate se_notrans t value Pr(>t) Lower Upper +#> parent_0 76.55425584 NA NA NA NA NA +#> k_parent 0.12081956 NA NA NA NA NA +#> k_M1 0.84258649 NA NA NA NA NA +#> k_M2 0.04210878 NA NA NA NA NA +#> k_M3 0.01122919 NA NA NA NA NA +#> f_parent_to_M1 0.32240199 NA NA NA NA NA +#> f_parent_to_M2 0.16099854 NA NA NA NA NA +#> f_M1_to_M3 0.27921501 NA NA NA NA NA +#> f_M2_to_M3 0.55641331 NA NA NA NA NA +#> sigma 1.14005399 NA NA NA NA NA
mkinerrmin(f_soil)
#> err.min n.optim df #> All data 0.09649963 9 20 #> parent 0.04721283 2 6 #> M1 0.26551209 2 5 -- cgit v1.2.1