作者DrRd (就這樣吧)
看板R_Language
標題[問題] type III anova
時間Thu May 2 13:39:09 2013
[問題類型]:
程式諮詢(我想用R 做某件事情，但是我不知道要怎麼用R 寫出來)
[軟體熟悉度]:
請把以下不需要的部份刪除
新手(沒寫過程式，R 是我的第一次)
[問題敘述]:
要進行type III的變異數分析
共有A、B兩個因子，其中A是固定效果，B是隨機效果
網路上有看到可以用car裡面的Anova(,type=3)
但它沒有讓我指定何者為固定何者為隨機
算出來的結果也不正確
請問除了手算之外，有其他的套件可以算嗎？
[程式範例]:
library(car)
mod = aov(y~a*b)
Anova(mod,type=3)
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.112.4.187
1F：→ andrew43:mixed model 的作法主要有: 05/02 17:07
2F：→ andrew43:1. 在 aov() 的 formula 中加入 Error() 05/02 17:07
3F：→ andrew43:2. 使用 nlme:lme() 或 lme4:lmer() 05/02 17:08
4F：→ andrew43:至於type3 anova，也因不同mixed model有不同實作法. 05/02 17:09
5F：→ DrRd:我可以請教error()的規則嗎？ 05/03 22:00
6F：→ andrew43:Error()就是用來設定 f 值的分母用的. 05/03 22:30
7F：→ andrew43:例如在 fixed A 與 random B 的完全隨機設計下, 05/03 22:33
8F：→ andrew43:aov(Y ~ A + B + Error(A:B)) 可檢驗A和unrestricted B 05/03 22:47
9F：→ andrew43:上述只限於平衡設計. 05/03 22:47
10F：→ andrew43:不過說老實話, 要這樣操作真的需要數統背景才不會出錯. 05/03 22:48
11F：→ andrew43:unbalanced mixed model 在 r 中遠不及商業分析軟體方便 05/03 22:48
13F：→ andrew43:mixed/random model 在不平衡設計往往需要藉由數個模型 05/03 22:53
14F：→ andrew43:相互比較或擷取部份資訊才能得到完整解答. 05/03 22:54
15F：→ andrew43:若手邊沒有教科書, 我自己操作也抖抖的, 很怕出錯. 05/03 22:55
16F：→ DrRd:感謝，我看說明建議unbalance用lme，不過也沒很好懂XD 05/03 22:57
17F：→ andrew43:或許你可以直接貼文你的實際需要, 我一起來解解看. 05/03 23:13
18F：→ andrew43:lme的作者有出專書解說他的作品怎麼用. 05/03 23:15
19F：→ andrew43:lmer適用性較廣但作者明確指出他故意讓lmer不直接求 f和p 05/03 23:16
20F：→ andrew43:所以對一般使用者仍有不便之處. 05/03 23:17
21F：→ andrew43:我的心得是目前mixed model在R中還蠻混沌的. 05/03 23:17
22F：→ andrew43:老實說, 遇到複雜一點mixed模型我都跳槽用SAS. sign. 05/03 23:19
23F：→ andrew43:不過, 唸過變方分析的人應該都明白 mixed model 本身就 05/03 23:25
24F：→ andrew43:很不好處理. 要在R中設計與操作自然困難一點. 05/03 23:25
25F：→ DrRd:感謝如此熱心的回應，我的實際需要就只是要算2 way anova 05/04 20:46
26F：→ DrRd:的mixed model，只是用R和SPSS跑出來的MS不一樣 XD 05/04 20:46
27F：→ andrew43:是 factorial unbalanced design 嗎? 05/04 21:33
28F：→ DrRd:對factorial unbalanced design 05/04 23:27
29F：→ andrew43:連 MS 都不一樣? 你可以貼文再說明你R怎麼做的嗎? 05/04 23:31