重庆分公司
重庆分公司，新征程启航
为企业提供网站建设、域名注册、服务器等服务
R语言中如何进行GO注释和富集分析
R语言中如何进行GO注释和富集分析，针对这个问题，这篇文章详细介绍了相对应的分析和解答，希望可以帮助更多想解决这个问题的小伙伴找到更简单易行的方法。
成都创新互联专注于企业营销型网站、网站重做改版、下冶网站定制设计、自适应品牌网站建设、H5技术、成都做商城网站、集团公司官网建设、外贸营销网站建设、高端网站制作、响应式网页设计等建站业务，价格优惠性价比高，为下冶等各大城市提供网站开发制作服务。
GO注释和富集分析使用TBtools完成
小编使用的数据是甜樱桃叶绿体蛋白编码基因做GO注释，然后挑部分基因做富集分析，挑选的基因是
rpoC1
rpoB
rpoA
rpoC2
atpI
atpF
atpE
atpH
atpB
atpA
accD
rbcL
rpl22
rpl23
rpl20
rps8
rps7
rps16
rps15
rps14
rps18
做完富集分析得到文件GOenrichmentOutput.txt..GO.Enrichment.final.xls
根据GOplot包的示例数据挑选出其中的5列
Class GO_Name GO_ID GenesOfSelectedSetInGOterm corrected p-value(BH method)
作为数据集1
数据集2包括
ID,logFC,AveExpr,t,P.Value,adj.P.Val,B
数据集2的列变量应该都是转录组数据分析的结果
比如logFC应该是倍数变化Fold change 然后取log
AveExpr应该是平均表达量等
然后模仿帮助文档的例子构造数据集
help(package="GOplot")
library(GOplot)
data(EC)
file1<-file.choose()
file2<-file.choose()
df1<-read.table(file1,sep="\t",header=T)
df2<-read.csv(file2,header=T)
colnames(df1)<-colnames(EC$david)
df<-circle_dat(df1,df2)
dim(df)
df
气泡图
GOBubble(df)
加一些参数的气泡图
GOBubble(df,table.legend = F,table.col = T,ID=T,
display = "multiple")+
scale_x_continuous(limits=c(0,5))+
scale_y_continuous(limits = c(0,10))
弦图
pdf("chordpra.pdf",height=15,width = 15)
chord<-chord_dat(df,sample(df$genes,6),sample(df$term,8))
GOChord(chord)
dev.off()
现在基本可以根据自己的数据来构造GOplot的输入文件
关于R语言中如何进行GO注释和富集分析问题的解答就分享到这里了，希望以上内容可以对大家有一定的帮助，如果你还有很多疑惑没有解开，可以关注创新互联行业资讯频道了解更多相关知识。
新闻标题：R语言中如何进行GO注释和富集分析
分享路径：http://cqcxhl.com/article/peccgg.html