R/analyze-auc_df.R

In wlandau/fbseqStudies: fbseqStudies

#' @title Function \code{auc_df}
#' @description make a auc data frame out of rds files extracted from simulation lists
#' @export
#' @param path to directory of extracted results files
auc_df = function(path){
  path = newdir(path)
  out = NULL
  for(f in list.files(path)){
    print(f)
    x = readRDS(paste0(path, f))
    for(n in names(x)){
      x[[n]]$fpr = x[[n]]$tpr = NULL
      x[[n]] = data.frame(x[[n]], heterosis = n)
    }
    mt = c(file = f, meta(f))
    m = as.data.frame(matrix(mt, nrow = nrow(x[[1]]), ncol = length(mt), byrow = T))
    colnames(m) = names(mt)
    for(n in names(x)) x[[n]] = cbind(x[[n]], m)
    for(n in names(x)) out = rbind(out, x[[n]])
  }
  out$analysis = gsub("normalnormal", "normal", out$analysis)
  out$analysis = gsub("normalLaplace", "Laplace", out$analysis)
  out$analysis = gsub("normalt", "t", out$analysis)
  out$analysis = gsub("normalhorseshoe", "horseshoe", out$analysis)
  out
}