Неправильный порядок букв Tukey в пакете R multcompView

Question

Я не понимаю порядок букв в функции multcompletters из multcompView. Согласно документации это должно быть в соответствии со средним значением группы. В следующем примере средняя группа получила c (из abc) и должна была получить b. Это ошибка?

require(multcompView)
# Data 
datacol <- c(21.1,20.2,21.8,20.9,23.3,21.1,20.2,21.8,20.9,23.3,19.8,16.4,
16.9,16.0,17.6,17.5,16.9,13.3,18.0,17.6,13.5,12.2,15.2,15.1,15.2,14.0)

# Group 
faccol <- c(rep(c(1,2),each=10),rep(3,6))

# Combined Dataframe 
tukeyset <- data.frame(datacol,as.factor(faccol))
colnames(tukeyset)[2] <- "faccol"

# Tukeytest 
tukeyres <- TukeyHSD(x=aov(lm(datacol~faccol,data=tukeyset)))
Tlevels <- tukeyres$faccol[,4]
multcompLetters(Tlevels) # WRONG ORDER, even reversed 
# Boxplot 
boxplot(tukeyset$datacol~tukeyset$faccol)
# adding the labels 
text(x=c(1,2,3),y=c(aggregate(data=tukeyset,datacol~faccol,mean)$datacol),
     labels=as.character(multcompLetters(Tlevels,reversed=TRUE)$Letters)[order(names(multcompLetters(Tlevels,reversed=TRUE)['Letters']$Letters))])

Answer 1

О нет, только что столкнулся с этой проблемой! 2 часа потратил на то, чтобы разобраться наконец!

Вы не можете вызывать multcompLetters. Это даст вам совершенно неправильный порядок.

Вы должны использовать multcompLetters2, multcompLetters3 или multcompLetters4.

Еще один очень важный момент - вам нужно преобразовать входной набор данных в фрейм данных, но не в тиббл!Тиббл для этого не работает.