实验报告八-SAS聚类分析与判别分析(8)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 2 1.0000 0.0000 0.9999 0.0001 1.0000 0.0000 0.9989 0.0011 1.0000 0.0000 1.0000 0.0000 0.9977 0.0023 1.0000 0.0000 1.0000 0.0000 1.0000 0.0000 1.0000 0.0000 1.0000 0.0000 1.0000 0.0000 1.0000 0.0000 0.9991 0.0009 0.9973 0.0027 0.9923 0.0077 1.0000 0.0000 1.0000 0.0000 0.9998 0.0002 1.0000 0.0000 1.0000 0.0000 0.9999 0.0001 1.0000 0.0000 0.9963 0.0037 0.9989 0.0011 1.0000 0.0000 0.9996 0.0004 0.0000 1.0000 0.0000 1.0000 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 . . . . . . . 2 2 2 2 2 2 2 2 2 1 2 2 2 2 2 2 2 1 2 2 2 2 2 1 1 1 2 1 1 1 * * * * * * * * * 0.0001 0.9999 0.0000 1.0000 0.0000 1.0000 0.0000 1.0000 0.0022 0.9978 0.0023 0.9977 0.0000 1.0000 0.0000 1.0000 0.0903 0.9097 0.7013 0.2987 0.0114 0.9886 0.3584 0.6416 0.0457 0.9543 0.0001 0.9999 0.0005 0.9995 0.0000 1.0000 0.0000 1.0000 0.8277 0.1723 0.0004 0.9996 0.0000 1.0000 0.0000 1.0000 0.0000 1.0000 0.0000 1.0000 0.9999 0.0001 1.0000 0.0000 1.0000 0.0000 0.0000 1.0000 0.9442 0.0558 0.9976 0.0024 0.9999 0.0001 61 62 63 64 65 66 67 . . . . . . . 2 2 1 2 2 1 2 * * * * * * * 0.2091 0.7909 0.0000 1.0000 0.9913 0.0087 0.0868 0.9132 0.0010 0.9990 1.0000 0.0000 0.1379 0.8621 * Misclassified observation 由表8.47组中成员的后验概率（ Posterior Probability of Membership in RaD）得知，第40个观测值被误判为第二类健将级，根据线性判别函数得出第40个观测值在第一、二类的后验概率分别为0.7013、0.2987；第48个观测值被误判为第二类健将级，根据线性判别函数得出第48个观测值在第一、二类的后验概率分别为0.8277、0.1723。并由线性判别函数得到，未定级的第54到67名运动员分别为一级、一级、一级、健将级、一级、一级、一级、健将级、健将级、一级、健将级、健将级、一级、健将级。 表8.48 Number of Observations and Percent Classified into Group Classification Summary for Calibration Data: LMF.P85 Cross-validation Summary using Linear Discriminant Function Number of Observations and Percent Classified into Group From Group . 1 8 57.14 2 6 Total 14 42.86 100.00 0 1 28 28 100.00 25 100.00 0.00 2 2 8.00 23 92.00 100.00 29 Total 38 56.72 67 43.28 100.00 0.5 Priors 0.5 利用交叉验证判别分类得出表8.48观察值频数及分类百分比（Number of Observations and Percent Classified into RaD），由表得知，第一类一级的样本被判别函数归入第一、二类的频数分别为28、0，百分比分别为100%、0%；已知第二类健将级的样本被判别函数归入第一、二类的频数分别为2、23，百分比分别为8%、92%。 综上，利用Bayes判别法得出未定级的第54到67名运动员分别为一级、一级、一级、健将级、一级、一级、一级、健将级、健将级、一级、健将级、健将级、一级、健将级。 ⑵当Σ1≠Σ2时： 利用proc discrim过程步实现Bayes判别分析： proc discrim data=Lmf.p85 pool=no crosslist; class Group; priors equal; var x1 x2 x3 x4 x5 x6; run; 结果：（由于简单信息与上述⑴过程中相等，因此不再列出。） 表8.49 Within Covariance Matrix Information Group Covariance Natural Log of the Matrix Rank Determinant of the Covariance Matrix 1 2 表8.49组内协方差矩阵信息（Within Covariance Matrix Information），协方差矩阵的秩（Covariance Matrix Rank）为6，协方差矩阵行列式的自然对数（Natural Log of the Determinant of the Covariance Matrix）分别为0.90637、6.42369。 表4.50 Generalized Squared Distance to Group From Group 1 2 1 0.90637 2 23.55603 6 6 0.90637 6.42369 122.36031 6.42369 表8.50广义组间距离平方（Generalized Squared Distance to RaD）两两配对的组间平方距离为 D2(1|1)?0.90637，D2(1|2)?23.55603， 2D2(2|1)?122.36031D，(2|2)?6.42369。 表8.51 Posterior Probability of Membership in Group Classification Results for Calibration Data: LMF.P85 Cross-validation Results using Quadratic Discriminant Function Obs From Group Classified into Group 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 1 1 1 1 1 1 1 1 1 1 * 1 2 1.0000 0.0000 0.9996 0.0004 0.9997 0.0003 0.9937 0.0063 1.0000 0.0000 1.0000 0.0000 0.9934 0.0066 1.0000 0.0000 1.0000 0.0000 0.9998 0.0002 1.0000 0.0000 1.0000 0.0000 1.0000 0.0000 0.9999 0.0001 0.9984 0.0016 0.9929 0.0071 0.0000 1.0000 1.0000 0.0000 1.0000 0.0000 0.9981 0.0019 1.0000 0.0000 0.9966 0.0034 0.9999 0.0001 0.9998 0.0002 0.9854 0.0146 0.9837 0.0163 1.0000 0.0000