Question: WGCNA soft-threshold power selection
0
gravatar for mahnazkiani
9 days ago by
mahnazkiani50
United States
mahnazkiani50 wrote:

Hello,

I have a question about the results of WGCNA for selecting soft-thresholding power, all power numbers are below the line in the scale independence graph, and I don't know which power I need to select for next step.I couldn't insert the picture to be more clear, but code that I used was as follows:

# choose soft-thresholding power
powers = c(c(1:10), seq(from = 12, to=20, by=2))
sft = pickSoftThreshold(datExpr, powerVector = powers, verbose = 5)
sizeGrWindow(9, 5)
par(mfrow = c(1,2))
cex1 = 0.9
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")
abline(h=0.9,col="red")

plot(sft$fitIndices[,1], sft$fitIndices[,5],
     xlab="Soft Threshold (power)",ylab="Mean Connectivity", type="n",
     main = paste("Mean connectivity"))
text(sft$fitIndices[,1], sft$fitIndices[,5], labels=powers, cex=cex1,col="red")

for (i in 1:dim(datExpr)[2]) {
  datExpr[,i] <- c(as.double(datExpr[[i]]))
}
rna-seq • 141 views
ADD COMMENTlink modified 9 days ago • written 9 days ago by mahnazkiani50

Where is your graph?

ADD REPLYlink written 9 days ago by Kevin Blighe45k

soft-T

ADD REPLYlink modified 9 days ago by genomax69k • written 9 days ago by mahnazkiani50

Choose 6, for now. I am also going to hypothesise that your dataset is rather 'flat'. If you plot a histogram (hist()), how does it appear? Which type of input data re you using? - please elaborate on the processing steps.

ADD REPLYlink modified 9 days ago • written 9 days ago by Kevin Blighe45k

I am using log2(x+1) normalized total counts RNa-seq data.

ADD REPLYlink written 6 days ago by mahnazkiani50

This is covered in the FAQ of the WCGNA webpage https://horvath.genetics.ucla.edu/html/CoexpressionNetwork/Rpackages/WGCNA/faq.html See #6: "I can't get a good scale-free topology index no matter how high I set the soft-thresholding power."

It kind of depends on your data and how many samples you have

ADD REPLYlink written 7 days ago by rmash0

mahnazkiani : How to add images to a Biostars post

ADD REPLYlink written 9 days ago by genomax69k

Thank you for sharing the link

ADD REPLYlink written 9 days ago by mahnazkiani50
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 949 users visited in the last hour