Question

我正在尝试使用JAGS来推断（随机）纯生育过程中的出生率。

在化学语言中，该模型等同于反应：X-> 2X，速率α* X（也可以看作链式反应的模型）

这是我用来生成进程的R代码（在固定时间）和用于推断参数alpha的jags代码。

library(rjags)

y <- 1; # Starting number of "individuals"
N <- 25 # number of time samplings
alpha <- 0.2 # per-capita birth rate
# Generate the time series
for(i in 2:N) {
  y<- c(y,y[i-1]+rpois(1,alpha*y[i-1]))
};

# The jags code
model_string <- "model{
  for(i in 2:N) {
    New[i] ~ dpois(alpha*y[i-1])
    y[i] <- y[i-1] + New[i]
  }
  alpha ~ dunif(0, 2)
}"

# Create and run the jags model
model <- jags.model(textConnection(model_string), data = list(y = y,N = N), n.chains = 3, n.adapt= 10000)
update(model, 5000); # Burnin for 10000 samples
mcmc_samples <- coda.samples(model, variable.names=c("alpha"), n.iter=5000)

当我运行代码时，出现以下错误：

Error in jags.model(textConnection(model_string), data = list(y = y, N = N),  : 
  RUNTIME ERROR:
Compilation error on line 4.
y[2] is a logical node and cannot be observed

我尝试了不同的东西，例如将alpha * y [i-1]放入一个新变量（比如lambda [i]）或者用New [i-1]改变对New [i]的调用，但没有任何效果。知道为什么会失败吗？另一种更聪明的方法吗？

提前谢谢。

用JAGS推断纯Borth过程

0 个答案: