错误：索引超出范围，取变量子集

问题描述

我对JAGS很陌生。我的数据包括15周内5条溪流的幼虫计数（每周两次观察和溪流）。为此，我使用具有150行，2列和15个条目（y）的数组。我的目标是借助开放的N混合物模型估算每周的欧米茄。我还尝试了未标记的程序包，但无法弄清楚pcountopen（）如何工作。

sink("model1.txt") 
cat("
model {
# Priors
  lambda ~ dunif(0,5)
  omega ~ dunif(0,1)
  p ~ dunif(0,1)
# Likelihood
  # Biological (ecological) model for true abundance (State model)
  for (i in 1:R) {                  #loop over R sites (150)
  N[i,1] ~ dpois(lambda)            #Abundance
  for (k in 2:15) {                 #loop over weeks (15)
   N[i,k] ~ dpois(omega[k] * N[i,k-1])
   # Observation model for replicated counts
   for (j in 1:T){                  #loop over temporal reps (2)
     y[i,j,k-1] ~ dbin(p,N[i,k-1]) #Detection
     } #j
   } #k
 } #i
}
",fill = TRUE)
sink()

# Bundle and summarize data set 
R = nrow(y) #sites = 150
T = ncol(y) #replications per observation day = 2
win.data <- list(y = y,R = R,T = T)

# Initial values 
Nst <- apply(y,c(1,3),max) +1
inits <- function() {list(N = Nst)}

# Parameters monitored
params1 <- c("omega","lambda","p")

# MCMC settings 
ni <- 20 ; nt <- 1 ; nb <- 0 ; nc <- 8

# Call JAGS (ART 1 min) and summarize posteriors 
library(jagsUI)
fm1 <- jags(win.data,inits,params1,"model1.txt",n.chains = nc,n.thin = nt,n.iter = ni,n.burnin = nb)

这给了我以下错误： jags.model中的错误（文件=模型文件，数据=数据，初始化=初始化，n.chains = n.chains 、：运行时错误：第12行出现编译错误。超出欧米茄范围的索引

感谢您的帮助！

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）

jags r r unmarked-package