Simulation problem

From: "MANOJ KHURANA" Subject: [NMusers] Simulation problem Date: Mon, November 22, 2004 1:42 pm Dear Group, Continuing my crawls I have another querry if someone can suggest me how to perform PKPD simulation using NONMEM (If possible w.r.to PDX-POP) and what changes need to be incorporated in the data set (an example of data set arrangement and Control Stream file would be of real help). Below is the ctl stream I wrote but its not running so obviously there is something wrong that I am doing, kindly correct me. Can we perform estimation as well as simulation simultaneously. Also please advice me on how to perform individual simulation versus population simulation using NONMEM? I tried FIXING THETAs, REMOVING ETAs as well as OMEGA BLOCK for individual simulation but run aborts with a msg that "sigma record found without omega record". Thanks in advance Manoj Khurana $PROB RUN# 001 $INPUT C ID EVID TIME DV TYPE AMT $DATA 001.CSV IGNORE=C $SUBROUTINES ADVAN4 TRANS4 $PK TVCL=THETA(1) CL=TVCL*EXP(ETA(1)) TVV2=THETA(2) V2=TVV2*EXP(ETA(2)) TVQ=THETA(3) Q=TVQ*EXP(ETA(3)) TVV3=THETA(4) V3=TVV3*EXP(ETA(4)) TVKA=THETA(5) KA=TVKA*EXP(ETA(5)) S2=V2/1000 $ERROR IPRED=F EMAX=THETA(6)+ETA(6) EC50=THETA(7)+ETA(7) E=(EMAX*F)/(EC50+F) Y=F*(1+ERR(1)) IF (TYPE.EQ.2) THEN Y=E+ERR(2) ENDIF $THETA (10, FIX) ;CL (100, FIX) ;V2 (12, FIX) ;Q (120, FIX) ;V3 (2, FIX) ;KA (100, FIX) ;EMAX (10, FIX) ;EC50 $OMEGA 0.03 FIX 0.08 FIX 0.3 FIX 0.13 FIX 0.99 FIX 0.7 FIX 0.01 FIX $SIGMA 0.02 FIX 1.6 FIX $SIMULATION (2674474) ONLYSIMULATION $TABLE ID TIME IPRED TYPE NOPRINT ONEHEADER FILE=001.TAB