R语言KEGG富集分析及结果解读：enrichKEGG函数使用指南

R语言KEGG富集分析及结果解读：enrichKEGG函数使用指南

本指南将带您了解如何在R语言中使用enrichKEGG函数进行KEGG通路富集分析，并将结果转换为易于理解和使用的格式。

代码示例：

edo <- enrichKEGG(foldchanges, organism = 'hsa', pvalueCutoff = 0.05, qvalueCutoff = 0.05)
edox <- setReadable(edox, 'org.Hs.eg.db', 'ENTREZID')
geneset <- cbind(edox@result[['Description']],edox@result[['geneID']])

代码解释：

1. enrichKEGG函数: - 该函数用于对差异表达基因列表 (foldchanges) 进行KEGG通路富集分析。 - organism = 'hsa' 指定物种为人类 (hsa)。 - pvalueCutoff = 0.05 和 qvalueCutoff = 0.05 设置显著性水平的p值和q值阈值，用于筛选显著富集的通路。
2. setReadable函数: - 该函数将富集分析结果 edox 中的基因ID转换为可读性更好的ENTREZID格式。
3. cbind函数: - 将富集结果中的通路描述 (edox@result[['Description']]) 和 ENTREZID格式的基因ID (edox@result[['geneID']]) 合并到 geneset 数据框中，方便后续分析和可视化。

总结：

通过以上代码，您可以利用R语言轻松地进行KEGG富集分析，并将结果转换为易于理解和使用的格式，为后续研究提供数据支持。