создание гистограммы ggplot в цикле ggplot

Question

этот вопрос является расширением предыдущей публикации: ggplot and loops

Я использовал приведенный выше пример для создания гистограмм в al oop. Я изменил приведенный выше пример, чтобы сгенерировать гистограмму и соответствующие столбцы ошибок. Я несколько преуспел, однако панели ошибок не заполняются в соответствии с отдельной переменной.

Буду очень признателен за помощь!

Мои модификации:

#make a dummy dataframe
D <- data.frame(
  x1 = runif(20),
  x2 = rnorm(20),
  x1_se = runif(20, 0.01, 0.09),
  x2_se = runif(20, -1, 1),
  treatment = rep(c("control","test"), each = 10)
)

# for reference later
p_names <- c("treatment","x1","x2")
se_names <- c("treatment","x1_se","x2_se")
trt <- rep(c("control","test"), each = 10)

# subset the standard error into its own dataframe
se <- D[,se_names]
names(se) <- str_remove(names(se), "_se")

plots <- list()

# the loop
for(nm in p_names) {

  #trt <- trt
  plots[[nm]] <- ggplot(data= D,  aes(x = trt, fill = trt))  + 
    geom_bar(aes_string(y = D[[nm]]), stat="identity", position = "dodge", color = "black") +
    geom_errorbar(aes(ymin= D[[nm]] - se[[nm]], 
                      ymax=   D[[nm]] + se[[nm]]), position=position_dodge(.9)) + ylab(nm)
}

print(plots[["x1"]])
print(plots[["x2"]])
<span class="math-container">```</span>

Answer 1

Не имеет смысла иметь один se для каждого наблюдения, если вы пытаетесь построить график среднего и se для каждого столбца, выполните следующие действия:

p_names = c("x1","x2")
for(nm in p_names) {
  plots[[nm]] <- ggplot(data= D,aes_string(x ="treatment",y=nm,fill="treatment"))+
    stat_summary(fun.y=mean,color = "black",geom="bar") +
    stat_summary(fun.data=mean_se,geom="errorbar",width=0.2)
}