有人可以帮我了解一下情况吗regr_slope
数据集中有行时返回 NULL?例如:
log=> select * from sb1 order by id, ts;
id | elapsed | ts
------+---------+----------------
317e | 86 | 1552861322.627
317e | 58 | 1552861324.747
317e | 52 | 1552861325.722
317e | 58 | 1552861326.647
317e | 82 | 1552861327.609
317e | 118 | 1552861328.514
317e | 58 | 1552861329.336
317e | 58 | 1552861330.317
317e | 54 | 1552861330.935
3441 | 68 | 1552861324.765
3441 | 84 | 1552861326.665
3441 | 56 | 1552861327.627
3441 | 50 | 1552861330.952
5fe6 | 42 | 1552993248.398
5fe6 | 44 | 1552993255.883
5fe6 | 44 | 1553166049.261
c742 | 62 | 1552861322.149
c742 | 68 | 1552861322.455
(18 rows)
log=> select id, regr_slope(elapsed, ts) as trend from sb1 group by id;
id | trend
------+----------------------
c742 |
317e |
5fe6 | 5.78750952760444e-06
3441 |
(4 rows)
有趣的是,相同的数据集和函数在 Oracle 11.2 中返回不同的结果:
SQL> select * from sb1 order by id, ts;
ID ELAPSED TS
---------- ---------- ----------------
317e 86 1552861322.627
317e 58 1552861324.747
317e 52 1552861325.722
317e 58 1552861326.647
317e 82 1552861327.609
317e 118 1552861328.514
317e 58 1552861329.336
317e 58 1552861330.317
317e 54 1552861330.935
3441 68 1552861324.765
3441 84 1552861326.665
3441 56 1552861327.627
3441 50 1552861330.952
5fe6 42 1552993248.398
5fe6 44 1552993255.883
5fe6 44 1553166049.261
c742 62 1552861322.149
c742 68 1552861322.455
18 rows selected.
SQL> select id, regr_slope(elapsed, ts) from sb1 group by id;
ID REGR_SLOPE(ELAPSED,TS)
---------- ----------------------
c742 19.6078431
5fe6 5.7875E-06
317e -1.0838511
3441 -3.8283951
我不知道这是否意味着 Postgres、Oracle 有问题,两者都没有问题,还是两者都有问题,尽管结果5fe6
是相同的。
深入研究代码后,我得到了答案:
问题在于,在这种情况下,PostgreSQL v12 之前的简单方法会导致不必要的大舍入误差。
让我们考虑一下id = 'c742'
:
公式为regr_slope
是这样的:
regr_slope
:= (N ⋅ Σ(Xi⋅Yi) - ΣXi ⋅ ΣYi) / (N ⋅ Σ(Xi2) - ΣXi ⋅ ΣXi)
问题出在除数上:
SELECT 2::float8 * (1552861322.149::float8 * 1552861322.149::float8 +
1552861322.455::float8 * 1552861322.455::float8) -
(1552861322.149::float8 + 1552861322.455::float8) *
(1552861322.149::float8 + 1552861322.455::float8);
?column?
----------
-2048
(1 row)
由于结果为负,PostgreSQL 返回 NULL 结果。
使用精确计算(使用numeric
):
SELECT 2 * (1552861322.149 * 1552861322.149 +
1552861322.455 * 1552861322.455) -
(1552861322.149 + 1552861322.455) *
(1552861322.149 + 1552861322.455);
?column?
----------
0.093636
(1 row)
自从 PostgreSQL 提交以来e954a727f0c8872bf5203186ad0f5312f6183746 https://git.postgresql.org/gitweb/?p=postgresql.git;a=commitdiff;h=e954a727f0c8872bf5203186ad0f5312f6183746,事情得到了改善,在 PostgreSQL v12 中,PostgreSQL 也返回了正确的结果:
select id, regr_slope(elapsed, ts) from sb1 group by id;
id | regr_slope
------+-----------------------
c742 | 19.607858781290517
317e | -1.0838511987808963
5fe6 | 5.787509483586743e-06
3441 | -3.828395463097356
(4 rows)
本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系:hwhale#tublm.com(使用前将#替换为@)