R/plotModels.R

In dosorio/ccCINTROP: Calculate and plot the CC50

#' @export plotModels
plotModels <- function(files, modelUmbral=0.85, join = FALSE, nCol=1, ...){
  set.seed(1)
  nFiles <- length(list.files(files))
  if(nFiles>0){
    file <- list.files(files,full.names = TRUE)
    mExperiments <- lapply(file,function(file){
      data <- read.table(file,header = TRUE)
      data <- data[,paste0("X",sort(as.numeric(gsub("X","",colnames(data)))))]
      data <- cData(data)
      data <- sExperiments(data)
      mExperiments <- suppressMessages(generateModels(data,modelUmbral = modelUmbral,join = TRUE))
    })
    names(mExperiments) <- sub(pattern = "(.*)\\..*$", replacement = "\\1", basename(file))
    mExperiments <- mExperiments[lengths(mExperiments)>0]
  } else {
    file <- files
    data <- read.table(file,header = TRUE)
    data <- data[,paste0("X",sort(as.numeric(gsub("X","",colnames(data)))))]
    data <- cData(data)
    data <- sExperiments(data)
    mExperiments <- suppressMessages(generateModels(data,modelUmbral = modelUmbral,join = join))
    if (join == TRUE){
      mExperiments <- list(mExperiments)
      names(mExperiments) <- sub(pattern = "(.*)\\..*$", replacement = "\\1", basename(file))
    }
  }
  if(length(mExperiments)==0){
    plot(c(0,max(unlist(rangeValues))),range(0,1.15),type = "n", bty = "n")
  } else{
    nplr::overlay(mExperiments,showLegend = FALSE, Cols = rainbow(length(mExperiments)),ylim=c(0,1.15),las=1,pch=16)
    legend("topright",legend = names(mExperiments),col = rainbow(length(mExperiments)),bty = "n",ncol = nCol, lwd = 2, pch = 19)
  }
}