R Markdown新窗口显示火山图时基因标签缺失问题求助

问题背景 #

使用Seurat筛选差异基因后，通过EnhancedVolcano在R Markdown中绘制火山图，在R Markdown内嵌预览中基因标签显示正常，但点击“在新窗口显示”按钮后，所有基因标签消失。

解决思路 #

• 调整图形设备尺寸
  新窗口默认画布尺寸可能过小，导致标签被裁剪或无法渲染。可以在绘图前手动指定窗口大小，或者在R Markdown代码块中设置输出尺寸：

  # 手动打开指定尺寸的新窗口
  dev.new(width = 12, height = 10)
  
  EnhancedVolcano(DEgenes,
      lab = rownames(DEgenes),
      x = 'avg_log2FC',
      y = 'p_val_adj',
      ylim = c(0,450),
      title = 'InEx vs ReA SF',
      xlab = bquote(~Log[2]~ '(Fold Change)'),
      ylab = bquote(~-Log[10]~ '(Adj. p)'),
      axisLabSize = 12,
      FCcutoff = 1.493, 
      pCutoff = 10e-40,
      pointSize = 1,
      drawConnectors = TRUE,
      widthConnectors = 0.75,
      colAlpha = 1,
      max.overlaps = 45)
  

  或者在R Markdown代码块中直接设置：

  # 你的EnhancedVolcano绘图代码
  

• 优化标签渲染参数
  EnhancedVolcano基于ggplot2，新窗口的渲染引擎可能对文本处理逻辑不同，调整以下参数可改善标签显示：

  • 增加labSize设置标签大小
  • 调大max.overlaps允许更多标签同时显示
  • 设置force增强标签的排斥力度，避免被挤出画布

  EnhancedVolcano(DEgenes,
      lab = rownames(DEgenes),
      x = 'avg_log2FC',
      y = 'p_val_adj',
      ylim = c(0,450),
      title = 'InEx vs ReA SF',
      xlab = bquote(~Log[2]~ '(Fold Change)'),
      ylab = bquote(~-Log[10]~ '(Adj. p)'),
      axisLabSize = 12,
      labSize = 3, # 增大标签尺寸
      FCcutoff = 1.493, 
      pCutoff = 10e-40,
      pointSize = 1,
      drawConnectors = TRUE,
      widthConnectors = 0.75,
      colAlpha = 1,
      max.overlaps = 100, # 放宽重叠限制
      force = 2 # 增强标签排斥力
  )
  

• 提升R Markdown输出分辨率
  如果是HTML输出，新窗口打开的可能是低分辨率缩略图，在R Markdown的YAML头部设置更高dpi：

  output:
    html_document:
      dpi: 300
  

• 直接导出高清图片
  绕过R Markdown的预览机制，直接用ggsave()导出高清图片到本地，打开后标签会正常显示：

  # 先将绘图对象保存为变量
  p <- EnhancedVolcano(DEgenes,
      lab = rownames(DEgenes),
      x = 'avg_log2FC',
      y = 'p_val_adj',
      ylim = c(0,450),
      title = 'InEx vs ReA SF',
      xlab = bquote(~Log[2]~ '(Fold Change)'),
      ylab = bquote(~-Log[10]~ '(Adj. p)'),
      axisLabSize = 12,
      FCcutoff = 1.493, 
      pCutoff = 10e-40,
      pointSize = 1,
      drawConnectors = TRUE,
      widthConnectors = 0.75,
      colAlpha = 1,
      max.overlaps = 45)
  
  # 导出高清图片
  ggsave("volcano_plot.png", p, width = 12, height = 10, dpi = 300)
  

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