聚类参考介绍.docVIP

聚类分析 PROC VARCLUS 语句用来调用VARCLUS过程. 其后的选项可以有输入输出的数据 集选项，有关分类个数的选项，聚类方式选项(包括主成分法、重心法)等等，VAR语句列出被分析的变量. 现有北京地区1951～1976年的气象资料，该资料见例9.1. 要用重心法对变量Dec,Jan,Feb,High7,High4,High8进行聚类可采用以下SAS程序： data tepmerat input year Dec Jan Feb high7 high4 high8; cards; 1951 1.0 -2.7 -4.3 4 -7 12 1952 -5.3 -5.9 -3.5 0 21 5 1953 -2.0 -3.4 -0.8 6 -9 5 1954 -5.7 -4.7 -1.1 10 17 6 1955 -0.9 -3.8 -3.1 1 5 11 1956 -5.7 -5.3 -5.9 -3 1 -12 1957 -2.1 -5.0 -1.6 -15 3 13 1958 0.6 -4.3 -0.2 10 -3 0 1959 -1.7 -5.7 2.0 -9 -5 -14 1960 -3.6 -3.6 -1.3 11 -3 18 1961 -3.6 -3.1 -0.8 5 -15 4 1962 -0.1 -3.9 -1.1 8 12 1 1963 -2.6 -3.0 -5.2 11 3 -3 1964 -1.4 -4.9 -1.7 -11 -8 7 1965 -3.9 -5.7 -2.5 -18 6 -6 1966 -4.7 -4.8 -3.3 -9 -6 15 1967 -6.0 -5.6 -4.0 4 0 20 1968 -1.7 -6.4 -5.1 -7 -2 -15 1969 -3.4 -5.6 -2.0 4 17 -23 1970 -3.1 -4.2 -2.9 9 -16 23 1971 -3.8 -4.9 -3.9 -13 5 -2 1972 -2.0 -4.1 -2.4 7 0 10 1973 -1.7 -4.2 -2.0 27 -11 4 1974 -3.6 -3.3 -2.0 17 -2 0 1975 -2.7 -3.7 0.1 -1 -13 10 1976 -2.4 -7.6 -2.2 5 9 -30 ； proc varclus data=temperat centroid; var Dec Jan Feb high7 high4 high8; run; 执行程序后，输出分成1类、2类、3类、4类的分法，仅考虑分为3类的表. Oblipuc Centroid Component Cluster Analysis Cluster Summary for 3 cluster(s) Cluster Members Cluster Variation Proportion variation Explained Explained 1 3 3.00000 1.80281 0.6009 2 2 1.00000 1.00000 1.0000 3 1 2.00000 1.26520 0.6326 Total variation explained=4.068004 Proportion=0.6780 R-spuared with VARIABLE Own Cluster Next Closeat 1-R* *2Ratio Cluster1 JAN 0.8172 0.2006 0.2287 HIGH7 0.4383 0.0231 0.5749 HIGH8 0.5769 0.1499 0.4977 Cluster2 HIGH4 1.0000 0.1653 0.0000 Cluster3 DEC 0.6326 0.0801 0.3994 FFB 0.6326 0.0441 0.3844 以上可见JAN，HIGH7，HIGH8为一组，HIGH4为一组，其余为一组. 2.判别分析 已知长江中游6月份降水分三级:偏少、偏多、正常，分别用1，2，3来表示这三个母体. 现利用1951～1975年共25年的观测资料(表15.6.1)建立四因子(指标)的三类判别函数.并将1976～1979年的观测资料(表15.6.2)作为检验资料(主要是为了示范这类问题的SAS程序). 表15.6.1 历年观测资料 年份 X1 X2 X3 X4 原分类 1951 0.58 82.0 44.0 40.6 1 1952 0.40 83.0 18.0 43.0 2 1953 0.55 85.0 36.0 30.7 2 1954 0.40 85.0