R/transformBounds.R

In CARNIVAL: A CAusal Reasoning tool for Network Identification (from gene expression data) using Integer VALue programming

## Vectorizing the bounds
##
## Enio Gjerga, 2020

transformBounds <- function(mt = mt, lpFile = lpFile){
  
  idx1 <- which(lpFile$`enter Problem`=="Bounds")
  idx2 <- which(lpFile$`enter Problem`=="Binaries")
  
  constraintSet <- lpFile$`enter Problem`[seq(from = idx1+1, 
                                              to = idx2-1, by = 1)]
  
  f.con <- matrix(data = 0, nrow = 2*length(constraintSet), ncol = nrow(mt))
  ## f.dir <- c()
  ## f.rhs <- c()
  f.dir <- rep("", length(constraintSet)*2)
  f.rhs <- rep("", length(constraintSet)*2)
  
  cnt <- 1
  for(ii in seq_len(length(constraintSet))){
    
    currConstraint <- strsplit(x = constraintSet[ii], split = " ", 
                               fixed = TRUE)[[1]]
    
    f.con[cnt, which(mt[, 1]==currConstraint[3])] <- 1
    ## f.dir <- c(f.dir, ">=")
    ## f.rhs <- c(f.rhs, currConstraint[1])
    f.dir[cnt] <- ">="
    f.rhs[cnt] <-  currConstraint[1]
    cnt <- cnt + 1
    f.con[cnt, which(mt[, 1]==currConstraint[3])] <- 1
    ## f.dir <- c(f.dir, "<=")
    ## f.rhs <- c(f.rhs, currConstraint[5])
    f.dir[cnt] <- "<="
    f.rhs[cnt] <-  currConstraint[5]
    cnt <- cnt + 1
    
  }
  
  ## ff <- list()
  ## ff[[length(ff)+1]] <- f.con
  ## ff[[length(ff)+1]] <- f.dir
  ## ff[[length(ff)+1]] <- f.rhs
  
  ff <- vector(mode = "list", length = 3)
  names(ff) <- c("con", "dir", "rhs")
  ff$con <- f.con
  ff$dir <- f.dir
  ff$rhs <- f.rhs
  
  return(ff)
  
}