animated barplot

#part_2: an animated barplot for the development of food price of some products in Uzbekistan over year and months using ggplot2 and gganimate
barplot
R
ggplot2
gganimate
geom_image
Author

Fazliddin Sultonov

Published

Feb 8, 2025

Datasource: Statistics Agency under the President of the Republic of Uzbekistan | Datastorage: GitHub

Data

Code 
data_raw <- read.csv2(
  "https://raw.githubusercontent.com/sultanovf/my_datasets/refs/heads/main/priceuz_202101-202411.csv")

  
price_uz <-
  data_raw |> 
  select(-runame) |> 
  add_column(runame, .after = "enname" )

datseq <- as.character(seq(as.Date("2021-01-01"),as.Date("2024-11-01"), by = "month"))
new_colnames <- c("uzname", "enname", "runame", unlist(strsplit(datseq, ",")))
colnames(price_uz) <- new_colnames


rmarkdown::paged_table(sample_n(price_uz, 10)) # show sample 10 rows
# oder {r, df_print = "paged"}
# sample_n(data_raw, 10

dataclean

Code 
# add column catecory
df_data <- price_uz |> 
  select(-enname) |> 
  filter(uzname %in% filtermap) |> 
  mutate(
    category = case_when(
      uzname %in% gosht ~ "Мясо",
      uzname %in% sabzavot ~ "Овощи",
      uzname %in% oziqovqat ~ "Питание",
      uzname %in% yog ~ "Mасло",
      TRUE ~ "Другой"
    )
  )

# pivot data
df_piv <-
  df_data |> 
  pivot_longer(
    cols = !c(uzname, runame, category),
    names_to = "date",
    values_to = "price"
  ) |> 
  mutate(
    date = as.Date(date)
  ) 

df_piv |> 
  sample_n(7)
#> # A tibble: 7 × 5
#>   uzname       runame          category date        price
#>   <chr>        <chr>           <chr>    <date>      <dbl>
#> 1 Non          Хлеб            Питание  2024-10-01  6114.
#> 2 Guruch       Рис             Питание  2024-04-01 20045.
#> 3 Tovuq goshti Курица          Мясо     2022-10-01 26676.
#> 4 Mol goshti   Говядина        Мясо     2021-01-01 59082.
#> 5 Sariyog      Масло сливочное Mасло    2021-05-01 51568 
#> 6 Tovuq goshti Курица          Мясо     2023-11-01 31566.
#> 7 Non          Хлеб            Питание  2023-08-01  5776.

Prepare Data for Plot

Code 
df2 <- 
  df_piv |>
  filter(
    !(runame %in% c("Баранина", "Говядина", "Масло сливочное", "Курица"))
  ) |> 
  arrange(date,price) |>
  group_by(date) |>
  # mutate(date = format(date, "%Y %b")) |>
  mutate(
    images = case_when(
      runame == "Пшеница" ~ paste0(img_path,uzname,".png"),
      runame == "Рис" ~ paste0(img_path,uzname,".png"),
      runame == "Хлеб" ~ paste0(img_path,uzname,".png"),
      runame == "Масло(хлопковое)" ~ paste0(img_path,uzname,".png"),
      runame == "Масло(подсолнечное)" ~ paste0(img_path,uzname,".png"),
      runame == "Морковь" ~ paste0(img_path,uzname,".png"),
      runame == "Чеснок" ~ paste0(img_path,uzname,".png"),
      runame == "Лук" ~ paste0(img_path,uzname,".png"),
      runame == "Картофель" ~ paste0(img_path,uzname,".png"),
      runame == "Сахар" ~ paste0(img_path,uzname,".png"),
      TRUE ~ uzname
      
    ),
    frame = row_number(),
    rank = 1:n(),
    #ord_runame = (fct_reorder(runame, price))
  ) |>  ungroup()

head(df2)
#> # A tibble: 6 × 8
#>   uzname    runame    category date       price images               frame  rank
#>   <chr>     <chr>     <chr>    <date>     <dbl> <chr>                <int> <int>
#> 1 Sabzi     Морковь   Овощи    2021-01-01 2073. gallery_img/product…     1     1
#> 2 Piyoz     Лук       Овощи    2021-01-01 2081. gallery_img/product…     2     2
#> 3 Non       Хлеб      Питание  2021-01-01 3137. gallery_img/product…     3     3
#> 4 Bugdoy    Пшеница   Питание  2021-01-01 3255. gallery_img/product…     4     4
#> 5 Kartoshka Картофель Овощи    2021-01-01 4033. gallery_img/product…     5     5
#> 6 Shakar    Сахар     Питание  2021-01-01 7360. gallery_img/product…     6     6

Plot animatet #1

Code 
# Title/subtitle/legendtitle
t2 <- paste("Цены на продукты питания")
p2 <- paste("категории: ") 


f_path <- "C:/Users/sultanov/Documents/fontawesome/otfs/"
font_add('Font Awesome 6 Brands', paste0(f_path, "Font Awesome 6 Brands-Regular-400.otf"))

showtext::showtext_auto()

p_df2 <-
  df2 |> 
  ggplot(aes(
    x= rank,
    y = price,
    fill = category))+ #    color = runame
  geom_col()+
  geom_text(aes(label = as.character(price)), hjust = 0, size = 4) +
  geom_image(aes(x = rank,  image = images),y = 0,   size = 0.1)+ 
  
  
  scale_y_continuous(
    limits = c(0,40000),breaks = seq(0, 40000, 5000),
    labels = function(x) paste0(x/1000, "k"),
    expand = expansion(c(0.09, 0.02))) +
  #expand_limits(y = 0)+
  #scale_y_continuous(labels = scales::comma) + # Format y-axis values
  #scale_x_reverse() +   # put rank 1 at the top/ Highest values on top
  #guides(color = FALSE, fill = FALSE)+
  
  coord_flip(clip = "off", expand = FALSE) +
  geom_label(aes(label = paste("Time:",format(date, "%Y %b"))),
             x = 1.6, y = 32000, size = 7, fill = "floralwhite", col = "gray")+
  labs(title = t2,
       subtitle = "DataSource: ©Агентство статистики при Президенте Республики Узбекистан",
       caption = social_caption_map, 
       x = "продукты", 
       y = "цены в тысячах(k = 1000)",
       fill = p2) +

  theme_minimal()+
  theme(plot.title = element_text(size = 16, hjust = 0.5 ),
        plot.subtitle = element_text(size = 14, color = "steelblue4", hjust = 0.5),
        plot.caption = ggtext::element_markdown(size = 12, hjust = 0.5),
        legend.position = "top",
        axis.text.y = element_blank(),
        axis.title.y = element_text(vjust = 11, size = 9),
        axis.title.x = element_text(size = 10),
        legend.box.background =element_rect("floralwhite")) +
  
  transition_states(date)+
  ease_aes("cubic-in-out") + #bounce-in cubic-in-out
  enter_fade() +
  exit_fade()

p_df2

# anim_1 <- p_df2 |>
#   animate(
#     duration = 20,
#     fps = 10,
#     start_pause = 5,
#     end_pause = 5,
#     #height = 530,
#     #width = 832,
#     #rewind = TRUE,
#     renderer = av_renderer()
#   )
# anim_1

Plot animated #2

Code 
p1_df2<- df2 |> 
  ggplot(aes(
    xmin= 0,
    xmax = price,
    y = rank,
    ymin = rank - 0.45,
    ymax = rank + 0.45,
    group = runame,
    fill = category))+
  geom_rect(alpha = 0.5) +
  
  facet_wrap(~date) +
  facet_null() +
  
  scale_x_continuous(
    limits = c(0,40000),breaks = seq(0, 40000, 5000),
    labels = function(x) paste0(x/1000, "k"),
    expand = expansion(c(0.05, 0.02)))+

  geom_image(aes(x=price, image = images), x= 0, size = 0.1)+
  geom_text(aes(label = paste("Time:", format(date, "%Y %b"))),
             x = 32000, y = 1.5, size = 6, color = "tomato4")+
  
  
  geom_label(aes(x = price, label = as.character(price)), colour = "white")+
  
  labs(title = t2, 
       subtitle = "DataSource: ©Агентство статистики при Президенте Республики Узбекистан",
       caption = social_caption_map,
       x = "цены в тысячах (k = 1000)", y = "продукты", 
       fill = p2)+   # legend title
  theme_minimal()+
  theme(
    plot.title = element_text(size = 16, color = "steelblue4", hjust = 0.5),
    plot.subtitle = element_text(size = 14, color = "tomato4", hjust = 0.5),
    plot.caption = element_textbox_simple(
      size = 12,
      lineheight = 1,
      padding = margin(5.5, 5.5, 5.5, 5.5),
      margin = margin(0, 0, 5.5, 0),
      hjust = -0.5,
      fill = "cornsilk"), #
    axis.ticks.y = element_blank(), 
    axis.text.y = element_blank(),
    # panel.grid.major.y = element_blank(),
    # panel.grid.minor.y = element_blank(),
    legend.position = "top",
    legend.box.background =element_rect("floralwhite")) +
  
  transition_reveal(date)+
  enter_fade() +
  exit_fade()

p1_df2

# anim_2 <- p1_df2 |>
#   animate(
#     duration = 20,
#     fps = 10,
#     start_pause = 5,
#     end_pause = 5,
#     # height = 530, #530
#     #width = 844, #633
#     #rewind = TRUE,
#     renderer = av_renderer()
#   )
# anim_2

Saved as Video

Video

Back to top