如何在Python中执行F检验

2023年 9月 9日 46.7k 0

统计学家使用 F 检验来检查两个数据集是否具有相同的方差。 F 检验以罗纳德·费舍尔爵士的名字命名。为了使用 F 检验,我们做出两个假设,一个原假设和一个备择假设。然后我们选择 F 检验认可的这两个假设中的任何一个。

方差是一种数据分布度量,用于说明数据与平均值的偏差。较高的值比较较小的值显示出更大的离散度。

在本文中,您将学习如何在Python编程语言中执行F-Test以及其使用案例。

F-测试过程

执行F-检验的过程如下:

  • 首先,定义原假设和备择假设。

    • 零假设或H0:σ12 = σ22(总体方差相等)

    • 替代假设或 H1:σ12 ≠ σ22(总体方差不相等)

  • 选择用于测试的统计数据。

  • 计算总体的自由度。例如,如果 m 和 n 是总体形状,则自由度分别表示为 (df1) = m–1 和 (df2) = n – 1。

  • 现在从 F 表中查找 F 值。

  • 最后,将双尾检验的 alpha 值除以 2,计算出临界值。

因此,我们使用总体的自由度定义 F 值。我们读取第一行中的 df1,而第一列中的 df2。

有各种用于独特自由度的 F 表。我们将步骤 2 中的 F 统计量与步骤 4 中计算的临界值进行比较。如果临界值小于 F 统计量,我们可以拒绝原假设。相反,当临界值在某个显着水平上大于F统计量时,我们可以接受原假设。

如何在Python中执行F检验

假设

在进行基于数据集的F-检验之前,我们做出了一些假设。

  • 数据总体服从正态分布,即符合钟形曲线。

  • 样本之间不相关,即人群中不存在多重共线性。

除了这些假设之外,在进行F检验时,我们还应考虑以下关键要点:

  • 最大方差值应该在分子中以执行右尾检验。

  • 在双尾检验中,将alpha除以2后确定临界值。

  • 检查是否存在方差或标准差。

  • 如果 F 表中没有自由度,则以最大值作为临界值。

F-Test在Python中的应用

语法

scipy stats.f()

登录后复制

参数

x : quantiles
q : lower or upper tail probability
dfn, dfd shape parameters
loc :location parameter
scale : scale parameter (default=1)
size : random variate shape
moments : [‘mvsk’] letters, specifying which moments to compute

登录后复制

Explanation

的中文翻译为:

解释

在这种方法中,用户必须将f_value和每个数组的可迭代长度传递给scipy.stats.f.cdf(),并将其减去1以执行F检验。

算法

  • 首先,导入NumPy和Scipy.stats库进行操作。

  • 然后创建两个具有两个不同变量名称的随机选择值列表,并将它们转换为 NumPy 数组,并使用 Numpy 计算每个数组的方差。

  • 定义一个函数来计算F分数,其中首先我们将数组的方差除以自由度为1。

  • 然后计算每个数组的可迭代长度,并将 f 值(方差比率)和长度传递到 CDF 函数中,并从 1 中减去该长度以计算 p 值。

  • 最后,函数返回 p_value 和 f_value。

示例

import numpy as np
import scipy.stats

# Create data
group1 = [0.28, 0.2, 0.26, 0.28, 0.5]
group2 = [0.2, 0.23, 0.26, 0.21, 0.23]

# Converting the list to an array
x = np.array(group1)
y = np.array(group2)

# Calculate the variance of each group
print(np.var(group1), np.var(group2))

def f_test(group1, group2):
f = np.var(group1, ddof=1)/np.var(group2, ddof=1)
nun = x.size-1
dun = y.size-1
p_value = 1-scipy.stats.f.cdf(f, nun, dun)
return f, p_value

# perform F-test
f_test(x, y)

登录后复制

输出

Variances: 0.010464 0.00042400000000000017

登录后复制

您可以观察到 F 检验值为 4.38712, 相应的 p 值为 0.019127。

由于p值小于0.05,我们将放弃零假设。因此,我们可以说这两个总体的方差不相等。

结论

读完本文后,您现在知道如何使用 F 检验来检查两个样本是否属于具有相同方差的总体。您已经了解了 F 测试过程、假设和 Python 实现。让我们用一些要点来总结这篇文章 -

  • F检验告诉你两个总体是否具有相等的方差。

  • 计算自由度并计算临界值。

  • 从F-表中找到F统计量,并将其与在前一步计算的关键值进行比较。

  • 根据临界值和 F 统计量比较接受或拒绝原假设。

以上就是如何在Python中执行F检验的详细内容,更多请关注每日运维网(www.mryunwei.com)其它相关文章!

相关文章

JavaScript2024新功能:Object.groupBy、正则表达式v标志
PHP trim 函数对多字节字符的使用和限制
新函数 json_validate() 、randomizer 类扩展…20 个PHP 8.3 新特性全面解析
使用HTMX为WordPress增效:如何在不使用复杂框架的情况下增强平台功能
为React 19做准备:WordPress 6.6用户指南
如何删除WordPress中的所有评论

发布评论