本文小编为大家详细介绍“WGCNA怎么根据KME值筛选hub基因”,内容详细,步骤清晰,细节处理妥当,希望这篇“WGCNA怎么根据KME值筛选hub基因”文章能帮助大家解决疑惑,下面跟着小编的思路慢慢深入,一起来学习新知识吧。
WGCNA分析中需要找出模块中也就共表达网络中的hub基因:
这里可以根据:KME (eigengene connectivity)值来筛选hub基因:Hub genes are those that show most connections in the network as indicated by their high KME (eigengene connectivity) value
代码部分很简单,输出所有的基因与模块的KME矩阵,然后筛选每个模块中kme最大的前几个基因就是hub基因:
#mergedMEs是经过相似模块合并之后的最终模块对应的模块特征基因:
MEs = mergedMEs
#输出结果:
datKME=signedKME(datExpr, MEs, outputColumnName="kME_MM.")
write.csv(datKME, "kME_MM_test.csv")
读到这里,这篇“WGCNA怎么根据KME值筛选hub基因”文章已经介绍完毕,想要掌握这篇文章的知识点还需要大家自己动手实践使用过才能领会,如果想了解更多相关内容的文章,欢迎关注天达云行业资讯频道。