Dseq2 gives the same results no matter the condition ???

0

Entering edit mode

Bioinformatics ▴ 30

@bioinformatics-10931

Last seen 3.9 years ago

United States

Hello,

I am trying to understand why Dseq2 gives me the same result if I set my condition differently ? I use Raw count values

First way to do

condition <- factor(c(rep("Normal", 10),rep("Tumor", 11)))
dds <- DESeqDataSetFromMatrix(dfm, DataFrame(condition), ~ condition)
keep <- rowSums(counts(dds)) >= 1
dds <- dds[keep,]
dds <- DESeq(dds)

Second way to do

condition <- factor(c(rep("Tumor", 11),rep("Normal", 10)))
dds <- DESeqDataSetFromMatrix(dfm, DataFrame(condition), ~ condition)
keep <- rowSums(counts(dds)) >= 1
dds <- dds[keep,]
dds <- DESeq(dds)

They both give me the same LogFold change???? how is that possible ?

dseq2 • 1.2k views

ADD COMMENT • link 7.5 years ago Bioinformatics ▴ 30

0

Entering edit mode

Probably because it is reordering it internally alphabetically or something like that. You could check the source code.

ADD REPLY • link 7.5 years ago chris86 ▴ 420

Login before adding your answer.