R语言chromoMap包报错及自定义基因颜色、添加图例问题

• 数据读取格式错误：你的输入文件为制表符分隔格式，read.csv默认以逗号作为分隔符，读取后所有字段会被合并到同一列，chromoMap无法识别所需的染色体、坐标字段。
• 变量未提前定义：代码中使用的col.set未做赋值操作，运行时会触发找不到对象的报错。
• 注释数据列顺序不符合要求：chromoMap要求注释数据的列顺序必须为染色体名称、起始位置、终止位置、注释属性，你提供的chr_anot第一列为基因名，顺序不匹配导致包无法解析坐标信息。

首先调整数据读取逻辑、列顺序，再指定分类配色模式即可，完整可运行代码如下：

# 若未安装依赖包先执行安装
# install.packages(c("chromoMap", "RColorBrewer"))
library(chromoMap)
library(RColorBrewer)

# 正确读取制表符分隔的输入文件
chr.data <- read.delim("chr_file.csv", header = TRUE, sep = "\t")
anno.data <- read.delim("chr_anot.csv", header = TRUE, sep = "\t")

# 调整注释数据列顺序，匹配chromoMap要求
anno.data <- anno.data[, c("chr", "start", "end", "gene")]

# 生成与基因数量匹配的配色集，也可自定义颜色向量
col.set <- brewer.pal(n = length(unique(anno.data$gene)), name = "Set2")

# 调用绘图函数
chromoMap(
  chrmList = list(chr.data),
  annoList = list(anno.data),
  labels = TRUE,
  data_based_color_map = TRUE,
  data_type = "categorical", # 指定数据为分类类型，按基因名分组配色
  data_color = list(col.set),
  legend = TRUE # 开启图例展示
)

说明：将data_type设置为categorical后，包会自动将注释数据第四列的基因名作为分类变量，每个分类对应分配一个颜色，开启legend=TRUE即可自动生成对应基因的颜色图例。

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