Soft-threshold in WGCNA
1
0
Entering edit mode
Arindam ▴ 80
@ag1805x-15215
Last seen 14 days ago
University of Eastern Finland

While choosing the soft-threshold, I observe that for my data (45 samples, case+control), the soft threshold does not reach R^2 ~ 0.8 for unsigned network + bicor but flattens around 0.7. Though a fit of 0.8 would have been definitely better, is it advisable to go with power at 0.7?

wgcna network bicor threshold • 2.6k views
ADD COMMENT
0
Entering edit mode
@peter-langfelder-4469
Last seen 14 days ago
United States

Please see WGCNA FAQ, point 6.

ADD COMMENT
0
Entering edit mode

I did came across that but still have some confusion.

the user should ensure that variables (probesets, genes etc.) have not been filtered by differential expression

Genes were only filtered to remove low counts.

If the scale-free topology fit index fails to reach values above 0.8 for reasonable powers ... chances are that the data exhibit a strong driver that makes a subset of the samples globally different from the rest.

The experimental design had batch effect which was accounted for. Hierarchical clustering after batch correction showed the sample grouping as desired based on the actual covariate.

the appropriate soft-thresholding power can be chosen based on the number of samples

This might be good when we at all do not reach any suitable power and for exploratory analysis. In my case, starting power = 14 the fit saturates around 0.7. Now, based on a sample size, if I choose power = 6, that means a much lower fit to scale-free topology. Which would be better in this case?

ADD REPLY

Login before adding your answer.

Traffic: 612 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6