搜索
查看: 2138|回复: 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");



本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有帐号?立即注册

x



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

使用道具 举报

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

本版积分规则

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

GMT+8, 2019-10-14 22:12 , Processed in 0.036828 second(s), 27 queries .

Powered by Discuz! X3.2

© 2001-2013 Comsenz Inc.