Advertisement
Guest User

Untitled

a guest
Jun 20th, 2019
68
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.83 KB | None | 0 0
  1. library(tidyverse)
  2. library(gganimate)
  3.  
  4. tbl <- tibble(x = seq(-8, 8, by = .01),
  5. A_1 = 4*x,
  6. B_1 = x^2,
  7. A_2 = 3*x,
  8. B_2 = x^2,
  9. A_3 = 2*x,
  10. B_3 = x^2,
  11. A_4 = x,
  12. B_4 = x^2,
  13. A_5 = 0*x,
  14. B_5 = x^2) %>%
  15. gather(group, density, A_1:B_5) %>%
  16. separate(group, c("group", "type"), sep = "_") %>%
  17. mutate(type = as.numeric(type)) %>%
  18. mutate(Title = case_when(
  19. type == 1 ~ "A = 0, B = 4",
  20. type == 2 ~ "A = 0, B = 3",
  21. type == 3 ~ "A = 0, B = 2",
  22. type == 4 ~ "A = 0, B = 1",
  23. TRUE ~ "A = B = 0"
  24. ))
  25.  
  26.  
  27. ggplot(tbl) + geom_line(mapping = aes(x = x, y = density, colour = group)) +
  28. transition_states(Title, transition_length = .5, state_length = 2, wrap = TRUE) +
  29. labs(title = '{closest_state}') + ylab("f(x)")
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement