jags: coda 包,如何 select 多个参数

jags: coda package, how to select multiple parameters

我使用 r2jags 到 运行 锯齿并进行采样。然后,我使用 coda 包来检查收敛性。由于我有很多主题和参数,我需要手动输入它们,如下面的代码。有没有更简单的方法,例如,为每个主题选择 gamma?或者,只 select gamma 从主题 1 到 10 而无需手动输入?

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)

一个可重现的例子会很好,但大概

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

会做你需要的吗?或者

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