Displaying 1 result from an estimated 1 matches for "rho_o_0".
Did you mean:
rho_d_0
2010 Oct 15
1
Problem using BRugs
...~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
data<-list("y","cap2","pol2","cap1","pol1","g","wo","wd","ww","mu","tau")
inits<-function() {list(beta=beta0, rho_o=rho_o_0, rho_d=rho_d_0,
rho_w=rho_w_0)}
parameters<-c("beta", "rho_o", "rho_d", "rho_w")
probit.sim<-BRugsFit(data,inits,parameters,modelFile="spatial.openbugs.txt",numChains=1,nIter=2000)
# my model
model {
for (i in 1:676) {
y[i] ~ dbern(p[i...