Entering edit mode
4.4 years ago
ovariohisterectomia
▴
40
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)