Как реализовать анализ точек изменения в R для нескольких категорий и построить графики с пользовательскими метками по осям x и y?

#r

Вопрос:

У меня есть фрейм данных, который выглядит примерно так:

 set.seed(42)  ## for sake of reproducibility
n <- 6
dat <- data.frame(date=seq.Date(as.Date("2020-12-26"), as.Date("2020-12-31"), "day"),
                  category=rep(LETTERS[1:4], n/2),
                  daily_count=sample(18:100, n, replace=TRUE)
                  )
dat
#         date    category  daily_count
#1  2020-12-26        A          60
#2  2020-12-27        B          32
#3  2020-12-28        C          39
#4  2020-12-29        D          75
#5  2020-12-30        A          25
#6  2020-12-31        B          53
#7  2020-12-26        C          60
#8  2020-12-27        D          32
#9  2020-12-28        A          39
#10 2020-12-29        B          75
#11 2020-12-30        C          25
#12 2020-12-31        D          53
 

Я пытаюсь реализовать анализ точек изменения для каждой категории и отобразить результаты на одном графике (фасете). Я хочу, чтобы дата на оси x и daily_count на оси y и название категории в верхнем центре каждого графика. Тем не менее, я использую следующий код для выполнения этой работы

 for (i in 1:length(unique(dat$category))){
  cat_name <- unique(dat$category)[i]
  df_change = dat[dat$category == cat_name,]
  ansvar = cpt.var(df_change$daily_count)
  plot(ansvar)
}
 

И это дает мне аналогичные графики, как показано ниже (они вычисляются на основе фактических данных):
введите описание изображения здесь