/* x3_2_2.sas */
/* Purpose:  */

title1 'Section 3.2.2. Note: special notation used in the GLM analysis.';

data a;
  input y a b;
  x=(b-5)/1;
  x2=x*x;
  cards;
0.606 1 2
0.660 1 4
0.984 1 6
0.908 1 8
0.761 0 2
0.933 0 4
1.072 0 6
0.979 0 8
;

proc print;
title2 'List of data.';

proc glm data=a;
  class a;
  model y=a x(a) x2(a) /noint xpx inverse solution;
title2 'Full model. "M".';

proc glm data=a;
  class a;
  model y=a x(a) /noint xpx inverse solution;
title2 'Reduced model. "H1".';

proc glm data=a;
  class a;
  model y=a x /noint xpx inverse solution;
title2 'Reduced model. "H2".';

proc glm data=a;
  model y=x /xpx inverse solution;
title2 'Reduced model. "H3".';

proc glm data=a;
  model y= /xpx inverse solution;
title2 'Reduced model. "H4".';

run;