2014-01-07 32 views
1

我對如何解析n = 0時for「for(i in 1:n)」中的for循環順序感興趣。1:n給出c(1,0),所以for循環運行索引1和0,還是for循環,然後跳過?當上界爲零時,如何解析JAGS中的循環?

一個例子:假設我有來自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/

+0

也許,看到這些的行爲:'對(我在1:0)打印( 「循環」)','爲(我在seq_len(0))打印(「循環」)''。 –

回答

-1

容易測試:

for (i in 1:0) print(i) 
[1] 1 
[1] 0 
+1

感謝但不是在R中,我的意思是JAGS! – jrnew