搜索
查看: 2494|回复: 0

请教WGCNA中计算软阈值的问题?

[复制链接]

1

主题

1

帖子

36

积分

新手上路

Rank: 1

积分
36
发表于 2017-7-25 14:34:42 | 显示全部楼层 |阅读模式
为什么随着软阈值增加, 反而Scale Free Topology Model Fit,signed R^2会降低这么多呢,求教?
# Choose a set of soft-thresholding powers
powers = c(c(1:10), seq(from = 12, to=20, by=2))
# Call the network topology analysis function
sft = pickSoftThreshold(datExpr, powerVector = powers, verbose = 5)#确定软阈值
# Plot the results:
sizeGrWindow(9, 5)
par(mfrow = c(1,2));
cex1 = 0.9;
# Scale-free topology fit index as a function of the soft-thresholding power
plot(sft$fitIndices[,1], -sign(sft$fitIndices[,3])*sft$fitIndices[,2],
     xlab="Soft Threshold (power)",ylab="Scale Free Topology Model Fit,signed R^2",type="n",
     main = paste("Scale independence"));
text(sft$fitIndices[,1], -sign(sft$fitIndices[,3])*sft$fitIndices[,2],
     labels=powers,cex=cex1,col="red");



8e4e84a50e3a9d8ef6f3391f4af5a0a.png
c35a1835b40dea1a31d365ffbef605b.png



上一篇:了解参考基因组及基因组注释信息
下一篇:(3)转录组之数据质控
回复

使用道具 举报

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

本版积分规则

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

GMT+8, 2020-2-20 02:18 , Processed in 0.022373 second(s), 28 queries .

Powered by Discuz! X3.2

© 2001-2013 Comsenz Inc.