jags: coda包，如何选择多个参数

Question

我使用r2jags运行jags并进行采样。然后，我使用coda包检查收敛性。由于我有许多主题和许多参数，我需要手动键入它们作为下面的代码。是否有更简单的方法，例如，为每个主题选择gamma？或者，只从主题1到10选择gamma，而不需要手动键入它们？

codaSamples = as.mcmc.list(samples$BUGSoutput)

acfplot(codaSamples[,c('gamma[1]','gamma[2]','gamma[3]','gamma[4]','gamma[5]','gamma[6]',
                         'gamma[7]','gamma[8]','gamma[9]','gamma[10]','gamma[11]','gamma[12]',
                         'gamma[13]','gamma[14]','gamma[15]')],lag.max=1000)

Answer 1

一个可复制的例子会很好，但想必是

cols <- sprintf("gamma[%d]",1:15); codaSamples[,cols]

会做你想做的吗？或

cols <- grep("^gamma",colnames(codaSamples)); codaSamples[,cols]