我试图了解混合设计研究的交互。
variable
是重复测量并且status
是主体间因子
library("multcomp", lib.loc="~/R/win-library/3.2")
library("nlme", lib.loc="C:/Program Files/R/R-3.2.0rc/library")
baseline<-lme(value~ 1, random = ~1|ID/variable/Status, data = SFAMixCon, method = "ML")
RepeatedAdded<-lme(value~ variable, random = ~1|ID/variable/Status, data = SFAMixCon, method = "ML")
BetweenAdded<-update(Variable, .~. + Status)
InteractionAdded<-update(VST, .~. + Status:variable)
有什么不同吗?:
anova(baseline,RepeatAdded, BetweenAdded, InteractionAdded)
发现任务有所不同,所以我想我会看看一些主要影响事后:
posthoc <-glht(Variable, linfct = mcp(variable = "Tukey"))
summary(posthoc)
confint(posthoc)
发现了一个重要的交互,那么为什么不看一些交互事后:
posthocX <-glht(VSTX, linfct = mcp(Status:variable ="Tukey"))
问题:
"**Error: unexpected '=' in "posthocX <-glht(VSTX, linfct = mcp(Status:variable ="**"
那么,有没有人需要查看使用这些工具的交互?