Стиль гистограммы Gnuplot

Question

У меня есть такой файл данных:

PS        105  169 169 169   68  0 0 255    105 169 169 169    0  169 169 169    0   169 169 169    0  169 169 169 
GMQ       43    0  255  0    39  0 0 255    40   34 139  34   105 169 169 169    105 169 169 169    0  169 169 169
GL        81    0  255  0    40  0 0 255    105 169 169 169    0  169 169 169    0   169 169 169    0  169 169 169   
CM_{1}    105  169 169 169   47  0 0 255    105 169 169 169    32 169 169 169    105 169 169 169    0  169 169 169
CM_{2}    105  169 169 169   44  0 0 255    105 169 169 169   105 169 169 169    105 169 169 169    0  169 169 169
P         105  169 169 169   105 0 0 255    105 169 169 169    0  169 169 169    0   169 169 169    0  169 169 169
MH        105  169 169 169   69  0 0 255    105 169 169 169    0  169 169 169    0   169 169 169    0  169 169 169

, в котором в столбцах 2 6 10 14 18 22 представлены значения, которые должны быть нанесены на график, а в остальных - код rgb.

Я хотел бы получить такую ​​графику,

пока с прикрепленным скриптом накладываются полосы:

Что я должен исправить в сценарии?

Это скрипт Gnuplot:

set term pngcairo enhanced size 2560,2048 font "Arial,50"
set output "Hist.png"
set encoding utf8

unset title

set yrange [0.0:150.]
set style data histogram
set style histogram cluster gap 1
set style fill solid
set boxwidth 1.
set xtics format ""
set ylabel "Training Set size [-]"
set grid ytics lw 2 
rgb(r,g,b) = int(r)*65536 + int(g)*256 + int(b)
unset key

plot "Data.out" u ($0):($2):(0.5):(rgb($3,$4,$5))   :xticlabels(1) w boxes lc rgb variable,\
             "" u ($0):($6):(0.5):(rgb($7,$8,$9))   :xticlabels(1) w boxes lc rgb variable,\
             "" u ($0):($10):(0.5):(rgb($11,$12,$13)):xticlabels(1) w boxes lc rgb variable,\
             "" u ($0):($14):(0.5):(rgb($15,$16,$17)):xticlabels(1) w boxes lc rgb variable,\
             "" u ($0):($18):(0.5):(rgb($19,$20,$21)):xticlabels(1) w boxes lc rgb variable

reset 

Спасибо за помощь, которую вы можете оказать мне.

С наилучшими пожеланиями.

Answer 1

Согласно демонстрационной странице gnuplot (http://gnuplot.sourceforge.net/demo/histograms.html) кластерные гистограммы должны быть построены так:

plot "Data.out" u 2:xtic(1), '' u 6, '' u 10, '' u 14, '' u 18

Однако я не нашел способа окрашивать столбцы гистограммы индивидуально в зависимости от данных из файла данных, как вы хотели бы. У меня сложилось впечатление, что ... lc rgb var не работает вместе со стилем гистограммы (но я могу ошибаться). Может быть, кто-то может сказать нам.

Поэтому мое предложение реализовать кластерную гистограмму «вручную» со стилем построения with boxes. Вы определяете последовательность своих столбцов, которую хотите построить ColSequence = "2 6 10 14 18", и еще несколько функций для сокращения команды построения.

Код:

### clustered histograms realized "manually" with boxes
reset session
unset title
unset key

$Data <<EOD
PS        105  169 169 169   68  0 0 255    105 169 169 169    0  169 169 169    0   169 169 169    0  169 169 169 
GMQ       43    0  255  0    39  0 0 255    40   34 139  34   105 169 169 169    105 169 169 169    0  169 169 169
GL        81    0  255  0    40  0 0 255    105 169 169 169    0  169 169 169    0   169 169 169    0  169 169 169   
CM_{1}    105  169 169 169   47  0 0 255    105 169 169 169    32 169 169 169    105 169 169 169    0  169 169 169
CM_{2}    105  169 169 169   44  0 0 255    105 169 169 169   105 169 169 169    105 169 169 169    0  169 169 169
P         105  169 169 169   105 0 0 255    105 169 169 169    0  169 169 169    0   169 169 169    0  169 169 169
MH        105  169 169 169   69  0 0 255    105 169 169 169    0  169 169 169    0   169 169 169    0  169 169 169
EOD

set ylabel "Training Set size [-]"
set yrange [0.0:150.]
set grid ytics lw 2 
set style fill solid 
set boxwidth 1.0
rgb(r,g,b) = int(r)*65536 + int(g)*256 + int(b)

ColSequence = "2 6 10 14 18"
myBoxWidth = 0.8
Gap = 2
ColCount = int(words(ColSequence))
ColNo(n) = int(word(ColSequence,n))
Color(n) = rgb(column(n+1),column(n+2),column(n+3))
BoxOffset(n) = 1./(ColCount+Gap)*n - 0.5
plot for [i=1:ColCount] $Data u ($0+BoxOffset(i)):(column(ColNo(i))): \
    (1./(ColCount+Gap)*myBoxWidth):(Color(ColNo(i))) w boxes lc rgb var, \
     $Data u 0:(0):xtic(1) w boxes
### end of code

Результат: