aboutsummaryrefslogtreecommitdiff
path: root/tests/testthat/summary_nlme_biphasic_s.txt
blob: 65aead6201ba0c9cd5b157ddd65c1139387c7c55 (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 3 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           -5.35           -0.08           -2.68           -4.37 
       g_qlogis 
          -0.33 

Fixed degradation parameter values:
     value  type
m1_0     0 state

Results:

   AIC  BIC logLik
  2683 2738  -1329

Optimised, transformed parameters with symmetric confidence intervals:
                lower   est.  upper
parent_0         99.6 100.43 101.26
log_k_m1         -5.5  -5.34  -5.18
f_parent_qlogis  -0.3  -0.08   0.09
log_k1           -3.2  -2.90  -2.60
log_k2           -4.6  -4.34  -4.07
g_qlogis         -0.5  -0.19   0.08

Correlation: 
                prnt_0 lg_k_1 f_prn_ log_k1 log_k2
log_k_m1        -0.177                            
f_parent_qlogis -0.164  0.385                     
log_k1           0.108 -0.017 -0.025              
log_k2           0.036  0.054  0.008  0.096       
g_qlogis        -0.068 -0.110 -0.030 -0.269 -0.267

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      0.1             0.3    0.6    0.5      0.3        3


Backtransformed parameters with asymmetric confidence intervals:
               lower  est. upper
parent_0       1e+02 1e+02 1e+02
k_m1           4e-03 5e-03 6e-03
f_parent_to_m1 4e-01 5e-01 5e-01
k1             4e-02 6e-02 7e-02
k2             1e-02 1e-02 2e-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   26  131       39      13      53
m1      144  479       NA      NA      NA

Contact - Imprint