兩組均數(shù)比較,可以用方差分析嗎?那些結(jié)果一樣的統(tǒng)計(jì)方法

2023年鄭老師多門科研統(tǒng)計(jì)課程:多次直播,含孟德爾隨機(jī)化方法

前段時(shí)間,有個(gè)學(xué)生就就論文修回中專家意見咨詢,專家對下面的表格兩組均數(shù)比較應(yīng)該采用t檢驗(yàn)而不是方差分析(F檢驗(yàn))。專家認(rèn)為,男性和女性的生命質(zhì)量比較應(yīng)該用t檢驗(yàn),而不是F檢驗(yàn)。真應(yīng)該這樣子的嗎?

氣的肝疼,我這里只能說,專家又扯蛋了。步子大了,真的很容易扯到蛋。
一般情況下,如果在同一張表格內(nèi)同時(shí)有兩組均數(shù)比較和多組均數(shù)比較,我的推薦無論兩組或者多組均數(shù)比較應(yīng)該都采用F檢驗(yàn)。因?yàn)閮山M均數(shù)比較時(shí),方差分析和t分析結(jié)果是一模一樣的。
1、t檢驗(yàn)和F檢驗(yàn)分析的結(jié)果
案例:
將出生28天的20只大白鼠隨機(jī)分成兩組,分別飼以高蛋白和低蛋白飼料,8周后觀察其體重(g)。問兩種不同飼料對大白鼠的體重影響有無差別?
高蛋白組:133,145,112,138,99,157,126,121,139,106,115
低蛋白組:118,75,106,87,94,110,102,124,130
t檢驗(yàn)的結(jié)果

F檢驗(yàn)也就是方差分析的結(jié)果

可以看出,兩個(gè)結(jié)果P=0.016,沒有任何差別。
2、兩樣本和多樣本秩和檢驗(yàn)分析的結(jié)果
兩樣本秩和檢驗(yàn)為Wilcoxon Mann-Whitney秩和檢驗(yàn),多樣本秩和檢驗(yàn)為Kruskal-Wallis H檢驗(yàn),一般情況下,2組數(shù)據(jù)秩和檢驗(yàn)都是為Wilcoxon Mann-Whitney秩和檢驗(yàn),但其實(shí),多樣本秩和檢驗(yàn)也是可以用的,結(jié)果也是一樣的。
我們就上述案例分別開展兩樣本和多樣本秩和檢驗(yàn)。
WilcoxonMann-Whitney檢驗(yàn)結(jié)果

Kruskal-Wallis H檢驗(yàn)結(jié)果

可以看出,兩個(gè)結(jié)果P=0.023,也沒有任何差別。
3、檢驗(yàn)和單因素線性回歸分析結(jié)果是完全一致的
我們在開展回歸分析的時(shí)候,經(jīng)常要先做單因素,再做多因素回歸。那么很多人舉棋不定,因?yàn)樗奈恼驴赡芤呀?jīng)做了,t檢驗(yàn)或者F檢驗(yàn),還要做單因素線性回歸嗎?答案是,不需要,因?yàn)榻Y(jié)果是一致的。
t檢驗(yàn)的結(jié)果

單因素線性回歸的結(jié)果

也可以看出,兩個(gè)結(jié)果P=0.016,沒有任何差別。
4、卡方檢驗(yàn)和單因素logistic回歸分析結(jié)果是一致的
卡方檢驗(yàn)的卡方值

卡方檢驗(yàn)的OR值

單因素logistic回歸分析的P值和OR值

雖然最后logistic的最終結(jié)果P值稍微有點(diǎn)差異,但主要指標(biāo)以及OR值完全一致,因此卡方檢驗(yàn)和單因素logistic回歸結(jié)果也是相同的。
5、統(tǒng)計(jì)小結(jié)
--t檢驗(yàn)和方差分析對兩樣本均數(shù)比較,t檢驗(yàn)和單因素線性回歸、兩樣本秩和和多樣本秩對兩組數(shù)據(jù)總體分布比較、卡方檢驗(yàn)和單因素logistic回歸分析,結(jié)果相同。
--如果統(tǒng)計(jì)表既要做t檢驗(yàn)也要做F檢驗(yàn),可以統(tǒng)一用F檢驗(yàn);秩和檢驗(yàn)也是如此,這樣可以統(tǒng)一表格風(fēng)格。
--回歸分析時(shí),當(dāng)自變量為二分類變量時(shí),單因素線性回歸和logistic回歸有時(shí)候可以用t檢驗(yàn)和卡方檢驗(yàn)代替。直接基于t檢驗(yàn)和卡方檢驗(yàn)做多因素線性回歸分析。

詳情請點(diǎn)擊下方:
https://mp.weixin.qq.com/s?__biz=MzAwOTYyMDY3OQ==&mid=2650407845&idx=1&sn=6121415ba521f6ae5ed8ef917ad9fd0e&chksm=83527e0db425f71b7ed6706dd2ae9f74f7bc64c93ce298321eba49cccbff09e7ab0d183c2808&token=549692846&lang=zh_CN#rd
