多元统计实验三方差分析.doc
文本预览下载声明
桂林电子科技大学
数学与计算科学学院实验报告
院系 数学与计算科学学院 学号 姓名 成绩 课程
名称 应用多元统计分析实验 实验项目
名 称 实验三 方差分析 一 ，实验目的
用SAS软件进行单个总体均值的推断 (16种统计量，P51,默认的只有5种)
也可以检验正态总体均值是否为零。
格式：
proc means probt; /* probt 用于检验正态总体均值是否为零*/
var x ;
run;
2. Univariate过程
功能：UNIVARIATE过程除了可以完成与MEANS过程相同的基本统计量外，还可以计算变量的极端值、分位数，生成频率表，并支持对数据进行正态性检验。
格式：
proc univariate normal ; /* normal 用于检验总体是否服从正态*/
var x ;
run;
3. IML 过程
由于在DATA数据步的过程中， SAS系统并没有提供向量或矩阵的输入，要想建立一个矩阵的数据，就要用到SAS/IML模块，它有着一套很强大、灵活的运算语言。在IML的环境下不但可以进行矩阵的各种运算，而且还可以与DATA步一起使用实现各种复杂的计算和操作功能
三，实验内容
中世纪英国的建筑风格，早期是罗马式，后来有哥特式。罗马式教堂平均长度为145.29m，中殿的平均高度为22.69m.n=16个哥特式教堂的样本数据如下：
序号
长度/m
中殿
高度/m
序号
长度/m
中殿
高度/m
1
158.19
30.49
9
124.05
21.95
2
68.58
22.86
10
89.92
26.82
3
91.44
15.85
11
83.21
16.76
4
127.41
18.90
12
126.49
20.42
5
124.66
20.73
13
55.47
13.72
6
129.54
26.21
14
161.54
31.39
7
112.78
17.37
15
186.23
31.39
8
154.28
24.99
16
144.17
25.6
试检验, .
四，实验过程原始记录（数据，图表，计算等）
data shiyan3;
input no x1 x2 @@;
long=x1-145.29;
haigh=x2-22.69;
cards;
1 158.19 30.49 9 124.05 21.95
2 68.58 22.86 10 89.92 26.82
3 91.44 15.85 11 83.21 16.76
4 127.41 18.90 12 126.49 20.42
5 124.66 20.73 13 55.47 13.72
6 129.54 26.21 14 161.54 31.39
7 112.78 17.37 15 186.23 31.39
8 154.28 24.99 16 144.17 25.6
;
proc means data=shiyan3 probt; /* probt 用于检验正态总体均值是否为零*/
var long haigh ;
run;
proc univariate data=shiyan3 normal ;
var long haigh;
run;
方法1.
从上表可以看出长度的检验结果p值为0.017，小于0.05.因此，拒绝原假设，即认为罗马式
教堂平均长度不是145.29m。
高度的检验结果的p值为0.3157，大于0.05，因此，接受原假设，即认为中殿的平均高度
为22.69m。
方法2：
1.长度
图2.1.1 长度正态性检验
图2.1.2 长度位置检验
由于样本量2000,因此，采用Shapiro-Wilk检验。由图2.1可知，p=0.88240.05,接受原假设，
即罗马式教堂平均长度服从正态分布。
并由Students t检验的p值为0.0170.05，所以y的均值不为0,即即罗马式教堂平均长度不是
145.29m。
2. 高度
图2.2.1 正态性检验
图2.2.2 位置检验
由于样本量2000,因此，采用Shapiro-Wilk检验。由图2.2.1可知，p=0.58790.05,接受原假设，
即罗马式教堂中殿的平均高度服从正态分布。
并由Students t检验的p值为0.91570.05，所以y的均值为0,即即罗马式教堂中殿的平均高度
为22.69m。
方法3：
data shiyan3;
input no x1 x2 @@;
cards;
1
显示全部