Question: DESEQ2 application on control vs treatment
2
gravatar for bioinfonerd
13 months ago by
bioinfonerd30
bioinfonerd30 wrote:

enter code hereHi everyone,

I have 5 control samples and 5 treatment samples for the same patients and I would like to extract the differentially expressed genes between control and treatment samples.

The method I am using now is:

dds <- DESeqDataSetFromMatrix(countData = counts_set,colData = annotation,design = ~ condition)
dds <- DESeq(dds)
res <- results(dds,contrast=c("condition","Treatment","Vehicle"))

where condition represent control and treatment as follows:

        condition
01_C   Vehicle
02_C   Vehicle
03_C   Vehicle
04_C   Vehicle
05_C   Vehicle
01_T Treatment
02_T Treatment
03_T Treatment
04_T Treatment
05_T Treatment

Could you please let me know if this is a correct way to do?

Thanks

nanostring deseq • 602 views
ADD COMMENTlink modified 13 months ago • written 13 months ago by bioinfonerd30

I've split your tag into two keywords: nanostring and deseq. That will make them more relevant.

ADD REPLYlink written 13 months ago by RamRS24k

Thanks very much!!!!!

ADD REPLYlink written 13 months ago by bioinfonerd30
3
gravatar for h.mon
13 months ago by
h.mon27k
Brazil
h.mon27k wrote:

It seems correct, except for the fact you say you have 5 control and 5 treatment samples, whereas your condition variable has 10 control and 10 treatment samples.

ADD COMMENTlink modified 13 months ago • written 13 months ago by h.mon27k

Thanks!! I have corrected it.

ADD REPLYlink written 13 months ago by bioinfonerd30

The code does indeed look fine - good work. Please also take time to look through the Vignette at the part where it discusses the results object.

ADD REPLYlink written 8 months ago by Kevin Blighe48k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1214 users visited in the last hour