找不到带有GWmodel对象的闪亮应用

时间:2020-09-05 12:57:57

标签: r shiny gwmodel

我是R的新手,我想基于一个名为GWModel的R包构建自己的Shiny应用程序。这是我的代码

library(shiny)
library(GWmodel)
library(sp)

ui<-shinyUI(fluidPage(
    titlePanel("Model"),
    
    
    sidebarLayout(
        sidebarPanel(
            textInput(inputId = "tanah",
                      label = "Luas Tanah",
                      value = 72),
            textInput(inputId = "bangunan",
                      label = "Luas Bangunan",
                      value = 36),
            textInput(inputId = "lat",
                      label = "Latitude",
                      value = -6.401150),
            textInput(inputId = "long",
                      label = "Longitude",
                      value = 106.770416),
            actionButton("Run_model", "Run model")
        ),
        
        mainPanel(
            
            tabsetPanel(
                tabPanel("model summary", tableOutput('summary'), verbatimTextOutput('summary2'))
            )
        )
    )))


server<- function(input,output,session){
    
    #Train data
    train.data <- read.csv("depoklagi.csv")
    
    
    
    #Train Data SPDF
    train.data.spdf <- SpatialPointsDataFrame(coords = train.data[,c("longlamudi","latlamudi")],
                                              data = train.data[,c("hrglamudi","jdllamudi","ltxlamudi","lbxlamudi","brlamudi")],
                                              proj4string = CRS("+proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0"))
    
    #Train Data DM
    train.data.DM <- gw.dist(dp.locat = coordinates(train.data.spdf))
    
    #Train Data BW (given)
    train.data.bw <-  bw.gwr(hrglamudi ~ ltxlamudi + lbxlamudi,
                             data = train.data.spdf,
                             approach = "CV",
                             kernel = "bisquare",
                             adaptive = T,
                             dMat = train.data.DM)
    
    ####################################################################################################################
    
    #Test data compile
    test.data <- reactive({
        X <- as.numeric(1)
        hrglamudi <- as.numeric(0)
        jdllamudi <- as.character("Prediction")
        ltxlamudi <- as.numeric(as.character(input$tanah))
        lbxlamudi <- as.numeric(as.character(input$bangunan))
        brlamudi <- as.numeric(2)
        latlamudi <- as.numeric(as.character(input$lat))
        longlamudi <- as.numeric(as.character(input$long))
        
        test.data <- cbind(X,hrglamudi,jdllamudi,ltxlamudi,lbxlamudi,brlamudi,latlamudi,longlamudi)
        test.data <- as.data.frame(test.data)
        
        test.data$X <- as.numeric(test.data$X)
        test.data$hrglamudi <- as.numeric(test.data$hrglamudi)
        test.data$jdllamudi <- as.character(test.data$jdllamudi)
        test.data$ltxlamudi <- as.numeric(as.character(test.data$ltxlamudi))
        test.data$lbxlamudi <- as.numeric(as.character(test.data$lbxlamudi))
        test.data$brlamudi <- as.numeric(test.data$brlamudi)
        test.data$latlamudi <- as.numeric(as.character(test.data$latlamudi))
        test.data$longlamudi <- as.numeric(as.character(test.data$longlamudi))
        
        test.data
    })
    
    #Test data SPDF
    test.data.spdf <- reactive({SpatialPointsDataFrame(coords = test.data()[,c("longlamudi","latlamudi")],
                                                       data = test.data()[,c("hrglamudi","jdllamudi","ltxlamudi","lbxlamudi","brlamudi")],
                                                       proj4string = CRS("+proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0"))
        
        test.data.spdf
    })
    
    #Test data DM
    test.data.DM <- reactive({gw.dist(dp.locat = coordinates(train.data.spdf),
                                      rp.locat = coordinates(test.data.spdf()))
        test.data.DM
    })
    
    #################################################################################################################### 
    
    
    #GWModel
    model <- eventReactive(input$Run_model,{gwr.predict(hrglamudi ~ ltxlamudi + lbxlamudi + brlamudi,
                                   data = train.data.spdf,
                                   predictdata = test.data.spdf(),
                                   bw = train.data.bw,
                                   kernel = "bisquare",
                                   adaptive = T,
                                   dMat1 = test.data.DM(),
                                   dMat2 = train.data.DM)
        model
    })
    
    output$summary <- renderTable(test.data())
    output$summary2 <- renderText(model()$SDF$prediction)
}

shinyApp(ui=ui, server=server)

当我运行该应用程序时,它不断返回“找不到对象'pd.locat'”,我已经读取了gwr.predict HereHere后面的代码,用过

我希望成为输出的是模型的预测数字。通过阅读本文和解决方案,我很感谢您的时间。谢谢

1 个答案:

答案 0 :(得分:0)

问题出在您的反应性元素上。在制作test.data的对象中,您是在反应堆中创建一个名为test.data的对象,然后对其进行修改,因此必须最后通过调用test.data返回该对象。作为函数的最后一个元素。在创建test.data.spdftest.data.DMmodel的反应堆中,情况并非如此。在上述每个函数中,您要调用一个产生输出的函数,然后让该函数产生直接返回输出的函数。您没有将其保存到随后需要返回的对象中。每个反应式的最后一行是对不存在的对象的调用。如果您在返回对象的每个反应式中都删除了行(即上面提到的三个中的每行的最后一行),则该应用程序应该可以工作。

library(shiny)
library(GWmodel)
library(sp)

ui<-shinyUI(fluidPage(
  titlePanel("Model"),
  
  
  sidebarLayout(
    sidebarPanel(
      actionButton("browser", "browser"),
#      tags$script("$('#browser').hide();"),
      
      textInput(inputId = "tanah",
                label = "Luas Tanah",
                value = 72),
      textInput(inputId = "bangunan",
                label = "Luas Bangunan",
                value = 36),
      textInput(inputId = "lat",
                label = "Latitude",
                value = -6.401150),
      textInput(inputId = "long",
                label = "Longitude",
                value = 106.770416),
      actionButton("Run_model", "Run model")
    ),
    
    mainPanel(
      
      tabsetPanel(
        tabPanel("model summary", tableOutput('summary'), verbatimTextOutput('summary2'))
      )
    )
  )))


server<- function(input,output,session){
  observeEvent(input$browser,{
      browser()
  })
  
  #Train data
  train.data <- read.csv("depoklagi.csv")
  
  
  
  #Train Data SPDF
  train.data.spdf <- SpatialPointsDataFrame(coords = train.data[,c("longlamudi","latlamudi")],
                                            data = train.data[,c("hrglamudi","jdllamudi","ltxlamudi","lbxlamudi","brlamudi")],
                                            proj4string = CRS("+proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0"))
  
  #Train Data DM
  train.data.DM <- gw.dist(dp.locat = coordinates(train.data.spdf))
  
  #Train Data BW (given)
  train.data.bw <-  bw.gwr(hrglamudi ~ ltxlamudi + lbxlamudi,
                           data = train.data.spdf,
                           approach = "CV",
                           kernel = "bisquare",
                           adaptive = T,
                           dMat = train.data.DM)
  
  ####################################################################################################################
  
  #Test data compile
  test.data <- reactive({
    X <- as.numeric(1)
    hrglamudi <- as.numeric(0)
    jdllamudi <- as.character("Prediction")
    ltxlamudi <- as.numeric(as.character(input$tanah))
    lbxlamudi <- as.numeric(as.character(input$bangunan))
    brlamudi <- as.numeric(2)
    latlamudi <- as.numeric(as.character(input$lat))
    longlamudi <- as.numeric(as.character(input$long))
    
    test.data <- cbind(X,hrglamudi,jdllamudi,ltxlamudi,lbxlamudi,brlamudi,latlamudi,longlamudi)
    test.data <- as.data.frame(test.data)
    
    test.data$X <- as.numeric(test.data$X)
    test.data$hrglamudi <- as.numeric(test.data$hrglamudi)
    test.data$jdllamudi <- as.character(test.data$jdllamudi)
    test.data$ltxlamudi <- as.numeric(as.character(test.data$ltxlamudi))
    test.data$lbxlamudi <- as.numeric(as.character(test.data$lbxlamudi))
    test.data$brlamudi <- as.numeric(test.data$brlamudi)
    test.data$latlamudi <- as.numeric(as.character(test.data$latlamudi))
    test.data$longlamudi <- as.numeric(as.character(test.data$longlamudi))
    
    test.data
  })
  
  #Test data SPDF
  test.data.spdf <- reactive({SpatialPointsDataFrame(coords = test.data()[,c("longlamudi","latlamudi")],
                                                     data = test.data()[,c("hrglamudi","jdllamudi","ltxlamudi","lbxlamudi","brlamudi")],
                                                     proj4string = CRS("+proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0"))
    
  })
  
  #Test data DM
  test.data.DM <- reactive({gw.dist(dp.locat = coordinates(train.data.spdf),
                                    rp.locat = coordinates(test.data.spdf()))
  })
  
  #################################################################################################################### 
  
  
  #GWModel
  model <- eventReactive(input$Run_model,{gwr.predict(hrglamudi ~ ltxlamudi + lbxlamudi + brlamudi,
                                                      data = train.data.spdf,
                                                      predictdata = test.data.spdf(),
                                                      bw = train.data.bw,
                                                      kernel = "bisquare",
                                                      adaptive = T,
                                                      dMat1 = test.data.DM(),
                                                      dMat2 = train.data.DM)
  })
  
  output$summary <- renderTable(test.data())
  output$summary2 <- renderText(model()$SDF$prediction)
}

shinyApp(ui=ui, server=server)