Почему последнее значение всегда отсутствует для последующего распределения вероятностей с использованием ‘bcp’ для определения точки изменения?

#r #time-series #bayesian

#r #временные ряды #байесовский

Вопрос:

Я использую пакет bcp для определения точки изменения и задаюсь вопросом, почему последнее значение распределения вероятностей всегда «NA». Мне интересно, возможно, алгоритму bcp требуется точка данных после определенной даты, чтобы предоставить значение. Мой учебный период заканчивается на неделе 12-29-2019, и мне интересно, нужна ли мне дополнительная неделя данных, чтобы действительно охватить учебный период в течение недели 12-29-2019.

Любое объяснение будет с благодарностью. Это происходило каждый раз, когда я использую функцию bcp, но я включаю воспроизводимый пример ниже:

Упаковка:

 library(gtrendsR)
library(tidyverse)
library(bcp)
  

Получение образца данных Google Trends:

 data = gtrends('stackoverflow', geo = "US-OR", time = "2018-04-01 2019-12-31")$interest_over_time %>% 
  select(-time:-category)
  

Использование пакета bcp:

 bcp = bcp(data$hits)
  
prob = bcp$posterior.prob
prob = as.data.frame(prob) 

bcp_dataframe = cbind(data, prob) %>% 
    select(date, prob)

  

Результирующий фрейм данных:

введите описание изображения здесь

Комментарии:

1. Я предлагаю отправить авторам электронное письмо. Я не смог подтвердить это в документации, но я предполагаю, что апостериорная вероятность, излучаемая для n временной точки, представляет вероятность изменения между n и n 1 ; следовательно, последняя запись не определена.