Hi,
I have been following the DTU workflow which is linked to the package rnaseqDTU and the paper: "Swimming downstream: statistical analysis of differential transcript usage following Salmon" https://f1000research.com/articles/7-952/v3
It is mentioned that, in favor of controlling the FDR of DRIMSeq one should apply a post-Hoc filter to exclude "transcripts with small SD of the per-sample proportions". This was not suggested to be done for results of DEXSeq and its workflow which is described in the same paper. The reason is that DEXSeq seemed to "be closer to controlling its FDR".
I wonder if this post-Hoc filter would be redundant if one is only interested in the intersection of DRIMSeq and DEXSeq results. As those FPs that were called by DRIMSeq would in the best case not be present in DEXSeq. (Having in mind that the intersection of results from both tools might lack some DTUs that were found by DRIMseq being more sensitve).
Thanks in advance
Fiona
Thanks for the answer and sorry for late answer. Is it ok to perform post-hoc filter on lfc transcripts after testing and reporting the unchanged pvalue? Is it ok to use
DEXSeq::estimateExonFoldChange
andDRIMSeq::coefficients
for this purpose if Id want to do the same post hoc filtering on both tool results?