[gnuplot] gnuplot을 사용하는 히스토그램?

매우주의하십시오 :이 페이지의 모든 답변은 비닝이 시작되는 위치 (원하는 경우 가장 왼쪽 빈의 왼쪽 가장자리)를 사용자의 손에서 내재적으로 결정합니다. 사용자가 데이터를 비닝하기 위해 이러한 기능 중 하나를 비닝 시작 위치에 대한 자신의 결정과 결합하면 (위의 링크 된 블로그에서 수행됨) 위의 기능이 모두 올바르지 않습니다. ‘Min’비닝을위한 임의의 시작점이있는 경우 올바른 기능은 다음과 같습니다.

bin(x) = width*(floor((x-Min)/width)+0.5) + Min

이것이 순차적으로 올바른 이유를 알 수 있습니다 (몇 개의 쓰레기통과 그 안에 어딘가에 점을 그리는 데 도움이 됨). 데이터 포인트에서 Min을 빼서 비닝 범위까지 얼마나되는지 확인합니다. 그런 다음 ‘폭’단위로 효과적으로 작업 할 수 있도록 이진으로 나눕니다. 그런 다음 결과를 ‘바닥’으로 만들어 해당 빈의 왼쪽 가장자리로 이동하고 0.5를 추가하여 빈의 중간으로 이동하여 너비를 곱하여 더 이상 빈 단위로 작업하지 않고 절대 스케일로 작업하십시오. 다시 마지막으로 시작 부분에서 뺀 최소 오프셋을 다시 추가하십시오.

이 기능을 실제로 고려하십시오.

Min = 0.25 # where binning starts
Max = 2.25 # where binning ends
n = 2 # the number of bins
width = (Max-Min)/n # binwidth; evaluates to 1.0
bin(x) = width*(floor((x-Min)/width)+0.5) + Min

예를 들어 1.1 값은 실제로 왼쪽 빈에 해당합니다.

• 이 함수는 왼쪽 빈 (0.75)의 중앙에 올바르게 매핑합니다.
• Born2Smile의 답변 bin (x) = width * floor (x / width)는 1에 잘못 매핑됩니다.
• mas90의 대답 bin (x) = width * floor (x / width) + binwidth / 2.0은 1.5에 잘못 매핑됩니다.

Binn 경계가 (n + 0.5) * binwidth (n이 정수를 넘음)에서 발생하는 경우 Born2Smile의 답은 정확합니다. bin 경계가 n * binwidth에서 발생하는 경우에만 mas90의 답이 정확합니다.

답변

n=100 #number of intervals
max=3. #max value
min=-3. #min value
width=(max-min)/n #interval width
#function used to map a value to the intervals
hist(x,width)=width*floor(x/width)+width/2.0
set boxwidth width*0.9
set style fill solid 0.5 # fill style

#count and plot
plot "data.dat" u (hist($1,width)):(1.0) smooth freq w boxes lc rgb"green" notitle

set style histogram rowstacked gap 0
set style fill solid 0.5 border lt -1
plot "./data.dat" smooth freq with boxes

bin(x,width,n)=x<=n*width? width*(n-1) + 0.5*binwidth:bin(x,width,n+1)
binwidth = 0.05
set boxwidth binwidth
plot "data.dat" u (bin($1,binwidth,1)):(1.0) smooth freq with boxes

rint(x)=(x-int(x)>0.9999)?int(x)+1:int(x)