Question: Correlating WGCNA modules to Diseased/Not-Diseased Trait
0
gravatar for RNAseqer
10 months ago by
RNAseqer 130
RNAseqer 130 wrote:

Hello all,

I have a dataset I am performing WGCNA analysis using the . The tutorial materials on Steve Hovarth's website are excellent but one question I have not found an answer to is this: When you have a trait that is not a continuous variable but rather a binary state such as diseased/not diseased and you are trying to correlate the modules to it, how should you encode the trait?

Should you assign a 0 and a 1 for non-diseased and diseased respectively? What about a 1 and a 2? Would both work equally well?

When so encoded how would the step where you correlate modules to traits work?

If you reversed this coding so that non-diseased was encoded as 2, and diseased was coded as 1 would you expect modules that negatively correlate to being diseased show up as a positive correlation?

Also, what if you had three conditions: Control, Traumatized, Diseased. Would encoding them as 1,2,3 make sense? Or would you want to limit your test to just those samples that were traumatized and diseased, if you were primarily interested in traum vs. diseased comparison.

I'd appreciate any advice you can offer!

modules wgcna traits • 371 views
ADD COMMENTlink modified 10 months ago by andres.firrincieli600 • written 10 months ago by RNAseqer 130
2
gravatar for andres.firrincieli
10 months ago by
andres.firrincieli600 wrote:

Should you assign a 0 and a 1 for non-diseased and diseased respectively? What about a 1 and a 2? Would both work equally well?

Binary traits should contain only 1 and 0

what if you had three conditions: Control, Traumatized, Diseased. Would encoding them as 1,2,3 make sense?

if you have three conditions is: 1 vs All (0)

When so encoded how would the step where you correlate modules to traits work?

Look at the section 3.a of this tutorial; the object datTraits is the binary matrix

ADD COMMENTlink written 10 months ago by andres.firrincieli600

For calculating the correlation p-value, the tutorial uses Student's t-test [corPvalueStudent()]. Should this be modified while using a categorical/binary variable?

ADD REPLYlink written 5 weeks ago by Arindam Ghosh280
1

Hi Arindam,

you can use the corPvalueStudent function

ADD REPLYlink written 5 weeks ago by andres.firrincieli600
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: 984 users visited in the last hour