\begin{Schunk}
\begin{Sinput}
 svn <- "OLS_reg.stan"
 rt1 <- stanc(file = file.path(sdir, sfn))
 sm1 <- stan_model(stanc_ret = rt1, verbose = FALSE)
 stanm0 <- sampling(sm1, data = list(N = nrow(dat), K = ncol(X), 
      y = y, X = X), pars = c("beta", "sigma"),
      chains = 4, iter = 2000, warmup = 1000)
\end{Sinput}
\end{Schunk}
