[問題] type III anova

[DrRd (就這樣吧)]

[問題類型]:

程式諮詢(我想用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
→ andrew43:mixed model 的作法主要有:
→ andrew43:1. 在 aov() 的 formula 中加入 Error()
→ andrew43:2. 使用 nlme:lme() 或 lme4:lmer()
→ andrew43:至於type3 anova，也因不同mixed model有不同實作法.
→ DrRd:我可以請教error()的規則嗎？
→ andrew43:Error()就是用來設定 f 值的分母用的.
→ andrew43:例如在 fixed A 與 random B 的完全隨機設計下,
→ andrew43:aov(Y ~ A + B + Error(A:B)) 可檢驗A和unrestricted B
→ andrew43:上述只限於平衡設計.
→ andrew43:不過說老實話, 要這樣操作真的需要數統背景才不會出錯.
→ andrew43:unbalanced mixed model 在 r 中遠不及商業分析軟體方便
→ andrew43:可以看看 http://preview.tinyurl.com/cws373a
→ andrew43:mixed/random model 在不平衡設計往往需要藉由數個模型
→ andrew43:相互比較或擷取部份資訊才能得到完整解答.
→ andrew43:若手邊沒有教科書, 我自己操作也抖抖的, 很怕出錯.
→ DrRd:感謝，我看說明建議unbalance用lme，不過也沒很好懂XD
→ andrew43:或許你可以直接貼文你的實際需要, 我一起來解解看.
→ andrew43:lme的作者有出專書解說他的作品怎麼用.
→ andrew43:lmer適用性較廣但作者明確指出他故意讓lmer不直接求 f和p
→ andrew43:所以對一般使用者仍有不便之處.
→ andrew43:我的心得是目前mixed model在R中還蠻混沌的.
→ andrew43:老實說, 遇到複雜一點mixed模型我都跳槽用SAS.  sign.
→ andrew43:不過, 唸過變方分析的人應該都明白 mixed model 本身就
→ andrew43:很不好處理. 要在R中設計與操作自然困難一點.
→ DrRd:感謝如此熱心的回應，我的實際需要就只是要算2 way anova
→ DrRd:的mixed model，只是用R和SPSS跑出來的MS不一樣 XD
→ andrew43:是 factorial unbalanced design 嗎?
→ DrRd:對factorial unbalanced design
→ andrew43:連 MS 都不一樣? 你可以貼文再說明你R怎麼做的嗎?