no gene names in deseq2 result file
1
0
Entering edit mode
Jitendra ▴ 10
@nabiyogesh-11718
Last seen 13 hours ago
United Kingdom
library('DESeq2')
countMatrix = read.table("RPR_count.txt",header=T,sep='\t',check.names=F)
dim(countMatrix)

> head (countMatrix)
                  HRPR_0D_V8_R1 LRPR_0D_V8_R1 HRPR_0D_V8_R2 LRPR_0D_V8_R2
1 ENSRNA049455608            27            15            13            20
2 ENSRNA049455955            50            31            66            20
3 ENSRNA049457359           149           325           146           203
4 ENSRNA049457570            52            25            30            22
5 ENSRNA049458758            16            19            44            32
6 ENSRNA049458986             6            14            17            20

colData <- read.table("Metadata_RPR.txt", check.names=F)
dim(colData)

# making first column as rownames in coldata
SampleInfo <- colData[,-1] 
rownames(SampleInfo) <- colData[,1] 

dds <- DESeqDataSetFromMatrix(countData = countMatrix, colData = SampleInfo, design = ~ geno_stage)

but in result file there is no gene names, could you please suggest how I can resolve it?

"","baseMean","log2FoldChange","lfcSE","stat","pvalue","padj"
"1",2133.39832753799,-1.48791540218224,0.133582135187263,-11.1385807697743,8.14062484586051e-29,2.12348199104271e-24
"2",518.774465513072,22.415307417241,2.21105823446886,10.1378186552493,3.75388225192356e-24,4.8960009270713e-20

Thanks

deseq2 • 1.9k views
ADD COMMENT
0
Entering edit mode

no, it is as column in countMatrix, should use these command to convert the first column as rownames

# making first column as rownames in countMatrix
count <- countMatrix[,-1] 
rownames(count) <- countMatrix[,1]
ADD REPLY
0
Entering edit mode
@mikelove
Last seen 8 hours ago
United States

DESeq2 will have gene names on the rownames of the results table if gene names are provided as input. Are you providing a counts matrix with gene names on the rows?

ADD COMMENT

Login before adding your answer.

Traffic: 533 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6