在JAGS中拟合GLMM时遭遇参数不匹配错误的技术求助

核心错误原因 #

触发Non-conforming parameters in function错误的直接原因是**Nre的定义错误**：你将Nre设为unique(temp$site)（索引值向量，如示例中的c(1,2,3)），但JAGS中for (i in 1:Nre)要求Nre是标量（即随机效应的总数量），向量无法被解析为有效的循环范围，导致参数维度不匹配。

此外代码还存在其他与模型目标不符的问题，一并整理修复：

问题梳理与修复步骤 #

1. 修正随机效应数量的定义 #

将Nre从索引值向量改为随机效应的数量，同时补充year随机效应的索引与数量（你目标模型包含site和year两个随机效应）：

count <- sample(0:20, 100, replace = TRUE)
rain <- rnorm(100, 0, 0.01)
ndvi <- rnorm(100, 0.2, 0.01)
burn <- rnorm(100, 0, 0.01)
site <- c(rep(1, 30), rep(2, 20), rep(3,50))
year <- c(rep(seq(2000:2009), 10))
temp <- as.data.frame(cbind(count, rain, ndvi, burn, site, year))

# 创建模型矩阵与变量
X <- model.matrix(~ rain + ndvi + burn, data = temp)
K <- ncol(X)
# site相关
re_site <- as.numeric(temp$site)
Nre_site <- length(unique(temp$site)) # 改为site的数量（标量）
# year相关：将原始年份转换为连续整数索引
re_year <- as.numeric(factor(temp$year))
Nre_year <- length(unique(temp$year)) # year的数量（标量）

win.data1 <- list(
  Y = temp$count, 
  X = X, 
  K = K, 
  N = nrow(temp), 
  Re_site = re_site, 
  Nre_site = Nre_site,
  Re_year = re_year,
  Nre_year = Nre_year
)

2. 修正JAGS模型代码 #

• 替换正态似然为Poisson似然（符合你拟合Poisson GLMM的目标），并使用log链接函数
• 正确拆分site和year的随机效应循环
• 删除无关的sigma.eps（Poisson模型不需要额外的残差标准差）

sink("lmm.text")
cat("
model {
  # 固定效应先验
  for (i in 1:K) {beta[i] ~ dnorm(0, 0.0001) }
  
  # site随机效应先验
  for (i in 1:Nre_site){
    a[i] ~ dnorm(0, tau.site)
  }
  # year随机效应先验
  for (i in 1:Nre_year){
    b[i] ~ dnorm(0, tau.year)
  }
  
  # 精度与标准差转换
  tau.site <- 1/(sigma.site*sigma.site)
  tau.year <- 1/(sigma.year*sigma.year)
  sigma.site ~ dunif(0.001, 10)
  sigma.year ~ dunif(0.001, 10)
  
  # Poisson似然（log链接）
  for (i in 1:N){
    Y[i] ~ dpois(mu[i])
    log(mu[i]) <- inprod(beta[], X[i,]) + a[Re_site[i]] + b[Re_year[i]]
  }
}
", fill = TRUE)
sink()

3. 修正初始值函数 #

确保初始值生成时使用的是随机效应的数量（标量），而非向量：

inits1 <- function(){
  list(beta = rnorm(K, 0, 0.01),
       a = rnorm(Nre_site, 0, 1),
       b = rnorm(Nre_year, 0, 1), 
       sigma.site = runif(1, 0.001, 10),
       sigma.year = runif(1, 0.001, 10))
}

4. 修正JAGS运行路径 #

确保模型文件路径正确（如果lmm.text在当前工作目录，无需src/前缀）：

params1 <- c("beta", "a", "b", "sigma.site", "sigma.year")

# 运行JAGS
J0 <- jags(data = win.data1,
          inits = inits1,
          parameters = params1,
          model.file = "lmm.text", # 修正路径
          n.thin = 10,
          n.chains = 3,
          n.burnin = 4000,
          n.iter = 5000)

J1 <- update(J0, n.iter = 10000, n.thin = 10)
out <- J1$BUGSoutput

后续转负二项模型的提示 #

若后续要替换为负二项模型，只需修改似然部分并添加离散参数的先验：

# JAGS模型中替换似然块
for (i in 1:N){
  Y[i] ~ dnegbin(p[i], r)
  mu[i] <- exp(inprod(beta[], X[i,]) + a[Re_site[i]] + b[Re_year[i]])
  p[i] <- r / (r + mu[i])
}
# 额外添加r的先验
r ~ dgamma(0.01, 0.01)

内容的提问来源于stack exchange，提问作者asmart2