在R中使用贝叶斯因子,DIC和W AIC选择模型

时间:2019-11-15 17:41:14

标签: r bayesian rjags

我目前正在研究一个自学问题,并且遇到了一些问题,这些问题能够使代码正常运行,甚至可以针对该事实进行设置。目的是在给定两个c值的情况下,能够获得每个模型的贝叶斯因子,DIC和W AIC

Problem involving the code

#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

0 个答案:

没有答案