question about disease progress and drop out joint model
Dear NMusers,
I am trying the joint modeling of disease progress and dropout model. NONMEM
did not indicate any error, but it looks like did not run it and there is no
results.
I am following Chuanpu Hu and Mak E. Sale's method. Our disease progress
model is y= base + (slope * time). Please find the code below and the
data file in attachment.
I would appreciate any feedback.
Kehua
$INPUT C ID TIME DV CMT EVID TYPE
$DATA 500.csv IGNORE=C
$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
Y=EXP(-CMHZ)
ENDIF
IF(TYPE.EQ.2.AND. DV.EQ.1) THEN
Y=EXP(-(CMHZ-HZLA))*(1-EXP(-HZLA))
ENDIF
$THETA
(0,70,1000); base
(0,.12,10); slope
(0,.0002,10); bshz
(0,.002,1); beta
(0,.002,1);
$OMEGA
0.03
0.03
$SIGMA 1 FIX
$ESTIM METHOD=1 LAPLACE MAXEVAL=9990
MSFO=008.MSF
$COV
500.csv
Description:
Binary data