|
别人回复了我的帖子,我发现 hourly std这么估算误差很大, 理论检验的话会发现展开后会多一些项数, 例子检验如下
session 1: 1, 20, -5
session 2: 3, 8, 19, 1
l = [1,20, -5, 3, 8, 19, 1]
正确的std为
In [4]: np.std([1,20, -5, 3, 8, 19, 1])
Out[4]: 8.8271381175859478
按老霍的公式呢:
In [20]: u = np.mean(l)
In [21]: l[3:]
Out[21]: [3, 8, 19, 1]
In [22]: np.sqrt((1/3.0*(np.sum(l[:3]) - 3*u)**2 + 1/4.0*(np.sum(l[3:]) - 4*u)**2)/2)
Out[22]: 2.2373985744503324
|
|