Re: DV Simulation
From:Toufigh Gordi
Subject: Re: [NMusers] DV Simulation
Date: Thu, 29 May 2003 14:03:07 -0400
Hi!
I am trying to simulate 100 subjects using the follwing coding:
$PROB ALL DOSES
; 2-COMP NLCL
$INPUT ID AMT RATE DUR DATE=DROP TIME DV EVID SEXF PKPD
$DATA simbase2.csv IGNORE=#
$SUBROUTINE ADVAN6 TRANS1 TOL=5
$MODEL NCOMP=2
COMP=CENT
COMP=PERIPH
$PK
"FIRST
" COMMON /PRCOMG/ IDUM1,IDUM2,IMAX,IDUM4,IDUM5
" INTEGER IDUM1,IDUM2,IMAX,IDUM4,IDUM5
" IMAX=1000000000
VM=THETA(1)*EXP(ETA(1))
KM=THETA(2)
V1=THETA(3)*EXP(ETA(2))
Q=THETA(4)
V2=THETA(5)
D1=DUR
S1=V1
$DES
CONC=A(1)/V1
DADT(1)=-A(1)*(VM/(KM+CONC))/V1-A(1)*Q/V1+A(2)*Q/V2
DADT(2)=A(1)*Q/V1-A(2)*Q/V2
CLN=VM/(KM+CONC)
K12=Q/V1
K21=Q/V2
K10=CLN/V1
BETA=0.5*((K12+K21+K10)-SQRT((K12+K21+K10)**2-4*K21*K10))
BTH=LOG(2)/BETA
VSS=V1+V2
$ERROR
Y=F*(1+EPS(1))
REPI=IREP
$THETA (2.92 FIX) ;1 VM
$THETA (2.87 FIX) ;2 KM
$THETA (42.6 FIX) ;3 V1
$THETA (0.421 FIX) ;4 Q
$THETA (48.1 FIX) ;5 V2
$OMEGA 0.105 FIX ;1 IIV VM
$OMEGA 0.0176 FIX ;2 IIV V1
$SIGMA 0.00987 FIX
$SIMULATION (20030521) ONLYSIM SUBPROBLEMS=100
$TABLE ID TIME DV REPI NOPRINT ONEHEADER FILE=SDTAB101
$TABLE ID VSS CLN K12 K21 K10 BETA BTH
NOPRINT ONEHEADER FILE=MYTAB101
The simulation works and the tables are generated. However, in the SDTAB101
file I get the results for only 1 subject and not 100, as written in the
model. Could anyone comment on what mistake I am making?
Thank you.
Toufigh Gordi