if (file.exists("betas.csv")) { d <- iaw$read.csv("betas.csv") system("gzip betas.csv") } else d <- iaw$read.csv("betas.csv.gz") s <- function( yearend ) { d21 <- subset(d, yyyymmdd == yearend*10000+1231 ) d21 <- d21[order(d21$tic),] d21 <- within(d21, bswa32 <- sprintf("%.02f", bswa32)) rownames(d21) <- NULL iaw$write.csv(d21, file <- paste0("betas-", yearend, ".csv")) } s(2021) s(2020)