I know if there is one control and one treatment group it is pretty straight forward to interpret log 2 fold change. But, I have time course experiment. I have infected cells with viruses and I harvested cells at different time points: o hr, 6 hour, 12 hr and 24 hour. Here o hr is my non-infected control. What is the best way to show log 2 fold change for each gene in this type of time course experiment. I ultimate aim to find highly expressed genes at different time points.
1 answer
You should analyze your data with a proper statistical framework such as DESeq2. The vignette offers a section about time series experiments plus there are plenty of questions about it on the Bioconductor support forum, just use google to find them. It might be worth looking at the shrinkage function in DESeq2 to get more appropriate log2 fold changes for lowly-expressed genes. I hope you have replicates. If not, please use google and the search function for posts on RNA-seq analysis without replicates. THere is a lot of material on this available both here on Biostars as well as on the BioC support forum and the web.
Log in to answer this question.
-
How to approach DESeq2 analysis on 2 plant lines infected with virus, with 2 time points (days post…
written by duplessisantzel •
-
DESeq2 contrasts and design formula for 2 conditions and 3 exposure times
written by pcar •
-
DEseq2: comparison the same when reversing factor order?
written by gt •
-
Comparing gene set enrichment across 3+ conditions?
written by penny.lane •
-
Deal with the time-course RNA-Seq datasets
written by wangdp123
-
What is the ideal FDR and fold change cutoff to consider in a RNA-seq experiment and which one shou…
written by aishu.jp •
-
Limma time-course: how to apply differential contrast if time zero is considered both as treatment …
written by harelarik
-
Best way to address different batches of RNA-seq
written by tud55122 •
-
Clustering Multiseries Data R
written by boczniak767
-
Heatmap Of Rna-Seq Data
written by Ying W
