[问题类型]:
程式咨询(我想用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
[关键字]:
选择性，也许未来有用
因为是超新手
所以折腾了好几个礼拜
跪求各位高手帮忙
感激不尽（跪）