-
Notifications
You must be signed in to change notification settings - Fork 0
/
EDA_KSI .sas
64 lines (49 loc) · 1.24 KB
/
EDA_KSI .sas
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
%let path=/home/u42959472/Capstone;
libname tps "/home/u42959472/Capstone";
/*Access Data */
options validvarname=v7;
proc import datafile="/home/u42959472/Capstone/KSI.csv"
dbms = csv
out = tps.ksi
replace;
guessingrows=max;
run;
/*Explore Data */
proc print data= tps.ksi (obs=30);
run;
proc contents data= tps.ksi varnum;
run;
proc sort data=tps.ksi
out=tps.ksi_nodups
nodupkey;
by _all_;
run;
proc freq data=tps.ksi_nodups nlevels;
table _char_ /noprint;
table _numeric_ /noprint;
run;
proc sort data=tps.ksi_nodups;
by Year Date Month Var6 Time;
run;
data tps.ksi_cleaned;
set tps.ksi_nodups;
run;
/* Analyze Data */
title "Number of KSI cases from 2008 to 2018";
proc freq data=tps.ksi_cleaned;
table Year*Fatal /nocum nopercent plots=freqplot;
run;
title;
title "Number of KSI cases due to aggressive driving";
proc freq data=tps.ksi_cleaned;
table Year*AG_DRIV /nocum nopercent plots=freqplot;
run;
title;
/* Fatal and Ag_Dri Cases */
%let Year=2018;
Title "Hazardous driving hot spot regions in the City of Toronto in &Year";
proc sgmap plotdata= tps.ksi_cleaned (where=(AG_DRIV=1 and Fatal=1 and Year=&Year));
openstreetmap;
scatter x=Longitude y=Latitude/ markerattrs=(symbol=star color=red size=10 px);
run;
title;