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)

发现任务有所不同，所以我想我会看一些主要的效果posthocs：

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

那么，是否有人必须使用这些工具来查看交互？