diff --git a/docs/articles/index.html b/docs/articles/index.html index 9d4be6c..f792288 100644 --- a/docs/articles/index.html +++ b/docs/articles/index.html @@ -76,7 +76,7 @@
mod_warfarin_nlmixr <- function() {
ini({
- #Fixed effects: population estimates
+ #Fixed effects: population estimates
THETA_ktr=0.106
THETA_ka=-0.087
THETA_cl=-2.03
@@ -110,7 +110,7 @@ warfarin PKPD modelTHETA_ec50=0.00724
THETA_kout=-2.9
THETA_e0=4.57
-
+
#Random effects: inter-individual variability
ETA_ktr ~ 1.024695
ETA_ka ~ 0.9518403
@@ -120,7 +120,7 @@ warfarin PKPD modelETA_ec50 ~ 0.7204165
ETA_kout ~ 0.3563706
ETA_e0 ~ 0.2660827
-
+
#Unexplained residual variability
cp.sd <- 0.144
cp.prop.sd <- 0.15
@@ -136,27 +136,29 @@ warfarin PKPD modelec50 = exp(THETA_ec50 + ETA_ec50)
kout = exp(THETA_kout + ETA_kout)
e0 = exp(THETA_e0 + ETA_e0)
-
+
#Structural model defined using ordinary differential equations (ODE)
DCP = center/v
PD=1-emax*DCP/(ec50+DCP)
-
+
effect(0) = e0
kin = e0*kout
-
+
d/dt(depot) = -ktr * depot
d/dt(gut) = ktr * depot -ka * gut
d/dt(center) = ka * gut - cl / v * center
d/dt(effect) = kin*PD -kout*effect
-
+
cp = center / v
pca = effect
-
+
#Model for unexplained residual variability
cp ~ add(cp.sd) + prop(cp.prop.sd)
pca ~ add(pca.sd)
})
- }
+}
+
+mod_warfarin_nlmixr <- mod_warfarin_nlmixr()
poso_estim_map()
map_warf_01 <- poso_estim_map(warf_01,mod_warfarin_nlmixr)
-#> using C compiler: ‘gcc (Ubuntu 11.4.0-1ubuntu1~22.04) 11.4.0’
map_warf_01
#> $eta
#> ETA_ktr ETA_ka ETA_cl ETA_v ETA_emax ETA_ec50
-#> -0.32052193 -0.54227235 0.79432182 -0.02986447 0.02383273 -0.28133882
+#> -0.32297013 -0.54003647 0.79433822 -0.02976759 0.02346678 -0.28140889
#> ETA_kout ETA_e0
-#> -0.30872226 -0.08386652
+#> -0.30867037 -0.08384263
#>
#> $model
#> ── Solved rxode2 object ──
@@ -201,11 +202,11 @@ data: first subject from t
#> THETA_ktr THETA_ka THETA_cl THETA_v THETA_emax THETA_ec50
#> 0.10600000 -0.08700000 -2.03000000 2.07000000 3.40000000 0.00724000
#> THETA_kout THETA_e0 cp.sd cp.prop.sd pca.sd ETA_ktr
-#> -2.90000000 4.57000000 0.14400000 0.15000000 3.91000000 -0.32052193
+#> -2.90000000 4.57000000 0.14400000 0.15000000 3.91000000 -0.32297013
#> ETA_ka ETA_cl ETA_v ETA_emax ETA_ec50 ETA_kout
-#> -0.54227235 0.79432182 -0.02986447 0.02383273 -0.28133882 -0.30872226
+#> -0.54003647 0.79433822 -0.02976759 0.02346678 -0.28140889 -0.30867037
#> ETA_e0
-#> -0.08386652
+#> -0.08384263
#> ── Initial Conditions ($inits): ──
#> depot gut center effect
#> 0 0 0 0
@@ -213,12 +214,12 @@ data: first subject from t
#> # A tibble: 1,451 × 18
#> time ktr ka cl v emax ec50 kout e0 DCP PD kin
#> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl>
-#> 1 0 0.807 0.533 0.291 7.69 0.968 0.760 0.0404 88.8 0 1 3.59
-#> 2 0.1 0.807 0.533 0.291 7.69 0.968 0.760 0.0404 88.8 0.0267 0.967 3.59
-#> 3 0.2 0.807 0.533 0.291 7.69 0.968 0.760 0.0404 88.8 0.102 0.885 3.59
-#> 4 0.3 0.807 0.533 0.291 7.69 0.968 0.760 0.0404 88.8 0.219 0.783 3.59
-#> 5 0.4 0.807 0.533 0.291 7.69 0.968 0.760 0.0404 88.8 0.373 0.681 3.59
-#> 6 0.5 0.807 0.533 0.291 7.69 0.968 0.760 0.0404 88.8 0.557 0.590 3.59
+#> 1 0 0.805 0.534 0.291 7.69 0.968 0.760 0.0404 88.8 0 1 3.59
+#> 2 0.1 0.805 0.534 0.291 7.69 0.968 0.760 0.0404 88.8 0.0267 0.967 3.59
+#> 3 0.2 0.805 0.534 0.291 7.69 0.968 0.760 0.0404 88.8 0.102 0.885 3.59
+#> 4 0.3 0.805 0.534 0.291 7.69 0.968 0.760 0.0404 88.8 0.219 0.783 3.59
+#> 5 0.4 0.805 0.534 0.291 7.69 0.968 0.760 0.0404 88.8 0.373 0.681 3.59
+#> 6 0.5 0.805 0.534 0.291 7.69 0.968 0.760 0.0404 88.8 0.557 0.590 3.59
#> # ℹ 1,445 more rows
#> # ℹ 6 more variables: cp <dbl>, pca <dbl>, depot <dbl>, gut <dbl>,
#> # center <dbl>, effect <dbl>
@@ -258,7 +259,7 @@ data: first subject from t