Langsung ke konten utama

Postingan

Menampilkan postingan dari Agustus, 2020

Script Pengaturan Dosis Berulang IV Bolus Dua Kompartemen

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