Hi R users,

I am trying to call openbugs from R.  And I got the following error message:

~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
model is syntactically correct
expected the collection operator c error pos 8 (error on line 1)
variable ww is not defined in model or in data set
[1] "C:\\DOCUME~1\\maomao\\LOCALS~1\\Temp\\RtmpqJk9R3/inits1.txt"
Initializing chain 1: model must be compiled before initial values loaded
model must be initialized before updating
model must be initialized before DIC an be monitored
Error in samplesSet(parametersToSave) :
  model must be initialized before monitors used
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

I did define variable "ww" in my data and model (they are listed below).  I
am not sure if this is due to some errors in my code (please see below) or
because openbugs cannot handle the model I am using.  In my model, y[i] also
depends on all other y[j]s.  Could you help me figure out the problem and
hopefully get the code to work?

Many thanks for your help.   --- Maomao

~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
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])

wwy[i]<- inprod(ww[i, 1:676] , y[])

woy[i]<- inprod(wo[i, 1:676] , y[])

wdy[i]<- inprod(wd[i, 1:676] , y[])

probit(p[i])<- rho_o * woy[i] + rho_d * wdy[i] + rho_w * wwy[i] + beta[1] +
beta[2] * cap2[i] + beta[3] * pol2[i] + beta[4] * cap1[i] + beta[5] *
pol1[i] + beta[6] * g[i]+ e[i]

                        }

# Priors

for (j in 1:6) {

beta[1:6] ~ dmnorm(mu[1:6], tau[1:6, 1:6])

                   }

rho_o ~ dunif(-1,1)

rho_d ~ dunif(-1,1)

rho_w ~ dunif(-1,1)

for (i in 1:676) {

e[i] ~ dnorm(0, 1)

                        }

                }

	[[alternative HTML version deleted]]

