[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