OPTION LS=76 NOCENTER;
DATA host;INPUT method score @@;
CARDS;
1 96 1 79 1 91 1 85 1 83 1 91 1 82 1 87 . .
2 77 2 76 2 74 2 73 2 78 2 71 2 80 . .
3 66 3 73 3 69 3 66 3 77 3 73 3 71 3 70 3 74 
;
PROC PRINT;
TITLE 'Print all the data';

PROC GLM;CLASS method;
MODEL score=method;
MEANS method/DUNCAN ALPHA=0.05; /* For multiple comparisons */
MEANS method/DUNCAN ALPHA=0.01;
TITLE 'One way ANOVA, hostility data';
 
PROC SORT DATA=host;BY method;
PROC UNIVARIATE DATA=host PLOT; 
VAR score;BY method;
TITLE 'Box plot to check normality';

PROC MEANS NOPRINT DATA=host;VAR score;
BY method;
OUTPUT OUT=meanscr N=size MEAN=mean STD=std_dev
STDERR=stderr;
PROC PRINT DATA=meanscr;
TITLE 'Output from PROC MEANS';