Entering edit mode
12 months ago
annaA
▴
10
Hello I need your help with my analysis.
So i have an bulk Rna-seq dataset which contains data from ips cells which develop into neurons in 5 days For each time point I have 7 replicates
So now with deseq2 I have made all the pair-wise comparisons and I want to continue by comparing each time point to all the others (e.g time point 1 vs all the other time points) in order to find specific genes for each developmental stage.
I have seen several ways online to do so but I am a bit confused.
Any ideas how t build these kind of contrasts ?
Thanks in advance, A
What confuses you, what did you try? Please understand that "it does not work"-like questions are hard to debug. Did you read the DESeq2 manual on contrasts?
So my code until now is as follows
my reference is 0_dpi
dds <- DESeq2::DESeq(deseq2Data)
pairwise comparisons
After I am filtering the results of each comparison with the following thresholds
Then I tried for example to compare 0_dpi to all otehr time points by this
this didnt work
Also I tried the LRT method with reduced model like
My question is which is the right way to go?
As said, this is not informative.
What
did not work, why, which errors, what is the problem?