model.file=file.path("model.txt")
n=58
inits=function(){list(tau=runif(4,0,10000),mu=runif(4,-10000,10000))}
params=c("mu","s2")
data=list("n","y")
output=bugs(data,inits,params,model.file,n.chains=3,n.iter = 3000,DIC = T)
#Diagnose
print(output)
#coda
output.coda=bugs(data,inits,params,model.file,n.chains=3,n.iter = 3000,codaPkg = T)
library(lattice)
library(coda)
result.coda=read.bugs(out.flat.coda)
summary(result.coda)
运行到上面那一步就出错了,后面的代码都没有运行,这个可能跟openbugs的存放路径有关吧,但是我也不太清楚


雷达卡


京公网安备 11010802022788号







