作者salutyz (写给城市的诗)
看板R_Language
标题[问题] 用lsmeans做ANOVA post hoc分析
时间Mon Dec 26 10:38:32 2016
[问题类型]:
程式谘询(我想用R 做某件事情，但是我不知道要怎麽用R 写出来)
[软体熟悉度]:
请把以下不需要的部份删除
新手(没写过程式，R 是我的第一次)
[问题叙述]:
请简略描述你所要做的事情，或是这个程式的目的
有A,B,C三个factors的repeated anova (type III)，已经跑出anova的结果，想要继续做post hoc
要用lsmean和constrast的方式（这是依据同领域朋友的建议，不过他是用SAS里proc glm 去跑底下 lsmeans ( for simple main effect) 和contrast(for differences)
查了很多资料用下面的程式跑之後
在ref.grid()和lsmeans()都出现Error message
Warning message:
In lsm.basis.aovlist(object, trms, xlev, grid, ...) :
Some predictors are correlated with the intercept - results are biased.
May help to re-fit with different contrasts, e.g. 'contr.sum'
不确定要怎样修改才对..... orz
[程式范例]:
fit.1 <- aov(value~A*B*C+Error(Subject/(A*B*C)),data=work.1)
(late.rg1 <- ref.grid(fit.1))
lsmeans(fit.1, "A")
[环境叙述]:
请提供 sessionInfo() 的输出结果，
里面含有所有你使用的作业系统、R 的版本和套件版本资讯，
让版友更容易找出错误
R version 3.1.3 (2015-03-09)
Platform: x86_64-apple-darwin13.4.0 (64-bit)
Running under: OS X 10.10.5 (Yosemite)
locale:
[1] C
attached base packages:
[1] stats graphics grDevices utils datasets methods base
other attached packages:
[1] lsmeans_2.20-23 estimability_1.1-1 car_2.0-25 nlme_3.1-126
[关键字]:
选择性，也许未来有用
因为是超新手
所以折腾了好几个礼拜
跪求各位高手帮忙
感激不尽（跪）
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 111.240.212.180
※ 文章网址: https://webptt.com/cn.aspx?n=bbs/R_Language/M.1482719915.A.EE7.html
1F：→ celestialgod: 你资料问题，要看资料 12/26 18:53
2F：→ salutyz: 啊！是EEG data的资料，一个变数是condition(4 levels) 12/27 00:15
3F：→ salutyz: 另外两个变数是收资料的电极，分前後位置（5 levels) 12/27 00:16
4F：→ salutyz: 和左右（5 levels） 12/27 00:16
6F：→ andrew43: 但A*B*C要先平摊成简单主效应 12/27 02:40
7F：→ andrew43: 或是你说明白你预期结果是什麽 12/27 02:42