Question: DESeq2 - How to study time course experiment taking subject variability into account?
0
gravatar for Marion
22 months ago by
Marion0
Marion0 wrote:

Hello,

I have been reading lots of posts about DESeq2 and time course/repeated measures experiment but I still cannot find the answer to my question. The fact that measures are repeated on the same set of samples is rarely mentionned but it seems crucial to me. My design is similar to the one in the table below : (3 Groups) x (3 Time points) x (more than 2 subjects per Group):

Subject Time Group
1 W2 Control
1 W3 Control
1 W4 Control
2 W2 Control
2 W3 Control
2 W4 Control
3 W2 Treat_1
3 W3 Treat_1
3 W4 Treat_1
4 W2 Treat_1
4 W3 Treat_1
4 W4 Treat_1
5 W2 Treat_2
5 W3 Treat_2
5 W4 Treat_2
6 W2 Treat_2
6 W3 Treat_2
6 W4 Treat_2

I would like to identify the differentially expressed genes at (at least) 1 time point between any pairs of Group while taking the intra-subject variability into account. Is it correct to use the code below?

dds <- DESeqDataSet(data, ~ Subject + Time + Group + Group:Time)
ddsLRT <- DESeq(dds, test="LRT", reduced = ~ Subject + Time + Group)
resLRT <- results(ddsLRT)

If I understood correctly, I will get a unique p-value per gene for the signifiance of the global test (differentially expressed at (at least) 1 time point between at least one pair of Group while taking the intra-subject variability into account).

Would it then be correct to perfom kind of post-hoc test using

res <- results(dds)

and changing the contrasts to get all desired comparisons?

I hope that I am being clear, don't hesitate to ask otherwise as it is my first post here.

Best regards,

Marion.

ADD COMMENTlink modified 22 months ago by Michael Love26k • written 22 months ago by Marion0
Answer: DESeq2 - How to study time course experiment taking subject variability into acc
0
gravatar for Michael Love
22 months ago by
Michael Love26k
United States
Michael Love26k wrote:

Check this section of the vignette. For comparing the time course difference across treatment 1, 2 and control, you can use the following approach:

https://bioconductor.org/packages/release/bioc/vignettes/DESeq2/inst/doc/DESeq2.html#group-specific-condition-effects-individuals-nested-within-groups

ADD COMMENTlink written 22 months ago by Michael Love26k
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 16.09
Traffic: 174 users visited in the last hour