metascape中聚类,然后对其中KEGG通路中基因进行定位
metascape是一个比较好的进行基因聚类分析的网站。
聚类后会发现有一些基因聚集到了KEGG通路中,但是这些基因具体在通路中的哪一个位置起作用,。
具体的基因可以通过点击图片中的Web得到对应的基因
KEGG通路可以通过点击ath04075得到kegg
当然是可以通过手动查询的。
但是一个一个的动手去搜索很累,就想到到了通过python 辅助查询
脚本如下:
from splinter import Browser
import re
b=Browser('chrome')
url1=input('请输入metascape上heatonecluster网页网址:')
number=input('请输入该聚类中需要查看的基因的个数:')
url2=input('请输入KEGG对应聚类的网页网址:')
b.visit(url1)
text=b.find_by_xpath('//*[@id="chart"]').text
ress=text.split('\n')[:int(number)]
res=[]
pa='\(.+?\)'
for gene in ress:r=re.findall(pa,gene)[0].lstrip('(').rstrip(')')res.append(r)
b.visit(url2)
i=1
for gene in res:b.fill('query',gene+'\n')print('已经查询%d/%s个基因'%(i,number))input1=input('继续查询请按enter键(终止请输入stop):')i+=1if input1!='stop':continueelse:break
然后按enter键就可以了。
metascape中聚类,然后对其中KEGG通路中基因进行定位相关推荐
- 如何获取kegg通路的基因列表
如何获取kegg通路的基因列表?相信很多人都有这个疑惑. 哪些时候需要知道通路的基因列表?举两个例子: 在找差异基因这一步,只想针对特定通路的基因去找.可以在找完差异基因后,和通路的基因列表取个交集. ...
- pvrect r语言 聚类_R语言实现KEGG通路富集可视化
用过KEGG的朋友应该都很熟悉里面的通路地图.你是否想过如果自己可以控制通路图将自己的基因绘制在一个通路图中,那么今天给大家介绍一个新推出的Bioconductor软件包pathview.这个包可以进 ...
- 20W+喜爱的Pathview网页版 | 整合表达谱数据KEGG通路可视化
Pathview网站简介 网址:https://pathview.uncc.edu/ 前段时间介绍了一个R包 - Pathview.它可以整合表达谱数据并可视化KEGG通路,操作是先自动下载KEGG官 ...
- R语者小case之——从KEGG原始网页批量生成通路的基因表格
我们经常要用到KEGG数据库来对基因做功能分析.经常长得好看的朋友问:如何获得整个通路的基因? 其实我们有多种方法可以获得通路中所有的基因情况,本文通过KEGG的原始网页生成某个通路的基因表格. 准备 ...
- Pathview包:整合表达谱数据可视化KEGG通路
Pathview是一个用于整合表达谱数据并用于可视化KEGG通路的一个R包,其会先下载KEGG官网上的通路图,然后整合输入数据对通路图进行再次渲染,从而对KEGG通路图进行一定程度上的个性化处理,并且 ...
- Pathview绘制KEGG通路图
今天小编和大家分享一个在KEGG通路分析中可能会经常用到的一个R包,名字叫Pathview.我将分五期进行,从快速上手,再到高级用法,最后在用实例对这个包做一个较全面的讲解.下面我们先简单了解一下这个 ...
- 靶中化合物设计/靶点及信号通路验证之酶抑制剂靶点预测
靶中化合物设计/靶点及信号通路验证之酶抑制剂靶点预测 酶抑制剂靶点预测,包括药物靶点评估.化合物库定制.AI虚拟筛选药物服务.分子对接筛选定制.药效团筛选定制.靶中化合物设计.靶点及信号通路验证,小编 ...
- 数据挖掘中聚类算法概述
1 聚类方法概述 聚类方法是将物理或抽象对象的集合组成为由类似的对象组成的多个类的过程被成为聚类.由聚类所组成的簇是一组数据对象的集合,这些对象与同一簇中的对象彼此类似,与其他簇中的对象相异.在许多应 ...
- linux下kegg注释软件,KEGG数据中全部代谢反应和代谢物注释信息的下载
# 加载函数与R包 ----------------------------------------------------------------- library(KEGGREST) librar ...
最新文章
- 有没有一种方法可以缓存GitHub凭证来推送提交?
- 通过数据挖掘组织营销潜力的三个重要途径
- Vue监听器与监听滥用
- 远程桌面连接出现身份验证错误。 要求的函数不受支持,这可能是由于 CredSSP 加密 Oracle 修正。...
- CF1039E-Summer Oenothera Exhibition【LCT,根号分治】
- mysql 深胡_Mysql胡说八道
- 解析DataWorks数据集成中测试连通性失败问题
- 分析go程序内存逃逸情况
- [IDDFS+背包] 洛谷P2744 [USACO5.3]量取牛奶Milk Measuring
- Translatium for Mac(多语言在线翻译工具)
- C++ string()什么意思
- Topaz Sharpen AI 人工智能锐化插件
- 小麦亩产一千八(kela)
- “0X????????”指令引用的“0x00000000”内存,该内存不能为“read”或“written”
- Java常用的开发软件下载地址以及问题解决
- Android 直接拨号和调用拨号盘
- java 新浪短网址生成器_新浪短网址链接生成器 2020最新新浪t.cn短链接在线生成接口推荐...
- Scratch迷宫之自动生成二 画迷宫
- 线上故障等级和处理规范,看阿里测试工程师如何拿捏
- 长沙银河计算机中等职业学校图片,长沙银河中等职业学校
热门文章
- 一篇文章带你入门zabbix监控系统
- 微信小程序上下固定,中间滚动,使用position:fixed遮挡正文
- 02java进阶03-异常、线程、同步、线程池、Lambda表达式、File类、递归
- 佛山植物组织培养实验室必备设备
- 如何在虚拟机的电脑里增加一块硬盘?
- mac 出现 Failed to create parent directories for [/data/applogs/xxl-job/xxl-job-admin.log]
- thinkpad T430 :Intel 7系主板在XP下装AHCI驱动的方法
- 字节数据-十六进制字符串-字符串转换工具类
- 逻辑思维在编程中的重要性!!!(回答问题有感)
- 2022-10-16