我有一个具有这种结构的数据框：
'data.frame': 39 obs. of 3 variables:
$ topic : Factor w/ 13 levels "Acido Folico",..: 1 2 3 4 5 6 7 8 9 10 ...
$ variable: Factor w/ 3 levels "Both","Preconception",..: 1 1 1 1 1 1 1 1 1 1 ...
$ value : int 14 1 36 17 5 9 19 9 19 25 ...
我想测试效果值〜变量，考虑到观察被分组在主题中。所以我想使用重复测量方差分析，其中＆＃34;变量＆＃34;被视为每个主题的重复措施。
电话是aov(value ~ variable + Error(topic/variable))。
到目前为止一切都好。
然后我想用glht(model, linfct= mcp(variable = 'Tukey'))执行事后测试，但我收到两个错误：
'glht'不支持类'aovlist'的对象 找不到“模型”的'model.matrix'方法！因为，取出错误术语解决错误，我认为这是问题。
那么，我怎样才能对重复测量anova进行事后测试？
谢谢！