READ;NVAR=60;NOBS=9092;NAMES=SEVERITY,PDO,NONIN,INCAP,FATAL,MASTER,
VEHAGE,HARLEY,HONDA,YAMAHA,SUZUKI,KAWASAKI,SPORTBIK,PASSNGR,ANIMAL,
GUARDRL,TREE,PEDBIKE,CURB,POLE,VEHICLE,WALL,CULVERT,DITCH,EMBANK,
TRAFCONT,AGE,FEMALE,EJECTED,PINNED,HELMET,MONTH,YEAR,LOCAL,STATEUS,
INTSTATE,RURAL,DAWNDUSK,DARK,LIGHTED,WETPAVE,INTSECT,HC,GRADE,CREST,
VEHFAIL,OFFRD,ALCOHOL,SPEEDING,NL,PO,MC,LM,LP,OP,CD,CH,BRCMONTH,
BRCPASS,SPDLIMIT
;file=D:\old_drive_d\new_laptop\CE697M\MC1.txt$

CREATE ;IF(BRCMONTH>24) BRC4=1; (ELSE) BRC4=0$
CREATE ;IF(MONTH=4) APRIL=1; (ELSE) APRIL=0$
CREATE ;IF(MONTH=7) JULY=1; (ELSE) JULY=0$
CREATE ;IF(SPDLIMIT>50) OVER50=1; (ELSE) OVER50=0$
CREATE ;IF(VEHAGE<5) NEW5=1; (ELSE) NEW5=0$


NLOGIT
;LHS=SEVERITY
;CHOICES = PDOPOSS,NONIN,INCAP,FATAL
;TREE = NOTINJ(PDOPOSS,NONIN),THREE(INCAP),FOUR(FATAL)
;MODEL:
U(PDOPOSS) = PDOPOSS*ONE+WETPAVE1*WETPAVE+HELMETN*HELMET+
INTSECT*INTSECT+NEW5*NEW5/
U(NONIN) = SPEEDN*SPEEDING+ALCOHOLN*ALCOHOL+PASSN*PASSNGR+
FEMALE1*FEMALE/
U(INCAP) = INCAP*ONE+TREEF*TREE+POLEF*POLE+DARKF*DARK/
U(FATAL) = FATAL*ONE+SPEEDF*SPEEDING+APRILF*APRIL+
JULYF*JULY+BRC4*BRC4/
U(NOTINJ) = HC*HC+AGE*AGE+ALCOHOL*ALCOHOL+
SPEED*SPEEDING+HELMET*HELMET+OVER50*OVER50
;IVSET: (THREE) = [1] / (FOUR) = [1]$