Updates to script. Added gene set enrichment analysis using topGO package. Changed RxNseq to also output intercept. For some reason, function is **much** slower than before...should optimize at some point