aboutsummaryrefslogtreecommitdiff
path: root/tests/testthat/summary_nlmixr_saem_biphasic.txt
blob: 144cbac7e4b0788fd49ed9963a712d4d7cc7b4b0 (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
nlmixr version used for fitting:    2.0.6 
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:
507 observations of 2 variable(s) grouped in 15 datasets

Degradation model predictions using RxODE

Fitted in test time 0 s

Variance model: Two-component variance function 

Mean of starting values for individual parameters:
       parent_0        log_k_m1 f_parent_qlogis          log_k1          log_k2 
         100.65           -5.38           -0.09           -2.74           -4.53 
       g_qlogis 
          -0.14 

Mean of starting values for error model parameters:
sigma_low_parent  rsd_high_parent     sigma_low_m1      rsd_high_m1 
            0.79             0.05             0.79             0.05 

Fixed degradation parameter values:
None

Results:

Likelihood calculated by gauss3_1.6  
   AIC  BIC logLik
  2400 2468  -1184

Optimised parameters:
                  est. lower upper
parent_0        100.35  99.0 101.7
log_k_m1         -5.35  -5.5  -5.2
f_parent_qlogis  -0.09  -0.3   0.1
log_k1           -2.76  -3.0  -2.6
log_k2           -4.49  -4.6  -4.4
g_qlogis         -0.17  -0.4   0.1

Correlation: 
                pr_0 l__1 f_p_ lg_1 lg_2
log_k_m1        -0.4                    
f_parent_qlogis -0.4  0.4               
log_k1           0.2 -0.2 -0.1          
log_k2           0.0  0.1  0.0  0.3     
g_qlogis         0.1 -0.1  0.0 -0.4 -0.5

Random effects (omega):
                    eta.parent_0 eta.log_k_m1 eta.f_parent_qlogis eta.log_k1
eta.parent_0                0.08         0.00                 0.0       0.00
eta.log_k_m1                0.00         0.02                 0.0       0.00
eta.f_parent_qlogis         0.00         0.00                 0.1       0.00
eta.log_k1                  0.00         0.00                 0.0       0.09
eta.log_k2                  0.00         0.00                 0.0       0.00
eta.g_qlogis                0.00         0.00                 0.0       0.00
                    eta.log_k2 eta.g_qlogis
eta.parent_0              0.00          0.0
eta.log_k_m1              0.00          0.0
eta.f_parent_qlogis       0.00          0.0
eta.log_k1                0.00          0.0
eta.log_k2                0.03          0.0
eta.g_qlogis              0.00          0.1

Variance model:
sigma_low_parent  rsd_high_parent     sigma_low_m1      rsd_high_m1 
            1.04             0.05             0.82             0.06 

Backtransformed parameters:
                est. lower upper
parent_0       1e+02 1e+02 1e+02
k_m1           5e-03 4e-03 6e-03
f_parent_to_m1 5e-01 4e-01 5e-01
k1             6e-02 5e-02 8e-02
k2             1e-02 1e-02 1e-02
g              5e-01 4e-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   25  150       45      11      61
m1      145  483       NA      NA      NA

Contact - Imprint