R/create_RRHO_file_intersect_then_rank.R

In graeberlab/general:

#' create RRHO file after intersecting
#' 

#' @param metric_file1 filepath of first metric file. usually 2 columns. should have no headers. first column is gene name. second column is metric. does not need to be sorted
#' @param metric_file2 filepath of second metric file. usually 2 columns. should have no headers. first column is gene name. second column is metric. does not need to be sorted
#' @param outfile_prefix prefix of the output file
#' @param metric.col.file1 # of the column that metric data is in for file 1
#' @param metric.col.file2 # of the column that metric data is in for file 2
#' @param rank1 default is descending rank (i.e. higher metric is rank 1)
#' @param rank2 default is descending rank (i.e. higher metric is rank 1)
#' 
#' @export
#' 

create_RRHO_file_intersect_then_rank<-function(metric_file1,metric_file2,outfile.prefix=paste0("rankrank.",metric_file1,"_vs_",metric_file2),metric.col.file1=2,metric.col.file2=2,rank1="desc",rank2="desc") {
  
  rankcol1= metric.col.file1-1
  rankcol2= metric.col.file2-1
  
  metric1 <- read.table(metric_file1, row.names=NULL,check.names=F,header=T,stringsAsFactors=F,sep="\t")
  metric2 <- read.table(metric_file2, row.names=NULL,check.names=F,header=T,stringsAsFactors=F,sep="\t")
  metric1 <- metric1[!duplicated(metric1[,1]),]
  metric2 <- metric2[!duplicated(metric2[,1]),]
  rownames(metric1)<-make.names(metric1[,1])
  rownames(metric2)<-make.names(metric2[,1])
  metric1 <- metric1[,-1,drop=FALSE]
  metric2 <- metric2[,-1,drop=FALSE]
  
  metric1 <- metric1[,rankcol1,drop=FALSE]
  metric2 <- metric2[,rankcol2,drop=FALSE]
  rownames(metric1) <-toupper(rownames(metric1))
  rownames(metric2)<-toupper(rownames(metric2))
  
  metric1.names = rownames(metric1)
  metric2.names = rownames(metric2)
  
  common_names<-intersect(metric1.names,metric2.names)
  print(paste0("File1 genes=",length(metric1.names)," File2 genes=",length(metric2.names)," Intersect genes=",length(common_names)))
  metric1.positions<- match(common_names,metric1.names)
  metric2.positions<- match(common_names,metric2.names)
  
  #matched 
  mm1 <- metric1[metric1.positions,,drop=FALSE]
  mm2 <- metric2[metric2.positions,,drop=FALSE] 
  
  #matched and ordered in descending order(i.e row 1 is highest value)
  if(rank1=="desc"){
    mmo1 <- mm1[order(-mm1[,1]),,drop=FALSE]
  } else {
    mmo1 <- mm1[order(mm1[,1]),,drop=FALSE]
  }
  if(rank2=="desc"){
    mmo2 <- mm2[order(-mm2[,1]),,drop=FALSE]
  } else {
    mmo2 <- mm2[order(mm2[,1]),,drop=FALSE]
  }
  
  
  #matched,ordered,and ranked
  mmr1 <- transform(mmo1, rank = 1:nrow(mmo1))
  mmr2 <- transform(mmo2, rank = 1:nrow(mmo2))
  
  #now use ranks from set 1 as reference
  ord<- match(rownames(mmr1),rownames(mmr2)) 
  mmr2<-mmr2[ord,,drop=FALSE]
  
  
  both <- data.frame(rownames(mmr1),rownames(mmr2),mmr1[,2],mmr2[,2],mmr1[,1],mmr2[,1]) 
  colnames(both) <- c("Unigene","Gene_Symbol",paste0(metric_file1,".rank"),paste0(metric_file2,".rank"),paste0(metric_file1,".metric"),paste0(metric_file2,".metric"))
  write.table(both, file = paste0(outfile.prefix, ".txt"), col.name=T,row.names=F, sep="\t", quote=F)
}