Question: Module- Triat relationship in WGCNA
0
gravatar for Prakash
7 months ago by
Prakash730
India
Prakash730 wrote:

Hello all,

I am working with 20 samples of RNAseq data ( 5 control and 5 KD with five different stimulation) and have identified several module using WGCNA. Next I would like to perform module- condition relationship analysis, basically i would like to identify condition specific module. I have designed the condition in binary format (see table below). can someone tell that the design is correct or not, if not, what would be best way to design trait in this scenario.

sample  EU  EI  EC  EP  EA  NU  NI  NC  NP  NA
E_Uns_R1    1   0   0   0   0   0   0   0   0   0
E_Uns_R2    1   0   0   0   0   0   0   0   0   0
E_IFNg_R1   0   1   0   0   0   0   0   0   0   0
E_IFNg_R2   0   1   0   0   0   0   0   0   0   0
E_CpG_R1    0   0   1   0   0   0   0   0   0   0
E_CpG_R2    0   0   1   0   0   0   0   0   0   0
E_pIC_R1    0   0   0   1   0   0   0   0   0   0
E_pIC_R2    0   0   0   1   0   0   0   0   0   0
E_all3_R1   0   0   0   0   1   0   0   0   0   0
E_all3_R2   0   0   0   0   1   0   0   0   0   0
N_Uns_R1    0   0   0   0   0   1   0   0   0   0
N_Uns_R2    0   0   0   0   0   1   0   0   0   0
N_IFNg_R1   0   0   0   0   0   0   1   0   0   0
N_IFNg_R2   0   0   0   0   0   0   1   0   0   0
N_CpG_R1    0   0   0   0   0   0   0   1   0   0
N_CpG_R2    0   0   0   0   0   0   0   1   0   0
N_pIC_R1    0   0   0   0   0   0   0   0   1   0
N_pIC_R2    0   0   0   0   0   0   0   0   1   0
N_all3_R1   0   0   0   0   0   0   0   0   0   1
N_all3_R2   0   0   0   0   0   0   0   0   0   1

Thanks

rna-seq forum • 417 views
ADD COMMENTlink modified 7 months ago by andrew.j.skelton735.5k • written 7 months ago by Prakash730
0
gravatar for Lluís R.
7 months ago by
Lluís R.810
Spain, Barcelona
Lluís R.810 wrote:

Your conditions seems like blocking factor. In the WGCNA it is recommended to make the correlation between the eigenvalue of the module with the phenotype variables. I would recommend however to do the correlation between the eigenvalue of the module with the eigenvalue of your condition.

However, if this is are really 5 different conditions it might be worthless, as the co-expressed genes might be few or not highly co-expressed...

ADD COMMENTlink written 7 months ago by Lluís R.810

but in this paper they have take three condiion in the last part of their analysis and did WGCNA can you explain ,how did they perform ?

ADD REPLYlink written 7 months ago by krushnach80460
2

Well I don't know anything from this paper beyond what it is described, but looking for the Figure 4 a. it seems like they made the correlation between the condition and the eigengene. If it is accurate or not is another thing... You can do as them, but note that they have more than 4 samples for each condition and up to 9.

Upon closer look to your matrix, it seems like you have several factors happening at the same time. If I understood correctly you have 2 conditions (N and E) and then (Uns, IFNg, pIC, CpG and all3) with two replicates for the second condition, which results in smaller groups for each condition, which are too small to do something with confidence.

ADD REPLYlink written 7 months ago by Lluís R.810

yes the sample size is a concern ..

ADD REPLYlink written 7 months ago by krushnach80460
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: 819 users visited in the last hour