如何在一个图上显示几个bootstrap置信区间

时间:2015-10-19 01:03:04

标签: r confidence-interval

如何在R中的一个图上绘制几个bootstrap置信区间? 我创建了一个Excel文档,其中包含置信区间的下限和上限。我想要一个情节,每年有两种颜色的置信区间(线条或方框都可以),用两种颜色表示A组或B组。

我尝试使用strcmp函数,但我尝试过的代码无效。

bwplot

数据:

bwplot(Lower+Upper~Year, data=yeargroupboot)

1 个答案:

答案 0 :(得分:0)

  • 我无法让bwplot以你想要的方式工作:bwplot(Lower+Upper~factor(Year)|Site,data=dd)是我最接近的
  • 您可以使用plotCI(来自plotrixgplots个套餐)手动执行此操作;
  • 或使用ggplot2,如下所示:

构建数据:

dd <- read.table(header=TRUE,text="
   Year     Upper      Lower   Site  
    2001     123        121      A  
    2001     115        113      B  
    2002     127        124      A  
    2002     114        113      B  
    2003     145        141      A  
    2003     100        99       B ")

library("ggplot2"); theme_set(theme_bw())
ggplot(dd)+
   geom_linerange(aes(Year,ymin=Lower,ymax=Upper,colour=Site),
                  position=position_dodge(width=0.25))+
   scale_x_continuous(breaks=2001:2003)