SAS中的聚类分析方法总结.docx
《SAS中的聚类分析方法总结.docx》由会员分享，可在线阅读，更多相关《SAS中的聚类分析方法总结.docx（13页珍藏版）》请在冰豆网上搜索。
SAS中的聚类分析方法总结
说起聚类分析，相信很多人并不陌生。
这篇原创博客我想简单说一下我所理解的聚类分析，欢迎各位咼手不吝赐教和拍砖。
按照正常的思路，我大概会说如下几个问题：
1.什么是聚类分析？
2.聚类分析有什么用？
3.聚类分析怎么做？
下面我将分聚类分析概述、聚类分析算法及sas实现、案例三部分来系统的回答这些问题。
聚类分析概述
1-聚类分析的定义
中国有句俗语叫“物以类聚，人以群分”一一剔除这句话的贬义色彩。
说白了就是物品根据物品的特征和功用可以分门别类，人和人会根据性格、偏好甚至利益结成不同的群体。
分门别类和结成群体之后，同类（同群）之间的物品（人）的特征尽可能相似，不同类（同群）之间的物品（人）的特征尽可能不同。
这个过程实际上就是聚类分析。
从这个过程我们可以知道如下几点：
1）聚类分析的对象是物（人），说的理论一点就是样本
2）聚类分析是根据物或者人的特征来进行聚集的，这里的特征说的理论一点就是变量。
当然特征选的不一样，聚类的结果也会不一样；
3）聚类分析中评判相似的标准非常尖键。
说的理论一点也就是相似性的度量非常矢键；
4）聚类分析结果的好坏没有统一的评判标准；
2.聚类分析到底有什么用？
1）说的官腔一点就是为了更好的认识事物和事情，比如我们可以把人按照地域划分为南方人和北方人，你会发现这种分法有时候也蛮有道理。
一般来说南方人习惯吃米饭，北方习惯吃面食；
2）说的实用一点，可以有效对用户进行细分，提供有针对性的产品和服务。
比如银行会将用户分成金卡用户、银卡用户和普通卡用户。
这种分法一方面能很好的节约银行的资源，另外一方面也能很好针对不同的用户实习分级服务，提高彼此的满意度。
再比如移动会开发全球通、神州行和动感地带三个套餐或者品牌，实际就是根据移动用户的行为习惯做了很好的用户细分一一聚类分析；
3）上升到理论层面，聚类分析是用户细分里面最为重要的工具，而用户细分则是整个精准营销里面的基础。
精准营销是目前普遍接纳而且被采用的一种营销手段和方式。
3.聚类分析的流程是怎样的？
比较简单的聚类分析往往只根据一个维度来进行，比如讲用户按照付费情况分成高端用户、中端用户和低端用户。
这个只需要根据商业目的统计一下相尖数据指定一个高端、中端和低端的分界点标准就可以。
如果是比较复杂的聚类分析，比如移动里面经常会基于用户的多种行为（通话、短信、gprs
流失扌屯等）来对用户进行细分，这个就是比较复杂的用户细分。
如果是这样的细分通常会作为一个比较标准的数据挖掘项目来执行，所以基本上会按照数据挖掘的流程来执行。
具体分如下几步：
1）业务理解主要是了解业务目标和数据挖掘的目标及执行计划
2）数据理解主要是弄清楚可已取哪些变量数据，具体怎么定义
3）数据整理
根据之前的定义提取需要的数据，并进行检测异常数据，并对变量进行挑选及探索，比如最终要用那些变量来执行聚类算法、那些变量是离散变量，需要做特殊处理、数据大概可以聚成几类、类别形状有不规则的情形吗？
4）建立模型矢键是选用什么样的距离（相似性度量）和算法：
I比如是样本比较小，形状也比较规则，可以选用层次聚类
I比如样本比较大，形状规则，各类的样本量基本相当，可以选用k-means算法
I比如形状规则，但是各类别之间的样本点的密度差异很大，可以选用基于密度的算法
5）模型评估
主要是评估聚类分析结果的好坏。
实际上聚类分析在机器学习里面被称之为无监督学习，是没有大家公认的评估方法的。
所以更多会从业务可解释性的角度去评估
聚类分析的好坏；
6）模型发布主要是根据聚类分析的结果根据不同的类的特诊去设计不同的产品、服务或者渠道策略，然后去实施营销
4.具体在sas里面如何执行？
通过前面的讲解我们已经知道，聚类分析涉及到如下6步，对应着6步SAS都会有相应的过程来执行。
1）距离的计算：
procdistance
2）数据标准化：
procstdize
3）聚类变量的选择：
procvarclus
4）初始类另U数的选择：
procmds禾口procprincomp
5）不规则形状的变换：
procaceclus
6）算法的选择：
层次聚类-proccluster划分型聚类・procfastclus（k-means）和
密度型聚类-procmodeclus
7）类另特征描述：
procmeans
以上四个部分就从是什么、为什么、怎么样三个角度对聚类分析做了简单的介绍。
接下来的帖子我会重点介绍SAS中各种聚类算法的差异、应用范围及实际的案例。
5.用procdistance做什么？
我们知道数据变量分四类：
名义变量、次序变量、interval变量和ritio变量。
但sas里面目
前的聚类算法都要求变量时ratio变量。
那想要对离散变量进行聚类怎么呢？
一种想法自然是讲所有的离散变量都转成0・1变量。
这会有如下几个问题：
1）变量的信息可能会有损失，比如次序型变量转成0・1变量后，次序信息就很难保留；
2）当离散变量的取值非常多时，转成0・1变量后生成的新变量也会非常多，这样也会造
成很多处理上的不便；
3）0-1变量也没法做标准化等等一些运算，因为这种运算其实是没有意义的那该如何处理离散变量的聚
类呢？
答案是用procdistance。
我们知道聚类过程中首先是从计
算距离或者相似度幵始的。
一个很自然的想法就是针对离散变量定义有意义的距离（对离散变量和连续变量混合类型的数据）。
Procdistance就是用来算这种距离的一个很好的过程。
距离或者相似度可以看成是连续数据，自然就可以用sas里面的聚类算法了。
6.用procstdize做什么？
前面说过聚类算法首先要算的距离，然后通过距离来执行后续的计算。
在距离计算的过程方差比较大的变量影响会更大，这个通常不是我们希望看到。
所以非常有必要讲参与聚类的变量转换成方差尽量相同。
Procstdize就能实现这种功能。
Procstdize不仅提供了将变量转换了均值为0，方差转换为1的标准化，还提供了很多其它类型的标准化。
比如，range标准化（变量减去最小值除以最大值和最小值得差）
7.用procvarclus做什么？
在做回归分析的时候，我们知道变量过多会有两个问题：
1）变量过多会影响预测的准确，尤其当无矢紧要的变量引入模型之后；
2）变量过多不可避免的会引起变量之前的共线性，这个会影响参数估计的精度聚类分析实际上也存在类似的问题，所以有必要先对变量做降维。
说到降维，马上有人会说这个可以用主成分啊，这个的确没
错。
但是主成分的解释性还是有点差。
尤其是第二主成分之后的主成分。
那用什么比较好呢？
答案是
procvarclus斜交主成分。
我们常说的主成分实际上正交主成分。
斜交主成分是在正交主成分的基础上再做了一些旋转。
这样得到的主成分不仅能保留主成分的优点（主成分变量相尖程度比较低）。
另外一方面又
能有很到的解释性，并且能达到对变量聚类的效果。
使同类别里面的变量尽可能相矢程度比较高，不同类别里面的变量相尖程度尽可能低。
这样根据一定的规则我们就可以在每个类别里面选取一些有代表性的变量，这样既能保证原始的数据信息不致损失太多，也能有效消除共线性。
有效提升聚类分析的精度。
8.用procmds禾口procprincomp做什么？
将原始数据降到两维，通过图形探测整个数据聚类后大致大类别数
9.用procaceclus做什么？
聚类算法尤其是k-means算法要求聚类数据是球形数据。
如果是细长型的数据或者非凸型数据，这些算法的表现就会相当差。
一个很自然的变通想法就是，能不能将非球形数据变换成球形数据呢？
答案是可以的。
这就要用到procaceclus。
10.标准化对聚类分析到底有什么影响？
1）在讲影响之前先罗列一下procstdize里面的标准化方法吧
PROC5TDIZE过程里面的标准化方法（表1）
METHOD
LOCATION
SCALE
MEAN
mean
1
MEDIAN
median
1
SUM
Q
Sum
EUCLEN
0
EuclideanLength
USTD
0
standarddeviationaboutorigin
STD
mean
standarddevj31k>n
RANGE
mirnrrwfn
T占nqe
MIDRANGE
midrange
MAXABS
0
n>axinium3bwls亡value
IQR
median
knterquanikerange
MAD
median
niedunabsolutedeviationfrommedian
ABW（c）
hAe：
aht1-steoM-estimate
biAlightAAesbnvite
AHUBER（c）
Hubert・step
TuberA*estim3te
AWAVE（c）
Wave1-stepM-esItmaie
WaveA-estimate
AGK（p）
mean
AGKesDFnate（ACECLUSJ
SPACINGfp）
midrnininKinrbspacinq
rnirirnumspacing
Up）
L（P）
L（p）（Mmkowskidtslan匚总s）
IN（dS）
readfromdatasei
readfromdataset
注：
具体计）1方法可查圜PROC寸D!
在过程详堀希助文榭
2）标准化对聚类分析的影响
从图1中不太容易看清楚标准化对于聚类分析的影响
标准化对于聚类的影响（图1）
从图2可以清晰的看到标准化对于聚类分析的影响
标准化对于聚类的影响（图2）
Plot1
注：
Plot2只是在的基就上对扇始埶据的）縮小了LOOO倍・y旷犬了1DD倍得到的
3）各种标准化方法的比较
一个模拟数据的例子，模拟数据有三个类别，每个类别有100个样本。
我们比较了各种标准
化方法之后再进行聚类的误判情况，可以大概看出各种标准化方法的差异。
但此例并不能说
明以下方法中误分类数小的方法就一定优与误分类数大的方法。
有时候还跟数据本身的分布
特征有矢。
这个例子也提醒我们有时候我们常用的std和range标准化并不见得是最好的选
择。
各种标准化方法得到的错分样本数（：
表2）
汇总麵握
聚英判足実邱蓋泯样*瞰st序
Obs
metihod
misc1ass•faud
cNiwq
1
EUCLEN
14
520.494
2
MEAN
14
520,494
3
MEDIAN
14
520.494
4
HONE
14
520.494
&
SPAC1NGC.3>
14
520.494
6
USTD
14
520.494
7
1QR
15
&15336
8
L<1）
15
515336
9
LC2）
IS
515.336
IO
M1DRANGE
15
515.336
11
RANGE
15
515.336
12
STD
16
515.336
13
SUM
16
515336
14
MAXABS
17
505.296
IS
AWAVE<25）
44
375.136
16
AGK-9）
46
228.364
17
AHUBER<・5>
54
247.639
18
ABW<.5〉
55
329672
IO
MAD
611
251.518
附：
本节相矢sas代码就作为回帖吧。
广1•模拟数据1;测试标准化方法对聚类的影响模拟数据，样本量相同，均值和方差不相同*/
datacompact;
keepxyc;
n=100;
scale=1;mx=0;my=0;c=1;linkgenerate;scale=2;mx=8;my=0;c=2;linkgenerate;scale=3;mx=4;my=8;c=3;linkgenerate;stop;
generate:
doi=1ton;
x=rannor
（1）*scale+mx;
y=rannor（1广scale+my;output;end;
return;
run;
title，模拟数据V;
procgplotdata=compact;ploty*x=c;
symbollc=blue;
symbol2c=black;
symbol3c=red;
run;
procstdizedata=compactmethod=stdout=scompacted2;
varxy;
run;
title'标准化后的模拟数据1'；procgplotdata=scompacted2;ploty*x=c;
symbollc=blue;
symbol2c=black;
symbol3c=red;
run;
**//*2.createresulttable*/
★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★it
*7dataresult;
lengthmethod$12;lengthmisclassified8;
lengthchisq8;
stop;
%letinputs=xy;
%letgroup=c;
%macrostandardize（dsn=,nc=,method=）;title"&methocT;
%if%bquote（%upcase（&method））=NONE%then%do;datatemp;
set&dsn;
run;
%end;
%else%do;
procstdizedata=&dsnmethod=&methodout=temp;var&inputs;
run;
%end;
procfastclusdata=tempmaxclusters=&ncleast=2out=clusoutnoprint;var&inputs;run;
procfreqdata=clusout;
tables&group*cluster/norownocolnopercentchisqout=freqout;outputout=statschisq;
run;
datatempsum;
setfreqoutend=eof;
by&group;retainmembersmodec;iffirst.&groupthendo;members=0;mode=0;end;members=members+count;
ifclusterNE.thendo;ifcount>modethendo;mode=count;
c=cluster;
end;
end;
iflast.&groupthendo;cum+（members-mode）;outputtemp;
end;
ifeofthenoutputsum;
run;
procprintdata=tempnoobs;var&groupcmembersmodecum;run;
dataresult;
mergesum（keep=cum）stats;if0thenmodifyresult;method=,,&methodn;misclassified=cum;chisq=_pchi_；
pchisq=p_pchi;
outputresult;
run;
%mendstandardize;
%standardize（dsn=compact5nc=33method=ABW（.5））;
%standardize（dsn=compact5nc=3,method=AGK（.9））;
%standardize（dsn=compact3nc=3,method=AHUBER（.5））;
%standardize（dsn=compact,nc=3,method=AWAVE（.25））;
%standardize（dsn=compact,nc=3,method=EUCLEN）;
%standardize（dsn=compact5nc=3,method=IQR）;
%standardize（dsn=compact,nc=3,method=：
L
（1））;
%standardize（dsn=compact,nc=3,method=L
（2））;
%standardize（dsn=compact,nc=3,method=MAD）;
%standardize（dsn=compact5nc=3,method=MAXABS）;
%standardize（dsn=compact,nc=3,methocl=MEAN）;
%standardize（dsn=compact5nc=3,method=MEDIAN）;%standardize（dsn=compact,nc=3,method=MIDRANGE）;
%standardize（dsn=compact,nc=3,methocl=：
NONE）;
%standardize（dsn=compact5nc=3,method=RANGE）;%standardize（dsn=compact,nc=3,method=SPACING（•3））；
%standardize（dsn=compact,nc=3,method=STD）;
%standardize（dsn=compact5nc=3,method=SUM）;
%standardize（dsn=compact5nc=3,method=USTD）;
procsortdata=result;
bymisclassified;
run;
title1汇总数据:
title2'聚类判定类别错误样本数排序：
procprintdata=result;
run;
升级会员 