火山图标记基因_火山地形图

火山图标记基因_火山地形图现在很多文章会要求火山图中显示某一个基因名称,或者是某一些基因的名称,这时候我们应该怎么做呢?很多人会选择使用AI或者PS来做,在这里我介绍的是,如何用R绘制这样的图加载R包library(ggpubr)library(ggthemes)如果没有安装,可以使用下面的代码进行安装install.packages(‘ggthemes’)install.packages(‘ggpubr’)…

大家好,又见面了,我是你们的朋友全栈君。如果您正在找激活码,请点击查看最新教程,关注关注公众号 “全栈程序员社区” 获取激活教程,可能之前旧版本教程已经失效.最新Idea2022.1教程亲测有效,一键激活。

Jetbrains全系列IDE稳定放心使用

现在很多文章开始出现这样的一种情况,在绘制火山图中,显示我们所关注的基因,那么如何去显示呢?很多人可能会这么做,在绘制普通的火山图之后,使用AI对图进行修改,添加部分基因,但是现在我要介绍的是如何用R绘制
library(ggpubr)
library(ggthemes)
data <- read.csv(“easy_input_limma.csv”, head=T,sep=’,’)
#绘制基本热图
data l o g p &lt; − − l o g 10 ( d a t a logp &lt;- -log10(data logp<log10(dataadj.P.Val)
ggscatter(data,x=“logFC”,y=“logp”)+theme_base()
在这里插入图片描述
#新增一列,选出上调和下调
data G r o u p = &quot; n o t − s i g n i f i c a n t &quot; d a t a Group=&quot;not-significant&quot; data Group=notsignificantdataGroup[which((dataKaTeX parse error: Expected ‘EOF’, got ‘&’ at position 16: adj.P.Val<0.05)&̲(datalogFC > 2))]=“up-regulated”
data G r o u p [ w h i c h ( ( d a t a Group[which((data Group[which((dataadj.P.Val<0.05)&(dataKaTeX parse error: Expected ‘EOF’, got ‘#’ at position 32: …own-regulated” #̲查看上调和下调的基因个数 ta…Group)
#绘制新的火山图
ggscatter(data,x=“logFC”,y=“logp”,color=“Group”)+theme_base()
在这里插入图片描述
#添加颜色和点的大小
ggscatter(data,x=“logFC”,y=“logp”,color=“Group”,palette=c(
“blue”,“gray”,“red”),size=1)+theme_base()
在这里插入图片描述
ggscatter(data,x=“logFC”,y=“logp”,color=“Group”,palette=c(
“blue”,“gray”,“red”),size=1)+theme_base()+geom_hline(yintercept=1.3,
linetype=“dashed”)+geom_vline(xintercept=c(-2,2),linetype=“dashed”)
在这里插入图片描述
#新增一列,这一列是需要展示的基因名在这里我展示前10
data l a b e l = &quot; &quot; d a t a &lt; − d a t a [ o r d e r ( d a t a label=&quot;&quot; data&lt;-data[order(data label=data<data[order(datalogp),]
#data S y m b o l &lt; − r o w n a m e s ( d a t a ) u p . g e n e &lt; − h e a d ( d a t a Symbol&lt;-rownames(data) up.gene&lt;-head(data Symbol<rownames(data)up.gene<head(dataSymbol[which(data G r o u p = = &quot; u p − r e g u l a t e d &quot; ) ] , 10 ) d o w n . g e n e &lt; − h e a d ( d a t a Group==&quot;up-regulated&quot;)],10) down.gene&lt;-head(data Group==upregulated)],10)down.gene<head(dataSymbol[which(dataKaTeX parse error: Expected ‘EOF’, got ‘#’ at position 31: …gulated”)],10) #̲将up.gene和down.g…label[match(data.top.genes,dataKaTeX parse error: Expected ‘EOF’, got ‘#’ at position 26: …data.top.genes #̲重新绘制 ggscatter(…label,font_label=8,repel=T)+theme_base()+geom_hline(yintercept=0.05,
linetype=“dashed”)+geom_vline(xintercept=c(-2,2),linetype=“dashed”)
在这里插入图片描述
很多人会说,这样的终究不太分明,我想让点的大小来反映FC,这样可以做到吗?当然可以了!如果有想了解的,可以试试,其实很简单的。

最近看到一个更加强大的包,绘制的火山图更加好看
EnhancedVolcano包
1、首先安装这个包
if (!requireNamespace(‘devtools’, quietly = TRUE))
install.packages(‘devtools’)
devtools::install_github(‘kevinblighe/EnhancedVolcano’)
加载包
library(EnhancedVolcano)
#res1是带有基因ID的列名,log2FoldChange是去过对数的差异倍数的列名,pvalue可以选择是校正前和校正后的,xlim是对log2FoldChange的范围进行限定
EnhancedVolcano(res1, lab = rownames(res1), x = ‘log2FoldChange’, y = ‘pvalue’,xlim = c(-5, 8))
绘制之后就是这样了
在这里插入图片描述
2、 还可以通过这个包直接修改点的大小和标签的大小
#title图片的名称, transcriptPointSize点的大小, transcriptLabSize标签的大小
EnhancedVolcano(res2,lab = rownames(res2), x = ‘log2FoldChange’,y = ‘pvalue’,
xlim = c(-8, 8),title = ‘N061011 versus N61311’,pCutoff = 10e-16,
FCcutoff = 1.5, transcriptPointSize = 1.5,transcriptLabSize = 3.0)
3、调整点的颜色和透明度
#colAlpha点的透明度,col是点的颜色,默认第一个是NS,第二个是log2FC,第三个是P,第四个是p&log2FC
EnhancedVolcano(res2, lab = rownames(res2),x = ‘log2FoldChange’,
y = ‘pvalue’,xlim = c(-8, 8),title = ‘N061011 versus N61311’, pCutoff = 10e-16,
FCcutoff = 1.5, transcriptPointSize = 1.5, transcriptLabSize = 3.0,
col=c(“grey30”, “forestgreen”, “royalblue”, “red2”), colAlpha = 1)
4、修改点的形状
shape点的形状,可以选择一个,也可以选择多个
选择一个形状
EnhancedVolcano(res2, lab = rownames(res2),x = ‘log2FoldChange’,
y = ‘pvalue’,xlim = c(-8, 8),title = ‘N061011 versus N61311’, pCutoff = 10e-16,
FCcutoff = 1.5, transcriptPointSize = 1.5, transcriptLabSize = 3.0,shape=4,
col=c(“grey30”, “forestgreen”, “royalblue”, “red2”), colAlpha = 1)
选择多个形状
EnhancedVolcano(res2, lab = rownames(res2),x = ‘log2FoldChange’,
y = ‘pvalue’,xlim = c(-8, 8),title = ‘N061011 versus N61311’, pCutoff = 10e-16,
FCcutoff = 1.5, transcriptPointSize = 1.5, transcriptLabSize = 3.0,shape=c(1,4,23,25),col=c(“grey30”, “forestgreen”, “royalblue”, “red2”), colAlpha = 1)
5、添加额外的阈值线,自定义阈值线,可以自己查看文档
hlineType线条的类型(平行于x轴),vline—是平行于Y轴的
cutoffLineCol线条的颜色,cutoffLineWidth线条的粗细
EnhancedVolcano(res2,lab = rownames(res2), x = ‘log2FoldChange’,y = ‘pvalue’,
xlim = c(-6, 6),title = ‘N061011 versus N61311’,pCutoff = 10e-12,FCcutoff = 1.5,transcriptPointSize = 1.5,transcriptLabSize = 3.0,colAlpha = 1,
cutoffLineType = ‘blank’,cutoffLineCol = ‘black’, cutoffLineWidth = 0.8,
hline = c(10e-12, 10e-36, 10e-60, 10e-84),
hlineCol = c(‘grey0’, ‘grey25’,‘grey50’,‘grey75’),
hlineType = ‘longdash’,hlineWidth = 0.8
,gridlines.major = FALSE,gridlines.minor = FALSE)
6、调整图例
#legend图例内容,legendPosition图例的位置,legendLabSize文字的大小
legendIconSize标签的大小
EnhancedVolcano(res2,lab = rownames(res2),x = ‘log2FoldChange’,y = ‘pvalue’,
xlim = c(-6, 6),pCutoff = 10e-12, FCcutoff = 1.5,cutoffLineType = ‘twodash’,cutoffLineWidth = 0.8,transcriptPointSize = 3.0,transcriptLabSize = 4.0,
colAlpha = 1,legend=c(‘NS’,‘Log (base 2) fold-change’,‘P value’, ‘P value & Log (base 2) fold-change’), legendPosition = ‘right’,legendLabSize = 16,legendIconSize = 5.0)
详细的请查看说明文档

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 举报,一经查实,本站将立刻删除。

发布者:全栈程序员-用户IM,转载请注明出处:https://javaforall.cn/179441.html原文链接:https://javaforall.cn

【正版授权,激活自己账号】: Jetbrains全家桶Ide使用,1年售后保障,每天仅需1毛

【官方授权 正版激活】: 官方授权 正版激活 支持Jetbrains家族下所有IDE 使用个人JB账号...

(0)
blank

相关推荐

  • Edge breaker记录

    Edge breaker记录和代码相关性非常强的是这篇文章,几乎就是伪代码了:http://citeseerx.ist.psu.edu/viewdoc/download?doi=10.1.1.58.7918&rep=rep1&type=pdf一、Compression阶段命令参数为umbrellla_times4.offB第一个参数就是输入网格,目前支持OFF和OVTable两种文件格式,这里…

  • SEVERE: Error filterStart 错误解决「建议收藏」

    SEVERE: Error filterStart 错误解决「建议收藏」使用jfinal写的项目,去掉jetty包,加入javaee.jar包,然后打包,部署到ubuntu的tomcat6下运行

  • stack overflow at line 9[通俗易懂]

    stack overflow at line 9[通俗易懂]web开发中,从设备得到一些信息显示在web页面中的文本框中,结果显示对话框1.起初认为可能是文本框的字符限制原因,有默认长度,故加上maxlength=128足够长,结果问题还是存在。2.因为这个问题在goole火狐 搜狗浏览器中都不存在,在IE6 360极速浏览器中存在,所以认为可能是浏览器的设置问题,从网上找了一下,解决办法是:打开IE的工具-Internet选项

  • 网站渗透测试

    网站渗透测试目录[TOC]公司的网站需要渗透测试,学习了一下渗透测试的方法,记录下,方便后期查阅。(1)暴力激活成功教程1.1风险分析:数据传输过程使用非加密的http协议,因此可对数据传输过程进行抓包分析;用户名、密码明文,且未设置验证码,导致可进行暴力激活成功教程,以获取身份凭证信息1.2加固建议1:使用https加密传输可以在apache或者tomcat下配置好证书,启用https就ok,网上很多

  • MySQL8.0正确修改密码的姿势[通俗易懂]

    MySQL8.0正确修改密码的姿势[通俗易懂]mysql更新完密码,总是拒绝连接、登录失败?MySQL8.0不能通过直接修改mysql.user表来更改密码。正确更改密码的方式备注:清空root密码MySQL8.0不能通过直接修改mysql.user表来更改密码。因为authentication_string字段下只能是MySQL加密后的43位字符串密码,其他的导致错误。错误不报出,但是无法再登录mysql,总是会提示无…

  • 深度学习中的batch(batch size,full batch,mini batch, online learning)、iterations与epoch

    概念引入    我们已知在梯度下降中需要对所有样本进行处理过后然后走一步,那么如果我们的样本规模的特别大的话效率就会比较低。假如有500万,甚至5000万个样本(在我们的业务场景中,一般有几千万行,有些大数据有10亿行)的话走一轮迭代就会非常的耗时。这个时候的梯度下降使用了全部的样本数据,所以叫做fullbatch。为了提高效率,我们可以把样本分成等量的子集。例如我们把500万样本分成1…

发表回复

您的电子邮箱地址不会被公开。

关注全栈程序员社区公众号