proc glm data=factexpt;
class A B;
model yield = A B A*B;
means A B / tukey lines;
title ‘Factorial Experiment Example’;
run;