Gnuplot - 使用直方图样式时如何将y值置于柱上方?

时间:2012-10-17 03:48:31

标签: gnuplot histogram

我目前正在使用脚本生成直方图,例如:

set style histogram cluster gap 4
plot for [COL=2:10] 'example.dat' u COL:xticlabels(1) title columnheader(COL)

现在我希望在直方图中的条形图上方添加y值(数字),但添加w labels会给出“此样式的列数不足”错误。

plot for [COL=2:10] 'example.dat' u COL:xticlabels(1) title columnheader(COL), \
    for [COL=2:10] 'example.dat' u COL title '' w labels

是否可以使用直方图样式添加y标签?

注意:我知道有一些示例用于绘制with boxes。如果可能的话,我希望使用直方图样式。

1 个答案:

答案 0 :(得分:9)

这是我提出的测试数据文件:

example.dat

 hi world foo bar baz qux
 1   2     3   4   5   6    
 4   5     7   3   6   5    

这是我用来绘制它的脚本:

set yrange [0:*]
GAPSIZE=4
set style histogram cluster gap 4
STARTCOL=2                 #Start plotting data in this column (2 for your example)
ENDCOL=6                   #Last column of data to plot (10 for your example)
NCOL=ENDCOL-STARTCOL+1     #Number of columns we're plotting 
BOXWIDTH=1./(GAPSIZE+NCOL) #Width of each box.
plot for [COL=STARTCOL:ENDCOL] 'example.dat' u COL:xtic(1) w histogram title columnheader(COL), \
    for [COL=STARTCOL:ENDCOL] 'example.dat' u (column(0)-1+BOXWIDTH*(COL-STARTCOL+GAPSIZE/2+1)-0.5):COL:COL notitle w labels

每个直方图簇在x轴上的总宽度为1个单位。我们知道我们需要多少宽度(盒子的数量+4,因为那是gapize)。我们可以计算每个框的宽度(1/(N+4))。然后我们将直方图绘制为正常。 (注意,我在绘图命令中添加了with histogram)。

根据内置帮助,标签需要3列数据(x y label)。在这种情况下,y位置和标签是相同的,可以直接从列COL中读取。第一个块的x位置以0为中心(总宽度为1)。因此,第一个块将位于x=-0.5+2*BOXWIDTH。这里的2是因为间隙是4个盒子宽度 - 左边两个,右边两个。下一个块将位于-0.5+3*BOXWIDTH等处。通常,(作为COL的函数)我们可以将其写为

-0.5+BOXSIZE*(COL-STARTCOL+1+GAPSIZE/2)

对于我们读取的每个附加块,我们需要将其向右移动1个单位。由于每个块对应于数据文件中的1行,因此我们可以使用伪列0(即column(0)$0),因为它会针对每个“记录/行”递增gnuplot读取。第0个记录保存标题,第一个记录保存第一个块。因为我们想要一个为第一条记录返回0的函数,所以我们使用column(0)-1。综上所述,我们发现x位置是:

(column(0)-1-0.5+BOXSIZE*(COL-STARTCOL+1+GAPSIZE/2))

这相当于我上面的内容。