R/annualize.R

In bocinsky/gutaker2020_rice_niche: Research Compendium for Gutaker et al. 2020 thermal niche model

#' Sum a GHCN record over a set of months to generate annual estimates.
#'
#' @param station GHCN station data
#' @param months The months over which to annualize
#'
#' @return A [data.frame] of annualized data.
#' @export
annualize <- function(station, months=1:12) {
  # first, remove unwanted months
  station <- station[station$MONTH %in% months, ]

  # then, split by year
  station.years <- split(station, station$YEAR)

  # then, sum over each year
  years.out <- lapply(station.years, function(year) {
    sum(year[, 3:33], na.rm = T)
  })

  return(as.data.frame(t(unlist(years.out))))
}