DM "output;clear;log;clear";
Options pagesize=45 linesize=80 PageNo=1 NoDate;

ODS RTF File="C:\teaching\stat530\p6_7.doc"; *this will make a output file. specify a folder and file name;
ODS Listing Close;								

Data p6_7;
    input A B C D Y;
	AB = A*B; AC = A*C; AD = A*D; BC = B*C; BD = B*D; CD = C*D; ABC=AB*C; ABD = AB*D; ACD = AC*D; BCD = BC*D; ABCD = ABC*D;
	Datalines;
	-1 -1 -1 -1 90
	-1 -1 -1 -1 93
	1 -1 -1 -1 74
	1 -1 -1 -1 78
	-1 1 -1 -1 81
	-1 1 -1 -1 85
	1 1 -1 -1 83
	1 1 -1 -1 80
	-1 -1 1 -1 77
	-1 -1 1 -1 78
	1 -1 1 -1 81
	1 -1 1 -1 80
   -1 1 1 -1 88
   -1 1 1 -1 82
	1 1 1 -1 73
	1 1 1 -1 70
	-1 -1 -1 1 98
	-1 -1 -1 1 95
	1 -1 -1 1 72
	1 -1 -1 1 76
	-1 1 -1 1 87
	-1 1 -1 1 83
	1 1 -1 1 85
	1 1 -1 1 86
	-1 -1 1 1 99
	-1 -1 1 1 90
	1 -1 1 1 79
	1 -1 1 1 75
	-1 1 1 1 87
	-1 1 1 1 84
	1 1 1 1 80
	1 1 1 1 80
	;


proc rank normal=vw; /* Computing ranked normal scores by residuals*/
var Y;
ranks nscore; run;
proc gplot ;
plot Y*nscore; /*plotting ranked residual vs. normal score*/
label nscore='Normal Score'; run;


proc glm data = p6_7; 
class A B C D AB AC AD BC BD CD ABC ABD ACD BCD ABCD;
model Y = A B C D AB AC AD BC BD CD ABC ABD ACD BCD ABCD;
estimate 'A' A 1 -1; 
run;

 
ODS Listing;
ODS RTF Close;