/* Homework 3: Cluster Analysis on Utility Companies                      */
/* WITH SELECTED OUTPUT													  */

OPTIONS PS=52;
OPTIONS LS=78;
OPTIONS NODATE;
OPTIONS PAGENO=1;

DATA UTILITIES;
   INFILE 'A:/UTILITY.PRN';
   INPUT COMPANY $ X1-X8;
   LABEL COMPANY = 'Company';
   LABEL X1 = 'Fixed-charge coverage ratio'
   LABEL X2 = 'Rate of return on capital';
   LABEL X3 = 'Cost per KW capacity in place';
   LABEL X4 = 'Annual load factor';
   LABEL X5 = 'Peak kWh demand growth from 1974 to 1975';
   LABEL X6 = 'Sales (kWh use per year)';
   LABEL X7 = 'Percent Nuclear';
   LABEL X8 = 'Total fuel costs (cents per kWh)';
RUN;

PROC STANDARD DATA=UTILITIES MEAN=0 STD=1 OUT=TRAN_UTIL;
   VAR X1-X8;
RUN;

PROC CLUSTER M=WARD NONORM OUT=SEEDS;
   VAR X1-X8;
RUN;

PROC TREE N=5;

PROC FASTCLUS DATA=TRAN_UTIL MAXC=4 MEAN=NEW OUT=FOURCLUS;
   VAR X1-X8;
RUN;

PROC FASTCLUS SEED=NEW MAXC=4 OUT=FOURCLUS;
   VAR X1-X8;
RUN;

PROC ANOVA;
   CLASS CLUSTER;
   MODEL X1-X8 = CLUSTER;
RUN;

PROC SORT;
   BY CLUSTER;
RUN;

PROC PRINT;
   VAR CLUSTER COMPANY;
RUN;