对于有序样品的聚类,在SPSS和Matlab上没有找到现成的程序。故只好自己来进行实现。
代码如下:
optimal_partition2.m
function [a,s,d]=optimal_partition2(x,method)
%% 有序样品的最优二分割法
%% x为m*n的矩阵,n为样品个数,m为指标个数
%% method指定义样品之间差异的方法,1为段内离差和,2为段内极差
%% a为最后所得到的分割点
%% s为最小损失函数值
%% d为段内离差和矩阵或段内极差矩阵
[m,n]=size(x);
d=zeros(n,n);
for j=2:n
for i=1:j
if method==1
tx=(sum((x(:,i:j))'))'/(j-i+1);
xa=zeros(m,1);
for k=i:j
xa=xa+(x(:,k)-tx).*(x(:,k)-tx);
end
d(i,j)=sum(xa);
else
d(i,j)=sum(max((x(:,i:j))')-min((x(:,(i:j)))'));
end
end
end
for i=1:n-1
S(i)=d(1,i)+d(i+1,n);
end
s=min(S);
a=find(S==s);
其中,这里假定多个变量时的段内极差为 ,
有序样品聚类matlab程序
- 格式:txt
- 大小:0.93 KB
- 文档页数:1