前面给大家介绍了
☞ 【实战】circleplot展示GO富集分析结果—附R代码
也用前面讲差异表达分析时,用到的TCGA数据库中的CHOL这套数据的差异表达分析结果来给大家演示了这个R包的使用方法。TCGA数据差异表达分析可以参考前面的文章。
当然circleplot这个函数,对于GEO数据库中的数据,也同样适用。只需要得到相应的差异表达分析结果就可以。GEO中数据差异表达分析可以参考前面的文章。
如果我们仔细研究circle_dat这个函数,你会发现。这个函数对差异表达分析的结果,要求其实很简单。只需要包含两列,第一列是基因的名字,第二列是差异表达分析的倍数改变logFC。☞circleplot展示GO富集分析结果文中,我们给出了一个例子
而对GO富集分析结果的要求包含category(BP,MF,CC),ID:GO号,term:GO条目的具体名称,富集分析矫正之后的p-value,还有富集到每个GO条目上的具体的基因名称。☞circleplot展示GO富集分析结果文中,我们也给出了一个例子。其实DAVID富集分析也可以得到类似的结果。
☞基因富集工具DAVID介绍(四)-GO富集分析及柱形图展示
我们知道KEGG富集分析的结果跟GO富集分析结果是有区别的。下面是KEGG富集分析的结果。
似乎circleplot这个函数是专门用来展示GO富集分析结果的。
但是我们在下面这篇文章中确实看到了circleplot展示的KEGG富集分析结果。
图A是GO富集分析结果,图B就是KEGG富集分析的结果。
经过小编不懈的努力,最终终于也成功画出了KEGG的circleplot。
关于GO和KEGG富集分析,以及柱形图和气泡图展示富集分析结果。可以参考往期内容。
完整KEGG富集分析+数据格式处理+circleplot的R代码+详细注释☟☟☟
【R语言】circleplot展示KEGG富集分析结果