Deleted:Different summary() and results() values in DeSeq2
2
0
Entering edit mode
Mafer • 0
@e7ba24a7
Last seen 13 months ago
Germany

I have run DeSeq2 in 2 different devices using the same count data and metadata tables, however when running the summary and results I get different values.

I made sure the files are read correctly, so decimals are separated by (.) in both devices and also after round the values I get the same count table values in both,


# sampleTable #
csvfile <- "metadata.csv"
sampleTable <- read.csv(csvfile, row.names=1, sep=",")
sampleTable$Stage <- as.factor( sampleTable$Stage )
summary(sampleTable)

# data matrix #
count_data_file <- "tpm.csv"
countdata <- read.csv(count_data_file, row.names=1, header=T, sep="\t",  dec ="."))
summary(countdata)


# construction of DESeqDataSet #
ddsMat <- DESeqDataSetFromMatrix( countData=round(countdata), colData=sampleTable, design= ~ Stage ) 
nrow(ddsMat)


# removal of not or low expressed genes #
dds <- ddsMat[ rowSums(counts(ddsMat)) > 1, ]
nrow(dds)

# differential expression analysis #

dds <- DESeq(dds)
res <- results(dds)
summary(res)

sessionInfo( )

These are the results obtained in each device, both running R v. 4.1.3

enter image description here

#deseq2 • 445 views
ADD COMMENT
This thread is not open. No new answers may be added
Traffic: 787 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