大家好,又见面了,我是你们的朋友全栈君。
#构建VAR模型
library(sandwich)
library(strucchange)
library(vars)
data.new<-data.frame(S1,S2)
VARselect(data.new,lag.max=20,type="trend") #选择最优的滞后阶数
var<-VAR(data.new,lag=1,ic="AIC")
summary(var)
coef(var)
plot(var)
sta=stability(var,type=c("OLS-CUSUM"),h=0.15,dynamic=FALSE,rescale=TRUE) #平稳性
plot(sta)
summary(sta)
#脉冲响应分析
var.irf=irf(var)
plot(var.irf)
serial.test(var,lags.pt =10,type="PT.asymptotic")
发布者:全栈程序员栈长,转载请注明出处:https://javaforall.cn/167150.html原文链接:https://javaforall.cn