如何在R Studio中使用ChromoMap为每个基因设置不同颜色并添加图例

解决方案 #

步骤1：调整注释数据结构

要实现单个基因独立配色，需要先在注释数据中新增一列作为分类分组标识，每个基因对应唯一的分组值，直接用基因名作为分组值即可：

anno.data$gene_group <- anno.data[,4] 
# 备注：anno.data[,4]为注释文件中存储基因名的列，可根据你实际的文件列结构调整索引

步骤2：修改绘图参数，开启图例并配置单基因配色

library(chromoMap)
# 配色数量需和你的基因总数保持一致，可自行调整色值避免重复
col.set = c("purple", "#4CBB17","#0096FF", "blue", "brown")
chr.data <- read.csv("chr_file.csv", header=T)
anno.data <- read.csv("chr_anot.csv", header = T)
# 新增基因分组列
anno.data$gene_group <- anno.data[,4]

chromoMap(
  chrlists = list(chr.data),
  annotations = list(anno.data),
  labels = T, 
  # 开启基于数据的配色规则
  data_based_color_map = T,
  # 指定配色数据为分类类型
  data_type = "categorical",
  # 指定注释数据中用作分组的列索引，此处对应新增的gene_group列
  data_cols = c(5),
  # 传入对应配色
  data_color = list(col.set),
  # 开启图例
  legend = T,
  # 自定义图例标题
  legend_title = "Gene Name"
)

补充说明

• 如果注释文件列结构不同，需调整data_cols的数值，对应你新增的基因分组列的实际位置
• 如果需要调整图例位置、字体大小，可额外添加legend_x、legend_y、legend_font_size参数自定义

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