为了找出影响价格波动的主要因素，我们使用逐步回归法来剔除一些对于应变量即价格影响很小的自变量剔除出我们的模型我们分别把WTI Price Field Production of Crude Oil (Thousand Barrels) 等自变量的名称改为x1，x2……最后的突发事件需要用到哑变量，哑变量只需要2个即可，我们将其作为X49，X50，X51,三个参数并将它们的值”正影响”，”无影响”，”负影响”分别改为-1,0,1.。
我们分别把WTI Price Field Production of Crude Oil (Thousand Barrels) 等自变量的名称改为x1，x2……最后的突发事件需要用到哑变量，哑变量只需要2个即可，我们将其作为X49，X50，X51,三个参数并将它们的值”正影响”，”无影响”，”负影响”分别改为-1,0,1.。
经过R语言处理以后我们得到模型
Y~x1 + x2 + x4 + x5 + x7 + x13 + x14 + x15 + x16 + x17 + x18 + x20 + x21 + x23 + x34 + x25 + x26 + x29 + x30 + x33 + x35 + x36 + x37 + x39 + x40 + x42 + x44 + x46 + x47 + x48 + x49 + x50
由此可得，影响较小的已经被剔除了。
Garch模型预测波动性
我们通过使用Garch模型来预测波动性，
先检验数据正态性，可以计算得出数据分布函数，QQ图，对数收益率序列折线图
> shapiro.test(rlogdiffdata)
Shapiro-Wilk normality test
data: rlogdiffdata
W = 0.94315, p-value = 1.458e-05由QQ图以及p值可见，数据大致上符合正态分布。
最后用VaR曲线来进行预警剧烈波动。
红色的点就是剧烈波动前的预警点。
强影响点分析
我们可以通过使用cook统计量来寻找强影响点，因此我们用R语言的influence.measures()函数来进行影响分析。
右侧带有*号的表示强影响点。
我们通过学生化残差来构造F检验，最终得到t检验，以此来检测异常点。通过
stdres<-rstudent(lm.sol)来得到学生化残差，然后通过公式
来计算Fj，并且最终转化为tj，
t=sqrt((144-51-1)*stdres^2/(144-51-stdres^2))
最后我们可以检查，如果
则它为异常点。
R语言执行
res<-t>abs(qt(.025, df=92))
可以直接得到大于对应t值的布尔值。
值为True的则可能为异常点。
预测
我们使用了HoltWinters来进行预测我们的油价区间
真实值基本都在预测的范围内，但是想要净赚预测还是比较困难的。
标签：变量,模型,Garch,x2,股票价格,影响,我们 来源： https://blog.csdn.net/qq_19600291/article/details/110931539
本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享； 2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关； 3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关； 4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除； 5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。