【问题标题】:Overlapping confidence interval with ggplot用ggplot重叠置信区间
【发布时间】:2018-09-06 04:07:10
【问题描述】:

我想用 geom_ribbon 在PSOPSOA 中绘制两个置信区间。我的数据集由

给出
id   Meta prob       mean    lowerci   upperci
1   PSO  0.1 6705423.00   9913.939  151671.3
2   PSO  0.2 6705423.00  24352.031  393418.7
3   PSO  0.3 6705423.00  64719.335  444035.0
4   PSO  0.4 6705423.00  85058.168  600026.5
5   PSO  0.5 6705423.00 140437.916  756819.1
6   PSO  0.6 6705423.00 179236.196  952494.7
7   PSO  0.7 6705423.00 211278.350  773605.9
8   PSO  0.8 6705423.00 169915.851 1078624.1
9   PSO  0.9 6705423.00 263216.389  936007.2
10  PSO  1.0 6705423.00 266200.032 1061063.0
11 PSOA  0.1   52460.43   9913.939  151671.3
12 PSOA  0.2  202813.18  24352.031  393418.7
13 PSOA  0.3  252836.56  64719.335  444035.0
14 PSOA  0.4  329661.97  85058.168  600026.5
15 PSOA  0.5  450833.52 140437.916  756819.1
16 PSOA  0.6  424932.84 179236.196  952494.7
17 PSOA  0.7  486794.40 211278.350  773605.9
18 PSOA  0.8  634493.58 169915.851 1078624.1
19 PSOA  0.9  521509.18 263216.389  936007.2
20 PSOA  1.0  648183.78 266200.032 1061063.0

我尝试使用上面的代码,但无法在 PSO 中绘制 geom_ribbon(图中的“HERE”)。

p <- ggplot(data=dat2, aes(x=prob, y=mean, colour=Meta)) + geom_point() + geom_line()
p <- p + geom_ribbon(aes(ymin=dat2$lowerci, ymax=dat2$upperci), linetype=2, alpha=0.1)

【问题讨论】:

  • 您无法在此处绘制功能区,因为功能区的 y 值与 PSOA 组相同。因此,看起来您没有在绘制 PSO 组,但它正在被绘制。它只是被PSOA 组所覆盖
  • 你真的想绘制类似ymax = mean + upperciymin = mean - lowerci 的东西吗?
  • PSO 的均值相等,但置信区间的界限不同。另外,我想ymax = upperciymax = mean + upperci很大ymax = upperci
  • 是的,PSO 组内的限制不同。但是,uppercilowerciPSOPSOA 组中是相同的,导致 PSO 被掩盖
  • 好的,但是你的CIsPSO 组没有意义,因为平均值甚至没有包含在其中

标签: r ggplot2 confidence-interval


【解决方案1】:

问题在于lowerciupperci 的数据与PSOPSOA 组的数据完全相同。绘制时,PSOA 组正在“掩盖”PSO 组,使其看起来没有被绘制。要查看它是如何工作的,您可以只使用数据的子集运行相同的代码:

ggplot(data=dat2[dat2$Meta == "PSO",], aes(x=prob, y=mean, colour=Meta)) + geom_point() + geom_line() +
        geom_ribbon(aes(ymin=lowerci, ymax=upperci), linetype=2, alpha=0.1)

注意PSO 的功能区看起来与您的原始图相同吗?

【讨论】:

  • 啊,是的,你是对的,没有看实际数字。将删除我的答案
猜你喜欢
  • 2015-09-14
  • 2021-02-05
  • 1970-01-01
  • 1970-01-01
  • 1970-01-01
  • 1970-01-01
  • 1970-01-01
  • 1970-01-01
  • 2020-08-04
相关资源
最近更新 更多