1

我正在尝试调整x并ymax在ggplot2包中的以下脚本中R调整 的绘图坐标errorbar，但是它返回错误。

  gplot <- function(prd) {

        ggplot() +
        geom_polygon(data=shp.t,aes(x=long,y=lat,group=group),
                      fill="white",colour="grey") +
        ## Plot errorbar
        geom_errorbar(data=te10.cent,size=2,colour="red",
                      alpha=.8,width=0,
                      aes_(x=quote(long.cent)-350,ymin=quote(lat.cent),　　　　　　　
                                 ymax=quote(lat.cent)+prd))
        }
gplot("Field Name") # Not number but field name of the data frame

（抱歉，我无法上传我正在使用的实际数据框。）

这些是我面临的错误：

Error in quote(lat.cent) + prd * .pt : non-numeric argument to binary Operator
Error in quote(lat.cent)+prd * .pt : non-numeric argument to binary Operator

-350如果和+prd从脚本中省略或在数据框中使用带有实际变量的“aes”，它在这两种情况下都有效。我尝试了其他脚本；"long.cent"-350而"lat.cent"+prd不是上面的脚本，但是它也返回相同的错误。

我搜索了解决方案，但所有解决方案都解释了如何在aes_不混合参数和算术计算的情况下使用。我需要将非标准表达式与算术计算混合来调整我的情节但是如何？

4

2 回答 2

2

我相信这将解决您的问题：

gplot <- function(prd) {

  ggplot() +
    geom_polygon(data = shp.t,
                 aes_(x = ~long,
                      y = ~lat,
                      group = ~group),
                 fill = "white",
                 colour = "grey") +
    ## Plot errorbar
    geom_errorbar(data = te10.cent,
                  size = 2,
                  colour = "red",
                  alpha = .8,
                  width = 0,
                  aes_(x = ~long.cent - 350,
                       ymin = ~lat.cent,　　　　　　　
                       ymax = ~lat.cent + prd))
}

可重现的例子：

一些数据： 

library(tidyverse)

data(iris)

iris %>%
  group_by(Species) %>%
  summarise_all(~mean(.)) -> summed_iris

gplot <- function(prd){
  ggplot(summed_iris) +
    geom_col(aes_(x = ~Species,
                  y = ~Sepal.Length))+
    geom_errorbar(aes_(x = ~Species,
                       ymin = ~Sepal.Length -prd,
                       ymax = ~Sepal.Length + prd))

}

gplot(0.5)

在此处输入图像描述

编辑：评论中的问题：

如果prd数据的列名是什么时候，最好预先计算这些值：

gplot <- function(prd){
  ymin <-  with(summed_iris, get("Sepal.Length") - get(prd))
  ymax <-  with(summed_iris, get("Sepal.Length") + get(prd))
  summed_iris <- data.frame(summed_iris, ymin, ymax)
  ggplot(summed_iris) +
    geom_col(aes_(x = ~Species,
                  y = ~Sepal.Length))+
    geom_errorbar(aes_(x = ~Species,
                       ymin = ~ymin,
                       ymax = ~ymax))

}
gplot("Petal.Length")

在此处输入图像描述

于 2018-06-11T13:54:55.823 回答
0

随着即将推出的 ggplot2 2.3.0（将于 2018 年 6 月下旬发布），aes_(...)已软弃用，建议使用 tidy eval。这意味着在此处的上下文中，!!prd当您想要变量的内容prd而不是符号时编写prd。

library(ggplot2)
library(dplyr)
#> 
#> Attaching package: 'dplyr'
#> The following objects are masked from 'package:stats':
#> 
#>     filter, lag
#> The following objects are masked from 'package:base':
#> 
#>     intersect, setdiff, setequal, union

iris %>%
  group_by(Species) %>%
  summarise_all(~mean(.)) -> summed_iris

gplot <- function(prd){
  ggplot(summed_iris) +
    geom_col(aes(x = Species, y = Sepal.Length)) +
    geom_errorbar(
      aes(x = Species, ymin = Sepal.Length - !!prd, ymax = Sepal.Length + !!prd)
    )
}

gplot(0.5)

如果您希望函数gplot()能够像 dos 一样接受符号aes()，那么您需要enquo()在函数定义的开头添加一个。它捕获提供给函数的表达式以及封闭环境。

gplot <- function(prd){
  prd <- enquo(prd)
  ggplot(summed_iris) +
    geom_col(aes(x = Species, y = Sepal.Length)) +
    geom_errorbar(
      aes(x = Species, ymin = Sepal.Length - !!prd, ymax = Sepal.Length + !!prd)
    )
}

gplot(Petal.Length)


scale <- 1/5
gplot(scale*Petal.Length)

于 2018-06-12T01:57:37.607 回答