BioMart

参考

R语言biomaRt工具包学习笔记

使用实例

if (!requireNamespace("BiocManager", quietly = TRUE))install.packages("BiocManager")
## 由于BioMart由R 3.X构建，而我本地R是4.X，因此这里进行R版本控制
BiocManager::install(version = '3.14')
BiocManager::install("biomaRt")
library(biomaRt)
### 设置查询数据库和数据集，以小鼠为例
ensembl <- useEnsembl(biomart = "genes",   dataset = "mmusculus_gene_ensembl")

运行时最后一行报错

Ensembl site unresponsive, trying useast mirror
Error in bmRequest(request = request, httr_config = httr_config, verbose = verbose) : Gateway Timeout (HTTP 504).

查询后似乎是ensembl.org网址不可用，该包也支持换其他数据库，但我这里没有再进行尝试。希望知道方法的大神帮忙指导。

clusterProfiler

参考

Ensemble ID格式及转换

使用实例

### 安装 clusterProfiler
### 先安装
if (!requireNamespace("BiocManager", quietly = TRUE))install.packages("BiocManager")BiocManager::install(version = '3.14')
BiocManager::install("clusterProfiler")
BiocManager::install("org.Mm.eg.db")library(clusterProfiler)
library(org.Mm.eg.db) # 小鼠数据库### 从rawcount.txt文件读取数据
data <- read.table("rawcounts.txt",head=TRUE)gene.ens.id<- data['ensemblid']
gene.ens.id<- c(t(gene.ens.id))## 有的ensembl id有版本号，直接转不行的, 这句代码是去除版本号的
gene.ens.id <- gsub("\\..*", "",  gene.ens.id)
gene.symbol <- bitr(geneID = gene.ens.id, fromType = "ENSEMBL",toType = c("SYMBOL"),OrgDb = org.Mm.eg.db)
gene.symbol
# 将转化结果保存到文件result.csv
write.table(gene.symbol,"result.csv",row.names=FALSE,col.names=TRUE,sep=",")

此处有33%Ensembl ID未找到对应Symbol

In bitr(geneID = gene.ens.id, fromType = "ENSEMBL", toType = c("SYMBOL"),  :33.01% of input gene IDs are fail to map...

基因表达数据处理-Ensembl ID转化为Gene symbol相关推荐

gene ID / Gene Symbol / Ensembl ID
1. 各种ID名称介绍 Gene ID 也称Entrez ID/EntrezGene ID ,是 NCBI 使用的能够对众多数据库进行联合搜索的搜索引擎, 其对不同的 Gene 进行了编号, 每个 g ...
FPKM值基因表达量的计算、基因ID转gene symbol的例子
高通量测序数据一般公司都会提供两种矩阵,一种是Row counts,前面说过的用于差异基因的筛选.第二种是FPKM值,可以理解为转录组基因的表达矩阵,可以用于做热图和基因表达变化的比较.但是数据挖掘中 ...
Entrez ID 和 gene symbol 和 Ensembl ID转换
简单粗暴网页版 bioDBnet - Biological Database Network db2db模块可以非常方便地实现各种ID之间的转换,当然 bioDBnet的功能有很多,可以试试. 手动迁 ...
uniport ID 转换为 gene symbol（ID转换）
网站 ID转换网址 input 选择 uniprot accession output 选择 Gene symbol 一般一次转换不能超过1万个基因ID,数据量太大就转换不出来,如果一次转换太多,需要 ...
获取gtf文件gene symbol ENSID gene_biotype
首先下载gtf文件,这里我们引用的是Ensembl的文件enensembl gtf文件下载这里面我们下载完文件后我们如何查看这个文件信息呢,首先我们用UEStudio 打开后我们看一下文件的数据结构 ...
android R.id.转化为view
LayoutInflater inflater=(LayoutInflater)getSystemService(Context.LAYOUT_INFLATER_SERVICE);View view= ...
在cygwin安装HOMER和最全使用说明
HOMER 基本功能 HOMER(Hypergeometric Optimization of Motif EnRichment)用于Motif identification和下一代测序分析,是一款很 ...
GEO芯片数据下载和探针ID转换（保姆级教程）
GEO芯片数据下载和探针ID转换(保姆级教程) 一.问题描述探针ID转换数据是否预处理过二.Rstudio的安装(建议阅读,避免后续转换时出错) 安装包的下载安装步骤三.(正文)芯片数据下载 ...
常用数据库的基因ID
1.Ensembl stable ID : Ensembl stable ID 的结构是根据不同物种设置的前缀, 加上数据所指的类型, 如基因蛋白质, 再加上一系列的数字. 有的时候可以有不同的版本, ...
差异基因 p log2foldchange_拟南芥的基因ID批量转换？差异基因，GO/KEGG数据库注释（转录组直接送你全套流程）...
新手遇到的问题都是类似的,比如批量ID转换虽然我写过大量的教程:ID转换大全不过都需要R基础,因为是大批量转换啊! 但热心肠的植物生物信息学教学大佬还是友善的给出了解决方案我也狗尾续貂制作了 ...

基因表达数据处理-Ensembl ID转化为Gene symbol

Ensembl ID转化为Gene symbol

BioMart

参考

使用实例

clusterProfiler

参考

使用实例

基因表达数据处理-Ensembl ID转化为Gene symbol相关推荐

最新文章

热门文章