R_非常简单的gene_denesity

R gene_d in chromsome

install.packages('RIdeogram')
require(RIdeogram)
file2 <- read.table('染色体大小文件.txt',header = T )
head(file2)
class(genne_d)
ideogram(karyotype = file2)
convertSVG("chromosome.svg", device = "png")
density <- read.table("ceshi2", sep = "\t", header = T,stringsAsFactors = T)
head(density)
ideogram(karyotype = file2, overlaid = density)
convertSVG("chromosome.svg", device = "png")

©著作权归作者所有,转载或内容合作请联系作者
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

推荐阅读更多精彩内容