PTT
Submit
Submit
选择语言
正體中文
简体中文
PTT
R_Language
[问题] 用lsmeans做ANOVA post hoc分析
楼主:
salutyz
(写给城市的诗)
2016-12-26 10:38:32
[问题类型]:
程式咨询(我想用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
[关键字]:
选择性,也许未来有用
因为是超新手
所以折腾了好几个礼拜
跪求各位高手帮忙
感激不尽(跪)
作者:
celestialgod
(天)
2016-12-26 18:53:00
你资料问题,要看资料
楼主:
salutyz
(写给城市的诗)
2016-12-27 00:15:00
啊!是EEG data的资料,一个变量是condition(4 levels)另外两个变量是收资料的电极,分前后位置(5 levels)和左右(5 levels)
作者:
andrew43
(讨厌有好心推文后删文者)
2016-12-27 02:09:00
参考本板
#1No0yAxH
但A*B*C要先平摊成简单主效应或是你说明白你预期结果是什么
继续阅读
Re: [问题] 不同资料集字段比对,再输出比对结果
celestialgod
[问题] 不同资料集字段比对,再输出比对结果
joson4921
Re: [问题] 多个变量的循环如何处理?
celestialgod
Re: [问题] 多个变量的循环如何处理?
abc2090614
Re: [问题] 多个变量的循环如何处理?
f496328mm
[问题] 多个变量的循环如何处理?
lovedmagic
[问题] 阶层资料整理
criky
[问题] ggplot个别图型合并后存盘问题
april6310
Re: [问题] Ordered Probit for panel data
v303ofsky
[问题] Ordered Probit for panel data
v303ofsky
Links
booklink
Contact Us: admin [ a t ] ucptt.com