使用jags遍历(子集)

Question

我有一个10000行12列的大型数据帧(discountdataset)。这些列包含不同的变量。前210行表示主题1(还有一列带有"subject1")，接下来的210行表示主题2，依此类推。

我想使用jags和一个循环函数来遍历数据帧中的所有52个主题，并为每个主题分配一个函数。我的代码如下所示：

#subsetting the dataframe by the variable subjectid
subsetdiscount <- split(discountdataset, as.factor(discountdataset$subjectid))

在这里，我的计划是循环并将以下jags函数分配给子集中的所有主题)，但是，它不起作用。我想我的错误是我想要传递给jags的变量"nt", "Choice"没有正确定义，或者没有更新。

library(rjags)

for (i in 1:length(subsetdiscount))
{

nt <- nrow (subsetdiscount)
Choice <- subsetdiscount$choice
amountSS <- subsetdiscount$val_basic
amountLL <- subsetdiscount$val_d
delayDIFF <- subsetdiscount$delay
con <- subsetdiscount$condition


data <- list("nt", "Choice", "amountSS", "amountLL", "delayDIFF", "con") #      to be passed on to JAGS

myinits <- list(
list(k = (c(0.01, 0.01))),
list(temp = (c(6, 6))))


parameters <- c("k", "temp")


samples <- jags(data, inits=myinits, parameters,
             model.file ="singlesubmodel_Ben_roundedchoice.txt", n.chains=2,         n.iter=20000, 
     n.burnin=1, n.thin=1, DIC=T)

Answer 1

尝试：

library(rjags)
library(R2jags)

subsetdiscount <- split(discountdataset, as.factor(discountdataset$subjectid))

output_models <- lapply(subsetdiscount, function(x) {
  nt <- nrow(x)
  Choice <- x$choice
  amountSS <- x$val_basic
  amountLL <- x$val_d
  delayDIFF <- x$delay
  con <- x$condition
  data <- list("nt", "Choice", "amountSS", "amountLL", "delayDIFF", "con") #      to be passed on to JAGS
  myinits <- list(list(k = (c(0.01, 0.01))),
    list(temp = (c(6, 6))))
  parameters <- c("k", "temp")
  samples <- jags(data, inits=myinits, parameters,
                  model.file ="singlesubmodel_Ben_roundedchoice.txt", 
                  n.chains=2, n.iter=20000, 
                  n.burnin=1, n.thin=1, DIC=T)
  return(samples)
})

output_models应该是一个列表，其中包含拆分主数据集所依据的每个因子的输出。请注意，如果没有任何提供的数据，测试这一点是相当困难的。因此，如果这不起作用，您可能需要提供一些数据用于测试。

我希望它能帮上忙。