TukeyHSD 错误

Question

我有一个如下所示的数据，我尝试执行方差分析并检查所有列之间的差异。它们彼此之间的差异有多大等等。

df<- structure(list(color = structure(c(3L, 4L, 3L, 4L, 4L, 4L, 4L, 
    4L, 4L, 4L), .Label = c("B", "G", "R", "W"), class = "factor"), 
        type = 1:10, X1 = c(0.006605138, 0.001165448, 0.006975109, 
        0.002207839, 0.00187902, 0.002208638, 0.001199808, 0.001162252, 
        0.001338847, 0.001106317), X2 = c(0.006041392, 0.001639298, 
        0.006140877, 0.002958169, 0.002744017, 0.003107995, 0.001729594, 
        0.001582564, 0.001971713, 0.001693236), X3 = c(0.024180351, 
        0.002189061, 0.027377442, 0.002886651, 0.002816333, 0.003527908, 
        0.00231891, 0.001695633, 0.00212034, 0.001962923)), .Names = c("color", 
    "type", "X1", "X2", "X3"), row.names = c(NA, 10L), class = "data.frame")

首先我使用以下命令执行方差分析

 anovar= aov(type~.,df)

然后总结输出如下：

summary(anovar)

到目前为止，它表现得很好，表现也很好。但是，当我尝试执行 TukeyHSD 时，我似乎遇到了结构问题。错误如下。我搜索了，我找不到任何类似的情况。任何评论将不胜感激

TukeyHSD(anovar)
# Error in rep.int(n, length(means)) : unimplemented type 'NULL' in 'rep3'
# In addition: Warning messages:
# 1: In replications(paste("~", xx), data = mf) : non-factors ignored: X1
# 2: In replications(paste("~", xx), data = mf) : non-factors ignored: X2
# 3: In replications(paste("~", xx), data = mf) : non-factors ignored: X3

Answer 1

As it says in the description of the TukeyHSD documentation, the function creates a set of confidence intervals on the differences between the means of levels of a factor with the specified family-wise probability of coverage.

This means that you need to have factors in your data set in order to run it. So if you select the factor as follows it works:

> TukeyHSD(anovar, which = 'color') #color is the only categorical data

  Tukey multiple comparisons of means
    95% family-wise confidence level

Fit: aov(formula = type ~ ., data = df)

$color
     diff       lwr      upr     p adj
W-R 4.375 -1.465325 10.21532 0.1121168

You also get a warning that non factors X1, X2, X3 are ignored.

In order to print the TukeyHSD object just save it and use plot. There is a plot method (as well as a print method) for class TukeyHSD objects.

forplot <- TukeyHSD(anovar, which = 'color')
plot(forplot)