Gnuplot(或matplotlib):创建带有分类的非柱状图

4 投票
2 回答
2698 浏览
提问于 2025-04-17 07:38

我想用gnuplot(如果需要的话也可以用matplotlib)来创建一个像这样的图表,但我不知道是否能做到或者该怎么做:

图表大致样子

当然,这只是一个粗略的草图。重要的是我需要绘制一对对的数值(在这个例子中,每对由一个红点和一个蓝点组成)。每对中的一个是单个值,另一个是表示一个范围的值(我想的是用平均值和误差条来表示范围的最大值和最小值,但我也乐意听听更好的主意)。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

它生成了下面的图,虽然方向是对的,但还不够理想(所有数据点的颜色都一样,而且对于单个值,你仍然能看到使用了误差条;另外,分组也不太好——如果组成一对的两个点能更靠近一点,图看起来会更舒服)。

我目前的进展

如果有人有任何建议(即使是创建一个看起来和我一开始给出的例子不完全一样的图表),我会非常感激。

非常感谢你的时间!

2 个回答

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 ')'
5

在Matplotlib中,你可以单独设置刻度。这个方法在Matplotlib的一个示例中有详细说明。

下面是一个开始的例子:

# 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...')

这是结果,你可以根据自己的需求进行定制:

在这里输入图片描述

撰写回答