library(rstan)
library(readxl)
can_tbill <- read_excel("~/can-tbill.xlsx")


dat1 = list(N = length(can_tbill$date),
            P = 4,
            y = as.matrix(can_tbill[,3:6]),
            tau = c(1,3,6,12))

mod1 = stan_model("dynamic-nelson-siegel-ar1.stan")

fit = sampling(mod1,
               dat1,
               iter = 5000)
print(fit)