Simulating Two PK Models Simultaneously.

[wayne200240]

Hi

I am working on antibiotic combination and would like to simulate two different PK models simultaneously to drive a PD model. I was wondering what is the best way to code this.

I have tried the following the code but it failed to run. Any advice would be much appreciated. Thank you very much.

Kind regards,

model <- new_ode_model(code = " dAdt[1] = -(CL/V)*A[1] - (Q/V)*A[1] + (Q/V2)*A[2]; dAdt[2] = -(Q/V2)*A[2] + (Q/V)*A[1]; dAdt[3] = - (CLR/VR)*A[3];") ; # A1 + A2 ==> two compartment model for Drug A ; # A3 --> One Compartment model for Drug B ; # PD model not included here.

r6 <- new_regimen(amt=c(rep(100,4),rep(2000,4)), time =c(0:3,0:3)*12, cmt = c(rep(1,4), rep(3,4)), type = c(rep("infusion",4), rep("infusion",4)), t_inf=c(rep(1,8), rep(1,8)))

dat3 <- sim_ode (ode = model, par = P_PKPD, regimen= R6 )

[melmelie]

I am trying to do the same, I think you can use a similar example to the PKPD example: http://pkpdsim.ronkeizer.com/example-pkpd-model.html You can use pk1 and pk2. You can specify the dosing compartment using cmt argument in the new_regimen function