statsmodel做时间序列分析,pacf结果不正确

yxnwh 2018-07-20 02:02:02
编译环境:
python3.6.6 win64

statsmodels库文件:
statsmodels-0.9.0-cp36-cp36m-win_amd64

运行时总是提示:
.\python36\lib\site-packages\statsmodels\regression\linear_model.py:1283: RuntimeWarning: invalid value encountered in sqrt
return rho, np.sqrt(sigmasq)


输出的pacf图形也与正确的相距甚远,代码如下,真确和错误的图如下

import pandas as pd
import numpy as np
import matplotlib.pyplot as plt
from statsmodels.graphics.tsaplots import plot_acf,plot_pacf

dta=[10930,10318,10595,10972,7706,6756,9092,10551,9722,10913,11151,8186,6422,
6337,11649,11652,10310,12043,7937,6476,9662,9570,9981,9331,9449,6773,6304,9355,
10477,10148,10395,11261,8713,7299,10424,10795,11069,11602,11427,9095,7707,10767,
12136,12812,12006,12528,10329,7818,11719,11683,12603,11495,13670,11337,10232,
13261,13230,15535,16837,19598,14823,11622,19391,18177,19994,14723,15694,13248,
9543,12872,13101,15053,12619,13749,10228,9725,14729,12518,14564,15085,14722,
11999,9390,13481,14795,15845,15271,14686,11054,10395]
dta_list = np.arange(2001,2091)
data=pd.DataFrame(dta,index=dta_list)
plt.figure
data.plot()
plt.show()

D_data=data.diff(1).dropna()
D_data.plot()
plt.show()

plot_acf(D_data).show()
plot_pacf(D_data).show()


...全文
533 3 打赏 收藏 转发到动态 举报
写回复
用AI写文章
3 条回复
切换为时间正序
请发表友善的回复…
发表回复
ZX1aoShua1 2020-01-27
  • 打赏
  • 举报
回复
请问一下 这个问题怎么解决的啊??迫切知道呢
蓝色_ 2018-11-07
  • 打赏
  • 举报
回复 7
兄弟我解决了,在函数plot_pacf中添加计算方法method='ywm'
蓝色_ 2018-11-06
  • 打赏
  • 举报
回复
请问楼主怎么解决的。我也遇到了相同的问题。

37,719

社区成员

发帖
与我相关
我的任务
社区描述
JavaScript,VBScript,AngleScript,ActionScript,Shell,Perl,Ruby,Lua,Tcl,Scala,MaxScript 等脚本语言交流。
社区管理员
  • 脚本语言(Perl/Python)社区
  • IT.BOB
加入社区
  • 近7日
  • 近30日
  • 至今

试试用AI创作助手写篇文章吧