# Prep data
myAreas <- c("British Columbia", "Alberta", "Saskatchewan",
             "Manitoba", "Ontario", "Quebec")
xx <- d1 %>% 
  filter(Year == 2016, Item == "Total area of farms", Unit == "Hectares")
x1 <- xx %>% filter(Area == "Canada")
xx <- xx %>% filter(Area != "Canada") %>% 
  mutate(Percent = round(100 * Value / x1$Value)) %>% 
  arrange(Percent) %>%
  mutate(Area = factor(Area, levels = rev(Area)),
         Area_Short = plyr::mapvalues(Area, d0$Area, d0$Area_Short),
         Label = ifelse(Area %in% myAreas, paste(Percent, "%"), NA))
# Plot
mp <- ggplot(xx, aes(x = Area_Short, y = Value / 1000000)) +
  geom_bar(aes(fill = Area), stat = "identity", 
           color = "black", alpha = 0.7) +
  geom_label(aes(label = Label), nudge_y = 1.1, size = 4) +
  scale_fill_manual(values = agData_Colors) +
  theme_agData(legend.position = "none", 
               axis.text.x = element_text(angle = 45, hjust = 1)) +
  labs(title = "Canadian Farmland (2016)", x = NULL, 
       y = "Million Hectares", caption = myCaption)
ggsave("farmland_canada_1_01.png", mp, width = 6, height = 4)

# Prep data
xx <- d1 %>% 
  filter(Year == 2016, Item == "Total area of farms", Unit == "Hectares")
x1 <- xx %>% filter(Area == "Canada")
xx <- xx %>% filter(Area != "Canada") %>% 
  mutate(Percent = round(100 * Value / x1$Value)) %>% 
  arrange(Percent) %>% filter(Percent > 0) %>%
  mutate(Label_y = cumsum(Percent),
         Label_y = Label_y - (Percent / 2),
         Label = ifelse(Percent > 6, paste(Percent, "%"), NA)) %>%
  arrange(desc(Percent)) %>%
  mutate(Area = factor(Area, levels = Area),
         Area_Short = plyr::mapvalues(Area, d0$Area, d0$Area_Short))
# Plot
mp <- ggplot(xx, aes(x = 1)) +
  geom_bar(aes(y = -Percent, fill = Area), stat = "identity", 
           color = "black", alpha = 0.7) +
  coord_polar("y", start = 0) +
  geom_label(aes(y = -Label_y, label = Label), 
             nudge_x = 0.275, size = 3, alpha = 0.9) +
  
  theme_agData_pie(plot.caption = element_text(size = 7)) +
  xlim(0.542, 1.45) +
  scale_fill_manual(name = NULL, values = agData_Colors) +
  labs(title = "Canadian Farmland (2016)", x = NULL, 
       y = "Million Hectares", caption = myCaption)
ggsave("farmland_canada_1_02.png", mp, width = 6, height = 4)

sum(xx$Value)

# Prep data
xx <- d1 %>%
  filter(Item == "Total area of farms",
         Unit == "Hectares", Area != "Canada")
# Plot
mp <- ggplot(xx, aes(x = Year, y = Value / 1000000)) +
  geom_col(alpha = 0) +
  geom_line(color = "darkgreen", alpha = 0.7, size = 1) + 
  geom_point(alpha = 0.3, size = 0.5) +
  facet_wrap(. ~ Area, scales = "free_y", ncol = 5) +
  scale_x_continuous(breaks = seq(1925, 2015, 20), expand = c(0.01,0)) +
  theme_agData(axis.text.x = element_text(angle = 45, hjust = 1)) +
  labs(title = "Canadian Farmland", x = NULL, 
       y = "Million Hectares", caption = myCaption)
ggsave("farmland_canada_1_03.png", mp, width = 12, height = 5)

# Prep data
xx <- d1 %>%
  filter(Item == "Total area of farms", Unit == "Hectares",
         Area %in% c("Ontario", "Quebec"))
# Plot
mp <- ggplot(xx, aes(x = Year, y = Value / 1000000)) +
  geom_col(position = "dodge", alpha = 0) +
  geom_line(aes(color = Area), size = 1.25, alpha = 0.7) +
  geom_point(alpha = 0.3) +
  scale_color_manual(name = NULL, values = c("darkblue", "steelblue")) +
  scale_x_continuous(breaks = seq(1925, 2015, 10), expand = c(0.01,0)) +
  theme_agData(legend.position = "bottom") +
  labs(title = "Ontario & Quebec Farmland", x = NULL, 
       y = "Million Hectares", caption = myCaption)
ggsave("farmland_canada_1_04.png", mp, width = 6, height = 4)

# Prep data
myItems <- c("Land in crops", "Tame or seeded pasture", "Summerfallow land")
xx <- d1 %>%
  filter(Item %in% myItems, Area == "Canada",
         Measurement == "Hectares")
# Plot
mp <- ggplot(xx, aes(x = Year, y = Value / 1000000)) +
  geom_col(alpha = 0) +
  geom_line(aes(color = Item), alpha = 0.7, size = 1) + 
  geom_point(alpha = 0.3, size = 0.5) +
  facet_wrap(Item ~ ., scales = "free") +
  scale_color_manual(values = c("darkred", "darkgreen", "darkblue")) +
  theme_agData(legend.position = "none") +
  labs(title = "Farmland in Canada", x = NULL, 
       y = "Million Hectares", caption = myCaption)
ggsave("farmland_canada_1_05.png", mp, width = 8, height = 3)

# Prep data
xx <- d2 %>% 
  filter(Item %in% c("Total area of farms", "Total number of farms"),
         Unit %in% c("Number", "Hectares")) %>%
  select(-Item, -Measurement) %>% 
  spread(Unit, Value) %>%
  mutate(Size = Hectares / Number)
# Plot
mp <- ggplot(xx %>% filter(Area == "Canada"), aes(y = Size, x = Year)) + 
  geom_col(alpha = 0) +
  geom_line(color = "darkgreen", alpha = 0.7, size = 1) + 
  geom_point(alpha = 0.3, size = 0.75) +
  scale_x_continuous(breaks = seq(1925, 2015, 10), expand = c(0.01,0)) +
  theme_agData() +
  labs(title = "Average Farm Size in Canada", x = NULL, 
       y = "Hectares / Farm", caption = myCaption)
ggsave("farmland_canada_2_01.png", mp, width = 6, height = 4)

# Plot
mp <- ggplot(xx %>% filter(Area != "Canada"), aes(y = Size, x = Year)) + 
  geom_col(alpha = 0) +
  geom_line(color = "darkgreen", alpha = 0.7, size = 1) + 
  geom_point(alpha = 0.3, size = 0.5) +
  facet_wrap(Area ~ ., ncol = 5) +
  scale_x_continuous(breaks = seq(1925, 2005, 20)) +
  theme_agData(axis.text.x = element_text(angle = 45, hjust = 1)) +
  labs(title = "Average Farm Size in Canada", x = NULL,
       y = "Hectares / Farm", caption = myCaption)
ggsave("farmland_canada_2_02.png", mp, width = 12, height = 5)

# Prep data
xx <- d2 %>% 
  filter(Area == "Canada", Unit != "Acres",
         Item %in% c("Total area of farms", "Total number of farms")) %>%
  mutate(Value = ifelse(Unit == "Hectares", Value / 100000, Value / 1000),
         Unit = plyr::mapvalues(Unit, 
            c("Hectares", "Number of farms reporting"), 
            c("Total area of farms (100,000 hectares)",
              "Number of farms (x1000)")))
# Plot
mp <- ggplot(xx, aes(y = Value, x = Year)) + 
  geom_col(position = "dodge", alpha = 0) +
  geom_line(aes(color = Unit), size = 1.5, alpha = 0.7) +
  geom_point(alpha = 0.3, size = 0.5) +
  scale_color_manual(name = NULL, values = c("steelblue", "darkgreen")) +
  scale_x_continuous(breaks = seq(1925, 2015, 10), expand = c(0.01,0)) +
  theme_agData(legend.position = "bottom") +
  labs(title = "Average Farm Size in Canada", 
       y = NULL, x = NULL, caption = myCaption)
ggsave("farmland_canada_2_03.png", mp, width = 6, height = 4)

# Prep data
xx <- d2 %>% 
  filter(Area != "Canada", Unit != "Acres",
         Item %in% c("Total area of farms", "Total number of farms")) %>%
  mutate(Value = ifelse(Unit == "Hectares", Value / 100000, Value / 1000),
         Unit = plyr::mapvalues(Unit, 
            c("Hectares", "Number of farms reporting"), 
            c("Total area of farms (100,000 hectares)",
              "Number of farms (x1000)")))
# Plot
mp <- ggplot(xx, aes(y = Value, x = Year)) + 
  geom_col(position = "dodge", alpha = 0) +
  geom_line(aes(color = Unit), size = 2, alpha = 0.7) +
  geom_point(alpha = 0.3, size = 0.5) +
  facet_wrap(Area ~ ., scales = "free_y", ncol = 5) +
  scale_color_manual(name = NULL, values = c("steelblue", "darkgreen")) +
  scale_x_continuous(breaks = seq(1925, 2005, 20)) +
  theme_agData(legend.position = "bottom",
               axis.text.x = element_text(angle = 45, hjust = 1)) +
  labs(title = "Average Farm Size in Canada", 
       y = NULL, x = NULL, caption = myCaption)
ggsave("farmland_canada_2_04.png", mp, width = 12, height = 5)

# Prep data
myColors <- alpha("darkgreen", seq(0.1, 1, by = 0.1))
myAreas <- c("British Columbia", "Alberta", "Saskatchewan",
             "Manitoba", "Ontario", "Quebec")
xx <- d3 %>% 
  filter(Area %in% myAreas, Measurement != "Total number of farms") %>%
  mutate(Year = factor(Year))
#
mp <- ggplot(xx, aes(y = Value / 1000, x = Measurement, fill = Year)) + 
  geom_col(position = "dodge", color = "black", lwd = 0.2) +
  facet_grid(Area ~ .) +
  scale_fill_manual(name = NULL, values = myColors) +
  theme_agData(axis.text.x = element_text(angle = 45, hjust = 1)) +
  labs(title = "Farm Size in Canada", x = NULL,
       y = "Thousand Farms", caption = myCaption)
ggsave("farmland_canada_3_07.png", mp, width = 9, height = 8)

# Prep data
myColors <- c("maroon4", "darkgoldenrod3", "darkgreen",
              "darkred", "steelblue4", "darkslategrey")
myAreas <- c("British Columbia", "Alberta", "Saskatchewan",
             "Manitoba", "Ontario", "Quebec")
xx <- d3 %>% 
  filter(Area %in% myAreas, Measurement != "Total number of farms") %>%
  mutate(Area = factor(Area, levels = myAreas))
#
mp <- ggplot(xx, aes(y = Value / 1000, x = Measurement, fill = Area)) + 
  geom_col(color = "black", lwd = 0.2) +
  facet_grid(Area ~ ., scales = "free_y") +
  scale_fill_manual(values = myColors) +
  theme_agData(legend.position = "none",
               axis.text.x = element_text(angle = 45, hjust = 1)) +
  labs(title = "Farm Size in Canada | {round(frame_time)}", 
       y = "Thousand Farms", x = NULL, caption = myCaption) +
  transition_time(Year)
anim_save("farland_canada_gif_01.gif", mp,
          nframes = 300, fps = 25, end_pause = 60, 
          width = 900, height = 1200, res = 150, units = "px")

# Prep data
xx <- d4 %>% filter(Area == "Canada")
# Plot
mp <- ggplot(xx, aes(x = Year, y = Value)) +
  geom_line(alpha = 0.7, size = 1.5, color = "darkgreen") +
  theme_agData() +
  labs(title = "Canadian Farm Prices", subtitle = "Farm land and buildings",
       y = "Dollars Per Acre", x = NULL, caption = myCaption)
ggsave("farmland_canada_4_01.png", mp, width = 6, height = 4)

# Prep data
xx <- d4 %>% filter(Area != "Canada")
# Plot
mp <- ggplot(xx, aes(x = Year, y = Value)) +
  geom_line(alpha = 0.7, size = 1.5, color = "darkgreen") +
  facet_wrap(Area ~ ., ncol = 5) +
  theme_agData() +
  labs(title = "Canadian Farm Prices", subtitle = "Farm land and buildings",
       y = "Dollars Per Acre", x = NULL, caption = myCaption)
ggsave("farmland_canada_4_02.png", mp, width = 12, height = 6)

# Prep data
xx <- d4 %>% filter(Year == 2022) %>% arrange(rev(Value)) %>%
  mutate(Area = factor(Area, levels = .$Area),
         Group = ifelse(Area == "Canada", "Can", "Provs"))
# Plot
mp <- ggplot(xx, aes(x = Area, y = Value, fill = Group)) +
  geom_col(alpha = 0.7, color = "black") +
  scale_fill_manual(values = c("steelblue","darkgreen")) +
  theme_agData(legend.position = "none",
               axis.text.x = element_text(angle = 45, hjust = 1)) +
  labs(title = "Canadian Farm Prices", subtitle = "Farm land and buildings",
       y = "Dollars Per Acre", x = NULL, caption = myCaption)
ggsave("farmland_canada_4_03.png", mp, width = 6, height = 4)