1

我试图了解混合设计研究的交互。

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 ="**"

那么,有没有人需要查看使用这些工具的交互?

4

0 回答 0