plot - gnuplot 计算多列统计数据

Question

我有一个简单的 9 列文件。我不想为每一列计算某些统计数据然后绘制它(使用 gnuplot)。

1) 这就是我计算除第一列之外的每一列的统计信息的方式。

stats 'data' every ::2 name "stats"

2) 在输出屏幕中我可以看到操作成功。注意列数/记录数为8

* FILE: 
  Records:      8
  Out of range: 0
  Invalid:      0
  Blank:        0
  Data Blocks:  1

* COLUMNS:
  Mean:          6.5000       491742.6625
  Std Dev:       2.2913          703.4865
  Sum:          52.0000       3.93394e+06
  Sum Sq.:     380.0000       1.93449e+12

  Minimum:       3.0000 [0]   490312.0000 [2]
  Maximum:      10.0000 [7]   492643.5000 [7]
  Quartile:      4.5000       491329.5000
  Median:        6.5000       491911.1500
  Quartile:      8.5000       492252.2500

  Linear Model: y = 121.8 x + 4.91e+05
  Correlation:  r = 0.3966
  Sum xy:       2.558e+07

3) 现在我可以通过像这样附加 _x 和 _y 来访问前两列的统计信息

print stats_median_x
print stats_median_y

我的问题是:

• 如何访问其余 6 列的统计信息(比如中位数)？
• 我怎样才能绘制一条横跨所有中位数的线并与某个 X 轴相对应？

我知道我可以简单地添加一个 python 脚本来预先计算所有这些，但如果有一种使用 gnuplot 本身的简单方法来完成它，我宁愿避免它。

谢谢!

Answer 1

简答

• “如何访问其他列的统计信息？”
  使用 stats 'data'using n，您将访问第 n 列...
• “例如，我如何绘制所有中位数？”
  例如set print 和 do for 循环可以创建可用于绘图的数据文件。

一个可行的解决方案

    set print "StatDat.dat" 
    do for [i=2:9] { # Here you will use i for the column.
      stats  'data.dat' u i nooutput ; 
      print i, STATS_median, STATS_mean , STATS_stddev # ...
    } 
    set print
    plot "StatDat.dat" us 1:2 # or whatever column you want...

再说几句
通过 help stats 向 gnuplot 本身寻求帮助，可以阅读很多有趣的东西 :-)。

Syntax:
stats 'filename' [using N[:M]] [name 'prefix'] [[no]output]]
This command prepares a statistical summary of the data in one or two columns of a file. The using specifier is interpreted in the same way as for plot commands. See plot for details on the index, every, and using directives.

• 从第一个突出显示的句子我们可以了解到它每次准备一列或最多两列的统计信息(很遗憾，让我们以后再看......)。
• 从第二个突出显示的句子可以看出它将遵循与 plot 命令相同的语法:
  所以 stats 'data'using 3 将为您提供 x
  第 3 列的统计信息和 stats 'data' 使用 4:5 在 x,y 中的第 4 和第 5...

关于你的解释的注释

1. 你说

   This is how I compute statistics for every column excluding the first one.
   stats 'data' every ::2 name "stats"

   这实际上不是前两列的统计数据不包括前两行，实际上它们的计数器是从 0 而不是 1 开始的。

2. 作为上述假设/解释的结果，当我们阅读

   Records: 8

   这意味着计算的行是8；您的文件有 10(可用)行，您指定 every::2 并跳过前两行，因此您有 8 条对统计有用的记录。
   的确如此，我们可以更好地理解何时在 help stats 中说

   STATS_records           # total number of in-range data records
   

   暗示“用于计算此统计数据”。

_{在 gnuplot 4.6 补丁级别 4 上测试
致力于 gnuplot 版本 5.0 补丁级别 1}