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() ```