R code for finding gene set