2021-12-24 KEGG

对于报错:
None of the keys entered are valid keys for 'SYMBOL'.


重新修改代码:

setwd("C:\\Users\\Administrator.DESKTOP-4UQ3Q0K\\Desktop")
library(clusterProfiler)
library("readxl")
d<-read.csv("mk_fm_grp2all.csv") 
names(geneList) <- as.character(d[,1])
geneList<-rownames(d) 
library(org.Mm.eg.db)
library(org.Hs.eg.db)
library(AnnotationHub)
library(dbplyr)
library(BiocFileCache)
geneList <- sort(geneList, decreasing = TRUE)
geneList = d$X

此时的geneList:


geneEntrezID <- bitr(geneList, fromType="SYMBOL", toType=c("ENTREZID","UNIPROT"), OrgDb="org.Hs.eg.db")

成功:
'select()' returned 1:many mapping between keys and columns
Warning message:
In bitr(geneList, fromType = "SYMBOL", toType = c("ENTREZID", "UNIPROT"), :
17.38% of input gene IDs are fail to map...

kk <- enrichKEGG(gene= geneEntrezID[,2],organism= 'human',pvalueCutoff = 0.05)
head(kk)
browseKEGG(kk, 'hsa04066')

第二次输入browseKEGG(kk, 'hsa04066'),就没有反应了,浏览器里也没有图,解决方式:



在网址里修改,改为所需的hsa04066


image.png
©著作权归作者所有,转载或内容合作请联系作者
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

推荐阅读更多精彩内容

  • 16宿命:用概率思维提高你的胜算 以前的我是风险厌恶者,不喜欢去冒险,但是人生放弃了冒险,也就放弃了无数的可能。 ...
    yichen大刀阅读 6,135评论 0 4
  • 公元:2019年11月28日19时42分农历:二零一九年 十一月 初三日 戌时干支:己亥乙亥己巳甲戌当月节气:立冬...
    石放阅读 6,948评论 0 2
  • 年纪越大,人的反应就越迟钝,脑子就越不好使,计划稍有变化,就容易手忙脚乱,乱了方寸。 “玩坏了”也是如此,不但会乱...
    玩坏了阅读 2,208评论 2 1
  • 感动 我在你的眼里的样子,就是你的样子。 相互内化 没有绝对的善恶 有因必有果 当你以自己的价值观幸福感去要求其他...
    周粥粥叭阅读 1,672评论 1 5
  • 昨天考过了阿里规范,心里舒坦了好多,敲代码也犹如神助。早早完成工作回家喽
    常亚星阅读 3,063评论 0 1