adding r script for rna-seq analysis