R/roc_selection.R

Defines functions roc_selection

Documented in roc_selection

roc_selection <- function(data, varname, outcomevar, low_bin, high_bin, step_bin){

  bin_range <- seq(low_bin,high_bin,step_bin)

  gk_gamma <- c()
  higham_gamma <- c()
  ag_all <- c()
  az_all <- c()
  sd_HR <- c()
  sd_FAR <- c()
  sd_criterion <- c()

  for (x in bin_range){

    binned_data <-bin_calculations(data, varname, outcomevar,x)
    gk_gamma <- c(gk_gamma, gamma_calc(data, varname, outcomevar, x))
    higham_gamma <- c(higham_gamma, higham_gamma(data, varname, outcomevar, x))
    sd_HR <- c(sd_HR, sd(binned_data$HR, na.rm = TRUE))
    sd_FAR <- c(sd_FAR, sd(binned_data$FAR, na.rm = TRUE))
    az_all <- c(az_all, az_calc(data, varname, outcomevar,x))
    ag_all <- c(ag_all, ag_calc(data, varname, outcomevar,x))
    sd_criterion <- c(sd_criterion, var(as.vector(criterion(data, varname, outcomevar,x))))

  }

  all_means <- data.frame(Bin = bin_range, gk_gamma, higham_gamma, ag_all, az_all, sd_HR, sd_FAR, sd_criterion)


  #plot(all_means$Bin,all_means$az_all, main = "Bin # by A_z Statistic", xlab = "Bin #", ylab = "A_g Statistic")
  return(all_means)
}
tmc2737/t2roc documentation built on May 29, 2019, 9:31 a.m.