$PROBLEM TIME COURSE MODEL RUN=405
$INPUT CCOM ID TIME DV EVID CTS8 BASE ECOG CENS
$DATA EST_JHQG_S273_S188_OverallSurvival.csv IGNORE=C

$SUBR ADVAN=6 TOL=6
$MODEL COMP=(HAZARD)

$PK
LAM  = THETA(1)*EXP(ETA(1))
SHP  = THETA(2)
BET1 = THETA(3)
BET2 = THETA(4)
BET3 = THETA(5)

$DES

DADT(1)= LAM*SHP*(T+0.00001)**(SHP-1)*EXP(BET1*CTS8+BET2*BASE+BET3*ECOG)

$ERROR
IF(NEWIND.NE.2) OLDCHZ=0
CHZ=A(1)-OLDCHZ
OLDCHZ=A(1)

SUR = EXP(-CHZ)
HAZNOW= LAM*SHP*(TIME+0.0001)**(SHP-1)*EXP(BET1*CTS8+BET2*BASE+BET3*ECOG)

IF(CENS.EQ.1) Y=SUR
IF(CENS.EQ.0) Y=SUR*HAZNOW

$THETA (0 0.15)	; BASE
$THETA 0.0001	; SHP
$THETA 0.0001	; BET1
$THETA 0.0001	; BET2
$THETA 0.0001	; BET3
$OMEGA 0 FIX

$ESTIMATION NOABORT MAXEVAL=9999 PRINT=5 METHOD=0 LIKE
$COV
$TABLE ID TIME DV SUR Y
FILE=sdtab1.tb NOPRINT ONEHEADER