Type I, Type II, and Type III ANOVAs

14.5 Type I, Type II, and Type III ANOVAs

It turns out that there is not just one way to calculate ANOVAs. In fact, there are three different types - called, Type 1, 2, and 3 (or Type I, II and III). These types differ in how they calculate variability (specifically the sums of of squares). If your data is relatively balanced, meaning that there are relatively equal numbers of observations in each group, then all three types will give you the same answer. However, if your data are unbalanced, meaning that some groups of data have many more observations than others, then you need to use Type II (2) or Type III (3).

事實證明,不只有一種方法可以計算方差撑刺。事實上,有三種不同的類型-稱為類型1芬膝、類型2和類型3(或類型I、II和III)形娇。這些類型在計算變異性的方式上有所不同(具體地說锰霜,就是“平方和”)。

均衡設計(balanced design)是指 每個實驗條件下的被試量相等桐早。否則就是非均衡設計(unbalanced design)癣缅。

如果你的數據相對比較平衡,也就是說每組的觀察次數相對相等哄酝,那么這三種類型的數據都會給你相同的答案友存。但是,如果您的數據是“不平衡的”炫七,即某些數據組的觀測值比其他組多得多爬立,那么您需要使用類型II(2)或類型III(3)。

  • 類型I方法(序貫型) 計算ANOVA效應:舉例如下:
    y~A+B+A:B
    A對y的影響万哪;A不做調整
    控制A時侠驯,B對y的影響;B根據A調整
    控制A和B的主效應時奕巍,A與B的交互效應吟策。A:B交互項根據A和B調整。
    R中的aov()默認使用的是類型I的止。
  • 類型II(分層型)
    效應根據同水平或低水平的效應做調整檩坚。
    A根據B調整,B根據A調整诅福,A:B交互項同時根據A和B調整匾委。
  • 類型III(邊界型)
    每個效應根據模型其他各效應做相應調整。
    A根據B和A:B做調整氓润,A:B交互項根據A和B調整赂乐。
    SPSS 和 SAS 默認使用的是類型III。如果要在R中得到相同的結果咖气,需要使用類型III挨措。aov()函數默認使用的是類型I的方法,若要使用類型III崩溪,需要用car包中的Anova()函數浅役,具體可參考 help(Anova,package="car")

樣本大小越不平衡伶唯,效應項的順序對結果的影響越大觉既。一般來說,越基礎性的效應越需要放在表達式前面。具體來講:
首先是協變量瞪讼,然后是主效應岭参,接著是雙因素的交互項,再接著是三因素的交互項尝艘,以此類推。
對于主效應姿染,越基礎性的效應越應放在表達式前面背亥。

The standard aov() function in base-R uses Type I sums of squares. Therefore, it is only appropriate when your data are balanced. If your data are unbalanced, you should conduct an ANOVA with Type II or Type III sums of squares. To do this, you can use the Anova() function in the car package. The Anova() function has an argument called type that allows you to specify the type of ANOVA you want to calculate.

R中的標準aov()函數使用第一類平方和。因此悬赏,只有當您的數據達到平衡時狡汉,它才是合適的。如果你的數據不平衡闽颇,你應該用II型平方和或III型平方和進行方差分析盾戴。要做到這一點,可以使用car包中的Anova()函數兵多。Anova()函數有一個名為type的參數尖啡,允許您指定要計算的方差分析的類型。

In the next code chunk, I’ll calculate 3 separate ANOVAs from the poopdeck data using the three different types. First, I’ll create a regression object with lm(). As you’ll see, the Anova() function requires you to enter a regression object as the main argument, and not a formula and dataset. That is, you need to first create a regression object from the data with lm() (or glm()), and then enter that object into the Anova() function. You can also do the same thing with the standard aov() function`.

在下一個代碼塊中剩膘,我將使用三種不同的類型從poopdeck數據計算3個獨立的ANOVA衅斩。首先,我將使用lm()創(chuàng)建一個回歸對象怠褐。正如您將看到的畏梆,Anova()函數要求您輸入回歸對象作為主參數,而不是輸入公式和數據集奈懒。也就是說奠涌,您需要先使用lm()函數(或glm())從數據創(chuàng)建一個回歸對象,然后將該對象輸入到Anova()函數中磷杏。標準的aov()函數也可以做同樣的事情溜畅。

# Step 1: Calculate regression object with lm()
time.lm <- lm(formula = time ~ type + cleaner,
              data = poopdeck)

Now that I’ve created the regression object time.lm, I can calculate the three different types of ANOVAs by entering the object as the main argument to either aov() for a Type I ANOVA, or Anova() in the car package for a Type II or Type III ANOVA:

現在我已經創(chuàng)建了回歸對象time.lm,我可以通過將該對象作為主參數輸入到aov(對于類型I ANOVA)或者car包(對于類型II或類型III ANOVA)中的Anova()參數來計算三種不同類型的ANOVA:

# Type I ANOVA - aov()
time.I.aov <- aov(time.lm)

# Type II ANOVA - Anova(type = 2)
time.II.aov <- car::Anova(time.lm, type = 2)

# Type III ANOVA - Anova(type = 3)
time.III.aov <- car::Anova(time.lm, type = 3)

As it happens, the data in the poopdeck dataframe are perfectly balanced (so we’ll get exactly the same result for each ANOVA type. However, if they were not balanced, then we should not use the Type I ANOVA calculated with the aov() function.

碰巧的是茴丰,poopdeck數據幀中的數據是完全平衡的(因此达皿,對于每種ANOVA類型,我們都會得到完全相同的結果贿肩。但是峦椰,如果它們不平衡,那么我們就不應該使用使用aov()函數計算的I型方差分析汰规。

要查看數據是否平衡汤功,可以使用以下函數:
To see if your data are balanced, you can use the function:

# Are observations in the poopdeck data balanced?
with(poopdeck,
     table(cleaner, type))
##        type
## cleaner parrot shark
##       a    100   100
##       b    100   100
##       c    100   100

As you can see, in the poopdeck data, the observations are perfectly balanced, so it doesn’t matter which type of ANOVA we use to analyse the data.

正如你所看到的,在poopdeck數據中溜哮,觀測數據是完全平衡的滔金,所以我們使用哪種類型的方差分析來分析數據并不重要色解。

有關不同類型的更多詳細信息,請查看https://mcfromnz.wordpress.com/2011/03/02/anova-type-iiiiii-ss-explained/.

For more detail on the different types, check out https://mcfromnz.wordpress.com/2011/03/02/anova-type-iiiiii-ss-explained/.

參考資料:https://bookdown.org/ndphillips/YaRrr/type-i-type-ii-and-type-iii-anovas.html

最后編輯于
?著作權歸作者所有,轉載或內容合作請聯系作者
  • 序言:七十年代末餐茵,一起剝皮案震驚了整個濱河市科阎,隨后出現的幾起案子,更是在濱河造成了極大的恐慌忿族,老刑警劉巖锣笨,帶你破解...
    沈念sama閱讀 206,126評論 6 481
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件,死亡現場離奇詭異道批,居然都是意外死亡错英,警方通過查閱死者的電腦和手機,發(fā)現死者居然都...
    沈念sama閱讀 88,254評論 2 382
  • 文/潘曉璐 我一進店門隆豹,熙熙樓的掌柜王于貴愁眉苦臉地迎上來椭岩,“玉大人,你說我怎么就攤上這事璃赡∨懈纾” “怎么了?”我有些...
    開封第一講書人閱讀 152,445評論 0 341
  • 文/不壞的土叔 我叫張陵鉴吹,是天一觀的道長姨伟。 經常有香客問我,道長豆励,這世上最難降的妖魔是什么夺荒? 我笑而不...
    開封第一講書人閱讀 55,185評論 1 278
  • 正文 為了忘掉前任,我火速辦了婚禮良蒸,結果婚禮上技扼,老公的妹妹穿的比我還像新娘。我一直安慰自己嫩痰,他們只是感情好剿吻,可當我...
    茶點故事閱讀 64,178評論 5 371
  • 文/花漫 我一把揭開白布。 她就那樣靜靜地躺著串纺,像睡著了一般丽旅。 火紅的嫁衣襯著肌膚如雪。 梳的紋絲不亂的頭發(fā)上纺棺,一...
    開封第一講書人閱讀 48,970評論 1 284
  • 那天榄笙,我揣著相機與錄音,去河邊找鬼祷蝌。 笑死茅撞,一個胖子當著我的面吹牛,可吹牛的內容都是我干的。 我是一名探鬼主播米丘,決...
    沈念sama閱讀 38,276評論 3 399
  • 文/蒼蘭香墨 我猛地睜開眼剑令,長吁一口氣:“原來是場噩夢啊……” “哼!你這毒婦竟也來了拄查?” 一聲冷哼從身側響起吁津,我...
    開封第一講書人閱讀 36,927評論 0 259
  • 序言:老撾萬榮一對情侶失蹤,失蹤者是張志新(化名)和其女友劉穎堕扶,沒想到半個月后腺毫,有當地人在樹林里發(fā)現了一具尸體,經...
    沈念sama閱讀 43,400評論 1 300
  • 正文 獨居荒郊野嶺守林人離奇死亡挣柬,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內容為張勛視角 年9月15日...
    茶點故事閱讀 35,883評論 2 323
  • 正文 我和宋清朗相戀三年,在試婚紗的時候發(fā)現自己被綠了睛挚。 大學時的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片邪蛔。...
    茶點故事閱讀 37,997評論 1 333
  • 序言:一個原本活蹦亂跳的男人離奇死亡,死狀恐怖扎狱,靈堂內的尸體忽然破棺而出侧到,到底是詐尸還是另有隱情,我是刑警寧澤淤击,帶...
    沈念sama閱讀 33,646評論 4 322
  • 正文 年R本政府宣布匠抗,位于F島的核電站,受9級特大地震影響污抬,放射性物質發(fā)生泄漏汞贸。R本人自食惡果不足惜,卻給世界環(huán)境...
    茶點故事閱讀 39,213評論 3 307
  • 文/蒙蒙 一印机、第九天 我趴在偏房一處隱蔽的房頂上張望矢腻。 院中可真熱鬧,春花似錦射赛、人聲如沸多柑。這莊子的主人今日做“春日...
    開封第一講書人閱讀 30,204評論 0 19
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽竣灌。三九已至,卻和暖如春秆麸,著一層夾襖步出監(jiān)牢的瞬間初嘹,已是汗流浹背。 一陣腳步聲響...
    開封第一講書人閱讀 31,423評論 1 260
  • 我被黑心中介騙來泰國打工蛔屹, 沒想到剛下飛機就差點兒被人妖公主榨干…… 1. 我叫王不留削樊,地道東北人。 一個月前我還...
    沈念sama閱讀 45,423評論 2 352
  • 正文 我出身青樓,卻偏偏與公主長得像漫贞,于是被迫代替她去往敵國和親甸箱。 傳聞我的和親對象是個殘疾皇子,可洞房花燭夜當晚...
    茶點故事閱讀 42,722評論 2 345

推薦閱讀更多精彩內容