
院 系:数学与统计学学院
专 业:__ _统计学
年 级: 2009级
课程名称: 统计分析
学 号:
姓 名:
指导教师:
2012年 5月 14 日
(一)实验名称
1.多元正态分布统计推断page115例题4.2.1
2.聚类分析
(二)实验目的
1.掌握简单多元正态分布统计推断问题。
2.掌握运用SAS编程对数据进行聚类分析。
(三)实验数据
第一题:
| 编号 | 身高 | 胸围 | 上半臂围 |
| 1 | 78 | 60.6 | 16.5 |
| 2 | 76 | 58.1 | 12.5 |
| 3 | 92 | 63.2 | 14.5 |
| 4 | 81 | 59 | 14.0 |
| 5 | 81 | 60.8 | 15.5 |
| 6 | 84 | 59.5 | 14.0 |
| city | x1 | x2 | x3 | x4 | x5 | x6 | x7 | x8 |
| 阿根廷 | 10.39 | 20.81 | 46.84 | 1.81 | 3.7 | 14.04 | 29.36 | 137.72 |
| 澳大利亚 | 10.31 | 20.06 | 44.84 | 1.74 | 3.57 | 13.28 | 27.66 | 128.3 |
| 奥地利 | 10.44 | 20.81 | 46.82 | 1.79 | 3.6 | 13.26 | 27.72 | 135.9 |
| 比利时 | 10.34 | 20.68 | 45.04 | 1.73 | 3.6 | 13.22 | 27.45 | 129.95 |
| 百莫大 | 10.28 | 20.58 | 45.91 | 1.8 | 3.75 | 14.68 | 30.55 | 146.62 |
| 巴西 | 10.22 | 20.43 | 45.21 | 1.73 | 3.66 | 13.62 | 28.62 | 133.3 |
| 缅甸 | 10. | 21.52 | 48.3 | 1.8 | 3.85 | 14.45 | 30.28 | 139.95 |
| 加拿大 | 10.17 | 20.22 | 45.68 | 1.76 | 3.63 | 13.55 | 28.09 | 130.15 |
| 智利 | 10.34 | 20.8 | 46.2 | 1.79 | 3.71 | 13.61 | 29.3 | 134.03 |
| 中国 | 10.51 | 21.04 | 47.3 | 1.81 | 3.73 | 13.9 | 29.13 | 133.53 |
| 哥伦比亚 | 10.43 | 21.05 | 46.1 | 1.82 | 3.74 | 13.49 | 27.88 | 131.35 |
| 库克群岛 | 12.81 | 23.2 | 52.94 | 2.02 | 4.24 | 16.7 | 35.38 | 1.7 |
| 哥斯达黎加 | 10.94 | 21.9 | 48.66 | 1.87 | 3.84 | 14.03 | 28.81 | 136.58 |
| 捷克斯洛伐克 | 10.35 | 20.65 | 45. | 1.76 | 3.58 | 13.42 | 28.19 | 134.32 |
| 丹麦 | 10.56 | 20.52 | 45. | 1.78 | 3.61 | 13.5 | 28.11 | 130.78 |
| 多米老家共和国 | 10.43 | 20.69 | 45.49 | 1.74 | 3.61 | 13.27 | 27.52 | 130.87 |
| 法国 | 10.11 | 20.38 | 45.28 | 1.73 | 3.57 | 13.34 | 27.97 | 132.3 |
| 德意志民主共和国 | 10.12 | 20.33 | 44.87 | 1.73 | 3.56 | 13.17 | 27.42 | 129.92 |
| 德意志联邦共和国 | 10.16 | 20.37 | 44.5 | 1.73 | 3.53 | 13.21 | 27.61 | 132.23 |
| 大不列颠及北爱尔兰 | 10.11 | 20.21 | 44.93 | 1.7 | 3.51 | 13.01 | 27.51 | 129.13 |
| 希腊 | 10.22 | 20.71 | 46.56 | 1.78 | 3. | 14.59 | 28.45 | 134.6 |
| 危地马拉 | 10.98 | 21.82 | 48.4 | 1. | 3.8 | 14.16 | 30.11 | 139.33 |
| 匈牙利 | 10.26 | 20.62 | 46.02 | 1.77 | 3.62 | 13.49 | 28.44 | 132.58 |
| 印度 | 10.6 | 21.42 | 45.73 | 1.76 | 3.73 | 13.77 | 28.81 | 131.98 |
| 印度来西亚 | 10.59 | 21.49 | 47.8 | 1.84 | 3.92 | 14.73 | 30.79 | 148.83 |
| 以色列 | 10.61 | 20.96 | 46.3 | 1.79 | 3.56 | 13.32 | 27.81 | 132.35 |
| 爱尔兰 | 10.71 | 21 | 47.8 | 1.77 | 3.72 | 13.66 | 28.93 | 137.55 |
| 意大利 | 10.01 | 19.72 | 45.26 | 1.73 | 3.6 | 13.23 | 27.52 | 131.08 |
| 日本 | 10.34 | 20.81 | 45.86 | 1.79 | 3. | 13.41 | 27.72 | 128.63 |
| 肯尼亚 | 10.46 | 20.66 | 44.92 | 1.73 | 3.55 | 13.1 | 27.38 | 129.75 |
| 韩国 | 10.34 | 20. | 46.9 | 1.79 | 3.77 | 13.86 | 29.23 | 136.25 |
| 朝鲜人民共和国 | 10.91 | 21.94 | 47.3 | 1.85 | 3.77 | 14.13 | 29.67 | 130.87 |
| 卢森堡 | 10.35 | 20.77 | 47.4 | 1.82 | 3.67 | 13. | 29.08 | 141.27 |
| 马来西亚 | 10.4 | 20.92 | 46.3 | 1.82 | 3.8 | 14. | 31.01 | 154.1 |
| 毛里求斯 | 11.19 | 22.45 | 47.7 | 1.88 | 3.83 | 15.06 | 31.77 | 152.23 |
| 墨西哥 | 10.42 | 21.3 | 46.1 | 1.8 | 3.65 | 13.46 | 27.95 | 129.2 |
| 荷兰 | 10.52 | 20.95 | 45.1 | 1.74 | 3.62 | 13.36 | 27.61 | 129.02 |
| 新西兰 | 10.51 | 20.88 | 46.1 | 1.74 | 3.54 | 13.21 | 27.7 | 128.98 |
| 挪威 | 10.55 | 21016 | 46.71 | 1.76 | 3.62 | 13.34 | 27.69 | 131.48 |
| 卢布亚新几内亚 | 10.96 | 27.78 | 47.9 | 1.9 | 4.01 | 14.72 | 31.36 | 148.22 |
| 菲律宾 | 10.78 | 21. | 46.24 | 1.81 | 3.83 | 14.74 | 30. | 145.27 |
| 波兰 | 10.16 | 20.24 | 45.36 | 1.76 | 3.6 | 13.29 | 27. | 131.58 |
| 葡萄牙 | 10.53 | 21.17 | 46.7 | 1.79 | 3.62 | 13.13 | 27.38 | 128.65 |
| 罗马尼亚 | 10.41 | 20.98 | 45.87 | 1.76 | 3. | 13.25 | 27.67 | 132.5 |
| 新加坡 | 10.38 | 21.28 | 47.4 | 1.88 | 3. | 15.11 | 31.32 | 157.77 |
| 西班牙 | 10.42 | 20.77 | 45.98 | 1.76 | 3.55 | 13.31 | 27.73 | 131.57 |
| 瑞士 | 10.25 | 20.61 | 45.63 | 1.77 | 3.61 | 13.29 | 27.94 | 130.63 |
| 瑞典 | 10.37 | 20.46 | 45.78 | 1.78 | 3.55 | 13.22 | 27.91 | 131.2 |
| 中国台北 | 10.59 | 21.29 | 46.8 | 1.79 | 3.77 | 14.07 | 30.07 | 139.27 |
| 泰国 | 10.39 | 21.09 | 47.91 | 1.83 | 3.84 | 15.23 | 32.56 | 149.9 |
| 土耳其 | 10.71 | 21.43 | 47.06 | 1.79 | 3.67 | 13.56 | 28.58 | 131.5 |
| 美国 | 9.93 | 19.75 | 43.86 | 1.73 | 3.53 | 13.2 | 27.43 | 128.22 |
| 苏联 | 10.07 | 20 | 44.6 | 1.75 | 3.59 | 13.2 | 27.53 | 130.55 |
| 西萨摩亚 | 10.82 | 21.86 | 49 | 2.02 | 4.24 | 16.28 | 34.71 | 161.83 |
1. page115例题4.2.1
2. 打开SAS软件并导入数据;
3.聚类分析;
4.做聚类树状图;
(五)分析实验结果
(六)实验体会(结论、评价、感谢与建议等)
(七)成绩评定
第一题:
程序如下:
data page45;
label no='序号' x1='身高' x2='胸围' x3='上半臂围';
input no class x1-x3 @@;
cards;
1 1 78 60.6 16.5
2 1 76 58.1 12.5
3 1 92 63.2 14.5
4 1 81 59 14.0
5 1 81 60.8 15.5
6 1 84 59.5 14.0
7 0 90 58 16
;
quit;
proc ttest data=page45 cochran;
class class;
var x1 x2 x3;
quit;
proc discrim method=normal pool=test manova;
class class;
var x1 x2 x3;
run;
proc glm data=page45;
class class ;
model x1 x2 x3=class/nouni;
manova h=class;
run;
结果如下:
方差检验结果可以看出,所有变量的p值都小于0.05,拒绝原假设,认为农村与城市的2周岁男婴上述三个指标的均值有显著差异。(a=0.05)
第二题:
程序如下(选取聚类分析1.txt的程序):
proc cluster data=eg method=ave standard outtree=beg
pseudo ccc;
var x1-x8;
id city;
run;
proc tree data=beg horizontal graphics nclusters=3 out=ceg;
run;
proc print data=ceg;
run;
结果如下:
聚类树状图
