运行RStudio,设置工作路径
查询工作路径
getwd()
设置工作路径
setwd("D:/RStudio/data")
包安装路径
.libPaths()
染色体基因位置图
#安装RIdeogram
install.packages('RIdeogram')
#调用RIdeogram
library(RIdeogram)
karyotype <- read.table("Ob_karyotype.txt", sep = "\t", header = T, stringsAsFactors = F)
density <- read.table("Ob_density.txt", sep = "\t", header = T, stringsAsFactors = F)
location <- read.table("Ob_location.txt", sep = "\t", header = T, stringsAsFactors = F)
ideogram(karyotype = karyotype, overlaid = density,label=location,label_type = "marker")
convertSVG("chromosome.svg", device = "tiff", dpi = 300)
convertSVG("chromosome.svg", device = "png")
Os.jpg
调整染色体颜色
ideogram(karyotype = karyotype, overlaid = density,label=location,label_type = "marker",colorset1 = c("#fc8d59", "#ffffbf", "#91bfdb"))
ideogram(karyotype = karyotype, overlaid = density,label=location,label_type = "marker",colorset1 = c("#4575b4", "#ffffbf", "#d73027"))#默认颜色
colorset1 = c("#f7f7f7", "#e34a33"), colorset2 = c("#f7f7f7", "#2c7fb8"),colorset1 = c("#e5f5f9", "#99d8c9", "#2ca25f"),colorset1 = c("#99FFFF","#66FFFF","#33FFFF")
Or.jpg
网友评论