3

我有一个输入需求模型，数据不平衡。因变量是yijt，其中i表示生产函数的投入(i=1,2,3)，j表示公司(j=1,..,21)，t表示时间(ti=1,..,Ti)。所以，因为面板是不平衡的，有三个方程，在 winbugs 中使用它的正确格式是什么？

我尝试以 STATA 长格式 3*107 的形式使用它，其中索引有点不同。行 ( i) 表示输入，列 ( j) 表示面板i*t。所以，没有明确的时间索引。我对列使用嵌套索引，将数据形成为面板。一小段代码：

for (i in 1:3){
     for(j in 1:107){
          log(mu[i,j]) <- a[i]+ u[firm[j]]+...

其中a[i]是一个特定于输入的自由参数，u是一个期望的特定于公司的误差分量，应该对公司进行估计（21）。为此，我使用嵌套索引，其中模型中 u 的表达式超过列u[firm[j]]，而先验 foru超过公司u[k]。

这是关于嵌套索引的问题。如果模型包含许多协变量，这些协变量也是时间相关的，我是否应该对所有这些协变量使用嵌套索引，否则 Winbugs 不会将其作为面板？

我尝试对所有多维变量使用多维数组。为了平衡很多 NA-s 被包括在内。例如，因变量 yijt 是一个 3 维数组y[ , , ]，（使用 list 命令），但由于某种原因，Winbugs 无法识别 NA-s？要点是后验分布不收敛，数据格式化可能是原因之一

4

1 回答 1

1

您最好将数据保留为长格式（y、i、j 和 t 列）和每个响应的行（如果缺少数据，y 将具有 NA）。

多级建模中心有一个很好的介绍性指南，关于使用 WinBUGS 以这种格式运行这些类型的模型[这里]

于 2013-03-25T13:54:09.940 回答