Модель VAR с экзогенной переменной не работает в Shiny - PullRequest
Новая
       23

Модель VAR с экзогенной переменной не работает в Shiny

0 голосов
/ 25 мая 2018

Давайте рассмотрим модель VAR с экзогенной переменной, чтобы различать два периода.Эта модель прекрасно работает следующим образом: 

library(shiny)
library(vars)

#--- Create Exogenous Variable 'periods'
data(Canada)
canTS <- Canada
periods <- as.matrix(data.frame(period=ifelse(index(canTS)>1996, 1, 0 ) ) )

#--- Fit the Model
fit1 <- VAR(Canada, p = 2, type = "none", exogen=periods)
coef(fit1)[[1]]

#-- Make Prediction
period2 <- as.matrix(data.frame(period = rep(1, 12)) ) # Future Exogen Values = 1 
predict(fit1, n.ahead=12, dumvar=period2)

rm(list=ls(all=TRUE)) # remove objects

Проблема возникает, когда используется приложение Shiny с ползунком ввода (который задает горизонт прогнозирования).Мы можем создать простое приложение Shiny как: 

ui <- fluidPage(
 sidebarLayout(
  sidebarPanel(
   sliderInput(inputId = "nAhead",
    label = "Forecast Period",
    min = 1,
    max = 12,
    value = 6)
  ),
  mainPanel(
   verbatimTextOutput("Model")
  )
 )
)

server <- function(input, output) {
 #--- Create Exogenous Variable 'periods'
 data(Canada)
 canTS <- Canada
 periods <- as.matrix(data.frame(period=ifelse(index(canTS)>1996, 1, 0 ) ) ) 

 #--- Fit the Model
 fit1 <- VAR(Canada, p = 2, type = "none", exogen=periods)

 #-- Make Prediction
 output$Model <- renderPrint({
 periods2 <- as.matrix(data.frame(period = rep(1, input$nAhead)) ) # forecasting window
 predict(fit1, n.ahead=input$nAhead, dumvar=periods2) 
 })
}

shinyApp(ui, server)

Я получаю следующее сообщение об ошибке: «Периоды объекта не найдены».Я не понимаю, почему функция predict не принимает новую матрицу данных, предоставленную dumvar.Есть идеи, как заставить это работать?Спасибо.

1 Ответ

0 голосов
/ 25 мая 2018

После строки 

periods <- as.matrix(data.frame(period=ifelse(index(canTS)>1996, 1, 0 ) ) )

просто добавьте следующую строку 

assign("periods",periods, envir = .GlobalEnv)

Тогда все заработает

...