buffalo_snow: Historic Snowfall by Month for Buffalo

In trinker/exampledata: Example Data for Teaching

Description

A dataset containing historical monthly snowfall records. Columns are in character because of T cell values. The user will want to determine what T means, replace it with an appropriate value, and convert the columns to numeric.

Usage

data(buffalo_snow)

Format

A data frame with 78 rows and 13 variables

Details

• SEASON. Year span for winter season

• JUL. Snow fall in this month

• AUG. Snow fall in this month

• SEP. Snow fall in this month

• OCT. Snow fall in this month

• NOV. Snow fall in this month

• DEC. Snow fall in this month

• JAN. Snow fall in this month

• FEB. Snow fall in this month

• MAR. Snow fall in this month

• APR. Snow fall in this month

• MAY. Snow fall in this month

• JUN. Snow fall in this month

References

https://www.weather.gov/buf/BuffaloSnow

Examples

## Not run: 
library(tidyverse)
library(exampledata)
library(viridis)

bufsnow <- buffalo_snow %>%
    gather(Month, Snow, - SEASON)  %>%
    extract(SEASON, c('Start'), '(\\d{4})-\\d{2}', remove = FALSE) %>%
    filter(!is.na(Snow)) %>%
    mutate(
        Month = gsub('(^.)(.+$)', '\\U\\1\\L\\2', Month, perl = TRUE)  %>%
            factor(levels = month.abb),
        Snow = replace_na(as.numeric(Snow), .01),
        Start = as.integer(Start),
        Year = case_when(as.integer(Month) %in% 1:6 ~ as.integer(Start + 1), TRUE ~ Start)
    ) %>%  
    rename(Season = SEASON) %>%
    arrange(Season, Year, Month) %>% 
    select(Season, Year, Month, Snow)


bufsnow %>%
    ggplot(aes(x = Month, y = Snow, color = as.factor(Year), group = Year)) +
        geom_line() +
        coord_polar()


bufsnow %>%
    mutate(Month = factor(Month, levels = c(month.abb[c(7:12, 1:6)]))) %>%
    ggplot(aes(x = Month, y = Season, fill= Snow)) +
    geom_tile() +
    scale_fill_viridis()


bufsnow %>%
    group_by(Season) %>%
    summarize(total = sum(Snow)) %>%
    ggplot(aes(y=total, x=Season, group = 1)) +
        geom_line(size = 1) +
        geom_point() +
        theme(axis.text.x = element_text(angle = 45, hjust = 1, vjust =1))


## End(Not run)

trinker/exampledata documentation built on May 25, 2019, 8:32 p.m.