Guest User

Untitled

a guest
Jan 19th, 2019
130
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 4.57 KB | None | 0 0
  1. require(ggplot2)
  2. require(scales)
  3. require(gridExtra)
  4.  
  5. mymelt <- structure(list(mydate = structure(c(15340, 15340, 15340, 15340,
  6. 15340, 15340, 15340, 15340, 15340, 15340, 15340, 15340, 15371,
  7. 15371, 15371, 15371, 15371, 15371, 15371, 15371, 15371, 15371,
  8. 15371, 15371, 15400, 15400, 15400, 15400, 15400, 15400, 15400,
  9. 15400, 15400, 15400, 15400, 15400, 15431, 15431, 15431, 15431,
  10. 15431, 15431, 15431, 15431, 15431, 15431, 15431, 15431, 15461,
  11. 15461, 15461, 15461, 15461, 15461, 15461, 15461, 15461, 15461,
  12. 15461, 15461, 15492, 15492, 15492, 15492, 15492, 15492, 15492,
  13. 15492, 15492, 15492, 15492, 15492, 15522, 15522, 15522, 15522,
  14. 15522, 15522, 15522, 15522, 15522, 15522, 15522, 15522, 15553,
  15. 15553, 15553, 15553, 15553, 15553, 15553, 15553, 15553, 15553,
  16. 15553, 15553), class = "Date"), variable = c("b", "bc", "f",
  17. "in", "it", "l", "of", "o", "pr", "s", "total", "tr", "b", "bc",
  18. "f", "in", "it", "l", "of", "o", "pr", "s", "total", "tr", "b",
  19. "bc", "f", "in", "it", "l", "of", "o", "pr", "s", "total", "tr",
  20. "b", "bc", "f", "in", "it", "l", "of", "o", "pr", "s", "total",
  21. "tr", "b", "bc", "f", "in", "it", "l", "of", "o", "pr", "s",
  22. "total", "tr", "b", "bc", "f", "in", "it", "l", "of", "o", "pr",
  23. "s", "total", "tr", "b", "bc", "f", "in", "it", "l", "of", "o",
  24. "pr", "s", "total", "tr", "b", "bc", "f", "in", "it", "l", "of",
  25. "o", "pr", "s", "total", "tr"), value = c(-23, 6.90000000000001,
  26. 459.799999999999, -403.6, -56.1, -95, -13.8, 32.6, 121.5, -15.7,
  27. 26.2000000000007, 12.5, -25.1, 238.3, 1047.2, -803.2, -151.5,
  28. -260.5, -59.6, -93.8, 461.5, -37.7, 26.7999999999993, -288.8,
  29. -46.4, 249, 1289.8, -783.2, -188.1, -414.9, -77.7, -61, 928.4,
  30. -36.8, 17.4000000000015, -841.7, -46.5, 276.2, 1384.8, -541.1,
  31. -71.8999999999999, -433.3, -61.3, -28.3, 494.699999999999, -23.4,
  32. -14.5999999999985, -964.5, -46.1, 376.2, 1020.1, -119.4, 56.8000000000001,
  33. -447.7, -9.50000000000001, 14.2, -9.20000000000164, 2.5, -42.7999999999993,
  34. -880.6, -52.9, 345.5, 892.599999999999, -241.8, 144.3, -428.2,
  35. -3.30000000000001, 91.9, -294.800000000002, -5.19999999999999,
  36. -42.1999999999971, -490.1, -64.5, 379.7, 679.299999999999, -143.1,
  37. 185.9, -419.8, -4.30000000000001, 182.4, -421.900000000002, 1.80000000000001,
  38. -59.8999999999978, -435.2, -80.2, 422.2, 645.499999999998, -391.4,
  39. 76.6000000000001, -387.4, -1.70000000000001, 211.2, -131.500000000002,
  40. -10.6, -40.8999999999978, -393.6), fill = c("#A4D3EE80", "#A478AB80",
  41. "#01AEF080", "#8DC73F80", "#F8931D80", "#FFAAAA80", "#8C8C8C",
  42. "#D38D5F80", "#23238E80", "#77B9B780", "#C8373780", "#EEDD8280",
  43. "#A4D3EE80", "#A478AB80", "#01AEF080", "#8DC73F80", "#F8931D80",
  44. "#FFAAAA80", "#8C8C8C", "#D38D5F80", "#23238E80", "#77B9B780",
  45. "#C8373780", "#EEDD8280", "#A4D3EE80", "#A478AB80", "#01AEF080",
  46. "#8DC73F80", "#F8931D80", "#FFAAAA80", "#8C8C8C", "#D38D5F80",
  47. "#23238E80", "#77B9B780", "#C8373780", "#EEDD8280", "#A4D3EE80",
  48. "#A478AB80", "#01AEF080", "#8DC73F80", "#F8931D80", "#FFAAAA80",
  49. "#8C8C8C", "#D38D5F80", "#23238E80", "#77B9B780", "#C8373780",
  50. "#EEDD8280", "#A4D3EE80", "#A478AB80", "#01AEF080", "#8DC73F80",
  51. "#F8931D80", "#FFAAAA80", "#8C8C8C", "#D38D5F80", "#23238E80",
  52. "#77B9B780", "#C8373780", "#EEDD8280", "#A4D3EE80", "#A478AB80",
  53. "#01AEF080", "#8DC73F80", "#F8931D80", "#FFAAAA80", "#8C8C8C",
  54. "#D38D5F80", "#23238E80", "#77B9B780", "#C8373780", "#EEDD8280",
  55. "#A4D3EE80", "#A478AB80", "#01AEF080", "#8DC73F80", "#F8931D80",
  56. "#FFAAAA80", "#8C8C8C", "#D38D5F80", "#23238E80", "#77B9B780",
  57. "#C8373780", "#EEDD8280", "#A4D3EE80", "#A478AB80", "#01AEF080",
  58. "#8DC73F80", "#F8931D80", "#FFAAAA80", "#8C8C8C", "#D38D5F80",
  59. "#23238E80", "#77B9B780", "#C8373780", "#EEDD8280")), .Names = c("mydate",
  60. "variable", "value", "fill"), row.names = c(NA, 96L), class = "data.frame")
  61.  
  62. myvals <- mymelt[mymelt$mydate == mymelt$mydate[nrow(mymelt)],] ## last date in mymelt should always be same as plotenddate as we subset earlier
  63. mymelt <- within(mymelt, variable <- factor(variable, as.character(myvals[order(myvals$value, decreasing = T),]$variable), ordered = TRUE))
  64.  
  65. p <- ggplot(mymelt, aes(x = mydate, y = value)) +
  66. geom_area(aes(fill = variable), position = "stack") +
  67. facet_wrap(~ variable, ncol = 4) +
  68. theme(axis.text.x = element_text(size = 8, angle = 90, colour = "grey50")) +
  69. theme()
  70. print(p)
  71.  
  72. p <- ggplot(mymelt, aes(x = mydate, y = value)) +
  73. geom_line(lwd=0.3) +
  74. facet_grid(. ~ variable) +
  75. theme_bw() +
  76. theme(axis.text.x = element_text(size = 5, angle = 90),
  77. axis.text.y = element_text(size = 8),
  78. axis.title.x = element_text(vjust = 0),
  79. axis.ticks = element_blank(),
  80. panel.grid.minor = element_blank())
  81. print(p)
  82. ggsave(plot=p, filename="plot.png", width = 8, height = 4)
Add Comment
Please, Sign In to add comment