R/generateModels.R

In dosorio/ccCINTROP: Calculate and plot the CC50

generateModels <- function(experiment, modelUmbral = 0.85, join = FALSE){
  models <- lapply(experiment,function(replicates){
    replicates <- as.data.frame.array(replicates)
    suppressWarnings(nplr::nplr(as.numeric(replicates[,2]),as.numeric(replicates[,1]),useLog = FALSE))
  })
  selectedModels <- (lapply(models,function(model){model@goodness$gof}) >= modelUmbral)
  models <- models[selectedModels]
  names(models) <- paste0("E",seq_len(length(models)))[selectedModels]
  if(join == TRUE & length(models)>0){
    experiment <- do.call(rbind.data.frame,experiment[selectedModels])
    models <- suppressWarnings(nplr::nplr(as.numeric(experiment[,2]),as.numeric(experiment[,1]),useLog = FALSE))
  }
  return (models)
}