WGCNA如何从module中挖掘关键基因

本文介绍了在WGCNA分析中,如何利用connectivity、module membership (MM) 和 gene significance (GS) 三个统计量在相关模块中筛选关键基因。connectivity衡量基因在网络中的连接度,MM值反映基因与模块主成分的相关性,GS值表示基因表达与表型的相关性。通过设定阈值,如GS>0.2且MM>0.8,可以找出与特定表型高度相关的hub gene,进而进行功能富集分析。
摘要由CSDN通过智能技术生成

欢迎关注微信公众号《生信修炼手册》!

识别到与表型数据相关的modules之后,还可以在该modules中进一步筛选基因,为了方便筛选,对于每个基因定义了以下三个统计量

1. connectivity

连接度,在之前的文章中,简单提过这个概念,类似于网络中节点的degree的概念,只不过在加权共表达网络中,由于每条边代表两个基因间的相关性的大小,对应一个数值,所以一个基因在共表达网络中的连接度定义为与该基因相连的所有边的数值之和。

另外,根据相连的基因是否和该基因位于同一个module, 又可以将边分为两类,和该基因位于同一个module内,定义为within, 位于不同的modules, 定义为out。在WGCNA中,可以通过intramodularConnectivity函数计算连接度,用法如下

# 计算基因间的邻接值
ADJ1=abs(cor(datExpr,use="p"))^6

# colorh1代表module的颜色
Alldegrees1=intramodularConnectivity(ADJ1, colorh1)

计算的结果如下

> head(Alldegrees1)
kTotal kWithin kOut kDiff
Gene1 31.80186 28.37595 3.425906 24.95005
Gene2 28.88249 26.47896 2.403522 24.07544
Gene3 25.38600 23.11852 2.267486 20.85103
<
  • 6
    点赞
  • 46
    收藏
    觉得还不错? 一键收藏
  • 6
    评论
评论 6
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值