大家好,我是邓飞,最近星球(飞哥的知识)有老师问了一个问题:
GAPIT软件,染色体的颜色是5个一循环,他有12个染色体,想每条染色体一个颜色绘制一条染色体:
我的回答:GAPIT大概率没有参数设置,但是可以把结果文件用CMplot进行可视化,这个肯定是没问题的,我回头写篇博客。那么,博客来了:
1,数据
library(CMplot)
data("pig60K")
# 挑选10000个位点
library(tidyverse)
set.seed(123)
dd = pig60K %>% sample_n(10000)
dd1 = dd %>% filter(Chromosome %in% 1:12)
dd1 %>% count(Chromosome)
head(dd1)
使用CMplot软件包的数据,提取12条染色体作为演示: