使用geom_segment为现有绘图上的选定数据点绘制箭头

时间:2015-10-12 21:24:28

标签: r plot ggplot2 upperbound

我使用以下数据在R中使用ggplot2包创建了一个图:

# Ind filt objs        flux  filt    objs err_flux#

01   4590 obj1 1005.448892  4590 errobj1   0.0401  
02   6220 obj1 1420.626789  6220 errobj1   0.0392  
03   7640 obj1 1855.581355  7640 errobj1   0.0432  
 *.......skipped lines...*    
20 21512 obj2  949.642188 21512 errobj2     0.1516    
21  2030 obj2    9.838299  2030 errobj2   **0.0000**  
22  2253 obj2   17.097003  2253 errobj2   **0.0000**  
23  2612 obj2   14.754347  2612 errobj2   **0.0000**  
24  3470 obj2   14.890868  3470 errobj2   **0.0000**  

我正在使用ggplot绘制 flux vs filt ,然后添加 err_flux 列中定义的y-error栏。 我想从 err_flux == 0。的数据点中绘制箭头。

到目前为止,这是我的计划:

g <- ggplot(finaldata[1:26, ], aes(x = filt, y = log10(flux), color = objs))
gf1 <- g + geom_point(size = 2, alpha = 0.8, pch = 15) + theme_bw()
gf <- gf1 + geom_errorbar(aes(ymax = log10(flux) + err_flux,
                              ymin = log10(flux) - err_flux),
                           size = 0.5, color = "darkgrey")
gff <-gf + xlab("wavelength") + ylim(0.3,4.5) +
          geom_smooth(method = "loess", se = FALSE)
gfinal <- gff + scale_color_hue(l = 10, c = 75)

结果图如下所示: enter image description here

以下是我试图添加箭头的内容(例如上图所示的'obj2'中最左边的4个点,但它不起作用:

for (i in 1:nrow(finaldata)) {
 if (finaldata$err_flux==0.){
   gfinal<-gfinal + geom_segment(x=finaldata$filt[i], y=finaldata$flux[i], xend=finaldata$filt[i], yend=finaldata$flux[i]-0.5,arrow=arrow(length=unit(0.25,"cm")))
  }
 }

此外,如果我将'geom_segment'应用于我的实际数据的子集,由于尺寸不同,它将不会添加到我现有的绘图( gfinal )。

任何建议/帮助都将不胜感激。

1 个答案:

答案 0 :(得分:2)

Annotate函数帮助我从选定的数据点向下绘制箭头(上界)。

我替换了我的程序的最后一部分,该部分使用了带有geom_segment的for和if循环到annotate,如下所示:

sub<-subset(finaldata[1:26,],finaldata[1:26,]$err_flux==0.)
gfinal+annotate("segment",x=sub$filt,xend=sub$filt,y=log10(sub$flux),yend=log10(sub$flux)-0.25,color="red",arrow = arrow(length = unit(0.2, "cm")))

导致了下图:See 4 <code>red</code> arrows on 4 leftmost data points for "obj2"