Question: Module- Triat relationship in WGCNA
0
gravatar for Prakash
21 months ago by
Prakash1.8k
India
Prakash1.8k 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 • 714 views
ADD COMMENTlink modified 21 months ago by andrew.j.skelton735.9k • written 21 months ago by Prakash1.8k
0
gravatar for Lluís R.
21 months ago by
Lluís R.890
Spain, Barcelona
Lluís R.890 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 21 months ago by Lluís R.890

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 21 months ago by krushnach80690
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 21 months ago by Lluís R.890

yes the sample size is a concern ..

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