runmod-3.1.R

library('rstan')

datafile <- 'concord-3-n20000.rds'
stanfile <- 'concord-2.stan'         # Stan model for model 3 is the same as for model 2.
outfile <- 'mod3.1.rds'

indata <- readRDS(datafile)

mod <- stan(stanfile, data=indata, cores=4,
            seed=867-5309, iter=2000, warmup=1000)

saveRDS(mod, outfile)

message('Outfile is ', outfile)