6-10今有6个铅弹头,用“中子活化”方法测得7种微量元素的含量数据(见表1)
(1)试用多种系统聚类法对6个弹头进行分类;并比较分类结果;(2)试用多种方法对7种微量元素进行分类
表1微量元素含量数据元素样品号Ag(银)(X1)Al(铝)Cu(铜)Ca(钙)Sb(锑)Bi(铋)Sn(锡)(X2)(X3)(X4)(X5)(X6)(X7)10
057985
5150347
9108586174261
084413
9700347
71079472000244030
072171
153054
05238601445949740
150101
7020307
030122901461638055
744002
8540229
657809912661252060
213000
7058240
910898028204135问题求解1对6个弹头进行分类对数据进行标准化变换,样品间距离定义为欧式距离,系统聚类的方法分别使用类平均法(AVE)、中间距离法(MID)、可变类平均法(FLE)和离差平方合法(WARD)
使用SAS软件CLUSTER过程对数据进行聚类分析(程序见附录1)
1类平均法图1类平均聚类法相关矩阵特征值图图2类平均聚类分析法聚类历史图由图2可知,NCL=1时半偏R2最大且伪F统计量在NCL=2,5时和伪t方统计量在NCL=1,4时较大
因此,将6个弹头分为两类(2)(2)121,2,4,6,3,5GG
SAS绘制的谱系聚类图如图3所示
图3类平均聚类分析法谱系聚类图1
2中间距离法图4中间距离聚类法相关矩阵特征值图图5中间距离聚类法聚类历史图由图5可知,中间距离法与类平均法结果一致
因此,也将6个弹头分为两类(2)(2)121,2,4,6,3,5GG
SAS绘制的谱系聚类图如图6所示