搜索
查看: 1258|回复: 1

[mRNA-seq] WGCNA模块太多怎么办?

[复制链接]

7

主题

9

帖子

114

积分

注册会员

Rank: 2

积分
114
发表于 2019-10-11 22:31:53 | 显示全部楼层 |阅读模式
WGCNA分析中(原始数据61000多个基因,处理掉7000多个缺失基因,最后保留53000多基因进行下一步分析),我在用一步法网络构建时模块太多:
1)第一次分析,按照默认参数分析的结果有96个模块
2)第二次分析,调整了参数,mergeCutHeight = 0.5,分析出的模块依然有71个,还是很多,代码如下
[AppleScript] 纯文本查看 复制代码
net = blockwiseModules(datExpr, power = 6, maxBlockSize = 6000,
                       TOMType = "unsigned", minModuleSize = 30,
                       reassignThreshold = 0, mergeCutHeight = 0.5,
                       numericLabels = TRUE, pamRespectsDendro = FALSE,
                       saveTOMs = TRUE,
                       saveTOMFileBase = "AS-green-FPKM-TOM",
                       verbose = 3)


我不敢再继续调整mergeCutHeight的数值了,怕结果不可信,现在我该解决这个问题?
回复

使用道具 举报

0

主题

4

帖子

265

积分

中级会员

Rank: 3Rank: 3

积分
265
发表于 2019-10-14 11:39:32 | 显示全部楼层
有必要用那么多基因吗。
过滤下不行吗 按表达量 差异  变异系数等
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

QQ|手机版|小黑屋|生信技能树 ( 粤ICP备15016384号  

GMT+8, 2020-7-13 10:12 , Processed in 0.028560 second(s), 28 queries .

Powered by Discuz! X3.2

© 2001-2013 Comsenz Inc.