首页
学习
活动
专区
圈层
工具
发布
社区首页 >问答首页 >,但进入了以下内容:错误解析模型文件:第5行"=“附近的语法错误

,但进入了以下内容:错误解析模型文件:第5行"=“附近的语法错误
EN

Stack Overflow用户
提问于 2022-03-27 20:05:36
回答 1查看 91关注 0票数 0

在运行jags模型时,我得到了以下错误消息:jags.model中的模块glm加载错误(model.file,data = data,inits = init.values,n.chains = n.chains,:

错误解析模型文件:第5行"=“附近的语法错误

这是我的密码:

代码语言:javascript
复制
install.packages('R2jags')
install.packages('rjags')
library(rjags)
library('R2jags')
library(lattice)
binomial.model.JAGS = function(){
  y ~ dbin(p,n)
  p=lambda*mu+rho*(1-mu)
  lambda ~ dunif(min = 0.2,max = 1.4)
  mu ~ dunif(min = 0,max = 1)
  rho ~ dunif(min = 0.1,max = 1.7)
}
n = 100000
y = 30000
data.JAGS = list(y = y, n = n)

inits.JAGS = list(list(lambda=0.8,mu=0.5,rho=0.9))
para.JAGS = c("p", "lambda", "mu", "rho")
fit.JAGS = jags(
                data=data.JAGS,inits=inits.JAGS,
                parameters.to.save=para.JAGS,
                n.chains=1,
                n.iter=9000,
                n.burnin = 1000,
                model.file = binomial.model.JAGS)

我真的不知道哪里出了问题。有人能帮我吗?非常感谢!!

EN

回答 1

Stack Overflow用户

发布于 2022-04-07 11:44:23

只要去掉你前科中所有的论点命名,你就可以走了。

代码语言:javascript
复制
binomial.model.JAGS = function(){
  y ~ dbin(p,n)
  p=lambda*mu+rho*(1-mu)
  lambda ~ dunif(0.2,1.4)
  mu ~ dunif(0,1)
  rho ~ dunif(0.1,1.7)
}
n = 100000
y = 30000
data.JAGS = list(y = y, n = n)

inits.JAGS = list(list(lambda=0.8,mu=0.5,rho=0.9))
para.JAGS = c("p", "lambda", "mu", "rho")
fit.JAGS = jags(
  data=data.JAGS,inits=inits.JAGS,
  parameters.to.save=para.JAGS,
  n.chains=1,
  n.iter=9000,
  n.burnin = 1000,
  model.file = binomial.model.JAGS)

         mu.vect sd.vect   2.5%    25%    50%    75%  97.5%
lambda     0.286   0.025  0.239  0.270  0.287  0.298  0.351
mu         0.936   0.064  0.758  0.930  0.957  0.977  0.996
p          0.300   0.001  0.297  0.299  0.300  0.301  0.303
rho        0.794   0.459  0.114  0.355  0.774  1.166  1.611
deviance  12.799   1.403 11.791 11.883 12.230 13.105 16.491

DIC info (using the rule, pD = var(deviance)/2)
pD = 1.0 and DIC = 13.8
DIC is an estimate of expected predictive error (lower deviance is better).
票数 0
EN
页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持
原文链接:

https://stackoverflow.com/questions/71639893

复制
相关文章

相似问题

领券
问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档