All,
Inter-occasion variability is typically coded like the example in NONMEM help:
$PROB prazosin data
$INPUT ID AGE RACE HT HCTZ AMT TIME DV EVID SS II OCC
$DATA praz21
$SUBROUTINE ADVAN2 TRANS2
$PK
OCC2=1-OCC
R=0
IF (RACE.EQ.2) R=1
TMP=(HT-160)*THETA(1)
IF (TMP.LE.0.) TMP=0.0
TMP2=(AGE-60)*THETA(2)
IF (AGE.LE.60.) TMP2=0.0
TVCL=THETA(3)+TMP-TMP2
TVCL=TVCL+THETA(4)*R+THETA(5)*HCTZ
IF (TVCL.LE.0.) EXIT 1 100
TVV=THETA(6)+HCTZ*THETA(7)
IF (TVV.LE.0.) EXIT 1 200
TVKA=THETA(8)
CL=TVCL*EXP(ETA(3)*OCC+ETA(5)*OCC2+ETA(1))
V =TVV *EXP(ETA(4)*OCC+ETA(6)*OCC2+ETA(2))
KA=TVKA*EXP(ETA(8)*OCC+ETA(9)*OCC2+ETA(7))
S2=V
$THETA (0,0.4) (0.,0.3) 20. 4 -4.
$THETA 80 -25
$THETA (0,1.)
$OMEGA BLOCK(2) 0.2 0.1 0.2
$OMEGA BLOCK(2) 0.2 0.1 0.2
$OMEGA BLOCK(2) SAME
$OMEGA 1.2
$OMEGA BLOCK(1) .5
$OMEGA BLOCK(1) SAME
$ERROR
Y=F*(1+EPS(1))
$SIGMA 0.1
$EST NOABORT SIG=3 MAX=3000 PRINT=10 POSTHOC
$COVARIANCE
$TABLE ID OCC ETA1 ETA3 ETA5 ETA2 ETA4 ETA6 NOPRINT FILE=praz32.tab
Data for one subject follow.
21 57 1 180 1 1000 0 0.62 1 1 12 0
21 57 1 180 1 0.0 1 12.06 0 0 0 0
21 57 1 180 1 0.0 3 6.81 0 0 0 0
21 57 1 180 1 0.0 4 4.89 0 0 0 0
21 57 1 180 1 0.0 5 4.04 0 0 0 0
21 57 1 180 1 0.0 6 2.82 0 0 0 0
21 57 1 180 1 0.0 7 2.72 0 0 0 0
21 57 1 180 1 0.0 8 1.78 0 0 0 0
21 57 1 180 1 0.0 10 1.07 0 0 0 0
21 57 1 180 1 0.0 12 0.75 0 0 0 0
21 57 1 180 1 0.0 0 2.83 3 0 0 0
21 57 1 180 1 1000 0 2.83 1 1 12 1
21 57 1 180 1 0.0 1 6.27 0 0 0 1
21 57 1 180 1 0.0 2 12.72 0 0 0 1
21 57 1 180 1 0.0 3 9.99 0 0 0 1
21 57 1 180 1 0.0 4 6.90 0 0 0 1
21 57 1 180 1 0.0 5 6.59 0 0 0 1
21 57 1 180 1 0.0 6 4.07 0 0 0 1
21 57 1 180 1 0.0 7 3.11 0 0 0 1
21 57 1 180 1 0.0 8 2.79 0 0 0 1
21 57 1 180 1 0.0 10 1.72 0 0 0 1
21 57 1 180 1 0.0 12 1.17 0 0 0 1
Does anyone know if changing the dosing records to be missing will affect the
results, so the following dataset is used:
$INPUT ID AGE RACE HT HCTZ AMT TIME DV EVID SS II OCC
Data for one subject follow.
21 57 1 180 1 1000 0 0.62 1 1 12 .
21 57 1 180 1 0.0 1 12.06 0 0 0 1
21 57 1 180 1 0.0 3 6.81 0 0 0 1
21 57 1 180 1 0.0 4 4.89 0 0 0 1
21 57 1 180 1 0.0 5 4.04 0 0 0 1
21 57 1 180 1 0.0 6 2.82 0 0 0 1
21 57 1 180 1 0.0 7 2.72 0 0 0 1
21 57 1 180 1 0.0 8 1.78 0 0 0 1
21 57 1 180 1 0.0 10 1.07 0 0 0 1
21 57 1 180 1 0.0 12 0.75 0 0 0 1
21 57 1 180 1 0.0 0 2.83 3 0 0 1
21 57 1 180 1 1000 0 2.83 1 1 12 .
21 57 1 180 1 0.0 1 6.27 0 0 0 2
21 57 1 180 1 0.0 2 12.72 0 0 0 2
21 57 1 180 1 0.0 3 9.99 0 0 0 2
21 57 1 180 1 0.0 4 6.90 0 0 0 2
21 57 1 180 1 0.0 5 6.59 0 0 0 2
21 57 1 180 1 0.0 6 4.07 0 0 0 2
21 57 1 180 1 0.0 7 3.11 0 0 0 2
21 57 1 180 1 0.0 8 2.79 0 0 0 2
21 57 1 180 1 0.0 10 1.72 0 0 0 2
21 57 1 180 1 0.0 12 1.17 0 0 0 2
With the modification of the control stream as follows:
$PK
OCC2=0
IF(OCC.EQ.2)OCC2=1
Matt.
________________________________
This e-mail (including any attachments) is confidential and may be legally
privileged. If you are not an intended recipient or an authorized
representative of an intended recipient, you are prohibited from using, copying
or distributing the information in this e-mail or its attachments. If you have
received this e-mail in error, please notify the sender immediately by return
e-mail and delete all copies of this message and any attachments.
Thank you.
