## ----setup, include = FALSE---------------------------------------------------
knitr::opts_chunk$set(
  collapse = TRUE,
  comment = "#>"
)
library(BayesURTrend)

## ----random-walk--------------------------------------------------------------
set.seed(123)
y_rw <- cumsum(rnorm(100))
res_rw <- bayes_ur_test(y_rw, p = 1, k = 1)
print(res_rw)
summary(res_rw)
plot(res_rw)

## ----stationary-ar1-----------------------------------------------------------
set.seed(456)
y_stat <- numeric(100)
for (t in 2:100) {
  y_stat[t] <- 0.6 * y_stat[t - 1] + rnorm(1)
}

res_stat <- bayes_ur_test(y_stat, p = 1, k = 1)
print(res_stat)
summary(res_stat)
plot(res_stat)