foc*_*foc 17 python time-series forecasting pandas statsmodels
我正在研究python中的时间序列.我觉得有用和有前途的图书馆是
也用于可视化:matplotlib
有没有人知道指数平滑的库?
不知何故,有些问题被合并或删除了,所以我会在这里发布我的答案.
在本地进行Python平滑处理.
'''
simple exponential smoothing
go back to last N values
y_t = a * y_t + a * (1-a)^1 * y_t-1 + a * (1-a)^2 * y_t-2 + ... + a*(1-a)^n * y_t-n
'''
from random import random,randint
def gen_weights(a,N):
ws = list()
for i in range(N):
w = a * ((1-a)**i)
ws.append(w)
return ws
def weighted(data,ws):
wt = list()
for i,x in enumerate(data):
wt.append(x*ws[i])
return wt
N = 10
a = 0.5
ws = gen_weights(a,N)
data = [randint(0,100) for r in xrange(N)]
weighted_data = weighted(data,ws)
print 'data: ',data
print 'weights: ',ws
print 'weighted data: ',weighted_data
print 'weighted avg: ',sum(weighted_data)
Run Code Online (Sandbox Code Playgroud)
您可以使用熊猫指数加权移动平均预测未来值http://pandas.pydata.org/pandas-docs/stable/generated/pandas.stats.moments.ewma.html 为
from pandas.stats.moments import ewma
import numpy as np
pred_period = 12
def predict(x,span,periods = pred_period):
x_predict = np.zeros((span+periods,))
x_predict[:span] = x[-span:]
pred = ewma(x_predict,span)[span:]
return pred
Run Code Online (Sandbox Code Playgroud)
| 归档时间: |
|
| 查看次数: |
26658 次 |
| 最近记录: |