如何在 sas 中执行 levene 检验
许多统计检验(例如单向方差分析)假设多个组之间的方差相等。
正式检验这一假设的一种方法是使用Levene 检验,该检验检验两个或多个组之间的方差是否相等。
该测试使用以下假设:
- 原假设 (H 0 ) :组间方差相等。
- 备择假设 ( HA ) :组间方差不相等。
如果检验的p 值小于所选的显着性水平,我们可以拒绝原假设并得出结论:我们有足够的证据表明组之间的方差不相等。
以下示例显示如何在 SAS 中执行 Levene 测试。
示例:SAS 中的 Levene 测试
假设我们在 SAS 中有以下数据集,显示了各种植物使用的肥料以及 18 种植物的生长情况(以英寸为单位):
/*create dataset*/ data my_data; input fertilizer $growth; datalines ; At 29 At 23 At 20 At 21 At 33 At 30 B 19 B 19 B17 B24 B25 B29 C 21 C22 C 30 C25 C24 C 33 ; run ; /*view dataset*/ proc print data =my_data;
假设我们想要执行单向方差分析来确定三种不同肥料之间的平均植物生长是否相等。
我们可以在 SAS 中使用以下语法,通过hovtest=levene(type=abs)语句执行单向方差分析,同时执行 Levene 测试以确定三个组是否具有相等的方差:
/*perform one-way ANOVA along with Levene's test*/
proc glm data = my_data;
class fertilizer;
model growth = fertilizer;
means fertilizer / hovtest =levene(type=abs);
run ;
第一个结果表显示单向方差分析的结果:
方差分析表中的 p 值为0.3358 。
由于该值不小于 0.05,因此我们可以得出结论,三种肥料之间的平均植物生长没有统计学上的显着差异。
为了检查这个结果是否可靠,我们必须检查方差齐性假设是否得到验证。
我们可以在输出的第二个表中看到 Levene 测试的结果:
从该表中,我们可以看到 Levene 检验的 p 值为0.6745 。
由于该值不小于 0.05,因此我们无法拒绝检验的原假设。
换句话说,我们可以假设三组的方差相等。
注意:我们在levene()函数中使用type=abs参数来指定在运行 Levene 测试时应使用残差的绝对值。这与R等其他统计软件使用的方法是一致的。
其他资源
以下教程解释了如何在 SAS 中执行其他常见统计测试:
如何在 SAS 中执行 Shapiro-Wilk 检验
如何在 SAS 中执行 Kolmogorov-Smirnov 检验
如何在 SAS 中使用 Proc Univariate 进行正态性检验