r - R：使用 nLme 和 multcomp 分解交互

Question

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

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

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