TukeyHSD 下标越界

Question

在运行进行看似简单的方差分析后，TukeyHSD 出现错误。

我的数据结构是下面的示例格式，我的实际数据中一共有5组：

data_frame:
A          B 
Group 1    2
Group 1    3
Group 1    5
Group 2    1
Group 2    7
Group 2    8

下面的结果来自 * 我的数据的实际输出，而不是上面的例子 *

aov(VA~as.factor(Etiologies),data_frame)
                          as.factor(Etiologies)     Residuals
Sum of Squares               37.85416               110.45051
Deg. of Freedom                     5                  81

Residual standard error: 1.167727
Estimated effects may be unbalanced

 summary(ANOVA_finalVA_all)
                               Df Sum Sq Mean Sq F value  Pr(>F)    
        as.factor(Etiologies)  5  37.85   7.571   5.552 0.00019 ***
         Residuals             81 110.45   1.364                    
        ---

然后当我运行 Tukey HSD 我得到以下错误

TukeyHSD(ANOVA_finalVA_all)
Error in FUN(X[[i]], ...) : subscript out of bounds

错误回溯：

8. lapply(args, "[[", "coefficients")
7. combine_mtables(...)
6. c.mtable(`Grand mean` = gmtable, tables)
5. c(`Grand mean` = gmtable, tables)
4. model.tables.aov(x, "means")
3. model.tables(x, "means")
2. TukeyHSD.aov(ANOVA_finalVA_all)
1. TukeyHSD(ANOVA_finalVA_all)

我猜 "coefficients" 有什么问题？这是方差分析 table

产生的系数

(Intercept)         Anova_VA_atFinal$EtiologiesA         Anova_VA_atFinal$EtiologiesB
1.73910734                           -0.78246714                            1.26089266 
Anova_VA_atFinal$EtiologiesC        Anova_VA_atFinal$EtiologiesD        Anova_VA_atFinal$EtiologiesE 
        0.07053282                            0.07662614                            1.09099566

据我所知，这似乎是正常的方差分析行为，并且我已确保我的组变量是因子。尽管 ANOVA 结果正常，但我似乎无法弄清楚为什么我会收到此错误。如果您能帮助解决此错误，我们将不胜感激！

Answer 1

问题出在包memisc上，使用前拆包即可避免此错误。

TukeyHSD 下标越界

TukeyHSD Subscript out of bounds

r

statistics

analysis

anova

tukey