Rで勉強2

Rで非線形回帰、2コンパートメントモデル。

p176 1コンパートメントモデル、IV投与
Rでやってみる
PKのパッケージにbiexp関数があるので
これを使う。
Estimation of initial and terminal half-life by fitting a biexponential model.

library(PK)
time <- c(0.1,0.2,0.4,0.7,1,1.5,2,3,4,6,8,12)
cp <- c(7.95,6.38,4.25,2.57,1.79,1.27,1.07,0.89,0.76,0.56,0.42,0.23)
res<-biexp(conc=cp, time=time, maxit=1000)
res
$parms
            initial  terminal
halflife  0.2581248 4.5789254
slope     2.6853184 0.1513777
intercept 8.5970653 1.3975529

傾きがα、βに相当するので
α=2.69h-1, β=0.15h-1,
切片の値からA=8.59 mg/L, B=1.40 mg/L
k21 = (Aβ+Bα) / (A+B) = 0.50
k10 = αβ/k21 = 0.80
k12 = α+β-k10-k21 = 1.54

AUC = A/a + B/b = 12.5 より
CL = S*F*D/AUC = 100 / 12.5 = 8 L/h (問題の設定にてdose = 100mg@iv)
k10 = CL /V1
V1 = CL / k10 = 10
k12 = CLd / V1より CLd = 15
Vb = SFD/(b*AUC)=100/(0.15*12.5)=53.3
となる。
全部をスクリプトに書いてconcとtimeだけ入れたら出るようにしたい。

Advertisements

Leave a Reply

Fill in your details below or click an icon to log in:

WordPress.com Logo

You are commenting using your WordPress.com account. Log Out / Change )

Twitter picture

You are commenting using your Twitter account. Log Out / Change )

Facebook photo

You are commenting using your Facebook account. Log Out / Change )

Google+ photo

You are commenting using your Google+ account. Log Out / Change )

Connecting to %s