/* x8_1.sas */
/* Purpose:  */

title1 'Example 8.1.';

data a;
  input x1 x2 x3 box$;
  x4=sqrt(x1*x1+x2*x2+x3*x3);
  x5=x4/x3;
  x6=(x1+x2)/x3;
  x7=(2*(x1*x2+x1*x3+x2*x3))/(x1*x2*x3);
  cards;
3.76    3.66    0.54   a
8.59    4.99    1.34   b
6.22    6.14    4.52   c
7.57    7.28    7.07   d
9.03    7.08    2.59   e
5.51    3.98    1.30   f
3.27    0.62    0.44   g
8.74    7.00    3.31   h
9.64    9.49    1.03   i
9.73    1.33    1.00   j
8.59    2.98    1.17   k
7.12    5.49    3.68   l
4.69    3.01    2.17   m
5.51    1.34    1.27   n
1.66    1.61    1.57   o
5.90    5.76    1.55   p
9.84    9.27    1.51   q
8.39    4.92    2.54   r
4.94    4.38    1.03   s
7.23    2.30    1.77   t
9.46    7.31    1.04   u
9.55    5.35    4.25   v
4.94    4.52    4.50   w
8.21    3.08    2.42   x
9.41    6.44    5.11   y
;

proc print data=a;
title2 'List of data.';

proc princomp data=a cov out=b;
  var x1-x7;
title2 'Principal component analysis.';

proc print data=b;
title2 'List of original and transformed data.';

proc plot data=b;
  plot prin2*prin1=box;
title2 'Plot of boxes using transformed variables.';

run;