生成脉冲响应置信区间。 Rs中的Tsdyn包

Question

我想使用R中的tsDyn包生成响应脉冲的置信区间。我尝试使用以下代码，但我认为这是不对的，如果有人提出任何建议，我将不胜感激。谢谢。

library(tsDyn)
library(boot)

set.seed(100)
data <- data.frame(var1 = rnorm(50), var2 = rnorm(50), rnorm(50))

fun_tvar <- function(data, sd =1){
  tvar <- TVAR(data, lag =2)
  irfs <- GIRF(tvar,shk = sd, replic = 10)
  list(ir_down = irfs[,3,1], ir_up = irfs[,3,2])
}

fc <- function(d, i){
  d2 <- d[i,]
  return(fun_tvar(d2))
}

boot1 <- boot(data3, fc, R=50)
boot1

有关GIRF函数here的信息。