autoplot转录本信息for循环画图并保存

话不多说，直接上代码

if (!requireNamespace("BiocManager", quietly = TRUE))install.packages("BiocManager")
BiocManager::install("ggbio",force = TRUE)
rm(list=ls())
library(ggbio)   #安装ggbio这个包#BiocManager::install("EnsDb.Mmusculus.v79")
BiocManager::install("EnsDb.Mmusculus.v75") #下载小鼠的注释ensDB包并加载#library(EnsDb.Mmusculus.v79)
library(EnsDb.Mmusculus.v75)

注释包的下载地址：
https://bioconductor.org/packages/3.15/data/annotation/

# hg38
#BiocManager::install("EnsDb.Hsapiens.v86")
#library(EnsDb.Hsapiens.v86)
#ensdb <- EnsDb.Hsapiens.v86
#ensdb <- EnsDb.Mmusculus.v79
ensdb <- EnsDb.Mmusculus.v75#基因名字的大小写，一定要和ensembl上的对应物种的对应起来，不然会报错
#如果用的是人的ensdb数据库，则基因名全部都是大写，如CDC45，在小鼠里面是Cdc45（小鼠里面，首字母大写）
#基因名一定不能错，不然会报错，和ensdb数据库里面的信息匹配不上#panel.border() = element_blank() 去掉图片的黑色边框
#panel.background = element_blank() 去掉背景颜色
#axis.line.x = element_line(colour = "black", size = 1, lineend = "square") 给横轴添加颜色和横线setwd("C:\\Users\\gaoshuang\\Desktop\\结果文件\\基因的转录本图\\results") #设置工作保存图的路径

开始画图

#画染色体的图，染色体X,Y，这个字母是需要大写的
#IRanges 里面放的是染色体位点信息
library(GenomicRanges)
p.ideo <- Ideogram(genome = "mm10")
p.ideo + xlim(GRanges("chrX", IRanges(161202097, 161209198)))  #X Y染色体都要大写


#画单个基因的转录本剪切信息
autoplot(ensdb,GeneNameFilter('Scml2'),gap.geom = "chevron",label.color = "black", color = "brown",fill = "brown") + theme(panel.border = element_blank(),panel.background = element_blank(),axis.line.x = element_line(colour = "black", size = 1, lineend = "square"))

#把各个转录本的id换为基因名
autoplot(ensdb, ~ symbol == "Scml2", names.expr="gene_name") + theme(panel.border = element_blank(),panel.background = element_blank(),axis.line.x = element_line(colour = "black", size = 1, lineend = "square"))

#画单个基因的转录本信息
autoplot(ensdb,GeneNameFilter("Scml2")) + theme(panel.border = element_blank(),panel.background = element_blank(),axis.line.x = element_line(colour = "black", size = 1, lineend = "square")) #画单个基因的转录本剪切信息(限制染色体位点信息的)
autoplot(ensdb,GeneNameFilter('Scml2'),label.color = "black",fill = "#1b9e77") + theme(panel.border = element_blank(),panel.background = element_blank(),axis.line.x = element_line(colour = "black", size = 1, lineend = "square"))+scale_x_continuous(limits = c(161234696,161237492),breaks = c(seq(161234696,161237492,by=1000)),position = "top")

使用R包ggbio画单个基因的 autoplot转录本信息图相关推荐

单个基因集富集分析泡泡图绘制
富集分析是生物信息分析中快速了解目标基因或目标区域功能倾向性的最重要方法之一.其中代表性的计算方式有两种: 一是基于筛选的差异基因,采用超几何检验判断上调或下调基因在哪些GO或KEGG或其它定义的通路 ...
RNA-seq流程学习笔记（14）-在windows10平台上利用R包合并表达矩阵、设置实验分组信息、列名及数据的导入导出
参考文章: RNA-seq(6): reads计数,合并矩阵并进行注释进入R学习和相关操作后,各种折腾已经快2周了,看了几个网站的教程,借了几本教科书,发现都是零零碎碎的知识.索性就按照之前的方法, ...
r语言echarts画箱线图_R语言学习 - 箱线图（小提琴图、抖动图、区域散点图）
箱线图箱线图是能同时反映数据统计量和整体分布,又很漂亮的展示图.在2014年的Nature Method上有2篇Correspondence论述了使用箱线图的好处和一个在线绘制箱线图的工具.就这样都 ...
R语言效用分析 ( 效能分析、Power analysis)、除了pwr包之外还有其它包、例如、基因研究中的效能分析、MBESS包可用于各种形式的效能分析和最少样本量确定、其他效用分析包的简要介绍
R语言效用分析 ( 效能分析.Power analysis).除了pwr包之外还有其它包.例如.基因研究中的效能分析(power analysis).MBESS包可用于各种形式的效能分析(power ...
r roc函数_画ROC曲线的R包总结
原标题:画ROC曲线的R包总结作者: Joseph Rickert 原文链接: https://rviews.rstudio.com/2019/03/01/some-r-packages-for-r ...
一个R包完成单细胞基因集富集分析（全代码）
singleseqgset | 单细胞RNA-Seq基因集富集分析 NGS系列文章包括NGS基础.转录组分析 (Nature重磅综述|关于RNA-seq你想知道的全在这).ChIP-seq分析 (Ch ...
人类基因id转换的R包：GeoTcgaData
一个非常全的人类基因id转换工具 GeoTcgaData包中有一个id转换的函数:id_conversion_vector() , 它可以对人类的各种基因id进行转换.2019年9月第一次发布,目 ...
R语言survival包coxph函数构建cox回归模型、ggrisk包ggrisk函数可视化Cox回归的风险评分图、使用风险得分的中位数计算最佳截断值cutoff（基于LIRI基因数据集）
R语言survival包coxph函数构建cox回归模型.ggrisk包ggrisk函数可视化Cox回归的风险评分图.使用风险得分的中位数计算最佳截断值cutoff(基于LIRI基因数据集) 目录
Microbiome：animalcules-交互式微生物组分析和可视化的R包
animalcules-交互式微生物组分析和可视化的R包 animalcules: interactive microbiome analytics and visualization in R Mi ...

使用R包ggbio画单个基因的 autoplot转录本信息图

autoplot转录本信息for循环画图并保存

开始画图

使用R包ggbio画单个基因的 autoplot转录本信息图相关推荐

最新文章

热门文章