方差齊次性檢驗Comparing Variances in R

在STHDA網(wǎng)站Comparing Variances in R 一文中哎榴,專門對正態(tài)性檢驗做了詳致的說明,翻譯并整理入下:

圖片.png

(一) F檢驗F-Test

F檢驗用于評估兩個總體(A和B)的方差是否相等。
F-Test: Compare Two Variances in R.

> var.test(len ~ supp, data = my_data)

    F test to compare two variances

data:  len by supp
F = 0.6386, num df = 29, denom df = 29, p-value = 0.2331
alternative hypothesis: true ratio of variances is not equal to 1
95 percent confidence interval:
 0.3039488 1.3416857
sample estimates:
ratio of variances 
         0.6385951 

(二) 比較方差的統(tǒng)計檢驗 Homogeneity of variances

有許多檢驗可以檢測不同組之間方差的均等性(均一性),包括:

  • F-test: Compare the variances of two samples. The data must be normally distributed.
  • Bartlett’s test: Compare the variances of k samples, where k can be more than two samples. The data must be normally distributed. The Levene test is an alternative to the Bartlett test that is less sensitive to departures from normality.
  • Levene’s test: Compare the variances of k samples, where k can be more than two samples. It’s an alternative to the Bartlett’s test that is less sensitive to departures from normality.
  • Fligner-Killeen test: a non-parametric test which is very robust against departures from normality.

Bartlett’s test用于測試k個樣本中方差的均勻性,其中k可以大于2叉袍。 適用于正態(tài)分布的數(shù)據(jù)。 當(dāng)數(shù)據(jù)分布為非正態(tài)分布時刽酱,下一部分將描述的Levene檢驗是Bartlett檢驗的更穩(wěn)健的替代方案喳逛。

2.1 Compute Bartlett’s test in R
# Bartlett’s test with one independent variable:
> bartlett.test(weight ~ group, data = PlantGrowth)

    Bartlett test of homogeneity of variances

data:  weight by group
Bartlett's K-squared = 2.8786, df = 2, p-value = 0.2371

# Bartlett’s test with multiple independent variables: 
> bartlett.test(len ~ interaction(supp,dose), data=ToothGrowth)

    Bartlett test of homogeneity of variances

data:  len by interaction(supp, dose)
Bartlett's K-squared = 6.9273, df = 5, p-value = 0.2261
2.2 Compute Levene’s test in R
library(car)
> # Levene's test with one independent variable
> leveneTest(weight ~ group, data = PlantGrowth)
Levene's Test for Homogeneity of Variance (center = median)
      Df F value Pr(>F)
group  2  1.1192 0.3412
      27               
> # Levene's test with multiple independent variables
> leveneTest(len ~ supp*dose, data = ToothGrowth)
Levene's Test for Homogeneity of Variance (center = median)
      Df F value Pr(>F)
group  5  1.7086 0.1484
      54               
2.3 Compute Fligner-Killeen test in R
> fligner.test(weight ~ group, data = PlantGrowth)

    Fligner-Killeen test of homogeneity of variances

data:  weight by group
Fligner-Killeen:med chi-squared = 2.3499, df = 2, p-value = 0.3088

參考資料:

  1. Comparing Variances in R
  2. 假設(shè)檢驗-方差齊性檢驗
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
  • 序言:七十年代末,一起剝皮案震驚了整個濱河市棵里,隨后出現(xiàn)的幾起案子润文,更是在濱河造成了極大的恐慌,老刑警劉巖殿怜,帶你破解...
    沈念sama閱讀 218,640評論 6 507
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件典蝌,死亡現(xiàn)場離奇詭異,居然都是意外死亡头谜,警方通過查閱死者的電腦和手機骏掀,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 93,254評論 3 395
  • 文/潘曉璐 我一進店門,熙熙樓的掌柜王于貴愁眉苦臉地迎上來柱告,“玉大人截驮,你說我怎么就攤上這事〖识龋” “怎么了葵袭?”我有些...
    開封第一講書人閱讀 165,011評論 0 355
  • 文/不壞的土叔 我叫張陵,是天一觀的道長乖菱。 經(jīng)常有香客問我坡锡,道長,這世上最難降的妖魔是什么窒所? 我笑而不...
    開封第一講書人閱讀 58,755評論 1 294
  • 正文 為了忘掉前任鹉勒,我火速辦了婚禮,結(jié)果婚禮上吵取,老公的妹妹穿的比我還像新娘禽额。我一直安慰自己,他們只是感情好海渊,可當(dāng)我...
    茶點故事閱讀 67,774評論 6 392
  • 文/花漫 我一把揭開白布。 她就那樣靜靜地躺著哲鸳,像睡著了一般臣疑。 火紅的嫁衣襯著肌膚如雪。 梳的紋絲不亂的頭發(fā)上徙菠,一...
    開封第一講書人閱讀 51,610評論 1 305
  • 那天讯沈,我揣著相機與錄音,去河邊找鬼。 笑死缺狠,一個胖子當(dāng)著我的面吹牛问慎,可吹牛的內(nèi)容都是我干的。 我是一名探鬼主播挤茄,決...
    沈念sama閱讀 40,352評論 3 418
  • 文/蒼蘭香墨 我猛地睜開眼如叼,長吁一口氣:“原來是場噩夢啊……” “哼!你這毒婦竟也來了穷劈?” 一聲冷哼從身側(cè)響起笼恰,我...
    開封第一講書人閱讀 39,257評論 0 276
  • 序言:老撾萬榮一對情侶失蹤,失蹤者是張志新(化名)和其女友劉穎歇终,沒想到半個月后社证,有當(dāng)?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體,經(jīng)...
    沈念sama閱讀 45,717評論 1 315
  • 正文 獨居荒郊野嶺守林人離奇死亡评凝,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點故事閱讀 37,894評論 3 336
  • 正文 我和宋清朗相戀三年追葡,在試婚紗的時候發(fā)現(xiàn)自己被綠了。 大學(xué)時的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片奕短。...
    茶點故事閱讀 40,021評論 1 350
  • 序言:一個原本活蹦亂跳的男人離奇死亡宜肉,死狀恐怖,靈堂內(nèi)的尸體忽然破棺而出篡诽,到底是詐尸還是另有隱情崖飘,我是刑警寧澤,帶...
    沈念sama閱讀 35,735評論 5 346
  • 正文 年R本政府宣布杈女,位于F島的核電站朱浴,受9級特大地震影響,放射性物質(zhì)發(fā)生泄漏达椰。R本人自食惡果不足惜翰蠢,卻給世界環(huán)境...
    茶點故事閱讀 41,354評論 3 330
  • 文/蒙蒙 一、第九天 我趴在偏房一處隱蔽的房頂上張望啰劲。 院中可真熱鬧梁沧,春花似錦、人聲如沸蝇裤。這莊子的主人今日做“春日...
    開封第一講書人閱讀 31,936評論 0 22
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽栓辜。三九已至恋拍,卻和暖如春,著一層夾襖步出監(jiān)牢的瞬間藕甩,已是汗流浹背施敢。 一陣腳步聲響...
    開封第一講書人閱讀 33,054評論 1 270
  • 我被黑心中介騙來泰國打工, 沒想到剛下飛機就差點兒被人妖公主榨干…… 1. 我叫王不留,地道東北人僵娃。 一個月前我還...
    沈念sama閱讀 48,224評論 3 371
  • 正文 我出身青樓概作,卻偏偏與公主長得像,于是被迫代替她去往敵國和親默怨。 傳聞我的和親對象是個殘疾皇子讯榕,可洞房花燭夜當(dāng)晚...
    茶點故事閱讀 44,974評論 2 355

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

  • 文章作者:Tyan博客:noahsnail.com | CSDN | 簡書 聲明:作者翻譯論文僅為學(xué)習(xí),如有侵權(quán)請...
    SnailTyan閱讀 2,432評論 0 5
  • 定義及分類 在進行 t 檢驗之前讓我們先看看它的定義:t 檢驗法就是在假設(shè)檢驗時利用 t 分布進行概率計算的檢驗方...
    研究僧小藍哥閱讀 26,490評論 1 13
  • pyspark.sql模塊 模塊上下文 Spark SQL和DataFrames的重要類: pyspark.sql...
    mpro閱讀 9,456評論 0 13
  • 好久沒來簡書了先壕,真有點對不起我的這十幾個粉絲瘩扼,因為最近我都要抖音上活動,所以忘了這里垃僚,更忘了要拍過程圖集绰,所以...
    唐小妞的水彩閱讀 224評論 0 3
  • 在生命歷程中,平凡的一天谆棺,也是我余生中栽燕,最寶貴的一天,最美麗的一天改淑。慵懶的睡個自然醒碍岔。八點多起床洗漱,喝水朵夏,吃...
    田輝_852a閱讀 324評論 0 0