使用如何计算趋势AWK？(How compute trend using awk?)

我有一个时间序列温度数据为：

ifile.txt
1921  25
1922  25.1
1923  24.2
1924  23.4
1925  24.4
1926  25.1
1927  23.6
1928  25.2
1929  23.9
1930  25.6

我想计算的趋势期间1923年至1929年。

我的算法来做到这一点是：

适合上述各点的线如

Y = mX + c ; where m is the trend

我能找到的整个时期的趋势，但我能不能选择一个特定的时期，要修改我的脚本。 请你检查我的脚本？

awk '/1923/,/1929/ {if (NR==FNR) x[NR] = $1;y[NR] = $2;
      sx += x[NR]; sy += y[NR];
      sxx += x[NR]*x[NR];
      sxy += x[NR]*y[NR];}
      END{
      det = NR*sxx - sx*sx;
      trend = (NR*sxy - sx*sy)/det;
      print trend;}' ifile.txt

它不打印正确的值。

 The correct value is 0.0679

Answer 1:

你可以把它简化一点，你就不算匹配的记录。 NR将在END块行的总数。

awk '/1923/,/1929/ {sx+=$1; sy+=$2; c++;
                    sxx+=$1*$1; sxy+=$1*$2}
           END     {det=c*sxx-sx*sx;
                    print det?(c*sxy-sx*sy)/det:"DIV0"}' file

还需要考虑分母为零。 这一趋势出现的

0.0678571

文章来源: How compute trend using awk?