+--------------------------------------------------------------------+
 |                                                                    |
 |                    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)