/* x6_1.sas */
/* Purpose:  */

title1 'Example 6.1.';

data a;
  input group$ si al number;
  cards;
earth 22.10 7.40 100000000
moon  19.4  5.9  1
moon  21.5  4.0  1
moon  19.2  4.0  1
moon  18.4  5.4  1
moon  20.6  6.2  1
moon  19.8  5.7  1
moon  18.7  6.0  1
;

proc print data=a;
title2 'List of data. Number of samples for earth set "high" (arbitrary).';

proc corr cov data=a;
  var si al;
  by group;
title2 'Simple statistics.';

proc glm data=a;
  class group;
  model si al=group;
  manova h=group /printe printh;
  weight number;
title2 'Test using GLM manova, and number of samples as weight. (Quick and dirty!)';

run;