diff options
Diffstat (limited to 'scripts/wgcna_analysis.R')
| -rw-r--r-- | scripts/wgcna_analysis.R | 10 |
1 files changed, 9 insertions, 1 deletions
diff --git a/scripts/wgcna_analysis.R b/scripts/wgcna_analysis.R index ee749e9..c3b1ac8 100644 --- a/scripts/wgcna_analysis.R +++ b/scripts/wgcna_analysis.R @@ -62,13 +62,21 @@ enableWGCNAThreads() powers <- c(c(1:10), seq(from = 12, to=20, by=2)) sft <- pickSoftThreshold(dataExpr, powerVector = powers, verbose = 5) +# check the power estimate + +if (is.na(sft$powerEstimate)){ + powerEst = 1 +}else{ + powerEst = sft$powerEstimate +} + # pass user options network <- blockwiseModules(dataExpr, #similarity matrix options corType = corType, #adjacency matrix options - power = 5, + power = powerEst, networkType = "unsigned", #TOM options TOMtype = TOMtype, |