OPTIONS PS=55 LS=70;
DATA ORTHO;
INFILE 'c:\hospital.txt' DELIMITER=',';

INPUT ZIP $  HID $ CITY $ STATE $ BEDS RBEDS OUTV ADM SIR SALESY 
	 SALES12 HIP95 KNEE95 TH TRAUMA REHAB HIP96 
	 KNEE96 FEMUR96;

IF STATE EQ 'FL' OR STATE EQ 'GA';
ARRAY X {12} BEDS RBEDS HIP95 KNEE95 HIP96 KNEE96 FEMUR96
         OUTV ADM SIR SALESY SALES12;
DO I=1 TO 7 ;
	 X{I} = SQRT(X{I});
END;
DO I=8 TO 12 ;
	 X{I} = LOG(1+X{I});
END;
RUN;

PROC CLUSTER METHOD=WARD;
VAR BEDS RBEDS OUTV ADM SIR  
	 HIP95 KNEE95 TH TRAUMA REHAB HIP96 
	 KNEE96 FEMUR96;
COPY BEDS RBEDS OUTV ADM SIR  
	 HIP95 KNEE95 TH TRAUMA REHAB HIP96 
	 KNEE96 FEMUR96 SALES12 SALESY  ;

PROC TREE NOPRINT NCL=7 OUT=TXCLUST;
COPY BEDS RBEDS OUTV ADM SIR  
	 HIP95 KNEE95 TH TRAUMA REHAB HIP96 
	 KNEE96 FEMUR96 SALES12 SALESY  ;
RUN;

PROC SORT DATA=TXCLUST; BY CLUSTER;

PROC MEANS DATA=TXCLUST;
VAR      BEDS RBEDS OUTV ADM SIR  
	 HIP95 KNEE95 TH TRAUMA REHAB HIP96 
	 KNEE96 FEMUR96 SALES12 SALESY  ;
BY CLUSTER;
OUTPUT   OUT = C 
         MEAN = MBEDS MRBEDS MOUTV MADM MSIR
         MHIP95 MKNEE95 MTH MTRAUMA MREHAB MHIP96
         MKNEE96 MFEMUR96 MSALES12 MSALESY  ;

PROC PRINT;
VAR CLUSTER MBEDS MRBEDS MOUTV MADM MSIR
         MHIP95 MKNEE95 MTH MTRAUMA MREHAB MHIP96
         MKNEE96 MFEMUR96 MSALES12 MSALESY  ;
RUN;

PROC UNIVARIATE FREQ DATA=TXCLUST;
VAR SALES12;
BY CLUSTER;

DATA C;
INPUT NCLUS R2CHANGE;
CARDS;
16 0.005127
15 0.006717
14 0.007246
13 0.007860
12 0.008437
11 0.008496
10 0.008555
9 0.009749
8 0.009836
7 0.014713
6 0.032362
5 0.037383
4 0.049615
3 0.063334
2 0.114961
;

PROC PLOT;
PLOT R2CHANGE*NCLUS;
RUN;