#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
; следовательно, последняя запись не определена.