Question: DESEQ2 application on control vs treatment
2
gravatar for bioinfonerd
22 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 • 1.1k views
ADD COMMENTlink modified 22 months ago • written 22 months ago by bioinfonerd30

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

ADD REPLYlink written 22 months ago by RamRS27k

Thanks very much!!!!!

ADD REPLYlink written 22 months ago by bioinfonerd30
3
gravatar for h.mon
22 months ago by
h.mon29k
Brazil
h.mon29k 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 22 months ago • written 22 months ago by h.mon29k

Thanks!! I have corrected it.

ADD REPLYlink written 22 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 16 months ago by Kevin Blighe60k
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: 1944 users visited in the last hour