Langkah-langkah simulasi farmakokinetika dengan R #Memanggil paket library(ggplot2) library(ggpubr) library(deSolve) #Mendefinisikan persamaan iv.2comp <-function(t, state, parameters) { with(as.list(c(state, parameters)),{ # rate of change dC1 <- k21*C2/V1 - (k12+ke)*C1 dC2 <- k12*C1/V2 - k21*C2 # return the rate of change list(c(dC1, dC2)) }) # end with(as.list ... } #dose = 100 #V1 = 10 #C1 pada waktu nol = dose/V1 state <- c(C1 = 10, C2 = 0) #pada waktu ke nol time <- seq(0, 100, .01) #Interval waktu pengamatan #mendefinisikan parameter parameters <- c(k12 = .23, k21 = .11, ke = 0.32, V1 = 10, V2 = .5) #Simulasi farmakokinetika untuk dosis tunggal out <- ode(y = state, times = time, func = iv.2comp, parms= parameters) dat1 <- as.data.frame(out) head(dat1) ## Plot simulasi farmakokinetika pada dosis tunggal c1 <- ggplot(dat1, aes(y = C1, x = time)) + geom_line(colour = 'red', size=1 + scale_y_log10() c2 <- ggplot(dat1, aes(y = C2, x