繪製直方圖條件

Question

我試圖做一個直方圖出的數據中，看起來像這樣的文件：

#Column 1 Column 2 
# 
0.0300  0.2126 
1.0000e-4 0.0104 
6.0000e-3 0.1299 
1.0000e-4 8.0600e-3 
1.0000e-4 0.0105 
0.0190  0.2204 
6.0000e-3 7.4900e-3 
1.0000e-4 0.0952 
6.0000e-3 7.4200e-3 
1.0000e-4 0.0131 
0.0190  0.3062 
0.0190  0.2561 
0.0300  0.9748 
0.0300  0.9406 
0.0300  0.0139 
1.0000e-4 0.0281 
0.0300  0.3625 
1.0000e-4 0.0945 
0.0300  0.5650 
1.0000e-4 0.1045 
6.0000e-3 0.2362 
1.0000e-4 0.0180 
1.0000e-4 0.1366 
1.0000e-4 0.0195 
0.0300  0.4652 
0.0190  0.3505 
0.0300  0.5146 
0.0190  0.4319 
6.0000e-3 0.2054 
6.0000e-3 0.2377 
0.0300  0.5281 
1.0000e-4 0.1128 
6.0000e-3 0.0623 

如果我使用的代碼：

n=20 #number of intervals 
max=0.03 #max value 
min=0 #min value 
width=(max-min)/n  #interval width 
hist(x,width)=width*floor(x/width)+width/2.0 

plot 'data' u (hist(\$1,width)):(1.0) smooth freq w boxes lc rgb "blue" lt 1 lw 0.5 notitle 

我得到正確的直方圖：

，但如果我使用條件線：

plot 'data' u (hist((\$2<=0.5?\$1:1/0),width)):(1.0) smooth freq w boxes lc rgb "blue" lt 1 lw 0.5 notitle 

我得到這個：

你可以看到，gnuplot沒有正確添加行，但他們密謀作爲單獨的列來代替。

有什麼辦法解決這個問題嗎？謝謝！

Answer 1

我懷疑這是gnuplot如何對待「缺失」數據的症狀。關於遺漏的數據，以下實際上略有不同：

plot 'data' u 1:2 w lines #connects lines across missing records 
plot 'data' u 1:($2) w lines #doesn't connect lines when a missing record is encountered 

我懷疑你看到的設計決定略有不同的症狀。 。

不幸的是它使典型gnuplot的數據過濾:-(沒用這裏幸運的是，你的條件很容易進入awk：

plot "< awk '{if ($2 <= 0.5) {print $0}}' test.dat " u (hist($1,width)):(1.0) smooth freq w boxes lc rgb "blue" lt 1 lw 0.5 notitle 

現在的gnuplot只看到你想讓它的數據（因爲它沒有看到任何NaN值，因此它不會創建任何新的計數器）。