R/calculateCC.R

In dosorio/ccCINTROP: Calculate and plot the CC50

calculateCC <- function(models, CC = 50, raw = FALSE ){
  sdCC <- lapply(models, function(replicates){
    suppressMessages(suppressWarnings(sd(unlist(nplr::getEstimates(replicates,1-(CC/100))[c(2,3,4)]))))
  })
  CC <- lapply(models, function(replicates){
    suppressMessages(suppressWarnings(nplr::getEstimates(replicates,1-(CC/100))[3]))
  })
  seCC <- lapply(models, function(replicates){
    replicates@stdErr[1]
  })
  aCC <- lapply(models,function(replicates){
    replicates@goodness$gof
  })
  if (length(models)==0){
    return (c(CC="NR", SD="NR", SE="NR",mADJ=paste("<")))
  }else{
  if (raw == TRUE){
    cc <- lapply(models, function(replicates){
      suppressMessages(suppressWarnings(nplr::getEstimates(replicates,1-(CC/100))))
    })
  } else{
    cc <- return(c(CC = ifelse(test = is.nan(round(mean(unlist(CC),na.rm = TRUE),6)),yes = ">",no = ifelse(round(mean(unlist(CC),na.rm = TRUE),6)<0,yes = "<",no = round(mean(unlist(CC),na.rm = TRUE),6))), SD = round(mean(unlist(sdCC),na.rm = TRUE),6), SE = round(mean(unlist(seCC),na.rm = TRUE),6), mADJ = round(mean(unlist(aCC),na.rm = TRUE),6)))
  }}
  return(cc)
}