## ----include = FALSE---------------------------------------------------------- knitr::opts_chunk$set( collapse = TRUE, comment = "", echo = FALSE, fig.height = 6, fig.width = 9, out.width = "100%" ) ## ----setup, include=FALSE----------------------------------------------------- library(knitr) library(Horsekicks) library(dplyr) library(tidyr) library(ggplot2) theme_set(theme_bw() + theme(plot.title = element_text(hjust = 0.5), strip.background = element_rect(fill = "antiquewhite"))) ## ----------------------------------------------------------------------------- kable(head(hkdeaths), align = "r") ## ----echo = TRUE-------------------------------------------------------------- hazards <- hkdeaths |> group_by(year) |> summarise(horse_kicks = sum(kick), falls = sum(fall), drownings = sum(drown), .groups = "drop") hazards_long <- hazards |> pivot_longer(cols = c(horse_kicks, falls, drownings), names_to = "cause", values_to = "deaths") ggplot(hazards_long) + aes(x = year, y = deaths, colour = cause) + geom_line(linewidth = 1.5) + scale_colour_viridis_d() + labs(title = "Annual Deaths from three Causes", x = "Year", y = "Deaths") ## ----echo=TRUE---------------------------------------------------------------- drown_1 <- glm(drownings ~ I(year - 1891), data = hazards, family = poisson(link = log)) exp(coef(drown_1)) |> rbind() |> as.data.frame(check.names = FALSE) |> kable(digits = 2) ## ----------------------------------------------------------------------------- b <- exp(coef(drown_1)) ## ----echo=TRUE---------------------------------------------------------------- ggplot(hkdeaths) + aes(x = year, y = drown) + geom_line(colour = "#21908CFF", linewidth=1.0) + facet_wrap(~ corps) + geom_smooth(method = "glm", formula = y ~ I(x-1891), method.args = list(family = poisson), linewidth = 1.0, colour = "#440154FF", se = FALSE) + labs(x = "Year", y = "Deaths", title = "Annual Drowning Deaths per Corps") ## ----eval=FALSE--------------------------------------------------------------- # # ggplot(hkdeaths) + aes(x = year, y = drown) + # # geom_line(colour = "#DF536B", linewidth=1.0) + facet_wrap(~ corps) + # # geom_smooth(method = "glm", formula = y ~ I(x-1891), method.args = list(family = poisson), # # linewidth = 2, colour = "#2297E6") + labs(y = "Deaths by drowning")