blob: f9171748078e196a3c7941bc8d37ef3949d991fb (
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
|
nlme version used for fitting: Dummy 0.0 for testing
mkin version used for pre-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
d_m1/dt = + f_parent_to_m1 * ((k1 * g * exp(-k1 * time) + k2 * (1 - g)
* exp(-k2 * time)) / (g * exp(-k1 * time) + (1 - g) *
exp(-k2 * time))) * parent - k_m1 * m1
Data:
509 observations of 2 variable(s) grouped in 15 datasets
Model predictions using solution type analytical
Fitted in test time 0 s using 2 iterations
Variance model: Constant variance
Mean of starting values for individual parameters:
parent_0 log_k_m1 f_parent_qlogis log_k1 log_k2
100.70 -6.30 -0.08 -3.09 -3.95
g_qlogis
0.03
Fixed degradation parameter values:
value type
m1_0 0 state
Results:
AIC BIC logLik
2711 2766 -1343
Optimised, transformed parameters with symmetric confidence intervals:
lower est. upper
parent_0 99.5 100.37 101.23
log_k_m1 -6.5 -6.23 -5.94
f_parent_qlogis -0.2 -0.08 0.08
log_k1 -3.6 -3.22 -2.84
log_k2 -4.5 -4.10 -3.68
g_qlogis -0.4 -0.10 0.17
Correlation:
prnt_0 lg_k_1 f_prn_ log_k1 log_k2
log_k_m1 -0.185
f_parent_qlogis -0.161 0.405
log_k1 0.056 -0.014 -0.016
log_k2 0.025 0.011 -0.004 0.026
g_qlogis -0.032 -0.046 -0.012 -0.109 -0.103
Random effects:
Formula: list(parent_0 ~ 1, log_k_m1 ~ 1, f_parent_qlogis ~ 1, log_k1 ~ 1, log_k2 ~ 1, g_qlogis ~ 1)
Level: ds
Structure: Diagonal
parent_0 log_k_m1 f_parent_qlogis log_k1 log_k2 g_qlogis Residual
StdDev: 1 2e-04 0.3 0.7 0.8 0.3 3
Backtransformed parameters with asymmetric confidence intervals:
lower est. upper
parent_0 1e+02 1e+02 1e+02
k_m1 1e-03 2e-03 3e-03
f_parent_to_m1 4e-01 5e-01 5e-01
k1 3e-02 4e-02 6e-02
k2 1e-02 2e-02 3e-02
g 4e-01 5e-01 5e-01
Resulting formation fractions:
ff
parent_m1 0.5
parent_sink 0.5
Estimated disappearance times:
DT50 DT90 DT50back DT50_k1 DT50_k2
parent 27 105 31 17 42
m1 352 1171 NA NA NA
|