import pandas as pd
import statsmodels.api as sm
import numpy as np# 其他包
#创建数据
t=300
data=pd.DataFrame({'Y1':np.random.normal(0,1,t),
'X1':-3*np.random.normal(0,1,t),
'X2':5*np.random.normal(0,1,t)},
index=pd.date_range(start='20200312',periods=t))
#拟合模型
lag_AR=3#AR滞后阶数
lag_MA=0#MA滞后阶数
model = sm.tsa.VARMAX(endog=data, order=(lag_AR, lag_MA)).fit(disp=False)#估计模型
attribute=dir(model) # 键入dir(model)以获得完整的属性列表
print(model.predict.__doc__)#查看帮助文档
model.summary()#统计信息汇总
#系数检验
params=model.params#系数
tvalues=model.tvalues#系数t值
bse=model.bse#系数标准误
pvalues=model.pvalues#系数p值
#模型检验
resid=model.resid#残差序列
model.aic#模型的AIC值
model.bic#模型的BIC值
model.df_model#模型自由度,自变量个数k
model.df_resid#残差自由度,n-k-1,样本个数-自变量个数-1
#预测
model.predict(exog=data[['X1','X2']])
python的statamodels模块的VAR模型实例
最新推荐文章于 2023-10-20 20:16:12 发布