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

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

ADD REPLYlink written 8 months ago by RamRS20k

Thanks very much!!!!!

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

Thanks!! I have corrected it.

ADD REPLYlink written 8 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 10 weeks ago by Kevin Blighe39k
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: 1139 users visited in the last hour