R ChromoMap散点图Y轴科学计数法显示及图例优化问题

方案1：预处理数据（推荐） #

你的p_wald数值极小（均在1e-4量级），直接显示原始值易导致图例精度丢失、Y轴刻度可读性差。最通用的解决方式是将p值转换为**-log₁₀(p值)**，数值范围会变为3~4左右，既解决图例显示异常，也让Y轴刻度更直观：

步骤1：转换数据 #

# 对p_wald取负对数，添加微小偏移避免log(0)报错
marker_ann_trans$neg_log_p <- -log10(marker_ann_trans$p_wald + 1e-10)

步骤2：修改绘图脚本 #

将数据列替换为转换后的结果，同步调整图例标签：

chromoMap(ch.file=list(chromosomes_sort), 
          data.files=list(marker_ann_trans[, c("rs", "chr", "ps", "ps_end", "neg_log_p")]),
          n_win.factor = 3,  
          data_based_color_map = T, 
          data_type = "numeric",
          canvas_width = 1400, 
          canvas_height = 1400, 
          left_margin = 80, 
          interactivity = F,
          legend = c(T), 
          lg_y = 750, 
          plot.legend.labels = c("-log10(Wald p-val)"), 
          data_colors = list(c("red2","blue")), 
          plots = "scatter", 
          ch_gap = 0.1, 
          ref_line = T,
          plot_color = c("black"))

该方案无需修改包底层代码，符合GWAS可视化常规做法，可彻底解决图例显示0.0的问题，Y轴刻度也会呈现清晰的整数/小数。

方案2：修改生成的HTML文件（保留原始p值） #

若必须显示原始p值，由于ChromoMap未提供直接设置轴刻度格式的参数，需手动修改生成的HTML文件：

• 运行原始绘图代码后，工作目录会生成默认名为chromoMap_plot.html的文件；
• 用文本编辑器打开该文件，找到包含d3.format的JavaScript代码段（负责刻度格式化）；
• 将Y轴刻度格式从类似d3.format(".2f")改为d3.format(".2e")（.2e表示保留两位小数的科学计数法）；
• 针对图例精度问题，同步修改图例刻度的格式化代码，比如改为d3.format(".4f")或科学计数法格式。

此方案适合必须保留原始p值的场景，但每次重新绘图都需重复操作。

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