基因本体(Gene Ontology)类别
在上一个练习中,您测试了生物通路的富集情况。现在,您将测试一组已知共同影响同一生物学过程的基因集合的富集情况,即基因本体(GO)类别。 来自白血病研究的差异表达基因中,哪些 GO 类别被过度代表?
本练习是课程的一部分
在 R 中使用 limma 进行差异表达分析
练习说明
第 2 章白血病研究的拟合模型对象 fit2 已加载到您的工作区。已加载 limma 包。
从数据框
fit2$genes中提取 Entrez 基因 ID。使用
goana测试富集的 GO 类别。将物种设为"Hs",表示 Homo sapiens(人类)。使用
topGO查看前 20 个富集的 GO 类别。将参数ontology设为"BP",以返回 Biological Processes(生物过程)。
交互式实操练习
通过完成这段示例代码来试试这个练习。
# Extract the entrez gene IDs
entrez <- ___
# Test for enriched GO categories
enrich_go <- ___(fit2[1:500, ], geneid = ___, species = ___)
# View the top 20 enriched GO Biological Processes
___(enrich_go, ontology = ___)