SAS的有序样品聚类程序,数据知道怎么填?但是后面的聚类图并没有出现结果

看着心好痛,找到这个程序的时候,很开心,但是数据输入进去后,结果并没有如期出现啊?至于类直径是对的,聚类图也没有出现,求指教?
如下:

现在把胡良平教授的教材《面向问题的统计学3》中的范例奉上,仅供参考。
本例一个分类指标,多个分类指标也是可以的。

Data mxwttjxt10_1;
input x@@;
cards;
9.3 1.8 1.9 1.7 1.5 1.3 1.4 2.0 1.9 2.3 2.1
;
run;
proc iml;
use mxwttjxt10_1;
read all into matrix_a;
close mxwttjxt10_1;
n=nrow(matrix_a);
call symputx(“n”,nrow(matrix_a));
call symputx(“m”,ncol(matrix_a));

/一、产生直径矩阵/
D=j(n,n,0); /建立n行n列的矩阵,元素初值均为0,以用来存储类的直径/
do i=1 to n;
do j=1 to n;
if i<j then D[i,j]=.; /用以得到下三角矩阵形式的直径矩阵/
else do;
temsum=matrix_a[i:j,]; /选取第i行到第j行作为子矩阵/
temsum2=temsum-shape(temsum[:,],i-j+1,&m); /用子矩阵中各元素减去本列元素的均值,得到离均差/
D[i,j]=round(ssq(temsum2),0.001); /

  • 0
    点赞
  • 6
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值