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だけ入れたら出るようにしたい。

広告

コメントを残す

以下に詳細を記入するか、アイコンをクリックしてログインしてください。

WordPress.com ロゴ

WordPress.com アカウントを使ってコメントしています。 ログアウト / 変更 )

Twitter 画像

Twitter アカウントを使ってコメントしています。 ログアウト / 変更 )

Facebook の写真

Facebook アカウントを使ってコメントしています。 ログアウト / 変更 )

Google+ フォト

Google+ アカウントを使ってコメントしています。 ログアウト / 変更 )

%s と連携中