Histogram of Bump Times

# hours/minutes to minutes
parse_time = function(h, m) {
    return(h + m/60)
}

# minutes to hours/minutes
format_time = function(m) {
    m = ifelse(m >= 13 * 60, m - 12 * 60, m)
    return(sprintf("%d:%02d", m %/% 60, m %% 60))
}

times = c(
    # Chapter 1
    parse_time(11, 46),
    parse_time(11, 18),
    parse_time(11, 30),
    parse_time(10, 30),
    parse_time(11, 58),
    parse_time(13, 12),
    parse_time(9, 57),
    parse_time(11, 17),
    parse_time(11, 14),
    parse_time(9, 13),
    parse_time(11, 13),
    parse_time(9, 13),
    parse_time(12, 56),
    parse_time(12, 46),
    parse_time(8, 57),

    # Chapter 2
    parse_time(11, 12),
    parse_time(10, 31),
    parse_time(9, 17),
    parse_time(11, 31),
    parse_time(11, 13),
    parse_time(11, 56),
    parse_time(10, 31)
)

min_hour = 8
max_hour = 14
slots_per_hour = 4
labels_per_hour = 4

breaks = seq(min_hour, max_hour, 1/slots_per_hour)
labels = seq(min_hour, max_hour, 1/labels_per_hour)

hist(times, main="Histogram of Bump Times", xlab="Time (PM ET)", ylab="Count Per 15 Minutes", col="#ddddff", xlim=c(parse_time(min_hour, 0), parse_time(max_hour, 0)), xaxt='n', breaks=breaks)

axis(side=1, at=labels, labels=format_time(labels * 60), las=3)

times_density = density(times, na.rm=T, from=min_hour, to=max_hour)
lines(times_density$x, times_density$y * 4, col="red")