Gnuplot(或matplotlib):在x轴上创建带有类别的非条形图

时间:2011-12-03 23:47:46

标签: python graph matplotlib gnuplot

我想使用gnuplot(或matplotlib,如果需要)创建一个这样的图形,但我不知道是否/如何完成:

rough outline of what the graph should look like

这当然只是一个粗略的草图。重要的是我需要绘制成对的值(在这个例子中,每对由红点和蓝点组成)。每对中的一个项目是单个值,另一个项目应该显示一系列值(我的想法是用误差条绘制平均值以指示范围的最大值和最小值,但我愿意接受更好的想法) 。除了给出各种类别的名称之外,x轴没有其他目的 - 所有重要的都是y值。

我很确定我可以使用直方图创建这样的东西(值对和x类别),但在这种情况下,盒子只会让我错误。

到目前为止我得到了什么:我有这个gnuplot命令:

plot 'TEST.out' using 0:2:3:xticlabel(1) w errorbars pt 7 notitle

与此数据文件一起使用(类别名称,y值,错误栏值):

cat1 15 0
cat1 18 3
cat2 13 0
cat2 10 4

它会产生下面的图,它朝向正确的方向,但是它还不理想(所有数据点都具有相同的颜色,对于单个值,您仍然可以看到使用了误差条;分组也不是'非常好 - 如果构成一对的两个点彼此更接近,这将使得情节更容易在眼睛上看。)

where I have got so far

如果有人有任何建议(即使创建的图表看起来与我在开头时给出的示例完全不同)。我会很感激。

非常感谢你的时间!

2 个答案:

答案 0 :(得分:5)

在Matplotlib中,您可以单独设置刻度。该方法在one of Matplotlib's examples中解释。

这是一个开始:

# Initializations:
from matplotlib import pyplot as pp
import numpy as np

pp.ion()  # "Interactive mode on": pyplot.* commands draw immediately

# Data:
series_red_y = [1.3, 1.4, 2.2]
series_blue_y = [1.6, 1.8, 1.8]
series_blue_err = [0.25, 0.25, 0.5]
names = ('Category 1', 'Category 2', 'Category 3')

# Where on the x-axis the data will go:    
series_red_x = np.arange(0, 3*len(series_red_y), 3)  # Step of 3: red dot, blue dot, empty space
series_blue_x = np.arange(1, 3*len(series_blue_y)+1, 3)  # Step of 3: red dot, blue dot, empty space

# Plotting:
pp.scatter(series_red_x, series_red_y, c='r', s=100)
pp.scatter(series_blue_x, series_blue_y, s=100)
pp.errorbar(series_blue_x, series_blue_y, yerr=series_blue_err, fmt=None,
            capsize=0)
pp.xticks((series_red_x+series_blue_x)/2., names)

# We wait until the user is ready to close the program:
raw_input('Press enter...')

以下是结果,您可以根据自己的特定需求进行自定义:

enter image description here

答案 1 :(得分:1)

这是一个gnuplot解决方案:

set xrange [-.25:2.75]
set xtics ('cat0' .25, 'cat1' 1.25, 'cat2' 2.25)
plot '< sed 1~2!d TEST.out' using 0:2:(.1) with circles, \
     '< sed 1~2d TEST.out' using (.5+$0):2:3 with errorbars

如果事先不知道类别的数量,您可能需要生成脚本的前两行。它应该很容易,例如在bash:

cat_count=$(wc -l < TEST.out)
let cat_count/=2
let cat_count--
echo xrange "[-.25:$cat_count.75]"
echo -n 'set xtics('
for i in $(seq 0 $cat_count) ; do
    echo -n "'cat$i'" $i.25
    ((i==cat_count)) || echo -n ,
done
echo ')'