+--------------------------------------------------------------------+
| |
| PK_PD_SIMULTANEOUS_1 EXAMPLE |
| |
+--------------------------------------------------------------------+
This is an example for simultaneous analysis of PK and direct PD.
Kinetics are one compartment. (See PK_PD_simultaneous_2, PK_PD_simul-
taneous_3).
$PROBLEM SIMULT PK/PD - ADVAN3
; THE DATA FILE CONTAINS BOTH CP AND EFFECT OBSERVATIONS.
; WHEN DV IS A CP OBSERVATION, CMT = 1 (OR 0)
; WHEN DV IS AN EFFECT OBSERVATION, CMT =2
$DATA data
$INPUT ID TIME DV AMT=DOSE CMT
$SUBROUTINES ADVAN3
$PK
K = THETA(1)*EXP(ETA(1))
K12= .001*K ; TRIVIAL LOSS TO EFFECT CMPT
K21= THETA(2)*EXP(ETA(2)) ; KEO
S1= THETA(3)*EXP(ETA(3)) ; V1 FOR DRUG
EMAX=THETA(4)*EXP(ETA(4))
C50= THETA(5)*EXP(ETA(5))
S2= S1*K12/K21 ; SO THAT CESS = CPSS
$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
REFERENCES: Guide VI Section IV.B.2, VII.C.3
Go to main index.
Created by nmhelp2html v. 1.0 written by Niclas Jonsson (Modified by AJB 5/2006,11/2007,10/2012)