MarcLavielle
commented
5 years ago Hi,
"unfortunately", I was not able to reproduce this behavior (see the attached R code) best,
Marc
Le sam. 1 sept. 2018 à 00:17, skurup789 notifications@github.com a écrit :
Hi Marc,
I am trying to simulate a model with a dual absorption profile using differential equations. I created 2 lists for administration each going into a different target compartment in the PK block which is mapped by type. However, I notice that when trying to do this, it only picks up the first list. I am also trying to test different administration routes in the same run by using the group function. However, similarly, it only picks up the first list specified in treatment.
This is an excerpt of the code I am using:
PK depot(target = A1, type = 1, p=Bio F) depot(target = A2, type = 2, p=Bio(1-F) , Tlag = alag) depot(target = A3, type = 3, p=1) Group 1 --> Subcutaneous
adm1 <- list(time=0, amount = 600000, type = 1) #1st depot adm1.2 <- list(time=0, amount = 600000, type = 2) #2nd depot adm1.3 <-list(time = 0, amount = 0, type = 3) # central
g1 = list(size = 10, level = 'population', treatment =c(adm1, adm1.2, adm1.3)) Group 2 --> Intravenous
adm2 <- list(time = 0, amount = 0 , type = 1) # 1st depot adm2.2 <- list(time = 0, amount = 0, type = 2) # 2nd depot adm2.3 <- list(time = 0, amount = 900000, type = 3) # central
g2 = list(size = 10, level = 'population', treatment =c(adm2, adm2.2, adm2.3))
simulx(model = pkmodel.pop, output = list(cp,pop.parm), group = list(g1,g2), settings = list(seed = 12345))
I was wondering whether I am doing something wrong here? Thanks in advance!
— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub https://github.com/MarcLavielle/mlxR/issues/18, or mute the thread https://github.com/notifications/unsubscribe-auth/AJ3yur5JaCnbApspYOu7zzv1YGSkM8pHks5uWbYEgaJpZM4WV66v .
-- https://soutienvenezuela.org/en/ https://soutienvenezuela.org/en/
Marc Lavielle
CMAP, Ecole Polytechnique91128 Palaiseau Cedex
Francehttps://soutienvenezuela.org/en/ Tel: (33) 1 69 33 46 00 (33) 6 37 31 93 05
Marc.Lavielle@inria.fr http://www.cmap.polytechnique.fr/~lavielle
"""""""""""""""""""""""""""""""""""""""""""""""""""""""""
skurup789
Hi Marc,
I am trying to simulate a model with a dual absorption profile using differential equations. I created 2 lists for administration each going into a different target compartment in the PK block which is mapped by type. However, I notice that when trying to do this, it only picks up the first list. I am also trying to test different administration routes in the same run by using the group function. However, similarly, it only picks up the first list specified in treatment.
This is an excerpt of the code I am using:
PK depot(target = A1, type = 1, p=BioF) depot(target = A2, type = 2, p=Bio(1-F) , Tlag = alag) depot(target = A3, type = 3, p=1)
Group 1 --> Subcutaneous
adm1 <- list(time=0, amount = 600000, type = 1) #1st depot adm1.2 <- list(time=0, amount = 600000, type = 2) #2nd depot adm1.3 <-list(time = 0, amount = 0, type = 3) # central
g1 = list(size = 10, level = 'population', treatment =c(adm1, adm1.2, adm1.3))
Group 2 --> Intravenous
adm2 <- list(time = 0, amount = 0 , type = 1) # 1st depot adm2.2 <- list(time = 0, amount = 0, type = 2) # 2nd depot adm2.3 <- list(time = 0, amount = 900000, type = 3) # central
g2 = list(size = 10, level = 'population', treatment =c(adm2, adm2.2, adm2.3))
simulx(model = pkmodel.pop, output = list(cp,pop.parm), group = list(g1,g2), settings = list(seed = 12345))
I was wondering whether I am doing something wrong here? Thanks in advance!