/* MARK 9650 Project 1: Cluster Analysis on Cereal data				  */
/* WITH SELECTED OUTPUT												  */

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

DATA BREAKFAST;
   INFILE 'A:/CEREAL.PRN';
   INPUT BRAND $ MANUFACTURER $ X1-X8;
RUN;
													 
PROC STANDARD DATA=BREAKFAST MEAN=0 STD=1 OUT=TRAN_BREAK;
   VAR X1-X8;
RUN;
														
PROC CLUSTER M=WARD NONORM OUT=SEEDS;
   VAR X1-X8;
RUN;

PROC TREE N=5;

PROC FASTCLUS DATA=TRAN_BREAK MAXC=3 MEAN=NEW OUT=THREECLUS;
   VAR X1-X8;
RUN;

PROC FASTCLUS SEED=NEW MAXC=3 OUT=THREECLUS;
   VAR X1-X8;
RUN;

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

PROC SORT;
   BY CLUSTER;
RUN;

PROC PRINT;
   VAR CLUSTER MANUFACTURER BRAND;
RUN;