options nocenter ps=100 ls=79;

data a;
   infile 'product.dat' firstobs=2;
   input A B C D E $ y;

proc glm;
   class A B C D E;
   model y = A B C D E / ss3;

proc glm;
   class A B C D E;
   model y = A B D E / ss3;

proc glm;
   class A B C D E;
   model y = B A E A*E / ss3;

proc glm;
   class A B C D E;
   model y = B E / ss3;
   output out=partial r=ry;

proc glm;
   class A B C D E;
   model ry = A | D / ss3;
   lsmeans A*D / out=lsm;

proc plot data=lsm;
   plot lsmean*A=D;

proc glm data=partial;
   class A B C D E;
   model ry = A | E / ss3;
   lsmeans A*E / out=lsm;

proc plot data=lsm;
   plot lsmean*A=E;