How to find DEG with deseq2 in R
0
0
Entering edit mode
16 months ago

I have downloaded a dataset(GSE98622) from GO databases. This data is best suited for my experiences, and I need to extract DEG from it, but this data is normalized and I can't use it. What should I do in this situation?

AKI <- read.table("GSE98622a.txt", sep = "\t")
x <-as.matrix(AKI)
rownames(x)<-x[,1]
x <- x[,-1]
AKI <- as.data.frame(x)
class(AKI)
gr <- c(rep("n", 17), rep("C", 16), rep("d", 16))
colData <- data.frame(group=gr, type="paired-end")
cds <- DESeqDataSetFromMatrix(data.matrix(AKI), colData, design = ~group)
resLFC.AKI <- data.frame(lfcShrink(data.matrix(AKI), coef="group_d_vs_n"))

Many thanks

GSE RNA-seq deseq2 • 667 views
ADD COMMENT
0
Entering edit mode

Did you get a matrix of already normalized read counts? What type of normalization has been used?

ADD REPLY
0
Entering edit mode

Matrix array but also we need numeric data type

The error is:

> cds <- DESeqDataSetFromMatrix(data.matrix(AKI), colData, design = ~group)
Warning message:
In DESeqDataSet(se, design = design, ignoreRank) :
  some variables in design formula are characters, converting to factors
> resLFC.AKI <- data.frame(lfcShrink(data.matrix(AKI), coef="group_d_vs_n"))
Error in lfcShrink(data.matrix(AKI), coef = "group_d_vs_n") : 
  is(dds, "DESeqDataSet") is not TRUE
ADD REPLY
1
Entering edit mode

Read the DESeq2 manual. Your entire code is wrong, manual will instruct you.

ADD REPLY

Login before adding your answer.

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