R/cleanData.R

In dosorio/ccCINTROP: Calculate and plot the CC50

cleanData <- function(file){
  set.seed(1)
  data <- read.table(file,header = TRUE)
  data <- data[,paste0("X",sort(as.numeric(gsub("X","",colnames(data)))))]
  b <- sum(is.na(data))
  while(length(boxplot(data,plot = FALSE)$out) != 0){
    for(i in seq_len(dim(data)[2])){
      data[data[,i] %in% boxplot.stats(data[,i])$out,i] <- NA
    }
  }
  a <- sum(is.na(data))
  message(paste(a-b, "values removed"))
  return(data)
}