可以使用R语言中的irf包来绘制pvar模型的方差分解图
library(vars)
library(irf)
# 导入数据
data("Canada")
# 建立pvar模型
model <- VAR(Canada, p = 2, type = "const")
# 进行方差分解
fevd <- fevd(model, n.ahead = 4)
# 绘制方差分解图
plot(fevd, main = "Variance Decomposition", xlab = "Horizon", ylab = "Percent")
代码使用了vars包中的VAR函数建立pvar模型,并使用fevd函数进行方差分解,然后使用plot函数绘制方差分解图。您可以根据需要修改图表的标题、横轴和纵轴标签等。