Code to produce AUC stats from Walker's data

walker_stats.R

library(farr)
library(tidyverse)

original <- read_csv("~/Downloads/WalkerDataCodeMar2021/original.csv")

original |>
  group_by(fyear) |>
  summarize(auc = auc(prob, aaer), .groups = "drop") |>
  mutate(avg_auc = mean(auc))

recode <- read_csv("~/Downloads/WalkerDataCodeMar2021/recode.csv")

recode |>
  group_by(fyear) |>
  summarize(auc = auc(prob, aaer), .groups = "drop") |>
  mutate(avg_auc = mean(auc))

rerun <- read_csv("~/Downloads/WalkerDataCodeMar2021/rerun.csv")

rerun |>
  group_by(fyear) |>
  summarize(auc = auc(prob, aaer), .groups = "drop") |>
  mutate(avg_auc = mean(auc))