Question

在运行看似简单的ANOVA后，我收到了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)

我认为＆＃34;系数＆＃34;？这是从ANOVA表中产生的系数

(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行为，我确保我的组变量是因子。尽管ANOVA结果正常，我似乎无法弄清楚为什么我会收到此错误。任何帮助解决此错误的任何帮助将不胜感激！

Answer 1

问题在于包装memisc，在使用前分离包装，可以避免这种错误。

TukeyHSD下标超出范围

1 个答案: