Jags有线,由于数据不同,偶尔会重新定义节点错误

时间:2017-07-18 04:38:42

标签: jags

我基于此创建了我的代码: http://users.aims.ac.za/~mackay/BUGS/Manual05/Examples1/node29.html

现在我使用不同的种子来模拟数据。奇怪的是,有些种子在第18行错误上给了我一个重新定义的节点dN [1,1],而其他种子却没有。有人可以帮忙吗?顺便说一句,为什么18号线上的dN [1,1]位于第一位? Jags如何计算线条?

错误消息:
“ RUNTIME错误:

第18行的编译错误。

尝试重新定义节点dN [1,1] “

  bugsmodel <- "
  # Set up data
  data{
    for(i in 1:N) 
    {
      for(j in 1:bigt) 
      {
      Y[i,j] <- step(obs.t[i] - t[j] + eps)
      dN[i, j] <- Y[i, j] * step(t[j + 1] - obs.t[i] - eps) * fail[i]
      }
    }
  }
  # Model
  model
  {
  for(i in 1:N){
    betax[i,1] <- 0
    for(k in 2:(p+1)){
    betax[i,k] <- betax[i,k-1] + beta[k-1]*x[i,k-1]
    }
  }
  for(j in 1:bigt) {
    for(i in 1:N) {
    dN[i, j] ~ dpois(Idt[i, j]) # Likelihood
    Idt[i, j] <- Y[i, j] * exp(betax[i,p+1]) * dL0[j] # Intensity
    }
    dL0[j] ~ dgamma(mu[j], c)
    mu[j] <- dL0.star[j] * c # prior mean hazard
  }
  c <- 0.001
  r <- 0.1
  for (j in 1 : bigt) {
    dL0.star[j] <- r * (t[j + 1] - t[j])
  } 
  for(k in 1:p){
    beta[k] ~ dnorm(0.0,0.000001)
  }
  }"

0 个答案:

没有答案