options ls=78;

data cls11;

infile 'class11.dat';

input sn region $ recklessness happiness;

thappiness=sqrt(happiness);

run;

proc print data=cls11;

run;

proc univariate data=cls11 plot normal;

var recklessness happiness thappiness;

title 'Univariate for All Variables';

run;

proc sort data=cls11;

by region;

run;

proc univariate data=cls11 plot normal;

by region;

var recklessness happiness thappiness;

title 'Univariates by Region for All Variables';

run;

proc anova data=cls11;

class region;

model recklessness=region;

means region/hovtest;

means region/t;

means region/bon;

means region/tukey;

means region/regwq;

title 'ANOVA for Recklessness by Region';

run;

proc anova data=cls11;

class region;

model happiness=region;

means region/hovtest;

means region/t;

means region/bon;

means region/tukey;

means region/regwq;

title 'ANOVA for Happiness by Region';

run;

proc anova data=cls11;

class region;

model thappiness=region;

means region/hovtest;

means region/t;

means region/bon;

means region/tukey;

means region/regwq;

title 'ANOVA for tHappiness by Region';

run;

quit;