方差分析是一种用于比较多个组或处理之间均值是否有显著差异的统计方法。在SciPy中,可以使用stats模块中的f_oneway函数来进行方差分析。
下面是一个简单的示例代码,演示如何使用SciPy进行方差分析:
from scipy import stats
# 定义三个样本数据
data1 = [1, 2, 3, 4, 5]
data2 = [2, 3, 4, 5, 6]
data3 = [3, 4, 5, 6, 7]
# 进行方差分析
f_statistic, p_value = stats.f_oneway(data1, data2, data3)
# 输出结果
print("F统计量:", f_statistic)
print("P值:", p_value)
# 判断是否拒绝零假设
alpha = 0.05
if p_value < alpha:
print("拒绝零假设,组间均值存在显著差异")
else:
print("接受零假设,组间均值不存在显著差异")
在上面的代码中,首先定义了三个样本数据data1、data2和data3。然后使用f_oneway函数进行方差分析,并得到F统计量和P值。最后根据P值与显著性水平alpha的比较,判断是否拒绝零假设。