不能在R中使用glht post hoc与重复测量方差分析？

我有一个具有这种结构的数据帧：

'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 ...

考虑到观察是按主题分组的，我想测试effect value~变量。所以我想使用重复测量方差分析，变量被认为是每个主题的重复测量。 调用是aovvalue~variable+Errortopic/variable

到现在为止一切都好

然后我想用glhtmodel执行一个事后测试，linfct=mcpvariable='Tukey'，但我收到两个错误：

“glht”不支持“aovlist”类的对象 找不到“model”的“model.matrix”方法！因为，去掉误差项，我想这就是问题所在

那么，如何对重复测量方差分析进行事后检验呢

---

谢谢

这个问题似乎离题了，因为它是关于选择一个合适的统计方法，而不是关于一个特定的编程问题。有趣的是，一个类似的问题在CrossValidated上被认为离题了，并移到了这里！无论如何，这与理论无关，但与glht的实现有关。我不清楚您的问题是否具体涉及glht函数的工作，或者是否有任何统计上有效的方法对重复测量方差分析进行事后测试。如果问题只是关于glht，那么请编辑您的问题以包含完整的。实际上，任何使用aov制作的模型，如果包含错误项，请在glht中生成此错误。