This is a test version of Biostars. For the public version, visit https://www.biostars.org.
WCGNA empty files for GO analysis

Hello all. I would appreciate your help with this. Im doing an WCGNA, everything is getting well but the only thing that I have a trouble is when I tried to generate a file with all my genes from each module (to perform a gene ontology analysis)

I've performed this command, it doesnt tell me there is an error, but when I open my .txt files those are all empty.

What would you could suggest for me?

Thank you so much

allLLIDs = pulmonData$gene_symbol 
intModules = unique(moduleColorsAutomatic) 
for (module in intModules)
{
  modGenes = (moduleColorsAutomatic==module) 
  modLLIDs = allLLIDs[modGenes];
  fileName = paste("LocusLinkIDs-", module, ".txt", sep="");
  cat(modLLIDs, sep="\n", file=fileName)
}

TotfileName = paste("LocusLinkIDs-all.txt", sep="");
cat(allLLIDs, sep="\n", file=TotfileName)
wgcna rna-seq
• 1,153 views • link
updated by Ram Ram • written by ovariohisterectomia • ovariohisterectomia

1 answer

I've already answered it.

The problem was that the pulmonData first column was named Gene ID, I changed to gene_symbol and now it works.

Thank you!

• 0 views • link
updated by Ram Ram • written by ovariohisterectomia • ovariohisterectomia

Log in to answer this question.

More posts like this