* Study on variability in air velocity through clean filters;
option nodate;
title "Air flow rate through air filters";

proc format;
value $forsite "A"="Site A" "B"="Site B";
Data flowrate;
infile "c:\csc423\week9\flowrate.txt";
input rate site $;
label rate = "flow rate";
format site $forsite.; 
dum_site=(site="B");

proc sort;
by site;

/*proc ttest;
class site;
var rate;
*/
/* draw side-by-side box plots*/
proc boxplot;
plot rate*site;
/* regression analysis - ANOVA model*/
proc reg;
model rate=dum_site;
plot student.*predicted.;
plot npp.*student.;
run; quit;