[BioC] paired time series data with two conditions

Simon Knott knott at usc.edu
Thu Nov 29 17:45:24 CET 2012


I have a short RNAseq time series dataset (two time pointsl T0 and TN)  
and I have two conditions Undrugged and Drugged. For each of the drug  
conditions the timepoints are paired (i.e one T0 corresponds to one  
T1). There are two T0/TN-paired replicates for each drug condition. So  
we have (Drugged_T0_A, Drugged_T1_A), (Drugged_T0_B, Drugged_T1_B),  
(Undrugged_T0_C, Undrugged_T1_C), (Undrugged_T0_D and Undrugged_T1_D).  
Is it possible to set up the test to identify genes that are altered  
more in the Drugged state (going from T0 to T1) than in the Undrugged  
state? I'm pretty sure this can be done with DESeqs glm functionality.

Essentially I want to find genes where (keeping the paired structure  
intact) log(Drugged_TN/Drugged_T0) << log(Undrugged_TN/Undrugged_T0) or
log(Drugged_TN/Drugged_T0) >> log(Undrugged_TN/Undrugged_T0).

Thanks So Much

Simon



More information about the Bioconductor mailing list