`combine_vars（）`：！至少一层必须包含所有刻度变量：fact_wrap（）＆amp; facet_grid（）

Question

我正在尝试在rthat上创建一个函数，并在我的WD上保存带有刻面的点图。以下功能的代码：

get_dotplot <- function(df, xvalue, avgvalue, sdvalue, svalue, gvalue, main, xaxis, yaxis, glegend, figure_title)
{
  dp <- ggplot(df, aes(x = xvalue, y = avgvalue, color = gvalue)) + 
    
    geom_point(stat = 'identity', aes(shape=svalue, color=gvalue))+
    
    geom_errorbar(aes(ymin=avgvalue-sdvalue, ymax=avgvalue+sdvalue))+
    
            facet_grid(cols = vars(svalue), scales = "fixed")+

        labs(x = xaxis, y = yaxis, title = main, color=glegend)+
 
        theme(axis.title.x.bottom = element_text(hjust = 0.5, vjust = 1),
              axis.title.y = element_text(hjust = 0.5, vjust = 1),
 
              axis.ticks.x = element_line(),
              axis.text.x = element_text(angle = 0, hjust = 1, vjust = 0.5, size = 7),
axis.ticks.x.bottom = element_line(colour = "grey", size = (0.5)),
axis.ticks.y.left = element_line(colour = "black", size = (0.4)),

panel.background=element_rect(colour = "black", size = 0.5, fill=NA),
panel.grid = element_blank())

    print(dp)
    ggsave(paste(figure_title, "png", sep = "."), plot = dp, scale = 1, dpi = 600)
}

get_dotplot(df, xvalue, avgvalue, sdvalue, svalue, gvalue, main, xaxis, yaxis, glegend, figure_title)

但是，我始终收到此错误消息：

Error in `combine_vars()`:
! At least one layer must contain all faceting variables: `svalue`.
* Plot is missing `svalue`
* Layer 1 is missing `svalue`
* Layer 2 is missing `svalue`
Backtrace:
  1. global get_dotplot_errorbar_yaxis(...)
  3. ggplot2:::print.ggplot(dp)
  5. ggplot2:::ggplot_build.ggplot(x)
  6. layout$setup(data, plot$data, plot$plot_env)
  7. ggplot2 f(..., self = self)
  8. self$facet$compute_layout(data, self$facet_params)
  9. ggplot2 f(...)
 10. ggplot2::combine_vars(data, params$plot_env, cols, drop = params$drop)

我怀疑这是因为faceTing，所以我在facet_wrap（）和facet_grid（）之间进行了播放，而没有结果。有人可以帮我吗？ 我检查了数据框中的Svalue变量，并且拼写正确。我还咨询了以前有关该主题的问题，但它们没有帮助。

数据集看起来像这样，但是有大量的个人和数天：

set.seed(108)
n <- 1:12
treatment <- factor(paste("trt", 1:2))
individuals <- sample(LETTERS, 2)
days <- c("12", "20", "25")
avg_var1 <- sample(1:100, 12) 
sd_var1 <- sample(1:50, 12) 
avg_var2 <- sample(1:100, 12) 
sd_var2 <- sample(1:50, 12) 
avg_var3 <- sample(1:100, 12) 
sd_var3 <- sample(1:50, 12)
test <- data.frame(n, treatment, individuals, days,avg_var1, sd_var1, avg_var2, sd_var2, avg_var3, sd_var3)

我定义了该函数的变量，如R：

df=test
xvalue=test$days
avgvalue=test$avg_var1
sdvalue = test$sd_var1
svalue=test$treatment
gvalue=test$individuals
main= "var1 in function of days"
xaxis="days"
yaxis="var1"
glegend="individuals"
figure_title ="var1_days"

Answer 1

如书面，您的代码将列传递到功能中，重复数据框中的数据。对于GGPLOT中使用的非标准评估，这似乎不是“很好”。从本质上讲，GGPLOT正在寻找DF中的列，称为“ svalue”用于构面（找不到它）。一旦修复了这一点，误差线就会发生相同的问题。

一定的一种方式就是仅传递列名，然后使用aes_string对变量。这对于刻面或计算值不起作用，因此在函数开始时计算这些值。这将给出：

get_dotplot <- function(df, xvalue, avgvalue, sdvalue, svalue, gvalue, main, xaxis, yaxis, glegend, figure_title)
{
  df$ymin <- df[[avgvalue]] - df[[sdvalue]]
  df$ymax <- df[[avgvalue]] + df[[sdvalue]]
  df$facets <- df[[svalue]]
  
  dp <- ggplot(df, aes_string(x = xvalue, y = avgvalue, color = gvalue)) + 
    
    geom_point(stat = 'identity', aes_string(shape=svalue, color=gvalue)) +
    
    geom_errorbar(aes(ymin=ymin, ymax=ymax))+
    
    facet_grid(cols = vars(facets), scales = "fixed")+
    
    labs(x = xaxis, y = yaxis, title = main, color=glegend)+
    
    theme(axis.title.x.bottom = element_text(hjust = 0.5, vjust = 1),
          axis.title.y = element_text(hjust = 0.5, vjust = 1),
          
          axis.ticks.x = element_line(),
          axis.text.x = element_text(angle = 0, hjust = 1, vjust = 0.5, size = 7),
          axis.ticks.x.bottom = element_line(colour = "grey", size = (0.5)),
          axis.ticks.y.left = element_line(colour = "black", size = (0.4)),
          
          panel.background=element_rect(colour = "black", size = 0.5, fill=NA),
          panel.grid = element_blank())
  
  print(dp)
  ggsave(paste(figure_title, "png", sep = "."), plot = dp, scale = 1, dpi = 600)
}

get_dotplot(df=test,
            xvalue="days",
            avgvalue="avg_var1",
            sdvalue = "sd_var1",
            svalue="treatment",
            gvalue="individuals",
            main= "var1 in function of days",
            xaxis="days",
            yaxis="var1",
            glegend="individuals",
            figure_title ="var1_days")