Question

我正在使用闪亮的应用来制作一些预测图。但是，我有一个问题，即应用程序只显示一些图表。我目前正在使用以下代码作为应用程序的服务器部分。

server <- function(input, output) {  

 output$graph <- renderPlot({
      row.number <- which(grepl(input$n1, Data$Kontotal))
      name <- Data [row.number, -(1:3)]
      name <- t(name)
      name <- ts(name, frequency=12, start=c(2007,1))
      n.row <- nrow(name)

      en<-max(time(name))
      ds<-as.data.frame(window(name,end=en))
      names(ds)<-'obs'
      ds$date<-as.Date(time(window(name,end=en)))     


       p1a <- ggplot(data=ds,aes(x=date,y=obs)) 
       p1a <- p1a+geom_line(col='red')
       p1a <- p1a+scale_x_date(name='',breaks='1 year',minor_breaks='1 month',labels=date_format("%b-%y"),expand=c(0,0))
       p1a <- p1a+scale_y_continuous(name='mio.kr')
       p1a <- p1a + theme_bw() + 
            theme(panel.border = element_blank(), 
                  panel.grid.major = element_blank(),
                  panel.grid.minor = element_blank(), 
                  axis.line = element_line(colour = "black")) })


 output$graph2 <- renderPlot({

  if (input$n3 == "HW") { 
           row.number <- which(grepl(input$n1, Data$Kontotal))
           name <- Data [row.number, -(1:3)]
           name <- t(name)
           name <- ts(name, frequency=12, start=c(2007,1))
           n.row <- nrow(name)


           fit <- auto.arima (name)
           fcast <- forecast(fit, input$n2)



           ### Plotting

           data.pred <- forc.ggplot (name, fcast)




           p1a <- ggplot(data=data.pred,aes(x=date,y=obs)) 
           p1a <- p1a+geom_line(col='red')
           p1a <- p1a+geom_line(aes(y=fitted),col='blue')
           p1a <- p1a+geom_line(aes(y=forecast))+geom_ribbon(aes(ymin=lo95,ymax=hi95),alpha=.25)
           p1a <- p1a+scale_x_date(name='',breaks='1 year',minor_breaks='1 month',labels=date_format("%b-%y"),expand=c(0,0))
           p1a <- p1a+scale_y_continuous(name='mio.kr')
           p1a <- p1a + theme_bw() + 
                theme(panel.border = element_blank(), 
                      panel.grid.major = element_blank(),
                      panel.grid.minor = element_blank(), 
                      axis.line = element_line(colour = "black"))
           p1a


      }  

  if (input$n3 == "arima") { 
        row.number <- which(grepl(input$n1, Data$Kontotal))
        name <- Data [row.number, -(1:3)]
        name <- t(name)
        name <- ts(name, frequency=12, start=c(2007,1))
        n.row <- nrow(name)


        fit <- auto.arima (name)
        fcast <- forecast(fit, input$n2)



        ### Plotting

        data.pred <- forc.ggplot (name, fcast)




        p1a <- ggplot(data=data.pred,aes(x=date,y=obs)) 
        p1a <- p1a+geom_line(col='red')
        p1a <- p1a+geom_line(aes(y=fitted),col='blue')
        p1a <- p1a+geom_line(aes(y=forecast))+geom_ribbon(aes(ymin=lo95,ymax=hi95),alpha=.25)
        p1a <- p1a+scale_x_date(name='',breaks='1 year',minor_breaks='1 month',labels=date_format("%b-%y"),expand=c(0,0))
        p1a <- p1a+scale_y_continuous(name='mio.kr')
        p1a <- p1a + theme_bw() + 
             theme(panel.border = element_blank(), 
                   panel.grid.major = element_blank(),
                   panel.grid.minor = element_blank(), 
                   axis.line = element_line(colour = "black"))
        p1a


   }  


 })

}

问题是只显示ARIMA图表。注意，但是，在这里发布的代码中我更改了代码的HW部分中的模型，因此它与ARIMA部分（包括模型）相同，只是为了检查它是因为ets函数出了问题（我用它估计了HW模型）。

它只是留下了一个空白区域，其中图形应该是HW的情况，但是绘制了ARIMA模型的右图。

有人试过类似的东西吗？

Answer 1

将p1a放在renderPlot()块之外的if块的末尾。 renderPlot尝试绘制最后一个表达式的块。如果是input$n3 =="HW"，则生成绘图对象p1a，然后继续评估if的下一个FALSE条件，因此不会将任何内容传递给plotOutput()。

另一种解决方案是使用switch()代替if。 switch()仅计算匹配元素的表达式。

闪亮的应用程序，只显示一些图表

1 个答案: