1

运行看似简单的方差分析后，我在 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)

我猜“系数”有什么问题吗？这是从方差分析表产生的系数

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

4

1 回答 1

2

问题出在包 memisc 上，使用前先分离包，这个错误是可以避免的。

于 2017-11-23T20:16:14.957 回答