У меня есть выход runjags
модели в R как mcmc.list
. Ниже приведен код для генерации 3 цепочек из 1000 образцов. Я хотел бы урезать все 12 цепей до последних 400 образцов. Я могу разделить цепочки и сохранить матрицы вывода цепочек в списке, но это больше не mcmc.list
, и я не могу понять, как превратить его обратно в mcmc.list.
Вот некоторые данные для запуска модели runjags
и преобразования выходных данных в mcmc.list
:
y <- rnorm(100)
jags.model ="
model {
#model
for (i in 1:N){
y[i] ~ dnorm(y.hat[i], tau)
y.hat[i] <- m0
}
#priors
m0 ~ dnorm(0, .0001)
tau <- pow(sigma, -2)
sigma ~ dunif(0, 100)
}
"
jags.data <- list(y = y, N = length(y))
jags.out <- runjags::run.jags(jags.model,
data = jags.data,
n.chains = 3,
adapt = 100,
burnin = 100,
sample = 1000,
monitor = c('m0'))
z <- coda::as.mcmc.list(jags.out)