I am analysing data which I have mapped to the reference genome using HISAT2. I estimated the expression using Stringtie.
Now I want to run DESeq2 on this data. Works well so far. But I am not sure how I have to "model" my experimental design.
I have:
- 5 x biol. replicates "Control (no treatment)"
- 5 x biol. replicates "Treatment low dosage (chemical compound)"
- 5 x biol. replicates "Treatment high dosage (chemical compound)"
Do I have to compare control<>low dosage and control<>high dosage individually or is there a better way to model the "gradient" of treatment using DESEq2?
1 answer
You can do a dose response model, though you might want to play around with the number values you use to represent zero, high and low treatment.
The results for this will say they are returning a log fold change, but they are not; they will instead be returning the slope of the line, and the p-value is how well the points line up on the line.
Log in to answer this question.
-
RNA-seq differential expression across multiple tissues, 2 trials, and 2 conditions
written by BioinfGuru
-
Linear models with limma: coefficients not estimatable, are the others OK?
written by fr
-
dUTP sequencing -- different strand-ratios for different samples
written by mschmid
-
Questions related to ChromHMM+using MACS2 peak calling output bed file as input for ChromHMM?
written by m.sadman.sakib
-
Comparison of samples of different read lengths in differential gene expression analysis
written by rthapa •
-
Majority of genes upregulated in RNA-seq study of toxic compound
written by mschmid
-
stringtie + ballgown for multivariate analysis of rna-seq data
written by nancy
-
ballgown for diff treatments+replicates+diff time points
written by bioinfo17 •
-
DESeq2 dataset setting with count table.
written by jhkim1972 •
-
DESeq2 with 2 factors: problem to fit the negative binomial and weird p-value distributions
written by Aurelie MLB
