如何在gnuplot中标记90百分位数

时间:2013-01-31 06:57:52

标签: gnuplot

如何在gnu情节中标记90百分位数。需要样品

plot [][] "fle.txt" using 1:2 with lines

这是我绘制图形的方式,我想在图表中标记90% 这是我的数据集

time(seconds)  frequency  cumulativeFrequency
10             2          2
11             6          8
12             8          16
13             10         26
14             7          33
15             5          38
16             5          43
17             4          47
18             2          49

2 个答案:

答案 0 :(得分:4)

我无法在gnuplot中看到任何方法,但使用python并不难 -

# percent90.py
import sys
def to_num(iterable):
    for line in iterable:
        columns = line.split()  # split into columns
        if not columns:  # empty line
            continue
        try:
            yield float(columns[1])  # This is a good number -- Yield it and keep going
        except ValueError:  # Not a number on the line
            pass            # just keep going -- do nothing
        except IndexError:
            print line

with open(sys.argv[1]) as fin:
     data = sorted(to_num(fin))

top_10 = data[int(len(data)*0.9):]  #top 10 percent of the data
print(top_10[0])

你可以称之为:

python percent90.py path/to/datafile

这将告诉您放置标记的位置。至于标记它,我可能会在gnuplot中做这样的事情:

YVAL = `python percent90.py path/to/datafile`
set arrow 1 from graph 0,first YVAL to graph 1,first YVAL ls 1 nohead front

答案 1 :(得分:3)

如果您只想计算数据的90%点,可以使用statsplot命令和gnuplot的内部变量来完成,然后绘制一条mgilson建议的行:

#!/usr/bin/env gnuplot

set terminal png
set output 'test.png'

# 'every ::1' skips header
stats 'fle.txt' every ::1 

mark = (STATS_max_y - STATS_min_y)*0.9 + STATS_min_y

set arrow 1 from graph 0,first mark to graph 1,first mark ls 1 nohead front

plot 'fle.txt' every ::1

此脚本生成此输出:

enter image description here