gnuplot - 如何在 gnuplot 中标记 90 个百分位

Question

如何在 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

Answer 1

我在 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

Answer 2

如果您只想计算数据的 90% 点，可以使用statsorplot命令和 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

此脚本产生以下输出：