Tukey HSD對於多個變量和單個變量返回不同的結果

Question

我試圖為多變量數據集運行Tukey HSD。 但是，當我對單個變量運行相同的測試時，結果卻完全相反。

在運行多個變量時，我在ANOVA輸出中觀察到以下錯誤：

87個效果中有8個是無法估計的

運行單個變量時，我在ANOVA輸出中觀察到以下錯誤：

估計效果可能不平衡

這與我收到的完全相反的Tukey HSD輸出有什么關系嗎？ 另外，如何繼續解決這個問題？

我使用aov（）並在我的數據集中擁有近500000個數據點。

更具體地說，以下代碼給了我不同的結果：

代碼1：

lm_test1 <- lm(y ~ x1+ x2, data=data)
glht(lm_test1, linfct = mcp(x1 = "Tukey"))

代碼2：

lm_test1 <- lm(y ~ x1, data=data)
glht(lm_test1, linfct = mcp(x1 = "Tukey"))

請告訴我這怎么可能...

Answer 1

經過一些調查研究，我發現這個問題的答案，所以想我應該張貼此。 R中的Anova是默認類型-I anova。 因此，這意味着我們輸入的第一個變量將在不控制任何其他因素的情況下考慮效果，另一方面，對於其他變量，將在控制其他變量的影響后顯示結果。 因此，由於我被輸入我的變量作為第二變量，示出的結果控制用於第一變量這是偶然的，在完全相反的方向看的直接作用之后。