library(vars)
data(Canada)
var_fit <- VAR(Canada, p = 1)
var_irf <- irf(var_fit, impulse = c("U", "rw"), response = "prod")


如何在图形中并排绘制两个脉冲响应

通常,我会使用par(mfrow = c(1,2)),但是它不能按预期工作。有什么帮助吗?

最佳答案

我发现了同样的问题。我“手动”解决了,这里是一个带有两个变量的VAR(1)模型的示例。

    impulse<-irf(model)
    irf1<-data.frame(impulse$irf$y1[,1],impulse$Lower$y1[,1],
                     impulse$Upper$y1[,1])
    irf2<-data.frame(impulse$irf$y1[,2],impulse$Lower$y1[,2],
                     impulse$Upper$y1[,2])
    par(mfrow=c(1,2), bg="azure2")
    matplot(irf1, type="l", lwd=2, col="blue2",
            ylab=expression(y[1]), lty=c(1,2,2))
    matplot(irf2, type="l", lwd=2, col="red2",
            ylab=expression(y[1]), lty=c(1,2,2))

关于r - 在R包“vars”中绘制并排脉冲响应,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/28870142/

10-09 01:47