Design in DESEQ2
0
0
Entering edit mode
5 months ago
hydrobeanz • 0

Greetings, I have RNA-seq data need to be analyzed, How to do? I have two groups, one disease one no disease. Within each group, they received treatments of treatment A, B, C, D. This is all in one column in MetaData table (DiseaseTreatment), like NoDiseaseA, DiseaseA, NoDiseaseC, DiseaseC etc etc. How do I compared within DESEQ2 design? I want to compare NoDiseaseA-D to DiseaseA-C, and NoDisease overall to Disease overall.

dds <- DESeqDataSetFromMatrix(countData=mydata, colData=MetaData, design=~DiseaseTreatment)

then do: compare1 <- results(dds, contrast=c("DiseaseTreatment","NoDiseaseA", "DiseaseA"), independentFiltering=TRUE, alpha=0.05)

Is that right?

Many thanks.

RNA-Seq R • 210 views
ADD COMMENT
1
Entering edit mode

I do not mean to sound harsh, but please go through the DESeq2 vignette:

Biostar

ADD REPLY

Login before adding your answer.

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