我正在尝试开发一个闪亮的应用程序,它将在浏览器上安装基于地图的热图,并允许您更改热图上显示的变量。地图是具有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$var
。 input$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 ...
答案 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()
不允许其中一个轴为手动轴,另一个轴为反应轴。