r效应估计图按组作为误差条

effect_data <- data.frame(
  Level = factor(rep(1:4, each = 12)),
  # Assuming 12 sets of analyses
  EffectEstimate = runif(48, min = -2, max = 2),
  # Replace with your actual effect estimates
  LowerCI = runif(48, min = -3, max = -1),
  # Replace with your actual lower CIs
  UpperCI = runif(48, min = 1, max = 3),
  # Replace with your actual upper CIs
  M1  = factor(rep(letters[1:12], times = 4))
)

从你的描述来看，我不确定你在看什么. 这里有两个变种.也许其中一个就是你想要的？

library(tidyverse)

# Variant A
ggplot(effect_data,
       aes(
         x = Level,
         y = EffectEstimate,
         group = M1,
         fill = M1
       )) +
  geom_errorbar(aes(ymin = LowerCI, ymax = UpperCI),
                position = position_dodge(0.7),
                width = 0.25) +
  labs(title = "Effect Estimates with Confidence Intervals",
       x = "Results from 12 different Models",
       y = "Effect Estimate") +
  geom_bar(stat = "identity",
           position = "dodge",
           width = 0.7) +
  theme_minimal()

# Variant B
ggplot(effect_data,
       aes(
         x = M1,
         y = EffectEstimate,
         group = Level,
         fill = Level
       )) +
  geom_errorbar(aes(ymin = LowerCI, ymax = UpperCI),
                position = position_dodge(0.7),
                width = 0.25) +
  labs(title = "Effect Estimates with Confidence Intervals",
       x = "Results from 12 different Models",
       y = "Effect Estimate") +
  geom_bar(stat = "identity",
           position = "dodge",
           width = 0.7) +
  theme_minimal()