title 'Three-Way ANOVA';
/* Three sets of instructions were used to assess the effects
	of anxiety on test performance: 
	(1) Instructions to reduce anxiety (Anxgp='Low')
	(2) Neutral instructions
	(3) Instructions designed to increase anxiety (Anxgp='High')
	Groups of n=10 men and women took either a verbal ability test or a
	math ability test under one of these conditions.
*/
data anxiety;
	n=10;
	input anxgp $ sex $ test $ mean var;
	std=sqrt(var);
	datalines;
Low     M  Verbal   12.7   28.1
Low     F  Verbal   14.0   25.
Low     M  Math     16.0   20.3
Low     F  Math     10.5   18.5
Neutral M  Verbal   13.0   22.1
Neutral F  Verbal   13.2   26
Neutral M  Math     12.6   23
Neutral F  Math     10.2   25
High    M  Verbal   12.1   27
High    F  Verbal   16.3   39.7
High    M  Math     12.2   43.6
High    F  Math     13.2   29.2
;
%stat2dat(data=anxiety, out=raw, 
	class=anxgp sex test,
	mean=mean, std=std, n=n);

proc glm data=raw order=data;
	class anxgp sex test;
	freq freq;
	model Y = anxgp|test|sex;
	run;