ODE Simulation
para <- c(Qecf = 0.2/1000,
Qcsf = 2.2/1000,
Vbr = 1.73,
Vlv = 50/1000,
Vtfv = 50/1000,
Vcm = 17/1000,
Vsas = 180/1000,
Vcen = 60.34,
CL = 15.41,
Q = 24.27,
CLin = 212.92/1000,
CLout = 211.7/1000,
CLpl = 1.89/1000,
CLlp = 3.48/1000,
CLpt = 29.01/1000,
CLtp = 68.12/1000,
Vper = 151.86)
- Setting initial condition
yini <- c(
Acen = 150.251^6, # Quantity in plasma compartment (ng)
Aper = 0.0, # ~ peripheral compartment
Aecf = 0.0, # ~ ECF
Alv = 0.0, # ~
Atfv = 0.0, # ~
Acm = 0.0, # ~
Asas = 0.0
)
PBPKmodAPAP <- function(t,state,parameters){
with(as.list(c(state,parameters)), {
Ccen <- Acen/Vcen
Cper <- Aper/Vper
Cecf <- Aecf/Vbr
Clv <- Alv/Vlv
Ctfv <- Atfv/Vtfv
Ccm <- Acm/Vcm
Csas <- Asas/Vsas
#ODE equations
dAcen = Q * Cper - Q * Ccen - CL * Ccen - CLin * Ccen + CLout * Cecf + Clv * CLlp - Ccen * CLpl + Ctfv * CLtp - Ccen * CLpt + Qcsf * Csas + 0.02/100*200000
dAper = Q * Ccen - Q * Cper
dAecf = CLin * Ccen -CLout * Cecf - Qecf * Cecf
dAlv = Qecf * Cecf - Clv * Qcsf - CLlp * Clv + CLpl * Ccen
dAtfv = Qcsf* Clv - Ctfv * Qcsf + CLpt *Ccen - CLtp * Ctfv
dAcm = Qcsf* Ctfv - Ccm * Qcsf
dAsas = Qcsf* Ccm - Csas * Qcsf;
list(c(dAcen, dAper, dAecf, dAlv, dAtfv, dAcm, dAsas),
"Ccen" = Ccen, "Cecf" = Cecf, "Clv" = Clv, "Ccm" = Ccm)
})
}
times <- c(0.1, 0.5, 1, 2, 4, 8, 12, 24, 48, 72, 96)
results <- ode(y=yini, times=times, func=PBPKmodAPAP, parms=para, method="lsoda")
head(results)
## time Acen Aper Aecf Alv Atfv Acm
## [1,] 0.1 1.150547e+13 0.000000e+00 0 0 0 0
## [2,] 0.5 8.879472e+12 1.576019e+12 13933501036 124243352 1472383420 14624629
## [3,] 1.0 6.539568e+12 2.927921e+12 26153780765 234625052 2026748977 52707681
## [4,] 2.0 3.826278e+12 4.346747e+12 39706398208 360908360 1721310901 126284652
## [5,] 4.0 1.883943e+12 4.936903e+12 47326838556 442457934 844178986 191000237
## [6,] 8.0 1.149371e+12 4.236764e+12 44401361031 437331405 431092901 189749335
## Asas Ccen Cecf Clv Ccm
## [1,] 0.0 190677268107 0 0 0
## [2,] 270876.3 147157313715 8054046842 2484867031 860272296
## [3,] 2391760.0 108378662368 15117792350 4692501037 3100451834
## [4,] 14146777.3 63411960897 22951675265 7218167204 7428508917
## [5,] 56408407.5 31222127295 27356554079 8849158678 11235308072
## [6,] 152596627.6 19048236042 25665526608 8746628104 11161725582
plot(results, select=c("Ccen", "Cecf" ,"Clv", "Ccm"))
Global Sensitivity Analysis
params <- c("Qecf", "Qcsf", "Vbr","Vlv", "Vtfv", "Vcm", "Vsas",
"Vcen", "Vper", "CL","Q","CLin","CLout","CLlp","CLpl","CLtp","CLpt")
omega_para <- data.frame(
omega_Vcen = 16.9,
omega_Vper = 4.14,
omega_CL = 4.58,
omega_Q = 15,
omega_CLin = 23,
omega_CLout = 5.67,
omega_CLlp = 15.6,
omega_CLpl = 17.9,
omega_CLtp = 23.4,
omega_CLpt = 7.76
)
- Calculate SD by omega NOTE: The given approach provide the extremely
large sd, which cannot be used in GSA (caused comuting error)
#sd_Vcen = sqrt(exp(omega_para$omega_Vcen^2)-1) # Extremely large value (1.045739e+62)
#sd_CL = sqrt(exp(omega_para$omega_CL^2)-1)
#sd_Vper = sqrt(exp(omega_para$omega_Vper^2)-1)
#sd_Q = sqrt(exp(omega_para$omega_Q^2)-1)
#sd_CLin = sqrt(exp(omega_para$omega_CLin^2)-1)
#sd_CLout = sqrt(exp(omega_para$omega_CLout^2)-1)
#sd_CLlp = sqrt(exp(omega_para$omega_CLlp^2)-1)
#sd_CLpl = sqrt(exp(omega_para$omega_CLpl^2)-1)
#sd_CLtp = sqrt(exp(omega_para$omega_CLtp^2)-1)
#sd_CLpt = sqrt(exp(omega_para$omega_CLpt^2)-1)
# The sd need to be checked
sd_Vcen = sqrt(log(omega_para$omega_Vcen^2)-1)
sd_Vper = sqrt(log(omega_para$omega_Vper^2)-1)
sd_CL = sqrt(log(omega_para$omega_CL^2)-1)
sd_Q = sqrt(log(omega_para$omega_Q^2)-1)
sd_CLin = sqrt(log(omega_para$omega_CLin^2)-1)
sd_CLout = sqrt(log(omega_para$omega_CLout^2)-1)
sd_CLlp = sqrt(log(omega_para$omega_CLlp^2)-1)
sd_CLpl = sqrt(log(omega_para$omega_CLpl^2)-1)
sd_CLtp = sqrt(log(omega_para$omega_CLtp^2)-1)
sd_CLpt = sqrt(log(omega_para$omega_CLpt^2)-1)
LL <- 0.1 # 10% lower limit
UL <- 1.1
q <- c(rep("qunif", 7), rep("qlnorm", 10))
q.arg <- list(
list(min = para["Qecf"]*LL, max= para["Qecf"]*UL), ##uniform with bounds corresponding to LL and UL of the nominal value
list(min = para["Qcsf"]*LL, max= para["Qcsf"]*UL),
list(min = para["Vbr"]*LL, max = para["Vbr"]*UL),
list(min = para["Vlv"]*LL, max = para["Vlv"]*UL),
list(min = para["Vtfv"]*LL, max = para["Vtfv"]*UL),
list(min = para["Vcm"]*LL, max = para["Vcm"]*UL),
list(min = para["Vsas"]*LL, max = para["Vsas"]*UL),
list(meanlog = log(max(para["Vcen"], 1e-10)), sdlog = log(sd_Vcen)), ##lognormal distribution with CV from Monolix
list(meanlog = log(max(para["Vper"], 1e-10)), sdlog = log(sd_Vper)),
list(meanlog = log(max(para["CL"], 1e-10)), sdlog = log(sd_CL)),
list(meanlog = log(max(para["Q"], 1e-10)), sdlog = log(sd_Q)),
list(meanlog = log(max(para["CLin"], 1e-10)), sdlog = log(sd_CLin)),
list(meanlog = log(max(para["CLout"], 1e-10)), sdlog = log(sd_CLout)),
list(meanlog = log(max(para["CLpl"], 1e-10)), sdlog = log(sd_CLpl)),
list(meanlog = log(max(para["CLlp"], 1e-10)), sdlog = log(sd_CLlp)),
list(meanlog = log(max(para["CLpt"], 1e-10)), sdlog = log(sd_CLpt)),
list(meanlog = log(max(para["CLtp"], 1e-10)), sdlog = log(sd_CLtp)))
- Pre-examination NOTE: The n = 20 generate th NA due to lack of
sample size. Increase the sample number can solve the issue.
length(params)==length(q)
## [1] TRUE
set.seed(1234)
x <- rfast99(params = params, n = 20, q = q, q.arg = q.arg, rep = 1) #
anyNA(x$a)
## [1] TRUE
x <- rfast99(params = params, n = 100, q = q, q.arg = q.arg, rep = 1) # replace n to 100
anyNA(x$a)
## [1] FALSE
outputs = c("Ccen","Cecf","Clv","Ccm")
out <- solve_fun(x, time=times, func = PBPKmodAPAP, initState = yini, outnames = outputs,
method="lsoda")
## Starting time: 2023-11-07 14:30:38
## Ending time: 2023-11-07 14:31:25
plot(out, vars=outputs[1])
out
##
## Call:
## rfast99(params = params, n = 100, q = q, q.arg = q.arg, replicate = 1)
##
## Model runs: 1700
##
##
## ==================================
## Sensitivity Indices
##
## first order:
## , , Ccen
##
## Qecf Qcsf Vbr Vlv Vtfv Vcm Vsas Vcen Vper CL
## 0.1 0.0023 0.0033 0.0017 0.0024 0.0014 0.0041 0.0015 0.9211 0.0076 0.0024
## 0.5 0.0023 0.0012 0.0008 0.0007 0.0044 0.0009 0.0105 0.7108 0.0048 0.0250
## 1 0.0021 0.0026 0.0009 0.0007 0.0030 0.0003 0.0040 0.5193 0.0043 0.1128
## 2 0.0022 0.0016 0.0016 0.0011 0.0114 0.0034 0.0029 0.4268 0.0445 0.3396
## 4 0.0068 0.0033 0.0022 0.0017 0.0060 0.0024 0.0031 0.1340 0.1880 0.5246
## 8 0.0034 0.0065 0.0014 0.0018 0.0003 0.0015 0.0024 0.0140 0.0881 0.3762
## 12 0.0007 0.0029 0.0009 0.0017 0.0004 0.0014 0.0007 0.0078 0.0519 0.4539
## 24 0.0001 0.0014 0.0007 0.0008 0.0019 0.0030 0.0002 0.0279 0.0043 0.6620
## 48 0.0008 0.0014 0.0017 0.0010 0.0066 0.0086 0.0027 0.0793 0.0626 0.6776
## 72 0.0027 0.0034 0.0037 0.0019 0.0150 0.0121 0.0082 0.1149 0.1672 0.6196
## 96 0.0055 0.0072 0.0060 0.0032 0.0260 0.0148 0.0166 0.1423 0.2536 0.5747
## Q CLin CLout CLlp CLpl CLtp CLpt
## 0.1 0.0138 0.0018 0.0035 0.0030 0.0017 0.0022 0.0007
## 0.5 0.1601 0.0049 0.0032 0.0055 0.0008 0.0022 0.0115
## 1 0.3318 0.0028 0.0017 0.0116 0.0012 0.0017 0.0031
## 2 0.4537 0.0023 0.0026 0.0064 0.0034 0.0017 0.0025
## 4 0.1231 0.0019 0.0055 0.0030 0.0054 0.0048 0.0033
## 8 0.2528 0.0010 0.0034 0.0003 0.0038 0.0029 0.0017
## 12 0.2884 0.0005 0.0009 0.0001 0.0022 0.0008 0.0008
## 24 0.1342 0.0007 0.0035 0.0005 0.0029 0.0006 0.0008
## 48 0.0122 0.0025 0.0103 0.0050 0.0039 0.0014 0.0011
## 72 0.0183 0.0044 0.0109 0.0139 0.0052 0.0026 0.0030
## 96 0.0529 0.0065 0.0109 0.0253 0.0077 0.0047 0.0061
##
## , , Cecf
##
## Qecf Qcsf Vbr Vlv Vtfv Vcm Vsas Vcen Vper CL
## 0.1 0.0101 0.0023 0.2268 0.0026 0.0063 0.0060 0.0094 0.1760 0.0197 0.0046
## 0.5 0.0095 0.0023 0.1823 0.0025 0.0056 0.0050 0.0066 0.1317 0.0073 0.0031
## 1 0.0093 0.0018 0.1462 0.0040 0.0072 0.0041 0.0069 0.0970 0.0232 0.0090
## 2 0.0089 0.0005 0.1016 0.0051 0.0184 0.0031 0.0074 0.0601 0.0256 0.0251
## 4 0.0071 0.0008 0.0500 0.0048 0.0216 0.0023 0.0085 0.0339 0.0250 0.0201
## 8 0.0057 0.0014 0.0081 0.0046 0.0199 0.0020 0.0126 0.0178 0.0234 0.0132
## 12 0.0053 0.0024 0.0010 0.0041 0.0182 0.0021 0.0176 0.0124 0.0223 0.0141
## 24 0.0045 0.0066 0.0018 0.0029 0.0152 0.0024 0.0031 0.0126 0.0212 0.0458
## 48 0.0039 0.0156 0.0019 0.0051 0.0174 0.0043 0.0040 0.0295 0.0197 0.1718
## 72 0.0081 0.0228 0.0034 0.0101 0.0249 0.0086 0.0110 0.0598 0.0187 0.2495
## 96 0.0195 0.0285 0.0054 0.0146 0.0345 0.0139 0.0211 0.0959 0.0220 0.2787
## Q CLin CLout CLlp CLpl CLtp CLpt
## 0.1 0.0281 0.4050 0.0030 0.0077 0.0105 0.0173 0.0040
## 0.5 0.0414 0.4964 0.0100 0.0053 0.0112 0.0125 0.0041
## 1 0.0507 0.4852 0.0231 0.0056 0.0103 0.0107 0.0041
## 2 0.0566 0.4859 0.0572 0.0075 0.0078 0.0097 0.0045
## 4 0.0475 0.5459 0.1062 0.0076 0.0064 0.0087 0.0056
## 8 0.0202 0.6064 0.1545 0.0059 0.0080 0.0082 0.0068
## 12 0.0160 0.6000 0.1941 0.0040 0.0077 0.0081 0.0074
## 24 0.0275 0.5417 0.2455 0.0011 0.0044 0.0081 0.0077
## 48 0.0272 0.3833 0.2002 0.0020 0.0021 0.0096 0.0076
## 72 0.0317 0.2705 0.1541 0.0046 0.0022 0.0126 0.0084
## 96 0.0379 0.2035 0.1255 0.0081 0.0024 0.0161 0.0101
##
## , , Clv
##
## Qecf Qcsf Vbr Vlv Vtfv Vcm Vsas Vcen Vper CL
## 0.1 0.0132 0.0065 0.0046 0.2256 0.0064 0.0030 0.0088 0.1254 0.0147 0.0052
## 0.5 0.0104 0.0060 0.0054 0.2183 0.0037 0.0017 0.0051 0.0908 0.0157 0.0105
## 1 0.0083 0.0067 0.0053 0.2102 0.0029 0.0009 0.0091 0.0640 0.0216 0.0165
## 2 0.0054 0.0113 0.0047 0.1955 0.0025 0.0004 0.0115 0.0332 0.0355 0.0253
## 4 0.0027 0.0256 0.0040 0.1469 0.0027 0.0010 0.0107 0.0109 0.0486 0.0361
## 8 0.0016 0.0557 0.0044 0.0623 0.0031 0.0022 0.0098 0.0024 0.0501 0.0480
## 12 0.0015 0.0840 0.0055 0.0354 0.0034 0.0026 0.0097 0.0010 0.0474 0.0572
## 24 0.0014 0.1602 0.0066 0.0677 0.0036 0.0028 0.0106 0.0055 0.0380 0.0849
## 48 0.0018 0.1806 0.0039 0.1727 0.0039 0.0026 0.0131 0.0299 0.0192 0.1356
## 72 0.0030 0.1447 0.0030 0.1864 0.0086 0.0023 0.0158 0.0602 0.0213 0.1677
## 96 0.0048 0.1167 0.0036 0.1322 0.0191 0.0023 0.0185 0.0920 0.0518 0.1811
## Q CLin CLout CLlp CLpl CLtp CLpt
## 0.1 0.0062 0.0079 0.0080 0.0064 0.2527 0.0078 0.0042
## 0.5 0.0098 0.0084 0.0073 0.0041 0.2591 0.0056 0.0016
## 1 0.0159 0.0092 0.0085 0.0027 0.2837 0.0047 0.0023
## 2 0.0214 0.0102 0.0102 0.0026 0.3601 0.0038 0.0038
## 4 0.0137 0.0112 0.0101 0.0222 0.5008 0.0032 0.0058
## 8 0.0101 0.0124 0.0078 0.0853 0.5875 0.0033 0.0063
## 12 0.0298 0.0134 0.0065 0.1292 0.5616 0.0035 0.0040
## 24 0.0678 0.0150 0.0074 0.1797 0.3983 0.0033 0.0010
## 48 0.0199 0.0140 0.0153 0.2106 0.2359 0.0024 0.0009
## 72 0.0129 0.0111 0.0241 0.2296 0.1793 0.0021 0.0014
## 96 0.0302 0.0095 0.0326 0.2432 0.1470 0.0026 0.0028
##
## , , Ccm
##
## Qecf Qcsf Vbr Vlv Vtfv Vcm Vsas Vcen Vper CL
## 0.1 0.0161 0.1876 0.0208 0.0027 0.1907 0.1001 0.0103 0.3311 0.0119 0.0021
## 0.5 0.0146 0.1824 0.0219 0.0071 0.1441 0.0911 0.0066 0.2669 0.0054 0.0022
## 1 0.0136 0.1700 0.0213 0.0075 0.1060 0.0860 0.0062 0.2104 0.0039 0.0064
## 2 0.0135 0.1374 0.0203 0.0059 0.0632 0.0766 0.0068 0.1384 0.0030 0.0210
## 4 0.0154 0.0837 0.0192 0.0051 0.0274 0.0487 0.0077 0.0638 0.0025 0.0526
## 8 0.0201 0.0364 0.0186 0.0048 0.0092 0.0106 0.0093 0.0139 0.0072 0.1167
## 12 0.0238 0.0178 0.0187 0.0047 0.0051 0.0135 0.0102 0.0097 0.0113 0.1968
## 24 0.0299 0.0040 0.0189 0.0043 0.0030 0.0453 0.0108 0.0114 0.0110 0.4348
## 48 0.0337 0.0170 0.0162 0.0030 0.0048 0.0706 0.0110 0.0387 0.0170 0.7395
## 72 0.0335 0.0410 0.0129 0.0028 0.0071 0.0734 0.0109 0.0652 0.0555 0.7219
## 96 0.0270 0.0709 0.0125 0.0033 0.0193 0.0763 0.0114 0.0815 0.1017 0.6327
## Q CLin CLout CLlp CLpl CLtp CLpt
## 0.1 0.0027 0.0020 0.0037 0.0128 0.0022 0.0036 0.3191
## 0.5 0.0115 0.0033 0.0049 0.0194 0.0021 0.0207 0.3479
## 1 0.0251 0.0044 0.0064 0.0198 0.0026 0.0461 0.3630
## 2 0.0417 0.0045 0.0085 0.0176 0.0039 0.1019 0.3779
## 4 0.0483 0.0055 0.0112 0.0166 0.0054 0.1933 0.4122
## 8 0.0288 0.0108 0.0133 0.0172 0.0064 0.2423 0.5045
## 12 0.0146 0.0137 0.0133 0.0173 0.0069 0.2815 0.6010
## 24 0.0055 0.0111 0.0094 0.0154 0.0067 0.3106 0.8293
## 48 0.0085 0.0068 0.0049 0.0121 0.0049 0.3320 0.4954
## 72 0.0194 0.0041 0.0045 0.0114 0.0049 0.3113 0.2506
## 96 0.0353 0.0030 0.0059 0.0187 0.0058 0.2362 0.1655
##
##
## interaction:
## , , Ccen
##
## Qecf Qcsf Vbr Vlv Vtfv Vcm Vsas Vcen Vper CL
## 0.1 0.0288 0.0488 0.0278 0.0330 0.0264 0.0619 0.0297 0.0627 0.0966 0.0447
## 0.5 0.0405 0.0116 0.0126 0.0229 0.0403 0.0202 0.0776 0.1355 0.0430 0.0451
## 1 0.0201 0.0241 0.0131 0.0196 0.0297 0.0060 0.0272 0.1471 0.0498 0.1821
## 2 0.0089 0.0103 0.0244 0.0234 0.1364 0.0401 0.0218 0.2486 0.0989 0.3160
## 4 0.0321 0.0278 0.0240 0.0298 0.0955 0.0405 0.0294 0.3048 0.2160 0.1509
## 8 0.0374 0.0700 0.0117 0.0251 0.0153 0.0284 0.0331 0.1108 0.0626 0.0715
## 12 0.0225 0.0506 0.0120 0.0189 0.0113 0.0312 0.0221 0.0528 0.0427 0.0840
## 24 0.0208 0.0245 0.0455 0.0211 0.0271 0.0563 0.0128 0.0244 0.0394 0.0797
## 48 0.0323 0.0654 0.0967 0.0425 0.0786 0.1118 0.0227 0.0654 0.0814 0.0925
## 72 0.0548 0.1430 0.1290 0.0704 0.1550 0.1554 0.0493 0.1439 0.1514 0.1471
## 96 0.0868 0.2318 0.1517 0.1030 0.2502 0.2010 0.0895 0.2317 0.2414 0.2068
## Q CLin CLout CLlp CLpl CLtp CLpt
## 0.1 0.0676 0.0266 0.0386 0.0516 0.0279 0.0291 0.0268
## 0.5 0.1583 0.0638 0.0448 0.0261 0.0114 0.0405 0.0517
## 1 0.1873 0.0276 0.0104 0.0543 0.0157 0.0182 0.0194
## 2 0.1541 0.0284 0.0122 0.0393 0.0386 0.0097 0.0330
## 4 0.4170 0.0160 0.0199 0.0428 0.0816 0.0311 0.0574
## 8 0.1891 0.0073 0.0166 0.0142 0.1090 0.0393 0.0639
## 12 0.0969 0.0065 0.0155 0.0112 0.0965 0.0245 0.0443
## 24 0.0873 0.0214 0.0550 0.0178 0.0561 0.0252 0.0200
## 48 0.0898 0.0601 0.1214 0.0565 0.0427 0.0534 0.0155
## 72 0.0975 0.0924 0.1456 0.1147 0.0595 0.0982 0.0209
## 96 0.1329 0.1310 0.1736 0.1848 0.0935 0.1541 0.0338
##
## , , Cecf
##
## Qecf Qcsf Vbr Vlv Vtfv Vcm Vsas Vcen Vper CL
## 0.1 0.0910 0.0321 0.4097 0.0530 0.1346 0.1524 0.1345 0.2911 0.2887 0.1651
## 0.5 0.1017 0.0253 0.4021 0.0467 0.2591 0.1356 0.1243 0.2367 0.1957 0.2103
## 1 0.1227 0.0278 0.3666 0.0585 0.4467 0.1188 0.1273 0.1939 0.3119 0.3431
## 2 0.1321 0.1180 0.2759 0.0780 0.4272 0.0942 0.1271 0.1611 0.3266 0.5836
## 4 0.1111 0.0806 0.1478 0.0920 0.2989 0.0646 0.1345 0.1485 0.3228 0.4429
## 8 0.0888 0.0794 0.0950 0.0993 0.2372 0.0388 0.1884 0.1318 0.3050 0.2721
## 12 0.0827 0.0939 0.1037 0.0951 0.2147 0.0307 0.3217 0.1128 0.2838 0.2020
## 24 0.0831 0.1330 0.1183 0.0676 0.1718 0.0365 0.4377 0.0989 0.2433 0.0902
## 48 0.1157 0.2322 0.1336 0.0854 0.1391 0.0825 0.0901 0.1199 0.1977 0.1500
## 72 0.1981 0.3433 0.1591 0.1268 0.1601 0.1501 0.0824 0.1693 0.1837 0.2746
## 96 0.2940 0.4437 0.1889 0.1521 0.2041 0.2321 0.1191 0.2400 0.2263 0.3649
## Q CLin CLout CLlp CLpl CLtp CLpt
## 0.1 0.2057 0.3393 0.0426 0.0465 0.1794 0.1493 0.0410
## 0.5 0.2474 0.2159 0.0609 0.0451 0.1641 0.1468 0.0335
## 1 0.2741 0.2040 0.0828 0.0627 0.1423 0.1442 0.0324
## 2 0.2931 0.1919 0.1083 0.1004 0.0988 0.1361 0.0427
## 4 0.2792 0.1575 0.1088 0.1233 0.0713 0.1236 0.0678
## 8 0.2198 0.1477 0.1509 0.1157 0.0946 0.1126 0.0822
## 12 0.1976 0.1681 0.2018 0.0880 0.0985 0.1095 0.0856
## 24 0.2168 0.2042 0.2669 0.0400 0.0719 0.1116 0.0883
## 48 0.2623 0.2865 0.3183 0.0380 0.0487 0.1337 0.0881
## 72 0.3006 0.3843 0.3861 0.0678 0.0540 0.1674 0.0940
## 96 0.3343 0.4700 0.4549 0.1068 0.0711 0.2068 0.1041
##
## , , Clv
##
## Qecf Qcsf Vbr Vlv Vtfv Vcm Vsas Vcen Vper CL
## 0.1 0.1431 0.0760 0.0782 0.2838 0.0324 0.0729 0.1200 0.2668 0.1305 0.1036
## 0.5 0.1277 0.0571 0.0773 0.2817 0.0270 0.0555 0.1541 0.2001 0.1211 0.0947
## 1 0.1114 0.0527 0.0766 0.2582 0.0296 0.0427 0.2261 0.1297 0.1810 0.0947
## 2 0.0840 0.0635 0.0768 0.2131 0.0330 0.0323 0.2029 0.0531 0.3179 0.1003
## 4 0.0525 0.1026 0.0810 0.1788 0.0468 0.0311 0.1663 0.0314 0.4011 0.1090
## 8 0.0410 0.1576 0.0925 0.2068 0.0571 0.0412 0.1514 0.0462 0.3955 0.1128
## 12 0.0372 0.1787 0.0985 0.2281 0.0581 0.0462 0.1464 0.0458 0.3859 0.1140
## 24 0.0226 0.1828 0.0946 0.2380 0.0569 0.0440 0.1457 0.0470 0.3478 0.1407
## 48 0.0202 0.2294 0.0767 0.2195 0.0579 0.0403 0.1646 0.0858 0.2314 0.2407
## 72 0.0333 0.2562 0.0799 0.1625 0.1028 0.0463 0.1900 0.1496 0.1398 0.3440
## 96 0.0514 0.2514 0.0955 0.1370 0.2041 0.0612 0.2173 0.2232 0.1498 0.4312
## Q CLin CLout CLlp CLpl CLtp CLpt
## 0.1 0.0421 0.1112 0.0563 0.1052 0.4122 0.0625 0.0675
## 0.5 0.0434 0.1147 0.0842 0.0833 0.4395 0.0518 0.0332
## 1 0.0525 0.1202 0.1245 0.0803 0.4369 0.0507 0.0471
## 2 0.0609 0.1275 0.1613 0.0743 0.4015 0.0488 0.0714
## 4 0.0726 0.1331 0.1747 0.0671 0.3309 0.0495 0.0925
## 8 0.1572 0.1404 0.1506 0.0579 0.2881 0.0530 0.1008
## 12 0.1686 0.1491 0.1257 0.0548 0.2926 0.0512 0.0757
## 24 0.1247 0.1719 0.0935 0.0866 0.3321 0.0458 0.0228
## 48 0.1116 0.1892 0.1317 0.1891 0.3713 0.0412 0.0113
## 72 0.1203 0.1874 0.2270 0.2834 0.3766 0.0468 0.0112
## 96 0.1486 0.1980 0.3088 0.3654 0.3714 0.0665 0.0165
##
## , , Ccm
##
## Qecf Qcsf Vbr Vlv Vtfv Vcm Vsas Vcen Vper CL
## 0.1 0.1694 0.2098 0.1316 0.1064 0.3517 0.2075 0.1544 0.2922 0.1312 0.0508
## 0.5 0.1408 0.2044 0.1748 0.1696 0.3413 0.2048 0.1674 0.2718 0.0704 0.0402
## 1 0.1341 0.1967 0.1951 0.1728 0.3135 0.1958 0.1569 0.2363 0.0559 0.0464
## 2 0.1412 0.1815 0.2058 0.1313 0.2507 0.1659 0.1382 0.1634 0.0504 0.0799
## 4 0.1709 0.1574 0.2089 0.1011 0.1678 0.0964 0.1239 0.0679 0.0531 0.1330
## 8 0.2280 0.1235 0.2113 0.0848 0.1165 0.0825 0.1079 0.0705 0.0999 0.1940
## 12 0.2689 0.1002 0.2147 0.0788 0.1022 0.1107 0.0902 0.1601 0.1246 0.2416
## 24 0.3327 0.0787 0.2169 0.0719 0.0868 0.1509 0.0578 0.0811 0.1688 0.2660
## 48 0.3754 0.0904 0.1960 0.0633 0.0825 0.1732 0.0470 0.1638 0.2381 0.1639
## 72 0.3828 0.1192 0.1762 0.0575 0.1326 0.1865 0.0488 0.2722 0.2689 0.1292
## 96 0.3301 0.1865 0.1843 0.0613 0.2442 0.2076 0.0575 0.3457 0.2946 0.1763
## Q CLin CLout CLlp CLpl CLtp CLpt
## 0.1 0.0490 0.1074 0.0980 0.2301 0.1445 0.0258 0.3054
## 0.5 0.0752 0.0988 0.0889 0.2285 0.1302 0.0757 0.2971
## 1 0.1322 0.0919 0.0843 0.1868 0.1268 0.1199 0.2908
## 2 0.1848 0.0808 0.0818 0.1550 0.1169 0.1586 0.2796
## 4 0.1652 0.0978 0.0811 0.1518 0.0952 0.2253 0.2549
## 8 0.1030 0.1495 0.0773 0.1799 0.0822 0.2919 0.2061
## 12 0.0770 0.1633 0.0738 0.1954 0.0809 0.2655 0.1656
## 24 0.0586 0.1266 0.0739 0.1785 0.0749 0.3307 0.0976
## 48 0.0822 0.0961 0.0907 0.1493 0.0621 0.4277 0.2563
## 72 0.1371 0.1235 0.1109 0.1437 0.0567 0.3765 0.3719
## 96 0.2094 0.1905 0.1422 0.1873 0.0562 0.3233 0.4211
##
##
## total order:
## , , Ccen
##
## Qecf Qcsf Vbr Vlv Vtfv Vcm Vsas Vcen Vper CL
## 0.1 0.0311 0.0521 0.0295 0.0355 0.0278 0.0660 0.0311 0.9838 0.1042 0.0470
## 0.5 0.0427 0.0129 0.0133 0.0236 0.0447 0.0211 0.0880 0.8464 0.0478 0.0701
## 1 0.0222 0.0267 0.0140 0.0203 0.0326 0.0063 0.0312 0.6663 0.0541 0.2949
## 2 0.0111 0.0119 0.0260 0.0246 0.1478 0.0434 0.0248 0.6754 0.1434 0.6556
## 4 0.0389 0.0311 0.0263 0.0315 0.1015 0.0429 0.0326 0.4388 0.4039 0.6755
## 8 0.0408 0.0764 0.0132 0.0270 0.0155 0.0299 0.0355 0.1249 0.1507 0.4477
## 12 0.0233 0.0535 0.0129 0.0206 0.0117 0.0326 0.0228 0.0606 0.0946 0.5379
## 24 0.0209 0.0260 0.0462 0.0219 0.0290 0.0593 0.0129 0.0522 0.0437 0.7417
## 48 0.0332 0.0667 0.0984 0.0435 0.0852 0.1204 0.0253 0.1447 0.1440 0.7700
## 72 0.0575 0.1464 0.1327 0.0724 0.1700 0.1675 0.0574 0.2588 0.3186 0.7667
## 96 0.0923 0.2391 0.1577 0.1062 0.2762 0.2158 0.1062 0.3740 0.4949 0.7815
## Q CLin CLout CLlp CLpl CLtp CLpt
## 0.1 0.0813 0.0284 0.0421 0.0546 0.0296 0.0313 0.0275
## 0.5 0.3184 0.0687 0.0479 0.0317 0.0122 0.0427 0.0632
## 1 0.5191 0.0304 0.0121 0.0659 0.0169 0.0199 0.0225
## 2 0.6079 0.0307 0.0149 0.0457 0.0421 0.0114 0.0355
## 4 0.5401 0.0179 0.0254 0.0459 0.0870 0.0359 0.0607
## 8 0.4419 0.0083 0.0201 0.0145 0.1127 0.0422 0.0656
## 12 0.3853 0.0070 0.0164 0.0114 0.0988 0.0253 0.0451
## 24 0.2215 0.0221 0.0585 0.0183 0.0590 0.0258 0.0209
## 48 0.1020 0.0626 0.1316 0.0616 0.0466 0.0548 0.0166
## 72 0.1157 0.0968 0.1565 0.1286 0.0647 0.1008 0.0239
## 96 0.1858 0.1375 0.1844 0.2101 0.1012 0.1588 0.0399
##
## , , Cecf
##
## Qecf Qcsf Vbr Vlv Vtfv Vcm Vsas Vcen Vper CL
## 0.1 0.1011 0.0344 0.6365 0.0556 0.1409 0.1584 0.1439 0.4672 0.3084 0.1697
## 0.5 0.1112 0.0276 0.5844 0.0492 0.2647 0.1406 0.1309 0.3684 0.2030 0.2134
## 1 0.1320 0.0296 0.5128 0.0625 0.4539 0.1229 0.1341 0.2910 0.3351 0.3521
## 2 0.1410 0.1185 0.3775 0.0830 0.4456 0.0973 0.1345 0.2212 0.3522 0.6087
## 4 0.1182 0.0815 0.1978 0.0969 0.3205 0.0669 0.1430 0.1824 0.3478 0.4630
## 8 0.0945 0.0808 0.1031 0.1039 0.2571 0.0408 0.2010 0.1495 0.3285 0.2854
## 12 0.0881 0.0963 0.1046 0.0991 0.2329 0.0328 0.3394 0.1252 0.3061 0.2161
## 24 0.0876 0.1396 0.1201 0.0705 0.1870 0.0389 0.4408 0.1115 0.2645 0.1361
## 48 0.1195 0.2478 0.1355 0.0905 0.1566 0.0868 0.0941 0.1494 0.2175 0.3218
## 72 0.2062 0.3660 0.1625 0.1369 0.1851 0.1587 0.0934 0.2292 0.2024 0.5241
## 96 0.3135 0.4722 0.1943 0.1667 0.2386 0.2460 0.1403 0.3359 0.2482 0.6436
## Q CLin CLout CLlp CLpl CLtp CLpt
## 0.1 0.2338 0.7443 0.0456 0.0542 0.1899 0.1667 0.0450
## 0.5 0.2888 0.7123 0.0709 0.0505 0.1753 0.1593 0.0376
## 1 0.3248 0.6892 0.1059 0.0683 0.1526 0.1549 0.0365
## 2 0.3497 0.6778 0.1655 0.1080 0.1066 0.1457 0.0472
## 4 0.3267 0.7034 0.2151 0.1309 0.0777 0.1323 0.0734
## 8 0.2400 0.7541 0.3054 0.1217 0.1026 0.1208 0.0890
## 12 0.2136 0.7681 0.3959 0.0920 0.1062 0.1175 0.0930
## 24 0.2444 0.7459 0.5123 0.0411 0.0763 0.1197 0.0960
## 48 0.2895 0.6698 0.5185 0.0401 0.0509 0.1433 0.0957
## 72 0.3324 0.6547 0.5402 0.0724 0.0562 0.1800 0.1025
## 96 0.3722 0.6735 0.5804 0.1149 0.0735 0.2229 0.1142
##
## , , Clv
##
## Qecf Qcsf Vbr Vlv Vtfv Vcm Vsas Vcen Vper CL
## 0.1 0.1563 0.0825 0.0828 0.5094 0.0388 0.0759 0.1288 0.3921 0.1452 0.1087
## 0.5 0.1381 0.0630 0.0827 0.5000 0.0307 0.0572 0.1593 0.2909 0.1367 0.1051
## 1 0.1196 0.0594 0.0818 0.4684 0.0325 0.0436 0.2352 0.1937 0.2026 0.1112
## 2 0.0894 0.0748 0.0815 0.4085 0.0355 0.0327 0.2144 0.0863 0.3534 0.1256
## 4 0.0552 0.1282 0.0850 0.3256 0.0495 0.0321 0.1770 0.0424 0.4497 0.1451
## 8 0.0427 0.2133 0.0969 0.2691 0.0602 0.0434 0.1612 0.0486 0.4457 0.1608
## 12 0.0387 0.2627 0.1039 0.2635 0.0615 0.0488 0.1561 0.0468 0.4333 0.1711
## 24 0.0241 0.3430 0.1011 0.3057 0.0605 0.0469 0.1563 0.0525 0.3858 0.2256
## 48 0.0221 0.4101 0.0806 0.3921 0.0618 0.0429 0.1776 0.1158 0.2506 0.3763
## 72 0.0364 0.4010 0.0828 0.3489 0.1114 0.0487 0.2058 0.2098 0.1611 0.5117
## 96 0.0563 0.3681 0.0991 0.2692 0.2232 0.0634 0.2358 0.3152 0.2016 0.6124
## Q CLin CLout CLlp CLpl CLtp CLpt
## 0.1 0.0484 0.1191 0.0643 0.1116 0.6650 0.0703 0.0717
## 0.5 0.0532 0.1231 0.0915 0.0874 0.6986 0.0574 0.0348
## 1 0.0684 0.1294 0.1330 0.0830 0.7206 0.0554 0.0494
## 2 0.0824 0.1377 0.1714 0.0770 0.7617 0.0526 0.0752
## 4 0.0863 0.1443 0.1848 0.0893 0.8317 0.0527 0.0984
## 8 0.1674 0.1528 0.1584 0.1431 0.8756 0.0563 0.1071
## 12 0.1983 0.1626 0.1323 0.1839 0.8542 0.0547 0.0796
## 24 0.1926 0.1870 0.1009 0.2663 0.7304 0.0491 0.0238
## 48 0.1315 0.2032 0.1470 0.3996 0.6071 0.0436 0.0122
## 72 0.1333 0.1985 0.2511 0.5130 0.5559 0.0489 0.0126
## 96 0.1788 0.2075 0.3414 0.6087 0.5184 0.0691 0.0193
##
## , , Ccm
##
## Qecf Qcsf Vbr Vlv Vtfv Vcm Vsas Vcen Vper CL
## 0.1 0.1855 0.3974 0.1523 0.1090 0.5424 0.3076 0.1647 0.6233 0.1431 0.0530
## 0.5 0.1553 0.3868 0.1966 0.1767 0.4854 0.2960 0.1740 0.5387 0.0757 0.0424
## 1 0.1477 0.3667 0.2163 0.1803 0.4195 0.2818 0.1631 0.4467 0.0597 0.0528
## 2 0.1547 0.3189 0.2260 0.1372 0.3139 0.2425 0.1450 0.3018 0.0533 0.1010
## 4 0.1863 0.2411 0.2281 0.1062 0.1952 0.1451 0.1317 0.1317 0.0556 0.1856
## 8 0.2481 0.1599 0.2299 0.0896 0.1257 0.0931 0.1172 0.0844 0.1071 0.3107
## 12 0.2927 0.1180 0.2334 0.0835 0.1073 0.1242 0.1004 0.1698 0.1359 0.4384
## 24 0.3626 0.0827 0.2358 0.0762 0.0899 0.1962 0.0686 0.0925 0.1798 0.7008
## 48 0.4091 0.1073 0.2121 0.0663 0.0873 0.2438 0.0580 0.2025 0.2552 0.9034
## 72 0.4162 0.1602 0.1891 0.0603 0.1397 0.2599 0.0597 0.3374 0.3244 0.8510
## 96 0.3571 0.2574 0.1967 0.0645 0.2634 0.2839 0.0688 0.4272 0.3963 0.8089
## Q CLin CLout CLlp CLpl CLtp CLpt
## 0.1 0.0517 0.1093 0.1017 0.2429 0.1467 0.0294 0.6245
## 0.5 0.0867 0.1021 0.0939 0.2478 0.1323 0.0964 0.6450
## 1 0.1573 0.0963 0.0906 0.2066 0.1294 0.1659 0.6538
## 2 0.2266 0.0853 0.0903 0.1726 0.1208 0.2605 0.6575
## 4 0.2135 0.1033 0.0923 0.1684 0.1006 0.4186 0.6670
## 8 0.1319 0.1603 0.0906 0.1971 0.0886 0.5342 0.7105
## 12 0.0916 0.1770 0.0870 0.2127 0.0878 0.5469 0.7667
## 24 0.0642 0.1377 0.0833 0.1939 0.0815 0.6413 0.9269
## 48 0.0907 0.1029 0.0956 0.1613 0.0670 0.7597 0.7517
## 72 0.1565 0.1277 0.1154 0.1551 0.0616 0.6878 0.6225
## 96 0.2448 0.1934 0.1481 0.2059 0.0619 0.5594 0.5866