blob: a19824ce5d10c06be268e71337f9f07c22868701 (
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
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
|
saemix 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:
510 observations of 2 variable(s) grouped in 15 datasets
Model predictions using solution type analytical
Fitted in test time 0 s
Using 300, 100 iterations and 4 chains
Variance model: Two-component variance function
Starting values for degradation parameters:
parent_0 k_m1 f_parent_to_m1 k1 k2
1e+02 7e-03 5e-01 1e-01 2e-02
g
5e-01
Fixed degradation parameter values:
None
Starting values for random effects (square root of initial entries in omega):
parent_0 k_m1 f_parent_to_m1 k1 k2 g
parent_0 101 0 0 0 0 0
k_m1 0 1 0 0 0 0
f_parent_to_m1 0 0 1 0 0 0
k1 0 0 0 1 0 0
k2 0 0 0 0 1 0
g 0 0 0 0 0 1
Starting values for error model parameters:
a.1 b.1
1 1
Results:
Likelihood computed by importance sampling
AIC BIC logLik
2334 2344 -1153
Optimised parameters:
est. lower upper
parent_0 1e+02 1e+02 1e+02
k_m1 7e-03 6e-03 7e-03
f_parent_to_m1 5e-01 4e-01 5e-01
k1 1e-01 9e-02 1e-01
k2 2e-02 2e-02 3e-02
g 5e-01 5e-01 5e-01
a.1 9e-01 8e-01 1e+00
b.1 5e-02 5e-02 6e-02
SD.parent_0 3e-02 -5e+01 5e+01
SD.k_m1 2e-01 1e-01 3e-01
SD.f_parent_to_m1 3e-01 2e-01 4e-01
SD.k1 4e-01 2e-01 5e-01
SD.k2 3e-01 2e-01 5e-01
SD.g 2e-01 6e-02 4e-01
Correlation:
pr_0 k_m1 f___ k1 k2
k_m1 -0.2
f_parent_to_m1 -0.3 0.1
k1 0.1 0.0 0.0
k2 0.0 0.0 0.0 0.1
g 0.1 -0.1 0.0 -0.2 -0.2
Random effects:
est. lower upper
SD.parent_0 0.03 -49.24 49.3
SD.k_m1 0.23 0.13 0.3
SD.f_parent_to_m1 0.30 0.19 0.4
SD.k1 0.40 0.25 0.5
SD.k2 0.34 0.21 0.5
SD.g 0.21 0.06 0.4
Variance model:
est. lower upper
a.1 0.93 0.80 1.06
b.1 0.05 0.05 0.06
Resulting formation fractions:
ff
parent_m1 0.5
parent_sink 0.5
Estimated disappearance times:
DT50 DT90 DT50back DT50_k1 DT50_k2
parent 13 73 22 6 32
m1 105 348 NA NA NA
|
