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

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

0 голосов
/ 27 сентября 2018

Я не понимаю порядок букв в функции 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))])
...