+--------------------------------------------------------------------+
| |
| PK_PD_SIMULTANEOUS_2 EXAMPLE |
| |
+--------------------------------------------------------------------+
This is an example for simultaneous analysis of PK and PD. Kinetics
are one-compartment. (See PK_PD_simultaneous_1, PK_PD_simultane-
ous_3).
$PROBLEM SIMULTANEOUS PK/PD - ADVAN6
; NOTE: EXAMPLE USES 1-COMPT PK;
; MORE COMPLEX COMPARTMENTAL MODEL CAN BE USED
; CMT = 1 (OR 0) FOR DV = CP; CMT = 2 FOR DV = EFFECT
$DATA data
$INPUT ID TIME DV AMT=DOSE CMT
$SUBROUTINE ADVAN6 TOL=5
$MODEL
COMP=(CENTRAL,DEFDOSE,DEFOBS)
COMP=EFFECT
$PK
K10= THETA(1)*EXP(ETA(1))
V1= THETA(2)*EXP(ETA(2))
S1= V1
K20= THETA(3)*EXP(ETA(3))
C50= THETA(4)*EXP(ETA(4))
EMAX=THETA(5)*EXP(ETA(5))
$DES
DADT(1)=-K10*A(1) ;NOTE NO LOSS TO A(2)
DADT(2)= K20*(A(1)/V1-A(2)) ;A(2) IS CE
$ERROR
IF (CMT.EQ.2) THEN
Y = EMAX*F/(C50+F)*EXP(ERR(2)) ; DV = EFFECT IF CMT=2
ELSE
Y = F*EXP(ERR(1)) ; DV = CP IF CMT=1 (OR 0)
ENDIF
; ADD APPROPRIATE $THETA, $OMEGA, $SIGMA, $ESTIMATION, ETC.
REFERENCES: Guide IV Section V.C.6, V.C.7
REFERENCES: Guide VI Section IV.B.2, VII.C.6
Go to main index.
Created by nmhelp2html v. 1.0 written by Niclas Jonsson (Modified by AJB 5/2006,11/2007,10/2012)