question about simulation based on joint disease progress and dropout model
Dear NMusers,
We are trying to do VPC based on a joint disease progress and dropout model.
But we did not get any results from the simulation. Does anybody have any
suggestions?
Thanks a lot.
Kehua
The control stream are below,
$SUBS ADVAN6 TOL=6
$MODEL COMP=(CUMHAZ)
COMP=(HZLAST)
$PK
BASE=THETA(1)*EXP(ETA(1))
SLOPE=THETA(2)+ETA(2)
BSHZ=THETA(3)
BETA=THETA(4)
$DES
SIZE=BASE+(SLOPE*TIME)
TEMP=BETA*SIZE
DADT(1)=EXP(TEMP)
DADT(2)=EXP(TEMP)
$ERROR
CMHZ=BSHZ*A(1)
HZLA=BSHZ*A(2)
IF (TYPE.EQ.1) THEN
IPRED=BASE+(SLOPE*TIME)
W=THETA(5)*IPRED
F_FLAG=0
Y=BASE+(SLOPE*TIME)+W*EPS(1)
ENDIF
IF(TYPE.EQ.2.AND.DV.EQ.0) THEN
F_FLAG=1
Y=EXP(-CMHZ)
ENDIF
IF(TYPE.EQ.2.AND. DV.EQ.1) THEN
F_FLAG=1
Y=EXP(-(CMHZ-HZLA))*(1-EXP(-HZLA))
ENDIF
$THETA 63.0 0.153 0.00811 0.00633 0.063
$OMEGA 0.507 0.0914
$SIGMA 1 FIX
$SIMULATION ONLYSIM SUBPROB=1 (1111111)
$TABLE ID TIME IPRED
NOPRINT ONEHEADER FILE=015.tab