我有一个具有这种结构的数据框：
'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 ...
我想测试效果值〜变量，考虑到观察是按主题分组的。所以我想使用重复测量方差分析，其中“变量”被认为是每个主题的重复测量。电话是aov(value ~ variable + Error(topic/variable))。
到此为止一切正常。
然后我想用 执行事后测试glht(model, linfct= mcp(variable = 'Tukey'))，但我收到两个错误：
'glht' 不支持 'aovlist' 类的对象 没有找到 'model' 的 'model.matrix' 方法！因为，取出错误项解决了错误，我认为这就是问题所在。
那么，如何对重复测量 anova 执行事后测试？
谢谢！