Question: WGCNA - Using real data for tutorial III
gravatar for hywel.williams
2.9 years ago by
hywel.williams0 wrote:

I have worked through the WGCNA tutorial I using my own data and I have used method 2b to construct my network and form modules and i have progressed to stage 3 where I've applied my traits and I have found some interesting findings.

I now want to look at hub genes from within those modules significantly associated with my trait of interest but when I move to the Tutorial III section I find that the input files bear no relationship to the files saved from the Tutorial I. Specifically, I am struggling to find a match from my data that corresponds to ModuleEigengeneNetwork1, GeneSummaryTutorial.csv, and TraitsTutorial.csv.


My aim is to get to Tutorial III part 7 and any guidance as to how I can get the data from Tutorial I to match the required input for Tutorial III would be greatly appreciated.


Many thanks - Hywel

ADD COMMENTlink modified 2.9 years ago • written 2.9 years ago by hywel.williams0
gravatar for Peter Langfelder
2.9 years ago by
United States
Peter Langfelder1.6k wrote:

I am a bit confused by your description. Are you working through tutorial I? If so, you should be looking at section 3, not Tutorial III, for an example of how to find important genes.

If you want to use Tutorial III, the best thing to do is to adapt the code in Section 2 ("Loading of expression data") so it loads and properly formats your own data.

For moduleEigengeneNetwork, use cbind() on the trait of interest and the module eigengenes (part of the output of blockwiseModules or you can use moduleEigengenes(datExpr, moduleColors) if you need to calculate them).

For the actual datExpr, you can use what you have from Tutorial I.

For traits, use a file containing your sample traits (conditions). The tutorial uses the variable 'y' to hold the trait of interest.

The tutorial is confusing (my fault) regarding GeneSummaryTutorial. This file is actually generated in Section 7, so we should not have been loading it in Section 2. The only thing it contains that's needed in the following is "truemodule" for which you can substitute your module colors (or labels). Just be aware that this is not really the 'true' module since in real data we don't really know 'true' modules.



ADD COMMENTlink written 2.9 years ago by Peter Langfelder1.6k
gravatar for hywel.williams
2.9 years ago by
hywel.williams0 wrote:

Dear Peter, that is very helpful, I will use your suggestions and repeat the analysis. Best wishes - Hywel

ADD COMMENTlink written 2.9 years ago by hywel.williams0
Please log in to add an answer.


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