read;nvar=32;nobs=825;names=
ID,INJFREQ,ROUTE,LENGTH,INCLANES,DECLANES,WIDTH,MIMEDSH,
MXMEDSH,SPEED,URB,FC,AADT,
SINGLE,DOUBLE,TRAIN,PEAKHR,GRADEBR,MIGRADE,MXGRADE,MXGRDIFF,
TANGENT,CURVES,MINRAD,ACCESS,MEDWIDTH,
FRICTION,ADTLANE,SLOPE,
INTECHAG,AVEPRE,AVESNOW;
FILE="D:\OLD_DRIVE_D\NEW_LAPTOP\CE697M\MIXED-LOGIT-Book.csv "$

/*ADT VARIABLE*/
create;laneadt=aadt/(inclanes+declanes)$
create;lnlanadt=log(laneadt)$
create;lnaadt=log(aadt)$
create;density=laneadt/length$

/*FRICTION VARIABLE*/
create;if(friction<=30)lowfri=1$
create;if(friction>30&friction<50)medfri=1$
create;if(friction>=50)hifri=1$

/*CURVE/LENGTH VARIABLE*/
/*dstat;rhs=curvmile$*/
create;curvmile=curves/length$
create;if(curvmile<=0.5)lowcvmil=1;(else)lowcvmil=0$
create;if(curvmile>0.5&curvmile<=2.5)medcvmil=1;(else)medcvmil=0$
create;if(curvmile>2.5)hicvmil =1;(else)hicvmil=0$
/*%TRUCK VARIABLE*/ 
create;truck=single+double+train$
create;pcttruck=truck/aadt$

/*MEDWIDTH VARIABLE*/
/*dstat;rhs=medwidth$
histogram;rhs=medwidth$*/
create;if(medwidth=1)med030=1$
create;if(medwidth=2)med3040=1$
create;if(medwidth=3)med4050=1$
create;if(medwidth=4)med5060=1$
create;if(medwidth=5)med60=1$

/*SPEED VARIABLE*/
/*dstat;rhs=speed$
histogram;rhs=speed$*/
create;if(speed<=50)speed1=1$
create;if(speed<=55)speed2=1$
create;if(speed>55)speed3=1$
create;if(speed>=55)speed4=1$

/*FUNCTIONAL CLASS VARIABLE*/
/*dstat;rhs=fc$
histogram;rhs=fc$*/
create;if(fc=1)local=1$
create;if(fc=5)intstate=1$

/*ACESS VARIABLE*/
create;if(access=0)none   =1$
create;if(access=1)partial=1$
create;if(access=2)full   =1$

/*SLOPE VARIABLE*/
create;if(slope=0)flat=1$
create;if(slope=1)slight=1$
create;if(slope=2)medium=1$
create;if(slope=0 |slope= 1)slpflat=1;(else)slpflat=0$
create;if(slope=2)slpmed=1;(else)slpmed=0$

/***WEATHER VARIABLE***/

/*histogram;rhs=avepre,avesnow$*/

create;if(avepre<=1.5)lowpre=1;(else)lowpre=0$
create;if(avepre>1.5&avepre<=2.5)medpre=1;(else)medpre=0$
create;if(avepre>2.5)hipre=1;(else)hipre=0$
/*dstat;rhs=avepre$
histogram;rhs=avepre$*/
create;if(avesnow<=1)norsnow=1$
create;if(avesnow>1)hisnow=1$

create;lanewid=(inclanes+declanes)/width$
dstat;rhs=lanewid$
create;if(lanewid<12)nlanwid=1;(else)nlanwid=0$
create;if(lanewid>12)wlanwid=1;(else)wlanwid=0$

create;intmi=intechag/length$
create;gbmile=gradebr/length$
nlogit;lhs=injfreq;
choices=pdo,pinj,inj;
model:
U(pdo)=a0+a1*laneadt+a3*minrad/
U(pinj)=b0+b2*truck/
U(inj)=c3*friction+c2*intmi+c1*gbmile
;fcn=a0(c),a1(c),a3(N),
b0(c),b2(N),c2(n),c3(c),c1(N);rpl;frequencies;parameter;pts=200,halton$