R/tbl_blister.R

In commfish/scallopr: Analyze ADF&G weathervane scallop survey data

#' tbl_blister
#'
#' @param scal_awl
#' @param tows
#'
#' @description  Returns a .csv table of scallop blister percentage by bed
#' @export tbl_blister
#'
#' @examples
#' tbl_blister(scal_awl, tows)
tbl_blister <- function(scal_awl, tows){

  tbl_names <- c("Bed", "N", "0%", "1-24%", "25-49%", "50-74%", "75-100%")

  scal_awl %>%
    left_join(tows) %>%
    filter(!is.na(mud)) %>%
    count(Bed, mud) %>%
    group_by(Bed) %>%
    mutate(prop = round(n / sum(n) * 100, 1)) %>%
    spread(mud, prop) %>%
    replace(is.na(.), 0) %>%
    mutate(`4` = 0) %>%
    summarise(N = sum(n),
              zero = round(max(`0`),1),
              one = round(max(`1`),1),
              two = round(max(`2`),1),
              three = round(max(`3`),1),
              four = round(max(`4`),1)) %>%
    rename_at(names(.), function(x) tbl_names) -> x

  write_csv(x, here::here(paste0("output/", YEAR, "/mud.csv")))
  x
}