Gnuplot顺畅的信心乐队

时间:2015-05-24 19:23:04

标签: gnuplot curve-fitting

根据这个问题给出的答案Gnuplot smooth confidence interval lines as opposed to error bars我能够得到相同的结果(y的误差是对称的,因此它是y加/减错误):

# x y errorY   
1   3   0.6 
2   5   0.4  
3   4   0.2
4   3.5 0.3

代码:

set style fill transparent solid 0.2 noborder
plot 'data.dat' using 1:($2-$3):($2+$3) with filledcurves title '95% confidence', \
     '' using 1:2 with lp lt 1 pt 7 ps 1.5 lw 3 title 'mean value'

现在通过连接每个y + errorY和y-errorY点来给出置信带。如果连接不仅仅是一条直线,而是一条平滑的线条,就像我可以用smooth csplines来平滑数据点那样我想要它。

1 个答案:

答案 0 :(得分:5)

这有点棘手,因为平滑仅适用于单个列,并且不能直接与filledcurves绘图样式结合使用。

因此,您必须首先生成两个临时数据文件,方法是将平滑的上下置信边界绘制为单独的数据文件

set table 'lower.dat'
plot 'data.dat' using 1:($2-$3) smooth cspline
set table 'upper.dat'
plot 'data.dat' using 1:($2+$3) smooth cspline
unset table

然后在绘制数据之前将这两个文件与paste lower.data upper.dat合并。如果您没有paste命令行程序,还可以使用paste.py之类的任何其他脚本来合并文件:

set terminal pngcairo
set output 'data.png'

set style fill transparent solid 0.2 noborder
plot '< paste lower.dat upper.dat' using 1:2:5 with filledcurves title '95% confidence', \
     'data.dat' using 1:2 with lines lt 1 smooth cspline title 'mean value',\
     '' using 1:2 with points lt 1 pt 7 ps 1.5 lw 3 title 'data points'

enter image description here