Entering edit mode
5.7 years ago
hirahameed222
•
0
I have generated the csv file by using following command
write.csv(results_genes, "chrX_gene_results.csv",
row.names=FALSE)
my output file lacks particular gene names in gene column the file appears as :
gene MSTRG.9464 1 0.456607798602799 0.46264129580405
gene MSTRG.61457 1 0.456608200007355 0.46264129580405
gene MSTRG.27219 0.999999999999999 0.456608285295604 0.46264129580405
Can some one tell me that how can I find and compose particular gene names in gene in column? is there any script ? or some thing is missing in my data which fails to fetch the gene name?
There is no need to SHOUT. I have removed the capital letters from your title.
It is furthermore unclear how you obtained "chrX_gene_results.csv", which is probably a crucial piece of information here.
My Bad i haven't checked that i have All-Caps in my post. I'll care for next time. I am following this protocol for differential expression analysis, I followed the same but my results have missing gene names.
Did you furnish a reference gtf file to stringtie/cufflinks pipeline? MSRTG are temporary names given by the workflow.
yes, I followed the protocol exactly.
What's the result if, in R, you type:
From your output, it seems you have reference gtf and your analysis included finding transcripts matching reference and also the novel transcripts. Wherever tool found transcripts, not in reference gtf, tool assigned MSRTG accessions to the transcripts.