R scripts for running DESeq analysis