$PROBLEM ESTIMATE DATA
$INPUT ID ET1A TIME OLDT TDOS P PAE DV
$DATA t1sim.txt IGNORE=#
$PRED
IND1A=0
IF (NEWIND.NE.2) THEN
   SW1A=0
  ; PAE=0
   FLG1=0
ENDIF 
PAE1=PAE
;TYPE=1 PART
;TIME TO FIRST AE STRUCTURAL PART
   ;WEIBULL FAILURE PARAMETER
   ALPH1=THETA(1)
   ALPH1C=THETA(2)
   GAM1=THETA(6)
   MED1=THETA(3)*(1-TDOS**GAM1*THETA(4)/(THETA(5)**GAM1+TDOS**GAM1))
   LAM1=0.693/MED1**ALPH1
    CP1=THETA(7)
   IF (TIME.LT.CP1) THEN
      H1ZT1=LAM1*TIME**ALPH1
      H1ZT0=LAM1*OLDT**ALPH1
   ELSE
      H1ZT1=LAM1*TIME**ALPH1C*CP1**(ALPH1-ALPH1C)
      H1ZT0=LAM1*OLDT**ALPH1C*CP1**(ALPH1-ALPH1C)
   ENDIF

   I1HAZ=H1ZT1-H1ZT0
   Q1=EXP(-I1HAZ)
   ;NEED PREVIOUS DV HERE
   ;MC PART
   ; baseline values for severity (logit)
   IF (PAE.EQ.0) THEN
      B11=THETA(8)
      B12=B11+THETA(9)
      B13=B12+THETA(10)
   ENDIF 

   IF(PAE.EQ.1) THEN
     B11=THETA(11)
     B12=B11+THETA(12)
     B13=B12+THETA(13)
     EMAX=0
   ENDIF
   IF(PAE.EQ.2) THEN
     B11=THETA(14)
     B12=B11+THETA(15)
     B13=B12+THETA(16)
     EMAX=0
   ENDIF

   IF(PAE.EQ.3) THEN
     B11=THETA(17)
     B12=B11+THETA(18)
     B13=B12+THETA(19)
     EMAX=0
   ENDIF

   
   RESP=0                
   ;-----------------------------------------------------
                  ; INTERSUBJECT VARIABILITY
  IF (ETA(1).GE.0) THEN
     ET1 = ((ETA(1)+1)**THETA(20)-1)/THETA(20)
  ELSE
     ET1 = -((1-ETA(1))**THETA(20)-1)/THETA(20)
  ENDIF

   
  
   
   ;-----------------------------------------------------
   ; logits for Y>=1, Y>=2, Y.=3
   A11 = B11 + ET1
   A12 = B12 + ET1
   A13 = B13 + ET1
   C11=EXP(A11)
   C12=EXP(A12)
   C13=EXP(A13)
   ;-----------------------------------------------------
   ; probabilities for Y>=1, Y>=2, Y>=3
   P11=C11/(1+C11)
   P12=C12/(1+C12)
   P13=C13/(1+C13)             
   ;IF (IND1A.EQ.0) THEN
      ;Probabilities for TRANSITION TO Y=0  Y=1, Y=2, Y=3 when it is beyond second transition
      PA1=1-P11
      PB1=P11-P12
      PC1=P12-P13
      PD1=P13
   ;ENDIF
   
IF (DV.EQ.0) THEN
   Y=PA1
ENDIF 
IF (DV.EQ.1) THEN
   Y=PB1
ENDIF 
IF (DV.EQ.2) THEN
   Y=PC1
ENDIF 
IF (DV.EQ.3) THEN
   Y=PD1
ENDIF 


   
   

   

$THETA

;AE1 STRUCTURAL PART
 (0.498 FIX)          ;1 ALPHA TTFAE1 SHAPE
 (.117 FIX)       ;2 SECOND ALPHA TTFAE1 SHAPE
 (1530 FIX)           ;3 MEDIAN TTFAE1 
 (0.99 FIX)       ;4 TDOS EMAX ON TTFAE1 MEDIAN
 (53.9 FIX)           ;5 TDOS ED50 ON TTFAE1 MEDIAN
 (6 FIX )         ;6 GAM1 ON TTFAE1 MEDIAN 
 (6 FIX)            ;7 CHANGE POINT FOR TTFAE1
 (-3.82)          ;8 B1 PAE1=0
 (-10,-0.897,0)         ;9 B2 PAE1=0
 (-10,-1.96,0)          ;10 B3 PAE1=0
 (2.78)           ;11 B1 PAE1=1
 (-10,-1.2,0)          ;12 B2 PAE1=1
 (-10,-6,0)        ;13 B3 PAE1=1
 (1.52)           ;14  B1 PAE1=2
 (-10,-0.136,0)         ;15  B2 PAE1=2
 (-10,-1.8,0)          ;16  B3 PAE1=2
 (2.48)           ;17 B1 PAE1=3
 (-10,-1.001,0)     ;18 B2 PAE1=3
 (-10,-2.108,0)         ;19 B3 PAE1=3
 (1 FIX)          ;20 JD ON ETA1
$OMEGA .163 
$EST METH=1 LAPLACE LIKE MAXEVALS=9999 PRINT =1
$COV PRINT=E MATRIX=R
$TABLE ID ETA1 NOHEADER NOPRINT FIRSTONLY FILE=T400A.TXT
;