15-假設(shè)檢驗之F檢驗

F檢驗(又稱為方差齊性檢驗)主要對于方差齊性或方差同質(zhì)性進(jìn)行檢驗殴蓬。


F分布

獨立樣本T檢驗前需要進(jìn)行方差齊性檢驗,F(xiàn)檢驗的功能就是進(jìn)行方差齊性檢驗蟋滴。
從兩個研究總體中隨機抽取樣本染厅,要對這兩個樣本進(jìn)行比較的時候,首先要判斷兩個總體方差是否相同津函,即方差齊性肖粮。若兩總體方差相等,則直接用t檢驗尔苦,若不等涩馆,可采用t’檢驗行施、變量變換或秩和檢驗等方法。
R中常用的三種F檢驗的方法魂那,bartlett.test方差齊性檢驗蛾号、var.test方差齊性檢驗、leveneTest方差齊性檢驗
前兩者是對原始數(shù)據(jù)的方差進(jìn)行檢驗涯雅,leveneTest是對方差模型的殘差進(jìn)行組間齊性檢驗须教。一般認(rèn)為是要求殘差的方差齊,所以一般的統(tǒng)計軟件都做的是leveneTest斩芭。

例:用高蛋白和低蛋白兩種飼料飼養(yǎng)1月齡的大白鼠轻腺,飼養(yǎng)3個月后,測定兩組大白鼠的增重量(g)划乖,兩組數(shù)據(jù)分別如下所示:
高蛋白組:134,146,106,119,124,161,107,83,113,129,97,123
低蛋白組:70,118,101,85,107,132,94
試問兩種飼料養(yǎng)殖的大白鼠增重量是否有顯著差異贬养?

x <- c(134,146,106,119,124,161,107,83,113,129,97,123)
y <- c(70,118,101,85,107,132,94)
#分組
group <- as.factor(c(rep(1, 12), rep(0, 7)))
z <- c(x, y)
df <- data.frame(z, group)

1、bartlett.test方差齊性檢驗

bartlett.test(z ~ group)

        Bartlett test of homogeneity of variances

data:  z by group
Bartlett's K-squared = 0.0066764, df = 1, p-value = 0.9349

2琴庵、var.test方差齊性檢驗

var.test(z ~ group)

F test to compare two variances

data:  z by group
F = 0.94107, num df = 6, denom df = 11, p-value = 0.9917
alternative hypothesis: true ratio of variances is not equal to 1
95 percent confidence interval:
 0.2425021 5.0909424
sample estimates:
ratio of variances 
          0.941066 

3误算、leveneTest方差齊性檢驗

library(car)
leveneTest(z, group)

Levene's Test for Homogeneity of Variance (center = median)
      Df F value Pr(>F)
group  1  0.0088 0.9264
      17    

三種檢驗結(jié)果都說明兩獨立樣本數(shù)據(jù)方差齊性,可以進(jìn)行獨立樣本T檢驗迷殿。

?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
  • 序言:七十年代末儿礼,一起剝皮案震驚了整個濱河市,隨后出現(xiàn)的幾起案子庆寺,更是在濱河造成了極大的恐慌蚊夫,老刑警劉巖,帶你破解...
    沈念sama閱讀 217,277評論 6 503
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件懦尝,死亡現(xiàn)場離奇詭異知纷,居然都是意外死亡,警方通過查閱死者的電腦和手機陵霉,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 92,689評論 3 393
  • 文/潘曉璐 我一進(jìn)店門琅轧,熙熙樓的掌柜王于貴愁眉苦臉地迎上來,“玉大人踊挠,你說我怎么就攤上這事乍桂。” “怎么了效床?”我有些...
    開封第一講書人閱讀 163,624評論 0 353
  • 文/不壞的土叔 我叫張陵睹酌,是天一觀的道長。 經(jīng)常有香客問我扁凛,道長忍疾,這世上最難降的妖魔是什么闯传? 我笑而不...
    開封第一講書人閱讀 58,356評論 1 293
  • 正文 為了忘掉前任谨朝,我火速辦了婚禮卤妒,結(jié)果婚禮上,老公的妹妹穿的比我還像新娘字币。我一直安慰自己则披,他們只是感情好,可當(dāng)我...
    茶點故事閱讀 67,402評論 6 392
  • 文/花漫 我一把揭開白布洗出。 她就那樣靜靜地躺著士复,像睡著了一般。 火紅的嫁衣襯著肌膚如雪翩活。 梳的紋絲不亂的頭發(fā)上阱洪,一...
    開封第一講書人閱讀 51,292評論 1 301
  • 那天,我揣著相機與錄音菠镇,去河邊找鬼冗荸。 笑死,一個胖子當(dāng)著我的面吹牛利耍,可吹牛的內(nèi)容都是我干的蚌本。 我是一名探鬼主播,決...
    沈念sama閱讀 40,135評論 3 418
  • 文/蒼蘭香墨 我猛地睜開眼隘梨,長吁一口氣:“原來是場噩夢啊……” “哼程癌!你這毒婦竟也來了?” 一聲冷哼從身側(cè)響起轴猎,我...
    開封第一講書人閱讀 38,992評論 0 275
  • 序言:老撾萬榮一對情侶失蹤嵌莉,失蹤者是張志新(化名)和其女友劉穎,沒想到半個月后捻脖,有當(dāng)?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體烦秩,經(jīng)...
    沈念sama閱讀 45,429評論 1 314
  • 正文 獨居荒郊野嶺守林人離奇死亡,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點故事閱讀 37,636評論 3 334
  • 正文 我和宋清朗相戀三年郎仆,在試婚紗的時候發(fā)現(xiàn)自己被綠了只祠。 大學(xué)時的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片。...
    茶點故事閱讀 39,785評論 1 348
  • 序言:一個原本活蹦亂跳的男人離奇死亡扰肌,死狀恐怖抛寝,靈堂內(nèi)的尸體忽然破棺而出,到底是詐尸還是另有隱情曙旭,我是刑警寧澤盗舰,帶...
    沈念sama閱讀 35,492評論 5 345
  • 正文 年R本政府宣布,位于F島的核電站桂躏,受9級特大地震影響钻趋,放射性物質(zhì)發(fā)生泄漏。R本人自食惡果不足惜剂习,卻給世界環(huán)境...
    茶點故事閱讀 41,092評論 3 328
  • 文/蒙蒙 一蛮位、第九天 我趴在偏房一處隱蔽的房頂上張望较沪。 院中可真熱鬧,春花似錦失仁、人聲如沸尸曼。這莊子的主人今日做“春日...
    開封第一講書人閱讀 31,723評論 0 22
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽控轿。三九已至,卻和暖如春拂封,著一層夾襖步出監(jiān)牢的瞬間茬射,已是汗流浹背。 一陣腳步聲響...
    開封第一講書人閱讀 32,858評論 1 269
  • 我被黑心中介騙來泰國打工冒签, 沒想到剛下飛機就差點兒被人妖公主榨干…… 1. 我叫王不留躲株,地道東北人。 一個月前我還...
    沈念sama閱讀 47,891評論 2 370
  • 正文 我出身青樓镣衡,卻偏偏與公主長得像霜定,于是被迫代替她去往敵國和親。 傳聞我的和親對象是個殘疾皇子廊鸥,可洞房花燭夜當(dāng)晚...
    茶點故事閱讀 44,713評論 2 354

推薦閱讀更多精彩內(nèi)容