runmod-8.7.R

library(cmdstanr)

datafile <- 'concord-8.rds'
stanfile <- 'concord-8.7.stan'
outfile <- 'mod8.7.rds'

cmdstan2rstan <- function(stanfile, indata) {
  mod <- cmdstan_model(stanfile)
  mx <- mod$sample(data=indata, seed=8675309, chains=4, parallel_chains=4)
  rstan::read_stan_csv(mx$output_files())
}

indata <- readRDS(datafile)
m8.7 <- cmdstan2rstan(stanfile, indata)
saveRDS(m8.7, outfile)