I hope I solved this by calling DESeq function before the plotDispEsts in line 132 ans I also loaded the readr library for write_tsv which gave me an error
https://git.mpi-cbg.de/bioinfo/ngs_tools/blob/master/dge_workflow/featcounts_deseq.R#L132
Status changed to closed