Time course RNA-seq data with replicates
1
0
Entering edit mode
@abhishek-singh-4725
Last seen 5 weeks ago
France

Hi,

I have RNA seq data for 7 time points with 3 replicates each, I would like to mark genes that are differentially expressed at least at one of the time points.

The DESeq manual details about dealing with time series data that have two conditions.

In general I would like to have a figure like this: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3462286/figure/F1/                        (figure B).

 

Kindly help.

 

Thank you

 

DEseq2 • 793 views
ADD COMMENT
1
Entering edit mode
@mikelove
Last seen 10 hours ago
United States

Have you seen the time course example in our workflow?

http://www.bioconductor.org/help/workflows/rnaseqGene/#time

Can you say more about your experimental design? Do you have a treated and control group or just the time series?

best,

Mike

ADD COMMENT
0
Entering edit mode

Hi,

I just have the time series, it is a developmental data set with three replicates for each stage.

Thank you for the link.

 

ADD REPLY
1
Entering edit mode

You can use a design of ~time, where time should be a factor, not a numeric variable:

class(dds$time)

And then you can use the following test:

dds <- DESeq(dds, reduced=~1, test="LRT")

res <- results(dds)

This will produce p-values testing if the mean level is constant across all times (the null hypothesis). See ?results and the section in the vignette about the LRT for more details.

ADD REPLY

Login before adding your answer.

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