作者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