/*dieser File dient der Clusteranalyse von Haushaltstätigkeiten. Zuvor muss mit
Praepdistanzsingle o.a. eine Distanzmatrix produziert werden. Diese steht dann in
SASUSER als SINGLEDISTANZ oder anderem Namen und wird zunächst übernommen*/

DATA SINGLEDISTANZ1;
SET SASUSER.SINGLEDISTANZ91_92;
RUN;

PROC IML;
USE SINGLEDISTANZ1;
READ ALL VAR _NUM_ INTO DIST;
DIST1 = DIST+1000*I(911);
HILF = {[911] 1}` * {[911] 1} -I(911);
DIST2 = (DIST-120*HILF);
NULL = {[911] 0}` * {[911] 0};
DIST3 = DIST2<>NULL;
DIST4 = DIST3##2;

MIN = DIST1[><,><];
PRINT MIN;

CREATE SINGLEDISTANZ2 FROM DIST4;
APPEND FROM DIST4;
QUIT;

DATA SINGLES;
SET SASUSER.SINGLETAG1TAKT;
RUN;

DATA IDENT;
SET SINGLES;
KEEP PERSNR;
RUN;

DATA DISTANZ (TYPE=DISTANCE);
MERGE IDENT SINGLEDISTANZ2;
/*Je nachdem ob das IML zugeschaltet ist, muss hier zwischen
Singledistanz2 und Singledistanz1 gewechselt werden*/
RUN;

PROC CLUSTER DATA=DISTANZ METHOD=Average OUTTREE=TREE PSEUDO /*NOPRINT*/ ;
ID PERSNR;
RUN;

PROC TREE DATA=TREE OUT=NEW NCLUSTERS=14;
ID PERSNR;
RUN;

PROC SORT;
BY PERSNR;
RUN;

DATA SINGLECLUS;
MERGE NEW SINGLES;
BY PERSNR;
DROP /*ZHA1-ZHA288*/ IDTAG IDHH WT1 WT2 ID1 IDTAG;
IF ALTER>100 then delete;
IF CLUSTER=1;
RUN;

PROC FREQ DATA=NEW;
TABLES CLUSTER;
RUN;