聚类分析 Q型 R型

对样本进行Q型聚类分析

对指标进行R型聚类分析

• 原理
• MATLAB代码

% PA_R.m
% 对指标进行R型聚类分析
clc,clear
a=load('data.txt');
b=zscore(a);%数据标准化
r=corrcoef(b);%计算相关系数矩阵
d=pdist(b','correlation');%计算相关系数导出的距离
z=linkage(d,'average');%按平均法聚类
h=dendrogram(z)%画聚类图
T=cluster(z,'maxclust',5);%把变量划分为5类
for i=1:5
    tm=find(T==i);%求第i类的对象
    tm=reshape(tm,1,length(tm));%变成行向量
    fprintf('第%d类的有%s\n',i,int2str(tm));%显示分类结果
end

• 结果：

第1类的有5
第2类的有1  2  3  4  6
第3类的有8
第4类的有9
第5类的有7

因此，可取1,5,7,8,9指标为代表

参考资料

• [1] 司守奎，孙兆亮.数学建模算法与应用(第2版)北京：国防工业出版社，2019：216-231.