Zeroclanzhang(讨论 | 贡献) (创建页面,内容为“{{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= | {{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] 是使用[[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的最新版本
节点状态 | / Win10及以上可用
在V1.0部署
|
---|---|
Welch检验 | |
节点开发者 | 决策链算法研发部 (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.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.0 2.1 2.2 Ruxton, G. D. (2006). "不等方差t检验是学生t检验和Mann–Whitney U检验的一种未被充分利用的替代方案". 行为生态学. 17 (4): 688–690. doi:10.1093/beheco/ark016.
- ↑ 3.0 3.1 Derrick, B; Toher, D; White, P (2016). "为什么韦尔奇检验在I型错误上是健壮的" (PDF). 心理学的定量方法. 12 (1): 30–38. doi:10.20982/tqmp.12.1.p030.
- ↑ [ 7.3.1. Do two processes have the same mean?], 工程统计手册, NIST. (2021-07-30访问的在线资源。)
- ↑ Allwood, Michael (2008). "双样本t检验中自由度的萨特思韦特公式" (PDF). p. 6.
- ↑ Yates; Moore; Starnes (2008). 统计实践 (第3版 ed.). 纽约: W.H. Freeman and Company. p. 792. ISBN 9780716773092.
- ↑ 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.
- ↑ 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.
- ↑ 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.
- ↑ 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.
查找其他类别的节点,请参考以下列表