#Load data
data(MCMC_data)
#Extract MCMC chains
ex <- MCMCchains(MCMC_data)
apply(ex, 2, mean)
#Extract MCMC chains for just 'beta' parameters
ex2 <- MCMCchains(MCMC_data, params = 'beta')
apply(ex2, 2, mean)
#Just 'beta[1]', 'beta[4]', and 'alpha[3]'
ex3 <- MCMCchains(MCMC_data, params = c('beta[1]', 'beta[4]', 'alpha[3]'),
ISB = FALSE, exact = TRUE)
apply(ex3, 2, sd)
Run the code above in your browser using DataLab