【问题标题】:How to define lower and upper bounds when averaging and plotting values?平均和绘制值时如何定义下限和上限?
【发布时间】:2018-08-28 14:02:44
【问题描述】:

我不明白如何定义平均值的下限和上限,以便在 ggplot 中绘制它们,如 herehere 所示。

我的数据框如下所示:

 > data
   scenario decade FutureRangeSize.FullDisp
38    rcp85   2060                    85847
1     rcp45   2020                   376997
75    rcp85   2100                     6328
42    rcp45   2060                    68878
54    rcp85   2060                    44859
37    rcp45   2060                    70499
13    rcp85   2020                   192015
80    rcp85   2060                     1571
48    rcp45   2100                    61264
26    rcp85   2020                   123585
49    rcp85   2100                    72115
58    rcp85   2060                    32498
52    rcp45   2060                    52573
19    rcp85   2020                   139007
61    rcp85   2060                    17571
36    rcp45   2100                    85076
27    rcp45   2060                   234809
84    rcp85   2100                      123
44    rcp85   2020                    52735
71    rcp45   2100                     8967

所以有两个场景(rcp45,rcp85)和从 2020 年到 2100 年的十年。我想分别为每个场景计算每个十年的平均值,并将下限和上限绘制为一个图中的阴影区域。

谁能提供一个工作示例如何做到这一点?

【问题讨论】:

  • 首先,请提供可用于帮助的示例数据。使用提供的数据是不可能的。再包括几年。将dput 的输出粘贴到链接问题中。
  • 我添加了数据集的一个随机子集。
  • 你的上限和下限是如何定义的?是标准差吗?
  • 我想做 sd 或 min/max 值

标签: r ggplot2 plot


【解决方案1】:

你可以从一个简单的tidyverse开始

library(tidyverse)

d %>% 
 group_by(scenario, decade) %>% 
 summarise_at(vars(FutureRangeSize.FullDisp), funs(Mean=mean, Sd=sd)) %>% 
  mutate(Sd=ifelse(is.na(Sd),0,Sd)) %>% 
  ggplot(aes(x=decade, y=Mean, color=scenario, fill=scenario)) + 
    geom_point() + 
    geom_ribbon(aes(ymin = Mean - Sd, ymax= Mean + Sd), alpha=0.2)

或者尝试stat_summary 描述的here 方法

ggplot(d, aes(x=decade,y=FutureRangeSize.FullDisp, color=scenario)) + 
 geom_point() +
 stat_summary(fun.data = 'mean_sdl', geom = 'smooth', alpha=0.2)

【讨论】:

    猜你喜欢
    • 1970-01-01
    • 1970-01-01
    • 1970-01-01
    • 1970-01-01
    • 2020-11-09
    • 1970-01-01
    • 2015-08-29
    • 2013-08-28
    • 2020-11-28
    相关资源
    最近更新 更多