R commands
# Loading the library
library("DESeq2")
# Reading and preparing data
c<-read.table("m.txt",head=TRUE,row.names = 1)
c<-round(c)
colData <-read.table("colData.txt",head=TRUE,row.names=1)
# Running DESEQ
dds<-DESeqDataSetFromMatrix(countData=c,colData=colData,design= ~ condition)
dds <-DESeq(dds)
res <- results(dds)
resOrdered <-res[order(res$pvalue),]
write.csv(as.data.frame(resOrdered),file="set0.csv")
# MA plot
plotMA(res,main="MF",ylim=c(-5,5))
# Printing differentially expressed genes
head(resOrdered)
resSig <- subset(resOrdered,padj<0.1)
References
Manuals
https://www.bioconductor.org/packages/release/bioc/vignettes/DESeq2/inst/doc/DESeq2.pdf https://www.bioconductor.org/packages/devel/bioc/vignettes/DESeq/inst/doc/DESeq.pdf
Papers
https://genomebiology.biomedcentral.com/articles/10.1186/s13059-014-0550-8 https://genomebiology.biomedcentral.com/articles/10.1186/gb-2010-11-10-r106