WGCNA剔除异常样本

WGCNA剔除异常样本

做WGCNA分析,需要对基因和样本进行一定的筛选。一般来说,会采用如下的筛选方案:

1.  筛选掉基因表达丰低, 在各个样本中变化不大的基因

2. 筛选掉异常的样本

基因的筛选,比较好处理。但是异常样本比较难办,可以参考WGCNA的一份代码,参考如下的代码:

# 计算相似性矩阵

A=adjacency(t(datExpr),type="signed")
# 计算网络的连接度
k=as.numeric(apply(A,2,sum))-1
# 连接度标准化
Z.k=scale(k)
# 设置连接度筛选阀值, 这部分是关键,筛选掉异常的样本
thresholdZ.k=-2.5
# 将异常样本进行标注
outlierColor=ifelse(Z.k<thresholdZ.k,"red","black")

# 删除异常样本
remove.samples= Z.k<thresholdZ.k | is.na(Z.k)
datExpr=datExpr[!remove.samples,]
datTraits=datTraits[!remove.samples,]


如果您对WGCNA数据挖掘感兴趣,请学习我们的WGCNA相关课程:

WGCNA加权基因共表达网络分析


  • 发表于 2018-05-25 10:53
  • 阅读 ( 6661 )
  • 分类:转录组

0 条评论

请先 登录 后评论
microRNA
microRNA

115 篇文章

作家榜 »

  1. omicsgene 698 文章
  2. 安生水 347 文章
  3. Daitoue 167 文章
  4. 生物女学霸 120 文章
  5. xun 82 文章
  6. 红橙子 78 文章
  7. rzx 74 文章
  8. CORNERSTONE 72 文章