DESeq2 nomalize value
1
0
Entering edit mode
史哲 • 0
@b0736200
Last seen 10 weeks ago
United States

Hello, guys, I have a problem getting the DEseq2 normalized value. I am using R4.3.2 and DESeq2 version 1.42.0. I have 32 samples for 4 groups. The raw count data is obtained from Hisat2-string tie pipeline. I would like to get DESeq2 normalized value, but only A1 sample value did not change from the raw count. Could anyone tell me how to resolve this? Do I make the wrong data. frame?

Thanks, Shitetsu

My code is below.

x<-read.table("Book14.txt",sep = "\t", header = T, row.names = 1, quote="")
is.data.frame(x)
x <- as.matrix(x)
coldata <- data.frame(sample = c("A1","A2","A3","A4","A5","A6","A7","A8","B1","B2","B3","B4","B5","B6","B7","B8","C1","C2","C3","C4","C5","C6","C7","C8","D1","D2","D3","D4","D5","D6","D7","D8"), condition = c("A","A","A","A","A","A","A","A","B","B","B","B","B","B","B","B","C","C","C","C","C","C","C","C","D","D","D","D","D","D","D","D"), row.names = "sample")
dds <- DESeqDataSetFromMatrix(countData = x, colData = coldata, design = ~ condition)
dds <- estimateSizeFactors(dds)
NC <- counts(dds, normalized=TRUE)
head(NC)

And the this is the capture of x and NC enter image description here enter image description here

DESeq2 • 337 views
ADD COMMENT
0
Entering edit mode
ATpoint ★ 4.0k
@atpoint-13662
Last seen 1 day ago
Germany

What is the output of sizeFactors(dds)?

ADD COMMENT
0
Entering edit mode

Hello, ATpoint, Thank you for your response. The output is this. enter image description here

ADD REPLY
0
Entering edit mode

There is your answer. The first sample has a size factor of exactly one, so its values stay exactly the same and the others are scaled to it. It's fine.

ADD REPLY
0
Entering edit mode

Thank you for your advice! I appreciate you!

ADD REPLY

Login before adding your answer.

Traffic: 488 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