Shiny - FILL值未在闪亮服务器中正确传递给ggplot - 未找到错误对象'input'

时间:2013-10-23 02:31:12

标签: r ggplot2 geospatial heatmap shiny

我正在尝试开发一个闪亮的应用程序,它将在浏览器上安装基于地图的热图,并允许您更改热图上显示的变量。地图是具有GIS形状文件的地理区域,所选的变量随后在地图上显示为热图。不幸的是我遇到的问题是变量没有正确传递给ggplot()而我的地图失败了。下面作为直接R脚本的server.r代码成功运行没有问题但是当适应Shiny时它失败了。

问题出现在server.r中的以下ggplot()代码中:

myplot1 <- myplot1 + aes(long, lat, group = group, fill = input$var) + ...

我收到错误:

  

eval(expr,envir,enclos)中的错误:找不到对象'input'

这与fill = input$var有关,它无法识别从ui.r传递的input$varinput$var是ui.r中选择的变量(var1,var2等),以便在热图上显示。这是我在代码中唯一无法识别input$var的实例。我在这一行之前使用了print(str(input$var)),它清楚地保存了所需变量的名称。如果我对其进行硬编码(例如fill=var1),那么ggplot()工作正常并且地图显示正确。

我还在environment = environment()中使用了ggplot(),但这会产生另一个错误:

  

错误:提供给连续刻度的离散值

我的解释是,它正在寻找变量描述的数据框,而是获取单个值。

我觉得这是我想念的简单事情 - 我应该宣布或重新分配。我将非常感谢有关此人的任何见解,指导或反馈。非常感谢!!

# server.R

library(shiny)

library(maps)
library(mapdata)
library(sp)
library(maptools)
library(scales)
library(RColorBrewer)
library(ggplot2)
library(rgeos)
library(plyr)
library(reshape)
library(mapproj)
library(rgdal)
library(grid)
library(gridExtra)

setwd("C:/Shiny")


# Step 1 Read/loading the target shapefile
gregion = readOGR(dsn="C:/Shiny", layer="duid")

# Step 2 Get row numbers from .dbf / explicitly identifies attribute rows by the .dbf offset.
gregion@data$id = rownames(gregion@data)

# Step 3 Makes centroid (point layer) from polygon "FORTIFY"
gregion.points = fortify(gregion, region="id")

# Step 4 Reading in .csv which will be joined to .dbf using "MERGE"
mydata <- read.csv("c:/Shiny/dataset.txt")

# Step 5 Joins the points to their corresponding attributes and finalizes the data preparation 
gregion.df = join(gregion.points, gregion@data, by="id")

# Step 6 Merge makes an inner join of the shapefile's data frame and the .csv on a common item (usually the spatial key)
mygeomdata <- merge(gregion.df, mydata, by.x="UID", by.y="UID")


# Define server logic required to plot various variables as heatmap
# Step 7 Create map
shinyServer(function(input, output) {

  # Compute the forumla text in a reactive expression since it is 
  # shared by the output$caption and output$mapPlot expressions

  formulaText <- reactive({
    paste("Variable:", input$var)
  })

  # Return the formula text for printing as a caption
  output$caption <- renderText({
    formulaText()
  })


  output$mapPlot <- renderPlot({

    myplot1 <- ggplot(mygeomdata)
    myplot1 <- myplot1 + aes(long, lat, group = group, fill = input$var) + labs(x = "Easting", y = "Northing") + scale_fill_gradient(low = "ghostwhite", high = "steelblue")
    myplot1 <- myplot1 + geom_polygon()
    myplot1 <- myplot1 + coord_equal()

    print(myplot1)

  })

})

#ui.R

library(shiny)

shinyUI(pageWithSidebar(

  # Application title
  headerPanel("Mapping"),

  # Sidebar with controls to select the variable to plot
  #

  sidebarPanel(
    selectInput("var", "Variable:",
                list("Variable 1" = "var1",
                     "Variable 2" = "var2"))
  ),

  # Show the caption and plot of the requested variable

  mainPanel(
    h3(textOutput("caption")),

    plotOutput("mapPlot")
  )
))

数据集的示例mydata <- read.csv("c:/Shiny/dataset.txt")是:

UID var1    var2    var3    var4    var5    var6    var7
1   0   0.001   0   0   0   0   0
2   0   0   0   0   1   0   0
3   0   0   0   0   0   0   0
4   0   0   0   0   1   0   0
5   0   0   0   0   1   0   0
6   0   0   0   0   1   0   0
7   0   0   0   0   0   0   0
8   0   0.004   0.026   0   0   0   0
9   0.499   0.014   0   0.499   1   0   0.033
10  0.573   0.002   0.015   0.573   1   0   0.427
11  1   0.003   0.01    1   1   0   0

mygeomdata具有以下结构:

 $ UID       : int  1 1 1 1 1 1 1 1 1 1 ...
 $ long      : num  393121 392895 392895 392840 392839 ...
 $ lat       : num  5501404 5502275 5502275 5502489 5502494 ...
 $ order     : int  1 2 3 4 5 6 7 8 9 10 ...
 $ hole      : logi  FALSE FALSE FALSE FALSE FALSE FALSE ...
 $ piece     : Factor w/ 4 levels "1","2","3","4": 1 1 1 1 1 1 1 1 1 1 ...
 $ group     : Factor w/ 5693 levels "0.1","1.1","10.1",..: 1 1 1 1 1 1 1 1 1 1 ...
 $ id        : chr  "0" "0" "0" "0" ...
 $ DUID      : Factor w/ 5656 levels "130023362","130023367",..: 1 1 1 1 1 1 1 1 1 1 ...
 $ PC        : Factor w/ 3617 levels "0","ZZZ0A3","ZZZ0A4",..: 3271 3271 3271 3271 3271 3271 3271 3271 3271 3271 ...
 $ DUIDAREA  : num  21687 21687 21687 21687 21687 ...
 $ ELEV      : num  14.8 14.8 14.8 14.8 14.8 ...
 $ GroupUp   : int  2 2 2 2 2 2 2 2 2 2 ...
 $ GroupUpT  : Factor w/ 2 levels "A","B": 2 2 2 2 2 2 2 2 2 2 ...
 $ var1      : num  0 0 0 0 0 0 0 0 0 0 ...
 $ var2      : num  0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 ...
 $ var3      : num  0 0 0 0 0 0 0 0 0 0 ...
 $ var4      : num  0 0 0 0 0 0 0 0 0 0 ...
 $ var5      : int  0 0 0 0 0 0 0 0 0 0 ...
 $ var6      : num  0 0 0 0 0 0 0 0 0 0 ...
 $ var7      : num  0 0 0 0 0 0 0 0 0 0 ...

4 个答案:

答案 0 :(得分:2)

即使您在尝试将此错误移植到Shiny时收到此错误,也与ggplot有关。简而言之,您正在使用aes的局部变量和函数参数,这比最初看起来更棘手。

请参阅this SO question以获取完整的讨论和一些不错的选择。

要记住一件事:对于ggplot,如果你把所有需要的东西作为数据框的一部分保存,你会发现它更容易。

agstudy's response in this question解释了发生了什么。

希望有所帮助。

答案 1 :(得分:1)

谢谢Ram Narasimhan。你向我指出了一些非常有用的帖子,我能够在一些挖掘后找到解决方案。以下是修改过的部分:

environment<-environment() 
  myplot1 <- ggplot(mygeomdata, aes(long, lat, group = group, fill = get(input$var)), environment = environment) 

上面的代码放在renderPlot调用之外。 非常感谢!

答案 2 :(得分:1)

实际上,有一个更简单的解决方案。像之前一样在renderPlot()函数中使用ggplot(),但使用aes_string()而不是aes():

output$mapPlot <- renderPlot({
    ggplot(mygeomdata, aes_string(x=long, y=lat, group = group, fill = input$var)) + geom_point()
})

答案 3 :(得分:1)

aes_string(x=input$x, y=input$y)也有效。只需确保两个输入都被视为无功输入(例如input$)。 ggplot()不允许其中一个轴为手动轴,另一个轴为反应轴。