R2OpenBUGS运行报错：expected right parenthesis error 求助

错误根源与修复要点 #

• dnorm参数误用：OpenBUGS中dnorm(mean, precision)的第二个参数是精度（方差的倒数），而非标准差。原代码直接传入Sig.Meas[i]（标准差），需替换为精度形式pow(Sig.Meas[i], -2)。
• 变量维度未声明：模型中Phi、M是对应每个测量样本（i=1_{I）和能量区间（j=1}J）的二维变量，必须明确维度后再赋值，否则OpenBUGS会解析失败。
• 循环逻辑错误：原代码在i循环内重复覆盖M、Phi，导致每个样本的计算结果被覆盖。需针对每个i和j定义元素，确保每个样本对应独立的能量区间计算。
• 初始值不符合先验：a_th的先验是dunif(0,1)，初始值不能取负数，且需生成与I=5匹配的5个初始值，而非1个。

修复后的完整代码 #

library(R2OpenBUGS)

# 修复后的BUGS模型
lineModel <- function() {
  # 声明二维变量维度，对应I个样本、J个能量区间
  for (i in 1:I) {
    for (j in 1:J) {
      # 计算每个样本-能量区间的Phi值
      Phi[i,j] <- a_th[i] * (En[j]/pow(T0,2)) * exp(-En[j]/T0)
      # 计算响应与Phi的乘积
      M[i,j] <- Resp[j] * Phi[i,j]
    }
    # 求和得到每个样本的预测测量值
    Meas_bar[i] <- sum(M[i,])
    # 使用精度作为dnorm的第二个参数
    Meas[i] ~ dnorm(Meas_bar[i], pow(Sig.Meas[i], -2))
  }
  
  # 先验分布定义
  for (i in 1:I) {     
    a_th[i] ~ dunif(0, 1) 
  }
}

# 数据部分（无需修改）
lineData <- list(I = 5, J = 21, T0 = 0.025, 
                 En = c(1.0000E-09, 1.2589E-09, 1.5849E-09, 1.9953E-09, 
                        2.5119E-09, 3.1623E-09, 3.9811E-09, 5.0119E-09, 
                        6.3096E-09, 7.9433E-09, 1.0000E-08, 1.2589E-08, 
                        1.5849E-08, 1.9953E-08, 2.5119E-08, 3.1623E-08, 
                        3.9811E-08, 5.0119E-08, 6.3096E-08, 7.9433E-08, 
                        1.0000E-07), 
                 
                 Meas = c(31, 84, 100, 95, 75), 
                 Sig.Meas = c(0.4, 0.4, 0.6, 0.8, 0.8),
                 
                 Resp = c(7.156E+00,6.971E+00,6.761E+00,6.527E+00,6.271E+00, 5.996E+00,
5.705E+00,5.402E+00,5.092E+00,4.777E+00,4.464E+00,4.153E+00,3.850E+00,3.555E+00,
3.273E+00,3.003E+00,2.749E+00,2.508E+00,2.284E+00,2.074E+00,1.882E+00))

# 修复后的初始值生成函数
lineInits <- function() {
  list(a_th = runif(5, 0, 1))  # 生成5个符合先验的初始值
}

# 运行模型
lineout <- bugs(lineData, 
                lineInits,
                parameters.to.save = c("a_th"),
                lineModel,
                n.burnin = 1500,
                n.chains = 1, 
                n.iter = 10000,
                working.directory = "F:/AAA_PROJECTs/Projet_DPR/Code_Scripts/Bayesian_Code/R_Codes/Proj02")

# 查看输出结果
print(lineout)

额外注意事项 #

• 建议将n.chains设置为2及以上，便于后续用收敛诊断工具（如coda包）判断模型是否收敛。
• 若使用model.file参数，需将模型代码保存为纯文本文件，注意OpenBUGS语法与R的细微差异（如幂运算pow()兼容性更好）。
• 确保working.directory指定的路径存在，OpenBUGS需要在该目录下生成临时日志和模型文件。

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