PTT數位生活區 / R_Language

[問題] type III anova

看板R_Language作者 (就這樣吧)時間11年前 (2013/05/02 13:39), 編輯推噓0(0推 0噓 29→)
留言29則, 2人參與, 最新討論串1/4 (看更多)
[問題類型]: 程式諮詢(我想用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
→
05/02 17:07, , 1F
mixed model 的作法主要有:
05/02 17:07, 1F
→
05/02 17:07, , 2F
1. 在 aov() 的 formula 中加入 Error()
05/02 17:07, 2F
→
05/02 17:08, , 3F
2. 使用 nlme:lme() 或 lme4:lmer()
05/02 17:08, 3F
→
05/02 17:09, , 4F
至於type3 anova，也因不同mixed model有不同實作法.
05/02 17:09, 4F
→
05/03 22:00, , 5F
我可以請教error()的規則嗎？
05/03 22:00, 5F
→
05/03 22:30, , 6F
Error()就是用來設定 f 值的分母用的.
05/03 22:30, 6F
→
05/03 22:33, , 7F
例如在 fixed A 與 random B 的完全隨機設計下,
05/03 22:33, 7F
→
05/03 22:47, , 8F
aov(Y ~ A + B + Error(A:B)) 可檢驗A和unrestricted B
05/03 22:47, 8F
→
05/03 22:47, , 9F
上述只限於平衡設計.
05/03 22:47, 9F
→
05/03 22:48, , 10F
不過說老實話, 要這樣操作真的需要數統背景才不會出錯.
05/03 22:48, 10F
→
05/03 22:48, , 11F
unbalanced mixed model 在 r 中遠不及商業分析軟體方便
05/03 22:48, 11F
→
05/03 22:51, , 12F
05/03 22:51, 12F
→
05/03 22:53, , 13F
mixed/random model 在不平衡設計往往需要藉由數個模型
05/03 22:53, 13F
→
05/03 22:54, , 14F
相互比較或擷取部份資訊才能得到完整解答.
05/03 22:54, 14F
→
05/03 22:55, , 15F
若手邊沒有教科書, 我自己操作也抖抖的, 很怕出錯.
05/03 22:55, 15F
→
05/03 22:57, , 16F
感謝，我看說明建議unbalance用lme，不過也沒很好懂XD
05/03 22:57, 16F
→
05/03 23:13, , 17F
或許你可以直接貼文你的實際需要, 我一起來解解看.
05/03 23:13, 17F
→
05/03 23:15, , 18F
lme的作者有出專書解說他的作品怎麼用.
05/03 23:15, 18F
→
05/03 23:16, , 19F
lmer適用性較廣但作者明確指出他故意讓lmer不直接求 f和p
05/03 23:16, 19F
→
05/03 23:17, , 20F
所以對一般使用者仍有不便之處.
05/03 23:17, 20F
→
05/03 23:17, , 21F
我的心得是目前mixed model在R中還蠻混沌的.
05/03 23:17, 21F
→
05/03 23:19, , 22F
老實說, 遇到複雜一點mixed模型我都跳槽用SAS. sign.
05/03 23:19, 22F
→
05/03 23:25, , 23F
不過, 唸過變方分析的人應該都明白 mixed model 本身就
05/03 23:25, 23F
→
05/03 23:25, , 24F
很不好處理. 要在R中設計與操作自然困難一點.
05/03 23:25, 24F
→
05/04 20:46, , 25F
感謝如此熱心的回應，我的實際需要就只是要算2 way anova
05/04 20:46, 25F
→
05/04 20:46, , 26F
的mixed model，只是用R和SPSS跑出來的MS不一樣 XD
05/04 20:46, 26F
→
05/04 21:33, , 27F
是 factorial unbalanced design 嗎?
05/04 21:33, 27F
→
05/04 23:27, , 28F
對factorial unbalanced design
05/04 23:27, 28F
→
05/04 23:31, , 29F
連 MS 都不一樣? 你可以貼文再說明你R怎麼做的嗎?
05/04 23:31, 29F
‣ 更多 DrRd 的文章
文章代碼(AID): #1HWVm0PP (R_Language)
R_Language 近期熱門文章
更多 近期熱門文章 >>
PTT數位生活區 即時熱門文章
更多 即時熱門文章 >>
‣ 更多 DrRd 的文章
文章代碼(AID): #1HWVm0PP (R_Language)