blob: 141f57a1955e308240ede73a15968ec3cfa6f85c (
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
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
|
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 = - k_parent_sink * parent - k_parent_m1 * parent
d_m1/dt = + k_parent_m1 * parent - k_m1_sink * m1
Model predictions using solution type eigen
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 100.7500 state
k_parent_sink 0.1000 deparm
k_parent_m1 0.1001 deparm
k_m1_sink 0.1002 deparm
Starting values for the transformed parameters actually optimised:
value lower upper
parent_0 100.750000 -Inf Inf
log_k_parent_sink -2.302585 -Inf Inf
log_k_parent_m1 -2.301586 -Inf Inf
log_k_m1_sink -2.300587 -Inf Inf
Fixed parameter values:
value type
m1_0 0 state
Results:
AIC BIC logLik
204.4486 212.6365 -97.22429
Optimised, transformed parameters with symmetric confidence intervals:
Estimate Std. Error Lower Upper
parent_0 99.600 1.57000 96.400 102.800
log_k_parent_sink -3.038 0.07626 -3.193 -2.883
log_k_parent_m1 -2.980 0.04033 -3.062 -2.898
log_k_m1_sink -5.248 0.13320 -5.518 -4.977
sigma 3.126 0.35850 2.396 3.855
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 99.600000 63.430 2.298e-36 96.400000 1.028e+02
k_parent_sink 0.047920 13.110 6.126e-15 0.041030 5.596e-02
k_parent_m1 0.050780 24.800 3.269e-23 0.046780 5.512e-02
k_m1_sink 0.005261 7.510 6.165e-09 0.004012 6.898e-03
sigma 3.126000 8.718 2.235e-10 2.396000 3.855e+00
FOCUS Chi2 error levels in percent:
err.min n.optim df
All data 6.398 4 15
parent 6.827 3 6
m1 4.490 1 9
Resulting formation fractions:
ff
parent_sink 0.4855
parent_m1 0.5145
m1_sink 1.0000
Estimated disappearance times:
DT50 DT90
parent 7.023 23.33
m1 131.761 437.70
Data:
time variable observed predicted residual
0 parent 99.46 99.59848 not tested
0 parent 102.04 99.59848 not tested
1 parent 93.50 90.23787 not tested
1 parent 92.50 90.23787 not tested
3 parent 63.23 74.07320 not tested
3 parent 68.99 74.07320 not tested
7 parent 52.32 49.91207 not tested
7 parent 55.13 49.91207 not tested
14 parent 27.27 25.01257 not tested
14 parent 26.64 25.01257 not tested
21 parent 11.50 12.53462 not tested
21 parent 11.64 12.53462 not tested
35 parent 2.85 3.14787 not tested
35 parent 2.91 3.14787 not tested
50 parent 0.69 0.71624 not tested
50 parent 0.63 0.71624 not tested
75 parent 0.05 0.06074 not tested
75 parent 0.06 0.06074 not tested
1 m1 4.84 4.80296 not tested
1 m1 5.64 4.80296 not tested
3 m1 12.91 13.02400 not tested
3 m1 12.96 13.02400 not tested
7 m1 22.97 25.04476 not tested
7 m1 24.47 25.04476 not tested
14 m1 41.69 36.69002 not tested
14 m1 33.21 36.69002 not tested
21 m1 44.37 41.65310 not tested
21 m1 46.44 41.65310 not tested
35 m1 41.22 43.31312 not tested
35 m1 37.95 43.31312 not tested
50 m1 41.19 41.21831 not tested
50 m1 40.01 41.21831 not tested
75 m1 40.09 36.44703 not tested
75 m1 33.85 36.44703 not tested
100 m1 31.04 31.98163 not tested
100 m1 33.13 31.98163 not tested
120 m1 25.15 28.78984 not tested
120 m1 33.31 28.78984 not tested
|
