r效应估计平铺图

发布于10月01日

我喜欢学习如何再现这个情节.我不知道这叫什么，但我有兴趣使用这个数据集重新创建它，其中X是我的预测变量，M1是我的分组变量.

所以Fladen=a，LandSort=b，Utlagan=c，等等

TM、RC、S2、S1、MM的级别为1、2、3、4

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(ggplot2)

ggplot(effect_data, aes(y = Level, x = EffectEstimate)) +
  geom_errorbarh(aes(xmin = LowerCI, xmax = UpperCI), 
                 height = 0,  #removing ticks
                 position = position_dodge(width = 0.3), 
                 color = "gray", alpha = 0.5,
                 size = 3) + 
  geom_point(color = "black", size = 3, position = position_dodge(width = 0.3)) +
  facet_wrap(~ M1, scales = "free_y", ncol = 4) +
  labs(title = "Effect Estimates with 95% CI", y = "Effect Estimate", x = "Level") +
  theme_minimal() +
  theme(
    legend.position = "none",
    strip.background = element_rect(fill="gray", color="black"), 
    strip.text = element_text(size=14, color="black"),  
    panel.border = element_rect(color = "black", fill = NA, linewidth = 1.5)
  )