Результаты Тьюки на boxplot в R

Question

Я работаю над размещением писем, представляющих результаты HSD Тьюки на блокпосте, который я сделал. Я нашел описание того, как это сделать, здесь: галерея r-графа, но я получаю сообщение об ошибке после запуска функции.

Вот мои данные:

dput(WaterConDryMass)
structure(list(ChillTime = structure(c(1L, 1L, 1L, 1L, 2L, 2L, 
2L, 2L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 5L, 5L, 5L, 5L), .Label = c("Pre-chill", 
"6", "13", "24", "Post-chill"), class = "factor"), dmass = c(0.22, 
0.19, 0.34, 0.12, 0.23, 0.33, 0.38, 0.15, 0.31, 0.34, 0.45, 0.48, 
0.59, 0.54, 0.73, 0.69, 0.53, 0.57, 0.39, 0.8)), .Names = c("ChillTime", 
"dmass"), row.names = c(NA, -20L), class = "data.frame")

Необходимые пакеты:

install.packages('ggplot2')
install.packages('multcompView')

С моим кодом для бокса:

WaterConDryMass$ChillTime <- factor(WaterConDryMass$ChillTime, levels=c("Pre-chill", "6", "13", "24", "Post-chill"))
ggplot(WaterConDryMass, aes(x = ChillTime, y = dmass)) +
  geom_blank() +
  theme_bw() +
  theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank()) +
  labs(x = 'Time (weeks)', y = 'Water Content (DM %)') +
  ggtitle(expression(atop(bold("Water Content"), atop(italic("(Dry Mass)"), "")))) +
  theme(plot.title = element_text(hjust = 0.5, face='bold')) +
  annotate(geom = "rect", xmin = 1.5, xmax = 4.5, ymin = -Inf, ymax = Inf, alpha = 0.6, fill = "grey90") +
  geom_boxplot(fill = 'green2') +
  geom_vline(aes(xintercept=1.5), linetype="dashed") +
  geom_vline(aes(xintercept=4.5), linetype="dashed")

Производит это:

Для ANOVA:

Model4 <- aov(dmass~ChillTime, data=WaterConDryMass)

Тьюки HSD:

TUKEY <- TukeyHSD(Model4)
plot(TUKEY , las=1 , col="brown" )

Код до функции (где есть ошибка)

    # I need to group the treatments that are not different each other together.
    generate_label_df <- function(TUKEY, variable){

      # Extract labels and factor levels from Tukey post-hoc 
      Tukey.levels <- TUKEY[[variable]][,4]
      Tukey.labels <- data.frame(multcompLetters(Tukey.levels)['Letters'])

      #I need to put the labels in the same order as in the boxplot :
      Tukey.labels$ChillTime=rownames(Tukey.labels)
      Tukey.labels=Tukey.labels[order(Tukey.labels$ChillTime) , ]
      return(Tukey.labels)
    }

    # Apply the function on my dataset
    LABELS=generate_label_df(TUKEY , "WaterConDryMass$ChillT")

Ошибка:

 Error in strsplit(x, sep) : non-character argument 

Код из остальной части скрипта (до которого я еще не дошел):

# A panel of colors to draw each group with the same color :
my_colors=c( rgb(143,199,74,maxColorValue = 255),rgb(242,104,34,maxColorValue = 255), rgb(111,145,202,maxColorValue = 255),rgb(254,188,18,maxColorValue = 255) , rgb(74,132,54,maxColorValue = 255),rgb(236,33,39,maxColorValue = 255),rgb(165,103,40,maxColorValue = 255))

# Draw the basic boxplot
a=boxplot(data$value ~ data$treatment , ylim=c(min(data$value) , 1.1*max(data$value)) , col=my_colors[as.numeric(LABELS[,1])] , ylab="value" , main="")

# I want to write the letter over each box. Over is how high I want to write it.
over=0.1*max( a$stats[nrow(a$stats),] )

#Add the labels
text( c(1:nlevels(WaterConDryMass$ChillTime)) , a$stats[nrow(a$stats),]+over , LABELS[,1]  , col=my_colors[as.numeric(LABELS[,1])] )

В конечном итоге я включу свой конкретный блок-график в "# draw basic box plot" в последнем фрагменте кода, но пока не зашел так далеко.

Цель состоит в том, чтобы получить прямоугольник с буквами выше верхнего усика, как это, но не обязательно с согласованием цветов, как в примере:

Answer 1

РЕДАКТИРОВАТЬ: Вот начало, чтобы закончить код, скопированный с вашего вопроса, чтобы получить ваш сюжет.

Я должен был изменить ярлыки ChillTime в структуре вашего dataframe с самого начала, поэтому они используют подчеркивания, а не дефисы. Аналогично, когда вы конвертируете ChillTime к фактору - уровни не могут иметь дефисы для multcompLetters работать. Наконец, вам просто нужно поместить имя переменной в вашу функцию (ChillTime) скорее, чем WaterConDryMass$ChillT,

library(ggplot2)
library(multcompView)

WaterConDryMass <- structure(list(ChillTime = structure(c(1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 5L, 5L, 5L, 5L), .Label = c("Pre_chill", "6", "13", "24", "Post_chill"), class = "factor"), dmass = c(0.22, 0.19, 0.34, 0.12, 0.23, 0.33, 0.38, 0.15, 0.31, 0.34, 0.45, 0.48, 0.59, 0.54, 0.73, 0.69, 0.53, 0.57, 0.39, 0.8)), .Names = c("ChillTime", "dmass"), row.names = c(NA, -20L), class = "data.frame")

WaterConDryMass$ChillTime <- factor(WaterConDryMass$ChillTime, levels=c("Pre_chill", "6", "13", "24", "Post_chill"))

ggplot(WaterConDryMass, aes(x = ChillTime, y = dmass)) +
    geom_blank() +
    theme_bw() +
    theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank()) +
    labs(x = 'Time (weeks)', y = 'Water Content (DM %)') +
    ggtitle(expression(atop(bold("Water Content"), atop(italic("(Dry Mass)"), "")))) +
    theme(plot.title = element_text(hjust = 0.5, face='bold')) +
    annotate(geom = "rect", xmin = 1.5, xmax = 4.5, ymin = -Inf, ymax = Inf, alpha = 0.6, fill = "grey90") +
    geom_boxplot(fill = 'green2') +
    geom_vline(aes(xintercept=1.5), linetype="dashed") +
    geom_vline(aes(xintercept=4.5), linetype="dashed")

Model4 <- aov(dmass~ChillTime, data=WaterConDryMass)
TUKEY <- TukeyHSD(Model4)
plot(TUKEY , las=1 , col="brown" )

generate_label_df <- function(TUKEY, variable){

    # Extract labels and factor levels from Tukey post-hoc 
    Tukey.levels <- TUKEY[[variable]][,4]
    Tukey.labels <- data.frame(multcompLetters(Tukey.levels)['Letters'])

    #I need to put the labels in the same order as in the boxplot :
    Tukey.labels$ChillTime=rownames(Tukey.labels)
    Tukey.labels=Tukey.labels[order(Tukey.labels$ChillTime) , ]
    return(Tukey.labels)
}

# Apply the function on my dataset
LABELS=generate_label_df(TUKEY , "ChillTime")

# A panel of colors to draw each group with the same color :
my_colors=c( rgb(143,199,74,maxColorValue = 255),rgb(242,104,34,maxColorValue = 255), rgb(111,145,202,maxColorValue = 255),rgb(254,188,18,maxColorValue = 255) , rgb(74,132,54,maxColorValue = 255),rgb(236,33,39,maxColorValue = 255),rgb(165,103,40,maxColorValue = 255))

# Draw the basic boxplot
a=boxplot(WaterConDryMass$dmass ~ WaterConDryMass$ChillTime , ylim=c(min(WaterConDryMass$dmass) , 1.1*max(WaterConDryMass$dmass)) , col=my_colors[as.numeric(LABELS[,1])] , ylab="value" , main="")

# I want to write the letter over each box. Over is how high I want to write it.
over=0.1*max( a$stats[nrow(a$stats),] )

#Add the labels
text( c(1:nlevels(WaterConDryMass$ChillTime)) , a$stats[nrow(a$stats),]+over , LABELS[,1]  , col=my_colors[as.numeric(LABELS[,1])]