Welch检验:修订间差异

来自决策链云智库
(创建页面,内容为“{{Infobox nodebasic|nodename=Welch检验|nodeimage=Analysis of Variance_Multiple Comparisons.png|developer=Dev.Team-DPS|productionstate=PC可用|productionstatedesc=在 V1.0部署|nodeenglishname=Has english name::Analysis of Variance_Multiple Comparisons|abbreviation=AVMC|funcmaincategory=数据分析|funcsubcategory=DataAGM Lv1 Cat::方差分析|nodecategory=数据挖掘|nodeinterpretor=R|nodeshortdescription=<p>Welch检验是指采…”)
 
无编辑摘要
标签手工回退
 
(未显示2个用户的11个中间版本)
第1行: 第1行:
{{Infobox nodebasic|nodename=Welch检验|nodeimage=Analysis of Variance_Multiple Comparisons.png|developer=Dev.Team-DPS|productionstate=PC可用|productionstatedesc=在[[DecisionTree | V1.0]]部署|nodeenglishname=[[Has english name::Analysis of Variance_Multiple Comparisons]]|abbreviation=AVMC|funcmaincategory=数据分析|funcsubcategory=[[DataAGM Lv1 Cat::方差分析]]|nodecategory=数据挖掘|nodeinterpretor=R|nodeshortdescription=<p>Welch检验是指采用Welch分布的统计量检验各组均值是否相等。Welch分布近似于F分布,采用Welch检验对方差齐性没有要求,所以当数据的分布不满足方差齐性的要求时,采用Welch检验比F检验更稳妥。\n用途:假设各组方差不等的情况下,用于检验两个或多个总体均值是否存在显著差异。\n参数:选择连续型数值变量和分类分组变量</p>|nodeinputnumber=4|nodeoutputnumber=3|nodeloopsupport=是|nodeifswitchsupport=否|nodeavailableplotlist=SplittingNephelogram|nodeavailabletablelist=Table_For_Downstream|nodeconfiguration=VariableList;DropManu|nodeinputports=WorkFlow-Control ;Transfer-Table ■|nodeoutputports=WorkFlow-Control ;Transfer-Table ■|statsapewikiurl=https://wiki.statsape.com/Welch检验_Plus|previousnode=[[方差齐性检验]]|nextnode=[[球形检验]]}}{{Navplate AlgorithmNodeList}}[[Category:方差分析]]
{{Infobox nodebasic  
|nodename=Welch检验
|nodeimage=Welch Test.png
|icon=Welch Test.svg
|simpleicon=Welch Test_Pure.svg
|developer=Dev.Team-DPS
|productionstate={{图标文件|Win}} / {{图标文件|W10}} Win10及以上可用
|productionstatedesc=在[[Update:DecisionLinnc 1.0.0.8|V1.0]]部署
|nodeenglishname=Welch Test
|abbreviation=WlcTest
|funcmaincategory=数据分析
|funcsubcategory=[[DataAGM Lv1 Cat::方差分析]]
|nodecategory=数据挖掘
|nodeinterpretor=R
|nodeshortdescription=<p>Welch检验是指采用Welch分布的统计量检验各组均值是否相等。Welch分布近似于F分布,采用Welch检验对方差齐性没有要求,所以当数据的分布不满足方差齐性的要求时,采用Welch检验比F检验更稳妥。</p><p>用途:假设各组方差不等的情况下,用于检验两个或多个总体均值是否存在显著差异。</p><p>参数:选择连续型数值变量和分类分组变量</p>
|nodeinputnumber=4
|nodeoutputnumber=3
|nodeloopsupport=是
|nodeifswitchsupport=否
|nodeavailableplotlist=SplittingNephelogram
|nodeavailabletablelist=F-Value;df;P-Value
|nodeconfiguration=VariableList;DropMenu
|nodeinputports=WorkFlow-Control ➤;Transfer-Variable ◆;Transfer-Table ■
|nodeoutputports=WorkFlow-Control ➤;Transfer-Variable ◆;Transfer-Table ■
|statsapewikiurl=https://wiki.statsape.com/Welch检验
|previousnode=[[方差齐性检验]]
|nextnode=[[球形检验]]
}}
 
在[[statistics|统计学]]中,'''韦尔奇的''t''检验''',或称'''不等方差''t''检验''',是一种两样本[[location test|位置检验]],用于检验两个[[population (statistics)|总体]]均值是否相等的(零)假设。它以其创造者[[Bernard Lewis Welch|伯纳德·刘易斯·韦尔奇]]的名字命名,是[[Student's t-test|学生''t''检验]]的一个适应性改进<ref name=Welch1947>{{Cite journal | last = Welch | first = B. L. | title = 当涉及到几个不同的总体方差时,“学生”的问题的推广 | journal = [[Biometrika]] | volume = 34 |issue=1–2 | pages = 28–35 | year = 1947 |doi =10.1093/biomet/34.1-2.28 | pmid = 20287819 | mr = 19277 }}</ref>,当两个样本具有不等方差甚至可能不等样本量时更为可靠。<ref name=Ruxton2006>{{Cite journal | last = Ruxton | first = G. D. | title = 不等方差t检验是学生t检验和Mann–Whitney U检验的一种未被充分利用的替代方案 |journal = [[Behavioral Ecology (journal)|行为生态学]] | volume = 17 | issue = 4 | pages = 688–690 | year = 2006 | doi = 10.1093/beheco/ark016| doi-access = free }}</ref><ref name="WhyWelch">{{cite journal|last1=Derrick|first1=B|last2=Toher|first2=D|last3=White|first3=P|title=为什么韦尔奇检验在I型错误上是健壮的|journal=心理学的定量方法|date=2016|volume=12|issue=1|pages=30–38|doi=10.20982/tqmp.12.1.p030|url=http://eprints.uwe.ac.uk/27232/27/p030.pdf|doi-access=free}}</ref> 这些检验通常被称为“非配对的”或“独立样本的”''t''检验,因为它们通常应用于两个被比较样本的统计单位不重叠的情况。鉴于韦尔奇的''t''检验不如学生的''t''检验流行<ref name=Ruxton2006/>,且可能对读者来说不太熟悉,一个更具信息性的名称是“韦尔奇的不等方差''t''检验”——或简称“不等方差''t''检验”。<ref name=WhyWelch/>
 
=='''假设'''==
 
学生的''t''检验假设两个总体的样本均值呈正态分布,并且这些总体具有相等的方差。韦尔奇的''t''检验则是为不等总体方差设计的,但仍维持正态分布的假设。<ref name=Welch1947/> 韦尔奇的''t''检验是[[Behrens–Fisher problem|贝伦斯-费舍尔问题]]的一个近似解决方案。
 
=='''计算'''==
 
韦尔奇的''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]个[[mean|样本均值]]及其[[standard error|标准误差]],[math]s_i[/math]表示[[Standard deviation#Corrected sample standard deviation|校正的样本标准差]],[[sample size|样本量]]为[math]N_i[/math]。与[[Student's t test|学生的''t''检验]]不同,分母并非基于[[pooled variance|合并方差]]估计。
 
与这种方差估计相关的[[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]\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''个方差估计相关的自由度。
 
由于我们有[[chi-square distribution|卡方分布]]的近似,所以该统计量大致符合[[Student's t-distribution|''t''分布]]。当[math]N_1[/math]和[math]N_2[/math]都大于5时,这种近似更为准确。<ref>{{cite web |url=https://secure-media.collegeboard.org/apc/ap05_stats_allwood_fin4prod.pdf#page=7 |title=双样本t检验中自由度的萨特思韦特公式 |page=6 |first=Michael |last=Allwood |date=2008 }}</ref><ref>{{cite book |last1=Yates |last2=Moore |last3=Starnes |title=统计实践 |edition=第3版 |page=792 |year=2008 |publisher=W.H. Freeman and Company |location=纽约 |isbn=9780716773092 }}</ref>
 
=='''统计检验'''==
 
一旦计算出''t''和''[math]\nu[/math]'',这些统计量可以与[[Student's t-distribution|''t''分布]]一起使用来检验两种可能的[[null hypothesis|零假设]]之一:
* 两个总体均值相等,此时应用[[two-tailed test|双尾检验]];或
* 其中一个总体均值大于或等于另一个,此时应用[[one-tailed test|单尾检验]]。
这些近似自由度是[[real number|实数]][math]\left(\nu\in\mathbb{R}^+\right)[/math],在面向统计的软件中就是这样使用的,而在电子表格中则向下舍入到最接近的整数。
 
=='''优势与局限性'''==
 
韦尔奇的''t''检验比学生''t''检验更为稳健,能够在方差不等和样本量不等的情况下,保持[[Type I and type II errors|第一类和第二类错误率]]接近名义水平。此外,即使在总体方差相等且样本量均衡的情况下,韦尔奇的''t''检验的[[Power (statistics)|功效]]也接近学生''t''检验的功效。<ref name=Ruxton2006/> 韦尔奇的''t''检验可以推广到2个以上的样本,<ref name="Welch1951">{{cite journal|last1=Welch|first1=B. L.|title=On the Comparison of Several Mean Values: An Alternative Approach|journal=Biometrika|date=1951|volume=38|issue=3/4|pages=330–336|doi=10.2307/2332579|jstor=2332579}}</ref> 这比[[one-way analysis of variance]] (ANOVA)更为稳健。
 
''不推荐''先进行方差相等的预检验,然后在学生''t''检验和韦尔奇''t''检验之间选择。<ref name=Zimmerman2004>{{Cite journal | last = Zimmerman | first = D. W. | title = A note on preliminary tests of equality of variances | journal = [[British Journal of Mathematical and Statistical Psychology]] | volume = 57 | pages = 173–181 | year = 2004 | issue = Pt 1 | doi = 10.1348/000711004849222| pmid = 15171807 }}</ref> 相反,如上所述,可以直接应用韦尔奇的''t''检验,而没有对学生''t''检验的任何实质性不利影响。韦尔奇的''t''检验对于偏态分布和大样本量仍然稳健。<ref name=Fagerland2012>{{Cite journal | last = Fagerland | first = M. W. | title = t-tests, non-parametric tests, and large studies—a paradox of statistical practice? | journal = BMC Medical Research Methodology| volume = 12 | page = 78 | year = 2012 | doi = 10.1186/1471-2288-12-78| pmc = 3445820 | pmid=22697476 | doi-access = free }}</ref> 对于偏态分布和较小的样本量,其可靠性会降低,在这种情况下可以考虑执行韦尔奇的''t''检验。<ref name=Fagerland2009>{{Cite journal | last1 = Fagerland | first1 = M. W. | last2 = Sandvik | first2 = L. | title = Performance of five two-sample location tests for skewed distributions with unequal variances | journal = [[Contemporary Clinical Trials]] | volume = 30 | issue = 5 | pages = 490–496 | year = 2009 | doi=10.1016/j.cct.2009.06.007| pmid = 19577012 }}</ref>
 
=='''示例'''==
 
以下三个示例比较了韦尔奇的''t''检验和学生''t''检验。样本来自使用[[R (programming language)|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值:
 
{| cellpadding="5" style="margin:auto;"  class="wikitable"
|-
!|
! colspan="3" style="text-align:center;"| Sample A1
! colspan="3" style="text-align:center;"| Sample A2
! colspan="4" style="text-align:center;"| Student's ''t''-test
! colspan="4" style="text-align:center;"| Welch's ''t''-test
|-
! 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_2[/math] || style="text-align:center;"| [math]\overline{X}_2[/math] || style="text-align:center;"| [math]s_2^2[/math]
! style="text-align:center;"| {{tmath|t}} || style="text-align:center;"| {{tmath|\nu}} || style="text-align:center;"| {{tmath|P}} || style="text-align:center;"| [math]P_\mathrm{sim}[/math]
! style="text-align:center;"| {{tmath|t}} || style="text-align:center;"| {{tmath|\nu}} || style="text-align:center;"| {{tmath|P}} || 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检验 ===
<syntaxhighlight lang="R">
 
oneway.test(formula, data, subset, na.action, var.equal = FALSE)
 
</syntaxhighlight>
 
== '''节点使用指南''' ==
* 用于确定两个独立样本是否有显著差异的统计检验
* 适用于当两个群体的方差不相等或样本大小不同的情况
* Welch检验不要求两个总体具有相同的方差
 
=== 方法选择 ===
* 无方法选择
 
=== 参数配置 ===
* 统计变量:选择连续型数值变量,每个变量会做一次Welch检验
* 分组变量:选择一个分类分组变量
* 筛选阈值:选择需要的P值阈值,节点会自动将满足阈值的变量筛选出,数据集也会同步筛选出满足的变量。
* 统计变量和分组变量要规避复用
* 此算法兼容空值
 
=== 注意事项 ===
* 确保两个样本是独立的,即一个样本中的观测并不影响另一个样本中的观测
* 对于极小的样本量,检验的效力会下降,这意味着检验发现实际差异的能力降低
 
== '''引用''' ==
{{Reflist}}
 
{{Navplate AlgorithmNodeList}}
 
[[Category:方差分析]]

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-支持逻辑判断
输入输出
可生成图片类型(推荐)
可生成数据表类型(推荐)
相关节点
上一节点方差齐性检验
下一节点球形检验



统计学中,韦尔奇的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.

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