OpenBugs返回

如何解决OpenBugs返回

我正在尝试编写Makridakis and Hyndman的书中的示例。 但是在某些时候，我们使用了重复性术语，但我不知道该如何解决。

当我尝试>>> level [i] -level [i-1]时，模型停止运行。 实际上，当我将“帐篷”与“水平”连接时，它不起作用。我不明白为什么，因为我可以将“水平”与“帐篷”连接起来，但不能反过来。

library(R2openbugs)

Dados <- c(362,385,432,341,382,409,498,387,473,513,582,474,544,681,557,628,707,773,592,627,725,854,661)

anos <- c(2015,2015,2016,2017,2018,2019,2020,2020)

meses <- c(1,2,3,4,1,4)

N       <- 24

s <- 4

data        <-  list(N = N,dados = Dados,meses = meses,anos = anos,s=s)

model1      <-  function()
{   
  for (i in 1:N) 
  {
    y[i] <- Consumo[i]
    y[i] ~ dnorm(mean[i],tau)
  }
  
  level[1] <- y[1]
  
  tent[1] <- y[1]

  mean[1] <- y[1]
    
  for (i in 2:N) 
  {
    level[i] <- (beta0*y[i])+((1-beta0)*(mean[i-1]+tent[i-1]))
    tent[i] <- (beta1*(level[i]-level[i-1]))+((1-beta1)*tent[i-1])
    mean[i] <- level[i]+tent[i]
   
  }
   
  beta0 ~ dnorm(3,20)
  beta1 ~ dnorm(3,20)
  tau ~ dgamma(1,1)
}