1

我对当 n = 0 时如何解析“for (i in 1:n)”中的 for 循环序列感兴趣。1:n 是否给出 c(1, 0)，因此 for 循环针对索引 1 和0，还是跳过for循环？

一个例子：假设我有来自 10 个组的观察，我有默认和特殊观察，我分别用正态分布和 t 分布建模。每个组可以有两种类型的观察，或者只有默认的，或者只有特殊的观察。是否可以通过这种方式设置模型：

model {
  for (g in 1:ngroups) {
    for (i in 1:ndefault.g[g]) { # ndefault.g[g] may be 0 for certain groups
      y.gi[g, i] ~ dnorm(yhat.gi[g, i], tauy.gi[g, i])
    }
    for (i in 1:nspecial.g[g]) { # nspecial.g[g] may be 0 for certain groups
      y.gi[g, i] ~ dt(yhat.gi[g, i], tauy.gi[g, i], dft)
    }
  }
...
}

谢谢！

交叉发布于：https ://sourceforge.net/p/mcmc-jags/discussion/610037/thread/d13fd9a2/

4

1 回答 1

-1

易于测试：

for (i in 1:0) print(i)
[1] 1
[1] 0
于 2014-01-07T06:45:45.280 回答