如何在R中比较带有每个区间的不同样本的置信区间和均值的数值数据来绘制图形？

Question

我想知道是否可以使用R自动绘制这样的图形：

我必须说以前的图像是使用Inkscape手动完成的。 但是，我想知道是否可以使用R自动构建此类图像。这将节省大量时间。

我有原始数据：

1   10.73   9.31    10.75   7.72    5.98    7.97
1   9.27    8.61    9.28    7.00    6.25    8.50
1   8.49    7.80    8.50    6.26    6.54    9.11
2   11.45   10.00   11.46   8.47    6.58    8.86
2   10.76   10.00   10.03   7.75    6.88    9.45
2   10.00   9.34    9.25    7.00    7.19    10.12

第一列是类别（在这种情况下，只有两个类别1和2）。 每个类别的上一行是较高的95％CL，每个类别的下一行是较低的95％CL。 其他行是平均值。

在这种情况下，我有6个不同的示例，这些示例显示在category列之后的6列中。

Answer 1

您可以通过以下方式在ggplot中执行此操作：首先重塑数据，使每一行都是一个置信区间，并为分组数据获取适当的x轴间距：

plot.dat <- data.frame(x=c(seq(.8, 1.2, length.out=6), seq(1.8, 2.2, length.out=6)),
                       lb=unlist(c(dat[3,-1], dat[6,-1])),
                       mid=unlist(c(dat[2,-1], dat[5,-1])),
                       ub=unlist(c(dat[1,-1], dat[4,-1])))
plot.dat
#       x    lb   mid    ub
# 1  0.80  8.49  9.27 10.73
# 2  0.88  7.80  8.61  9.31
# 3  0.96  8.50  9.28 10.75
# 4  1.04  6.26  7.00  7.72
# 5  1.12  6.54  6.25  5.98
# 6  1.20  9.11  8.50  7.97
# 7  1.80 10.00 10.76 11.45
# 8  1.88  9.34 10.00 10.00
# 9  1.96  9.25 10.03 11.46
# 10 2.04  7.00  7.75  8.47
# 11 2.12  7.19  6.88  6.58
# 12 2.20 10.12  9.45  8.86

然后，您可以使用geom_rect和geom_segment进行geom_segment以获取形状：

library(ggplot2)
ggplot(plot.dat, aes(xmin=x-.03, xmax=x+.03, ymin=lb, ymax=ub)) +
  geom_rect(fill=NA, color="black") +
  geom_segment(aes(x=x-.03, xend=x+.03, y=mid, yend=mid)) +
  theme_bw() +
  scale_x_continuous(breaks=c(1, 2)) +
  ylim(0, max(plot.dat$ub)) +
  xlab("Group") +
  ylab("Value")