WGCNA Trait File Issues

WGCNA Trait File Issues

1

I have RNASeq data that I am trying to use for WGCNA. However I am running into downstream issues that I believe stem from my trait file. I have 24 samples, 3 replicates of 8 different treatment groups. I am trying to follow the WGCNA tutorial but having to make minor adjustments because my treatment groups are qualitative, not quantitative. I did make each treatment group into a number, 1-8, instead of the treatment name. My input data file is “mat2” and my trait file is “Condition”. Here is my code for my trait file upload:

CData = read.csv(“~/Documents/Conditiondata.csv”)
dim(CData)
names(CData)
CData2 = CData[, -c(1)] #removing the first column
names(CData2)
dim(CData2)

Samples = rownames(mat2)
conditionRows = match(Samples, CData2$Sample)
Condition= CData2[conditionRows, ]
rownames(Condition) = CData2[conditionRows, 1]
collectGarbage()
rownames(Condition) == rownames(mat2) #confirming that row names match, and this comes back as TRUE.

Once I get to the point of making the module-trait heatmap, I end up with coerced NAs for the samples. I’ve attached an image of what my resulting heatmap. Any help is greatly appreciated! Thank you!


correlation


Deseq2


Gene


WGCNA

• 35 views


Login
before adding your answer.

Traffic: 1838 users visited in the last hour

Read more here: Source link