安装
1
2
3
|
install.package("RIdeogram")
# 直接安装出现报错,需要依赖rsvg包
conda install -c conda-forge librsvg r-rsvg
|
使用
准备输入文件
- karyotype.csv 核型文件,必须。包含5列
- 第1列 染色体名
- 第2, 3列 染色体起止位置(0-base)
- 第4, 5列 着丝粒起止位置(可省略)
- heatmap.csv 染色体heatmap,可选。包含4列
- 第1列 染色体名
- 第2, 3列 染色体指定位置
- 第4列 该位置对应的数值
- mark.csv 染色体旁形状标记,可选。包含6列
- 第1列 染色体名
- 第2列 形状(box, triangle, circle)
- 第3列 染色体名
- 第4, 5列 染色体指定位置
- 第6列 颜色
绘制
1
2
3
4
5
6
7
8
9
10
11
12
|
ideogram(karyotype = karyotype,
#如果没提供easy_input_heatmap.csv,就在下面两行前面加#
overlaid = gene_density,
#染色体上heatmap的配色
colorset1 = c("forestgreen", "gold", "deeppink3"),
#如果没提供easy_input_mark.csv,就在下面这行前面加#
label = Random_RNAs_500,
#染色体宽度,默认170
width = 200,
#图例的位置
Lx = 160, #图例左上角跟左边的距离
Ly = 20) #图例左上角跟顶端的距离
|
格式转换
1
2
3
4
|
#svg2tiff, svg2jpg, svg2png函数设定类似
svg2pdf("chromosome.svg",
width = 12, height = 8, #设定画布大小
dpi = 300)
|
参考来源
https://cran.r-project.org/web/packages/RIdeogram/vignettes/RIdeogram.html