dblogr.com/

Sugarbeet Production in USA

Graphs of sugarbeet production using USDA data

Prepare Data

# devtools::install_github("derekmichaelwright/agData")
library(agData)
# Prep data
dd <- agData_USDA_Sugarbeet
myCaption <- "www.dblogr.com/ or derekmichaelwright.github.io/dblogr/ | Data: USDA"

Production

2017

# Prep data
xx <- dd %>% 
  filter(Measurement == "Production", Year == 2017) %>% 
  group_by(State) %>%
  summarise(Value = sum(Value)) %>%
  ungroup() %>% arrange(desc(Value)) %>%
  mutate(State = factor(State, levels = .$State))
# Plot
mp <- ggplot(xx, aes(x = State, y = Value / 1000000)) + 
  geom_bar(stat = "identity", color = "black", 
           fill = "darkgreen", alpha = 0.7) + 
  theme_agData(axis.text.x = element_text(angle = 45, hjust = 1)) +
  labs(title = "Sugarbeet Production - 2017", x = NULL,
       y = "Million Tonnes", caption = myCaption)
ggsave("sugarbeet_usa_01.png", mp, width = 6, height = 4)

All Data

# Prep data
xx <- agData_USDA_Sugarbeet %>% 
  filter(Measurement == "Production") %>% 
  group_by(Year, State) %>%
  summarise(Value = sum(Value))
# Plot
mp <- ggplot(xx, aes(x = Year, y = Value / 1000000)) + 
  geom_line(color = "darkgreen", alpha = 0.7) + 
  facet_wrap(State ~ ., ncol = 6) +
  theme_agData(axis.text.x = element_text(angle = 45, hjust = 1)) +
  labs(title = "Sugarbeet Production", x = NULL,
       y = "Million Tonnes", caption = myCaption)
ggsave("sugarbeet_usa_02.png", mp, width = 10, height = 6)

Yield

Minnesota + Idaho + North Dakota

# Prep data
myColors <- c("darkblue", "darkgreen", "steelblue")
myStates <- c("Minnesota", "Idaho", "North Dakota")
xx <- agData_USDA_Sugarbeet %>% 
  filter(Measurement == "Yield", State %in% myStates) %>% 
  group_by(Year, State) %>%
  summarise(Value = mean(Value)) %>%
  ungroup() %>% arrange(desc(Value)) %>%
  mutate(State = factor(State, levels = myStates),
         Group = ifelse(Year < 2004, "1", "2"),
         Group = paste(State, Group) )
# Plot
mp <- ggplot(xx, aes(x = Year, y = Value / 1000, 
                     color = State, shape = State)) + 
  geom_point(alpha = 0.7) +
  geom_line(alpha = 0.2) +
  stat_smooth(aes(group = Group), geom = "line", method = "lm", 
              se = F, size = 1.25, alpha = 0.7) +
  scale_color_manual(name = NULL, values = myColors) +
  scale_shape_manual(name = NULL, values = c(15:17)) +
  theme_agData(legend.position = "bottom") +
  labs(title = "Sugarbeet Yield", x = NULL,
       y = "tonnes / ha", caption = myCaption)
ggsave("sugarbeet_usa_03.png", mp, width = 6, height = 4)

Wyoming + Colorado + Nebraska

# Prep data
myColors <- c("darkred", "antiquewhite4", "darkgoldenrod2")
myStates <- c("Wyoming", "Colorado", "Nebraska")
xx <- agData_USDA_Sugarbeet %>% 
  filter(Measurement == "Yield", State %in% myStates) %>% 
  group_by(Year, State) %>%
  summarise(Value = mean(Value)) %>%
  ungroup() %>% arrange(desc(Value)) %>%
  mutate(State = factor(State, levels = myStates),
         Group = ifelse(Year < 2004, "1", "2"),
         Group = paste(State, Group) )
# Plot
mp <- ggplot(xx, aes(x = Year, y = Value / 1000, 
                     color = State, shape = State)) + 
  geom_point(alpha = 0.7) +
  geom_line(alpha = 0.2) +
  stat_smooth(aes(group = Group), geom = "line", method = "lm", 
              se = F, size = 1.25, alpha = 0.7) +
  scale_color_manual(name = NULL, values = myColors) +
  scale_shape_manual(name = NULL, values = c(15:17)) +
  theme_agData(legend.position = "bottom") +
  labs(title = "Sugarbeet Yield", x = NULL,
       y = "tonnes / ha", caption = myCaption)
ggsave("sugarbeet_usa_04.png", mp, width = 6, height = 4)

dblogr.com/


© Derek Michael Wright