Welch检验：修订间差异

来自决策链云智库
无编辑摘要
无编辑摘要
标签：手工回退
 
（未显示同一用户的3个中间版本）
第46行： 第46行：
与这种方差估计相关的[[degrees of freedom (statistics)|自由度]][math]\nu[/math]&nbsp;是使用[[Welch–Satterthwaite equation|韦尔奇-萨特思韦特方程]]近似计算的：<ref>[ 7.3.1. Do two processes have the same mean?], 工程统计手册, [[NIST]].  (2021-07-30访问的在线资源。)</ref>
与这种方差估计相关的[[degrees of freedom (statistics)|自由度]][math]\nu[/math]&nbsp;是使用[[Welch–Satterthwaite equation|韦尔奇-萨特思韦特方程]]近似计算的：<ref>[ 7.3.1. Do two processes have the same mean?], 工程统计手册, [[NIST]].  (2021-07-30访问的在线资源。)</ref>


:[math]
[math]\nu \approx \frac{\left(\frac{s_{1}^{2}}{N_{1}}+\frac{s_{2}^{2}}{N_{2}}\right)^{2}}{\frac{s_{1}^{4}}{N_{1}^{2} \nu_{1}}+\frac{s_{2}^{4}}{N_{2}^{2} \nu_{2}}} .[/math]
\nu \quad  \approx \quad
\frac{\left( \; \frac{s_1^2}{N_1} \; + \; \frac{s_2^2}{N_2} \; \right)^2 }
{ \quad \frac{s_1^4}{N_1^2 \nu_1} \; + \; \frac{s_2^4}{N_2^2 \nu_2 } \quad }.
[/math]


当[math]N_1 = N_2[/math]时，这个表达式可以简化为：
当[math]N_1 = N_2[/math]时，这个表达式可以简化为：
:[math]
 
\nu \approx \frac
[math]\nu \approx \frac{s_{\Delta \bar{X}}^{4}}{\nu_{1}^{-1} s_{\bar{X}_{1}}^{4}+\nu_{2}^{-1} s_{\bar{X}_{2}}^{4}}[/math]
{s_{\Delta\bar{X}}^4}
{\nu_1^{-1} s_{\bar{X}_1}^4 + \nu_2^{-1} s_{\bar{X}_2}^4}.
[/math]


此处，[math]\nu_i = N_i-1[/math]是与第''i''个方差估计相关的自由度。
此处，[math]\nu_i = N_i-1[/math]是与第''i''个方差估计相关的自由度。
第89行： 第82行：
第二个示例是对于不等的方差([math]\sigma_1^2 = 9.0[/math]，[math]\sigma_2^2 = 0.9[/math])和不等的样本量([math]N_1 = 10[/math]，[math]N_2 = 20[/math])。较小的样本具有更大的方差：
第二个示例是对于不等的方差([math]\sigma_1^2 = 9.0[/math]，[math]\sigma_2^2 = 0.9[/math])和不等的样本量([math]N_1 = 10[/math]，[math]N_2 = 20[/math])。较小的样本具有更大的方差：


: [math]\begin{align}
[math]\begin{align}
A_1 &= \{17.2, 20.9, 22.6, 18.1, 21.7, 21.4, 23.5, 24.2, 14.7, 21.8\}
A_1 &= \{17.2, 20.9, 22.6, 18.1, 21.7, 21.4, 23.5, 24.2, 14.7, 21.8\}
\\
\\
第97行： 第90行：
第三个示例是对于不等的方差([math]\sigma_1^2 = 1.4[/math]，[math]\sigma_2^2 = 17.1[/math])和不等的样本量([math]N_1 = 10[/math]，[math]N_2 = 20[/math])。较大的样本具有更大的方差：
第三个示例是对于不等的方差([math]\sigma_1^2 = 1.4[/math]，[math]\sigma_2^2 = 17.1[/math])和不等的样本量([math]N_1 = 10[/math]，[math]N_2 = 20[/math])。较大的样本具有更大的方差：


: [math]\begin{align}
[math]\begin{align}
A_1 &= \{19.8, 20.4, 19.6, 17.8, 18.5, 18.9, 18.3, 18.9, 19.5, 22.0\}
A_1 &= \{19.8, 20.4, 19.6, 17.8, 18.5, 18.9, 18.3, 18.9, 19.5, 22.0\}
\\
\\
第108行： 第101行：
|-
|-
!|
!|
! colspan="3" style="text-align:center;"| 样本A1
! colspan="3" style="text-align:center;"| Sample A1
! colspan="3" style="text-align:center;"| 样本A2
! colspan="3" style="text-align:center;"| Sample A2
! colspan="4" style="text-align:center;"| 学生''t''检验
! colspan="4" style="text-align:center;"| Student's ''t''-test
! colspan="4" style="text-align:center;"| 韦尔奇''t''检验
! colspan="4" style="text-align:center;"| Welch's ''t''-test
|-
|-
! style="text-align:center;"| 示例
! style="text-align:center;"| Example
! style="text-align:center;"| [math]N_1[/math] || style="text-align:center;"| [math]\overline{X}_1[/math] || style="text-align:center;"| [math]s_1^2[/math]
! style="text-align:center;"| [math]N_1[/math] || style="text-align:center;"| [math]\overline{X}_1[/math] || style="text-align:center;"| [math]s_1^2[/math]
! style="text-align:center;"| [math]N_2[/math] || style="text-align:center;"| [math]\overline{X}_2[/math] || style="text-align:center;"| [math]s_2^2[/math]
! style="text-align:center;"| [math]N_2[/math] || style="text-align:center;"| [math]\overline{X}_2[/math] || style="text-align:center;"| [math]s_2^2[/math]

2024年2月9日 (五) 09:52的最新版本

Welch Test.png
节点状态
Windows / Windows 10 Win10及以上可用
在V1.0部署
Welch检验Welch Test.svg
节点开发者决策链算法研发部 (Dev.Team-DPS)
节点英文名Welch Test
功能主类别数据分析
英文缩写WlcTest
功能亚类别方差分析
节点类型数据挖掘
开发语言R
节点简介

Welch检验是指采用Welch分布的统计量检验各组均值是否相等。Welch分布近似于F分布，采用Welch检验对方差齐性没有要求，所以当数据的分布不满足方差齐性的要求时，采用Welch检验比F检验更稳妥。

用途：假设各组方差不等的情况下，用于检验两个或多个总体均值是否存在显著差异。

参数：选择连续型数值变量和分类分组变量

端口数量与逻辑控制(PC)
Input-入口4个
Output-出口3个
Loop-支持循环是
If/Switch-支持逻辑判断否
输入输出
可生成图片类型（推荐）
可生成数据表类型（推荐）
相关节点
上一节点方差齐性检验
下一节点球形检验
WikimediaUI-Globe.svg相关网站




在统计学中，韦尔奇的t检验，或称不等方差t检验，是一种两样本位置检验，用于检验两个总体均值是否相等的（零）假设。它以其创造者伯纳德·刘易斯·韦尔奇的名字命名，是学生t检验的一个适应性改进[1]，当两个样本具有不等方差甚至可能不等样本量时更为可靠。[2][3] 这些检验通常被称为“非配对的”或“独立样本的”t检验，因为它们通常应用于两个被比较样本的统计单位不重叠的情况。鉴于韦尔奇的t检验不如学生的t检验流行[2]，且可能对读者来说不太熟悉，一个更具信息性的名称是“韦尔奇的不等方差t检验”——或简称“不等方差t检验”。[3]

假设

学生的t检验假设两个总体的样本均值呈正态分布，并且这些总体具有相等的方差。韦尔奇的t检验则是为不等总体方差设计的，但仍维持正态分布的假设。[1] 韦尔奇的t检验是贝伦斯-费舍尔问题的一个近似解决方案。

计算

韦尔奇的t检验通过以下公式定义统计量t：

[math]t = \frac{\Delta\overline{X}}{s_{\Delta\bar{X}}} = \frac{\overline{X}_1 - \overline{X}_2}{\sqrt{ {s_{\bar{X}_1}^2} + {s_{\bar{X}_2}^2} }}\,[/math]
[math]s_{\bar{X}_i} = {s_i \over \sqrt{N_i}} \,[/math]

其中[math]\overline{X}_i[/math]和[math]s_{\bar{X}_i}[/math]分别是第[math]i^\text{th}[/math]个样本均值及其标准误差，[math]s_i[/math]表示校正的样本标准差，样本量为[math]N_i[/math]。与学生的t检验不同，分母并非基于合并方差估计。

与这种方差估计相关的自由度[math]\nu[/math] 是使用韦尔奇-萨特思韦特方程近似计算的：[4]

[math]\nu \approx \frac{\left(\frac{s_{1}^{2}}{N_{1}}+\frac{s_{2}^{2}}{N_{2}}\right)^{2}}{\frac{s_{1}^{4}}{N_{1}^{2} \nu_{1}}+\frac{s_{2}^{4}}{N_{2}^{2} \nu_{2}}} .[/math]

当[math]N_1 = N_2[/math]时，这个表达式可以简化为：

[math]\nu \approx \frac{s_{\Delta \bar{X}}^{4}}{\nu_{1}^{-1} s_{\bar{X}_{1}}^{4}+\nu_{2}^{-1} s_{\bar{X}_{2}}^{4}}[/math]

此处，[math]\nu_i = N_i-1[/math]是与第i个方差估计相关的自由度。

由于我们有卡方分布的近似，所以该统计量大致符合t分布。当[math]N_1[/math]和[math]N_2[/math]都大于5时，这种近似更为准确。[5][6]

统计检验

一旦计算出t和[math]\nu[/math]，这些统计量可以与t分布一起使用来检验两种可能的零假设之一：

  • 两个总体均值相等，此时应用双尾检验；或
  • 其中一个总体均值大于或等于另一个，此时应用单尾检验。

这些近似自由度是实数[math]\left(\nu\in\mathbb{R}^+\right)[/math]，在面向统计的软件中就是这样使用的，而在电子表格中则向下舍入到最接近的整数。

优势与局限性

韦尔奇的t检验比学生t检验更为稳健，能够在方差不等和样本量不等的情况下，保持第一类和第二类错误率接近名义水平。此外，即使在总体方差相等且样本量均衡的情况下，韦尔奇的t检验的功效也接近学生t检验的功效。[2] 韦尔奇的t检验可以推广到2个以上的样本，[7] 这比one-way analysis of variance (ANOVA)更为稳健。

不推荐先进行方差相等的预检验，然后在学生t检验和韦尔奇t检验之间选择。[8] 相反，如上所述，可以直接应用韦尔奇的t检验，而没有对学生t检验的任何实质性不利影响。韦尔奇的t检验对于偏态分布和大样本量仍然稳健。[9] 对于偏态分布和较小的样本量，其可靠性会降低，在这种情况下可以考虑执行韦尔奇的t检验。[10]

示例

以下三个示例比较了韦尔奇的t检验和学生t检验。样本来自使用R编程语言的随机正态分布。

对于所有三个示例，总体均值为[math]\mu_1 = 20[/math]和[math]\mu_2 = 22[/math]。

第一个示例是对于不等但接近的方差([math]\sigma_1^2 = 7.9[/math]，[math]\sigma_2^2 = 3.8[/math])和相等的样本量([math]N_1 = N_2 = 15[/math])。让A1和A2表示两个随机样本：

[math]A_1 = \{27.5, 21.0, 19.0, 23.6, 17.0, 17.9, 16.9, 20.1, 21.9, 22.6, 23.1, 19.6, 19.0, 21.7, 21.4\}[/math]
[math]A_2 = \{27.1, 22.0, 20.8, 23.4, 23.4, 23.5, 25.8, 22.0, 24.8, 20.2, 21.9, 22.1, 22.9, 20.5, 24.4\}[/math]

第二个示例是对于不等的方差([math]\sigma_1^2 = 9.0[/math]，[math]\sigma_2^2 = 0.9[/math])和不等的样本量([math]N_1 = 10[/math]，[math]N_2 = 20[/math])。较小的样本具有更大的方差：

[math]\begin{align} A_1 &= \{17.2, 20.9, 22.6, 18.1, 21.7, 21.4, 23.5, 24.2, 14.7, 21.8\} \\ A_2 &= \{21.5, 22.8, 21.0, 23.0, 21.6, 23.6, 22.5, 20.7, 23.4, 21.8, 20.7, 21.7, 21.5, 22.5, 23.6, 21.5, 22.5, 23.5, 21.5, 21.8\} \end{align}[/math]

第三个示例是对于不等的方差([math]\sigma_1^2 = 1.4[/math]，[math]\sigma_2^2 = 17.1[/math])和不等的样本量([math]N_1 = 10[/math]，[math]N_2 = 20[/math])。较大的样本具有更大的方差：

[math]\begin{align} A_1 &= \{19.8, 20.4, 19.6, 17.8, 18.5, 18.9, 18.3, 18.9, 19.5, 22.0\} \\ A_2 &= \{28.2, 26.6, 20.1, 23.3, 25.2, 22.1, 17.7, 27.6, 20.6, 13.7, 23.2, 17.5, 20.6, 18.0, 23.9, 21.6, 24.3, 20.4, 24.0, 13.2\} \end{align}[/math]

通过模拟t统计量的分布来获得参考p值，用于检验总体均值相等的零假设([math]\mu_1 - \mu_2 =0[/math])。结果在下表中总结，显示了双尾p值：

Sample A1 Sample A2 Student's t-test Welch's t-test
Example [math]N_1[/math] [math]\overline{X}_1[/math] [math]s_1^2[/math] [math]N_2[/math] [math]\overline{X}_2[/math] [math]s_2^2[/math] [mathmath>[/math] style="text-align:center;"| [mathmath>[/math] style="text-align:center;"| [mathmath>[/math] style="text-align:center;"| [math]P_\mathrm{sim}[/math] [mathmath>[/math] style="text-align:center;"| [mathmath>[/math] style="text-align:center;"| [mathmath>[/math] style="text-align:center;"| [math]P_\mathrm{sim}[/math]
1 15 20.8 7.9 15 23.0 3.8 −2.46 28 0.021 0.021 −2.46 24.9 0.021 0.017
2 10 20.6 9.0 20 22.1 0.9 −2.10 28 0.045 0.150 −1.57 9.9 0.149 0.144
3 10 19.4 1.4 20 21.6 17.1 −1.64 28 0.110 0.036 −2.22 24.5 0.036 0.042

当两个样本具有相似的方差和样本量时（示例1），韦尔奇的t检验和学生t检验给出了相同的结果。但请注意，即使您从具有相同方差的总体中采样，样本方差也会不同，两个t检验的结果也会有所不同。因此，使用实际数据时，两个检验几乎总是会给出略有不同的结果。

对于不等方差，当较小样本具有较大方差时（示例2），学生t检验给出了低p值；当较大样本具有较大方差时（示例3），学生t检验给出了高p值。对于不等方差，韦尔奇t检验给出的p值接近模拟p值。

节点使用的R语言示例代码

Welch检验

oneway.test(formula, data, subset, na.action, var.equal = FALSE)

节点使用指南

  • 用于确定两个独立样本是否有显著差异的统计检验
  • 适用于当两个群体的方差不相等或样本大小不同的情况
  • Welch检验不要求两个总体具有相同的方差

方法选择

  • 无方法选择

参数配置

  • 统计变量：选择连续型数值变量，每个变量会做一次Welch检验
  • 分组变量：选择一个分类分组变量
  • 筛选阈值：选择需要的P值阈值，节点会自动将满足阈值的变量筛选出，数据集也会同步筛选出满足的变量。
  • 统计变量和分组变量要规避复用
  • 此算法兼容空值

注意事项

  • 确保两个样本是独立的，即一个样本中的观测并不影响另一个样本中的观测
  • 对于极小的样本量，检验的效力会下降，这意味着检验发现实际差异的能力降低

引用

  1. ↑ 1.0 1.1 Welch, B. L. (1947). "当涉及到几个不同的总体方差时，"学生"的问题的推广". Biometrika. 34 (1–2): 28–35. doi:10.1093/biomet/34.1-2.28. MR 0019277. PMID 20287819.
  2. ↑ 2.0 2.1 2.2 Ruxton, G. D. (2006). "不等方差t检验是学生t检验和Mann–Whitney U检验的一种未被充分利用的替代方案". 行为生态学. 17 (4): 688–690. doi:10.1093/beheco/ark016.
  3. ↑ 3.0 3.1 Derrick, B; Toher, D; White, P (2016). "为什么韦尔奇检验在I型错误上是健壮的" (PDF). 心理学的定量方法. 12 (1): 30–38. doi:10.20982/tqmp.12.1.p030.
  4. ↑ [ 7.3.1. Do two processes have the same mean?], 工程统计手册, NIST. (2021-07-30访问的在线资源。)
  5. ↑ Allwood, Michael (2008). "双样本t检验中自由度的萨特思韦特公式" (PDF). p. 6.
  6. ↑ Yates; Moore; Starnes (2008). 统计实践 (第3版 ed.). 纽约: W.H. Freeman and Company. p. 792. ISBN 9780716773092.
  7. ↑ Welch, B. L. (1951). "On the Comparison of Several Mean Values: An Alternative Approach". Biometrika. 38 (3/4): 330–336. doi:10.2307/2332579. JSTOR 2332579.
  8. ↑ Zimmerman, D. W. (2004). "A note on preliminary tests of equality of variances". British Journal of Mathematical and Statistical Psychology. 57 (Pt 1): 173–181. doi:10.1348/000711004849222. PMID 15171807.
  9. ↑ Fagerland, M. W. (2012). "t-tests, non-parametric tests, and large studies—a paradox of statistical practice?". BMC Medical Research Methodology. 12: 78. doi:10.1186/1471-2288-12-78. PMC 3445820. PMID 22697476.
  10. ↑ Fagerland, M. W.; Sandvik, L. (2009). "Performance of five two-sample location tests for skewed distributions with unequal variances". Contemporary Clinical Trials. 30 (5): 490–496. doi:10.1016/j.cct.2009.06.007. PMID 19577012.

查找其他类别的节点，请参考以下列表

取自“https://wiki.statsape.com/index.php?title=Welch检验&oldid=9189”
最近修改
此页面最后编辑于2024年2月9日 (星期五) 09:52。
版权
DecisionLinnc Wiki (决策链云智库) © 2024 by Statsape is licensed under Attribution-NonCommercial-ShareAlike 4.0 International. To view a copy of this license, visit [http://creativecommons.org/licenses/by-nc-sa/4.0/| CreativeCommons-Licenses-4.0].