我目前正在研究一个自学问题,并且遇到了一些问题,这些问题能够使代码正常运行,甚至可以针对该事实进行设置。目的是在给定两个c值的情况下,能够获得每个模型的贝叶斯因子,DIC和W AIC
#Load the data
Y1=563
N1=2820
Y2=10
N2=27
library('rjags')
#Model 1 with c=1
Y1~dpois(N1*lambda1)
lambda1=dunif(0,1)
Y1~dpois(N1*lambda2)
lambda2=dunif(0,1)
#Model 2 with c=1
lambda0=dunif(0,1)
Y1~dpois(N1*lambda0)
Y2~dpois(N2*lambda0)
#Model 1 with c=10
Y1~dpois(N1*lambda1)
lambda1=dunif(0,10)
Y1~dpois(N1*lambda2)
lambda2=dunif(0,10)
#Model 2 with c=10
lambda0=dunif(0,10)
Y1~dpois(N1*lambda0)
Y2~dpois(N2*lambda0)
但是,当我输入此代码时,出现以下错误,
> lambda1=dunif(0,1)
Warning message:
In dunif(0, 1) : NaNs produced