GWAS分析中曼哈顿图如何添加基因信息

大家好,我是邓飞,虽然我知道GWAS分析后一定是要有曼哈顿图的,没有曼哈顿图的GWAS没有灵魂,但是,谁能想到,需要在曼哈顿图上添加上基因,怎么不在上面画蒙娜丽莎呢???

既然有需求,那就研究一下吧,毕竟,解决问题就意味着能力提升,做困难的事情,能力提升得更快,像我这种凡事依赖ChatGPT的人,如果ChatGPT能解决,我是肯定不会看说明文档的。

之前,我看过高手和小白的区别,小白都去知乎和csdn上找答案了,然后关注了一堆公众号。高手直接去看说明文档了。我充分体验到从小白到高手的过程,要是能当小白,我肯定不愿意看说明文档的。但是网上真的找不到,ChatGPT也是错误百出,于是,我就翻开说明文档,一下就搞定了,很快的。有点不真实,赶紧写篇博客,介绍一下如何对曼哈顿图中添加基因信息。

应用场景:

GWAS分析后,对结果进行曼哈顿图绘制,对显著性位点基因型基因注释

曼哈顿图的数据哪里来?当然是GWAS分析而来,GWAS有不同的模型,得到的结果都是类似的:SNP、染色体、物理位置、P值,这些就够了:

read-normal-img

read-normal-img

read-normal-img

read-normal-img

下面就是把这个基因名,添加到上面的曼哈顿图上面了。

关键参数介绍:

  • highlight = xx$SNP,这个是指明需要高亮的位点,SNP名称对应的又染色体和物理位置的信息,所以只需要指定SNP的名称就行了,一个向量,有几个基因,就高亮几个SNP,一个也行,八个也行。我曾经准备了所有的SNP信息,事实证明我想复杂了。

  highlight.text = xx$gene,基因名和SNP需要一一对应。

read-normal-img

read-normal-img

如果你还没有学会,那就点击文末链接系统学习吧!

GWAS分析先做后学

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值