R语言GO富集分析报错求助
R studio运行别人共享的关于差异基因GO富集分析作图的代码,报错如图。检索网页后未找到解决办法,更新R包到最新版后依然报错。请前辈们支招,谢谢!
报错:
Error in has_pairsim(x) :
Term similarity matrix not available. Please use pairwise_termsim function to deal with the results of enrichment analysis.
出错部分:
> for(i in c("BP","CC","MF")){
+ #GO???????
+ kk=enrichGO(gene = gene,OrgDb = org.Hs.eg.db, pvalueCutoff =1, qvalueCutoff = 1, ont=i, readable =T)
+ GO=as.data.frame(kk)
+ GO=GO[(GO$pvalue<pvalueFilter & GO$qvalue<qvalueFilter),]
+ write.table(GO,file=paste0(i,".txt"),sep="\t",quote=F,row.names = F) #???渻????
+ showNum=30
+ if(nrow(GO)<30){
+ showNum=nrow(GO)
+ }
+ if(nrow(GO)!=0){
+ #?״ͼ
+ pdf(file=paste0(i,".barplot.pdf"),width = 11,height = 7)
+ bar=barplot(kk, drop = TRUE, showCategory =showNum,color = colorSel)
+ print(bar)
+ dev.off()
+
+ #??ͼ
+ pdf(file=paste0(i,".bubble.pdf"),width = 11,height = 7)
+ bub=dotplot(kk,showCategory = showNum, orderBy = "GeneRatio", color = colorSel)
+ print(bub)
+ dev.off()
+
+ #Ȧͼ
+ pdf(file=paste0(i,".circos.pdf"),width = 8,height = 5.5)
+ cnet=cnetplot(kk, foldChange=geneFC, circular = TRUE, colorEdge = TRUE)
+ print(cnet)
+ dev.off()
+
+ #Enrichment Map
+ pdf(file=paste0(i,".emapplot.pdf"),width = 7,height = 6)
+ emap=emapplot(kk, showCategory = showNum, color = colorSel)
+ print(emap)
+ dev.off()
+ }
+ }
R语言GO富集分析报错求助相关推荐
- R语言-处理异常值或报错的三个示例
R语言-处理异常值或报错的三个示例 参考文章: (1)R语言-处理异常值或报错的三个示例 (2)https://www.cnblogs.com/awishfullyway/p/6611690.html ...
- R语言ggplot绘制地图-报错汇总(一)
R语言ggplot绘制地图-报错汇总 报错两例 报错1: 报错2: 报错两例 在用ggplot绘制地图时出现了两个报错,网上搜索了没有相关说明,虽然解决方式很蠢,但是可能对于出现同样报错的人会有帮助, ...
- R语言 Rstudio 安装 包 报错installation of package ‘*’ had non-zero exit status
在安装r语言向量自回归包(vars)时候报错了 记录一下供其他朋友参考,有同类相似其他问题也可以联系我一起解决. 先贴报错信息 > install.packages("v ...
- R语言加载xlsx报错错误: JAVA_HOME cannot be determined from the Registry解决方法
错误: package or namespace load failed for 'xlsx': loadNamespace()里算'rJava'时.onLoad失败了,详细内容: 调用: fu ...
- R语言下载CLL包报错arning in install.packages : package ‘CLL’ is not available for this version of RA ve
错误: Warning in install.packages :package 'CLL' is not available for this version of R.A version of t ...
- R语言学习 - 富集分析泡泡图
功能富集泡泡图 功能富集分析用来展示某一组基因(一般是单个样品上调或下调的基因)倾向参与哪些功能调控通路,对从整体理解变化了的基因的功能和潜在的调控意义具有指导作用,也是文章发表中一个有意义的美图.通 ...
- R语言armasubsets使用及报错问题处理
时间序列分析 在时间序列分析的建模过程中,我们常用armasubsets给我们的模型进行定阶,优化. armasubsets() in TSA package,判别准则BIC. armasubsets ...
- R语言处理FASTQ文件报错相关
setwd("E:/CTQ/plw/3w3d")##设置路径 install.packages("Biostrings")##下载包 > seqdata ...
- R语言|GO富集分析
GO富集分析 -- 差异圈图.聚类圈图 一.数据链接 链接: https://pan.baidu.com/s/1PPUW5YyJHjwxvkjmMi0Vtw?pwd=c9s8 提取码: c9s8 二. ...
最新文章
- java 对象被回收的例子_Java对象的后事处理——垃圾回收(二)
- 二维字符数组按长度排序_字符串长度 字符数组长度
- javascript面向对象包装类Class的类库解析
- 虚拟化系列-Windows server 2012 虚拟机管理
- 深山红叶PE工具箱嫦娥一号纪念版 V30[1115]
- Symbol MC1000 扫描 冲突问题 把下面文件做成scanwedge.reg的注册表文件,放在Application重起
- Python SQLAlchemy --3
- mysql innodb文件_MySQL的InnoDB文件介绍
- 【今日CV 视觉论文速览】16 Nov 2018
- js-this作用域
- .Net Core 1.1打包发布到Linux
- 事务消息和普通消息的区别
- LinuxShell笔记3 条件控制和流程控制
- CSS概述 CSS声明
- PADS logic 和PADS layout 连接
- cmd命令快速重启IIS服务
- vs2003网站项目加载失败解决方案
- 双精度,单精度和半精度
- monkey命令——压力测试——转载参考09
- 康得新董事长是谁?_您人生董事会中的谁?