Generate and save dispersion plot only once per dataset, now its done for each comparison:
# save dispersion plot
output_fname = paste(input_folder, "/results/edgeR/dispersion/", feature_type, "_dispersion.png", sep="")
png(output_fname)
plotBCV(dge)
dev.off()
Generate and save dispersion plot only once per dataset, now its done for each comparison: