Nonlinear Mixed Models: PK, One-compartment model (oral dose)

[안재형]

Theoph.csv

proc nlmixed data=Theoph;
  parms beta1=0.5, beta2=1.5, beta3=.1, s2b1=.1, s2b2=.1, cb12=0, s2e=2;
  V=beta1+u1;
  ka=beta2+u2;
  ke=beta3;
  eta= Dose*ka/V/(ka-ke)*(exp(-ke*time)-exp(-ka*time));
  model conc ~ normal(eta, s2e);
  random u1 u2 ~ normal([0,0],[s2b1,cb12,s2b2]) subject=subject;

  predict eta out=theophout;
run; 

                                Standard
Parameter   Estimate      Error    DF   t Value   Pr > |t|   Alpha   Lower     Upper   Gradient

beta1         0.4709      0.02454     10     19.19     <.0001     0.05     0.4162     0.5256   -0.00608
beta2         1.8987      0.3890      10      4.88      0.0006     0.05     1.0319     2.7656   0.001329
beta3        0.08530    0.004441     10     19.20     <.0001     0.05    0.07540    0.09519   -0.03315
s2b1        0.005732   0.002542     10      2.25      0.0478     0.05    0.000068   0.01140   -0.06633
s2b2          1.5846     0.9006       10      1.76      0.1090     0.05    -0.4221     3.5913   -0.00012
cb12        -0.00315    0.03026     10     -0.10      0.9191     0.05   -0.07059    0.06428   -0.03607
s2e           0.5184      0.07231     10      7.17      <.0001     0.05     0.3573     0.6795   0.000078 

proc sgplot data=theophout noautolegend;
  series x=time y=pred / group=subject;
  scatter x=time y=conc;
run;