Documentation Example 3 for PROC MODECLUS
/****************************************************************/
/* S A S S A M P L E L I B R A R Y */
/* */
/* NAME: modecex3 */
/* TITLE: Documentation Example 3 for PROC MODECLUS */
/* PRODUCT: STAT */
/* SYSTEM: ALL */
/* KEYS: Cluster Analysis, Test on Number of Cluster */
/* PROCS: MODECLUS, SGPLOT */
/* DATA: */
/* */
/* SUPPORT: saswfk UPDATE: July 25, 2010 */
/* REF: PROC MODECLUS, Example 3 */
/* MISC: */
/* */
/****************************************************************/
title 'Modeclus Analysis with the JOIN= option';
title2 'A Normal Cluster Surrounded by a Ring Cluster';
data circle; keep x y;
c=1;
do n=1 to 30;
x=rannor(5);
y=rannor(5);
output;
end;
c=2;
do n=1 to 300;
x=rannor(5);
y=rannor(5);
z=rannor(5)+8;
l=z/sqrt(x**2+y**2);
x=x*l;
y=y*l;
output;
end;
run;
proc modeclus data=circle m=1 r=1 to 3.5 by .25 join=20 short;
run;
proc modeclus data=circle m=1 r=2.5 join out=out;
run;
proc sgplot data=out noautolegend;
yaxis values=(-10 to 10 by 5);
xaxis values=(-15 to 15 by 5);
scatter y=y x=x / group=cluster Markerchar=cluster;
by _NJOIN_;
run;