#creating a rank file for GSEA
setwd("C:/Users/thood/Documents/R/TCGA")
x<-read.csv("res001_recurrent_fusions_w_gene_names_9691.csv")
x$fcSign=sign(x$log2FoldChange)
x$logP=-log10(x$padj)
x$metric=x$logP/x$fcSign

y<-x[,c("hgnc_symbol", "metric")]
y <- y[order(-y$metric),]
write.table(y,file="expression.rnk",quote=F,sep="\t",row.names=F)