R函数中ggtitle无法识别SNP变量名的问题求助

问题原因 #

你遇到的报错是因为{{SNP}}在ggplot的美学映射（aes）中用于引用数据框的列，但在ggtitle()里，它会尝试解析成列的值而非列的名称字符串。函数运行时会去寻找名为chrD3.54068099的对象（即列的某个取值），但该对象不存在，因此触发报错。

解决方案 #

需要提取参数对应的变量名称字符串，替代直接使用{{SNP}}作为标题。以下提供两种可行方案：

方案1：使用rlang包（tidyverse风格） #

利用rlang包的ensym()将参数转换为符号，再用as_label()转为字符串，以此获取SNP列的名称作为标题。

修改后的函数代码：

boxplot_function <- function(Group,SNP) {
  # 获取SNP参数对应的列名
  snp_title <- rlang::as_label(rlang::ensym(SNP))
  
  Group %>%
    group_by(`Pheno`,{{SNP}}) %>%
    count({{SNP}}) %>%
    group_by(`Pheno`) %>%
    mutate(`Proportion` = prop.table(`n`)) %>%
    ggplot(aes(x = `Pheno`, y = `Proportion`, fill = factor({{SNP}}, levels = c("HOM_ALT", "HET", "HOM_REF")))) + 
    geom_bar(position = "stack", stat = "identity") +
    labs(x = "Phenotype", fill = "Genotype") +
    geom_text(aes(label = round(`Proportion`, 2)), position = "stack", vjust = 2) +
    ggtitle(snp_title)
}

方案2：基础R实现（无需额外包） #

使用基础R的deparse(substitute(SNP))直接获取变量名称字符串：

修改后的函数代码：

boxplot_function <- function(Group,SNP) {
  # 获取SNP参数对应的列名
  snp_title <- deparse(substitute(SNP))
  
  Group %>%
    group_by(`Pheno`,{{SNP}}) %>%
    count({{SNP}}) %>%
    group_by(`Pheno`) %>%
    mutate(`Proportion` = prop.table(`n`)) %>%
    ggplot(aes(x = `Pheno`, y = `Proportion`, fill = factor({{SNP}}, levels = c("HOM_ALT", "HET", "HOM_REF")))) + 
    geom_bar(position = "stack", stat = "identity") +
    labs(x = "Phenotype", fill = "Genotype") +
    geom_text(aes(label = round(`Proportion`, 2)), position = "stack", vjust = 2) +
    ggtitle(snp_title)
}

验证调用 #

使用你提供的示例数据调用函数：

boxplot_function(Superelaxed_vcf_topvariants, chrF2.9280651)

此时图表标题会正确显示为chrF2.9280651，函数可正常运行。

内容的提问来源于stack exchange，提问作者Balegde