2

我正在尝试用贝叶斯拟合分层负二项式模型。尽管我的数据是专有的，但我能够使用人造黄油数据集重新创建相同的错误。我得到的错误如下：

> look <- rhierNegbinRw(Data = list(regdata = dat1), Mcmc = list(R = 1000, 
nprint = 100))
Z not specified - using a column of ones instead
Error in alpha <= 0 : 
  comparison (4) is possible only for atomic and list types

我按如下方式设置了模拟数据（回归完全是荒谬的——只是想让它工作）：

data(margarine)
chpr <- margarine$choicePrice
chpr$hhid <- as.factor(chpr$hhid)
N <- nlevels(chpr$hhid)
dat1 <- vector(mode = "list", length = N)
for (i in 1:N) {
  dat1[[i]]$y <- chpr[chpr$hhid==levels(chpr$hhid)[i], "PPk_Stk"]
  dat1[[i]]$X <- model.matrix( ~ choice + PBB_Stk, 
                          data = chpr[chpr$hhid == levels(chpr$hhid)[i], ])  

}

我将非常感谢对此问题的任何见解。

4

0 回答 0