Make DESeq2 or edgeR QC plots
diffQCplots(method, fit, y)
method
one of "DESeq2" or "edgeR"
fit
output of fitDiffTSS (if method = "edgeR")
y
output of fitDiffTSS (if method = "DESeq2")