本文winbugs小白,最近开始初学这方面内容,希望大家能够帮忙解惑><!
想请教,如果想要load inits,初始值代码应该如何写?
下面是所写代码:
model{
for(i in 1:regions){
Y[i]~dpois(mu[i])
log(mu[i])<-log(E[i])+beta0+x[i]*beta1+u[i]+v[i]
SMR[i]<-Y[i]/E[i]
lambda[i]<-mu[i]/E[i]
v[i]~dnorm(0.0,tau.h)
}
u[1:regions]~car.normal(adj[],weights[],num[],tau.c)
for(k in 1:sumNumNeigh){
weights[k]<-1
}
beta0~dnorm(0.0,1.0E-5)
beta1~dnorm(0.0,1.0E-5)
tau.h ~dgamma(1.0E-1,1.0E-3)
tau.c ~dgamma(1.0E-1,1.0E-1)
}
#data
list(regions=214,sumNumNeigh=1186,num=c(4,5,5,5,7,6,6,6,4,4,5,7,6,6,2,6,4,3,9,7,4,7,5,6,5,5,6,8,6,6,7,6,3,4,7,6,5,4,7,7,6,4,7,8,3,9,5,3,6,7,6,5,2,4,10,0,5,5,6,3,0,5,4,5,7,6,6,6,5,6,7,6,5,5,6,5,5,6,6,6,6,8,7,7,8,6,7,7,4,5,6,5,5,6,3,6,7,5,7,3,8,4,6,8,6,5,5,6,10,5,1,6,5,4,6,5,5,6,5,4,6,4,5,8,7,7,5,3,5,9,3,5,8,7,5,5,7,9,4,5,5,3,5,4,5,3,2,5,3,7,5,5,7,10,5,5,7,7,6,5,7,6,9,5,6,7,6,6,4……
后续数据过多就没有贴出,如有需要可以私发。