utils 1.2 KB

123456789101112131415161718192021222324252627282930313233
  1. #' Get layers to overlay significant changepoints on prophet forecast plot.
  2. #'
  3. #' @param m Prophet model object.
  4. #' @param threshold Numeric, changepoints where abs(delta) >= threshold are significant. (Default 0.01)
  5. #' @param cp_color Character, line color. (Default "red")
  6. #' @param cp_linetype Character or integer, line type. (Default "dashed")
  7. #' @param trend Logical, if FALSE, do not draw trend line. (Default TRUE)
  8. #' @param ... Other arguments passed on to layers.
  9. #'
  10. #' @return A list of ggplot2 layer.
  11. #'
  12. #' @import ggplot2
  13. #'
  14. #' @examples
  15. #' \dontrun{
  16. #' plot(m, fcst) + layer_changepoints(m)
  17. #' }
  18. #'
  19. #' @export
  20. layer_changepoints <- function(m, threshold = 0.01, cp_color = "red",
  21. cp_linetype = "dashed", trend = TRUE, ...) {
  22. layers <- list()
  23. if (trend) {
  24. trend_layer <- geom_line(aes_string("ds", "trend"), color = cp_color, ...)
  25. layers <- append(layers, trend_layer)
  26. }
  27. signif_changepoints <- m$changepoints[abs(m$params$delta) >= threshold]
  28. cp_layer <- geom_vline(xintercept = as.integer(signif_changepoints),
  29. color = cp_color, linetype = cp_linetype, ...)
  30. layers <- append(layers, cp_layer)
  31. layers
  32. }