我估计R中的SVAR,但A-B表格结果与Eviews非常不同,我不确定为什么会这样。另外,我认为正确的选项给了我错误信息。谁能帮助我? 这是我使用的R代码:
resA <- matrix(NA, nrow = 5, ncol = 5)
resA[2,4]=resA[2,5]=resA[3,4]=resA[3,5]=resA[4,2]=resA[4,3]=resA[4,5]=0
resA[5,2]=resA[5,3]=resA[5,4]=0
resA[1,1]=resA[2,2]=resA[3,3]=resA[4,4]=resA[5,5]=1
resA
model=VAR(vardata, p=2, type="const")
summary(model)
stt=matrix(0.1, nrow = 1, ncol = 10)
model1=SVAR(model, Amat=resA, lrtest=TRUE, estmethod="scoring", start=stt, conv.crit=0.0001, max.iter=500)
summary(model1)
irf.gap=irf(model1, impulse="gap", boot=FALSE, n.ahead=15, runs=100)
plot(irf.gap)
问题是最后一个命令IRF。它给我的反转形状比Eviews。我认为,因为Eviews只提到它正在使用Cholesky Decomposition和df调整(这个东西应该与CI相关)和#34;对Cholesky的回应S.D.创新+/- 2 S.E。&#34;,我想问题应该来自一个SD和2SE,仍然不确定R命令&#34; irf&#34;确实...
顺便说一下,R I&m; m使用的包是library(vars)
,而对于Eviews,我使用了IRF的默认设置。
更新: 问题发生是因为命令irf计算的结构脉冲响应函数不同于Eviews&#39; Cholesky分解。 任何人都会分享任何链接以及手动计算Eview版本的IRF的步骤非常感谢!