I am using a random intercept model with SITEID as random and NAUSEA as outcome.
Thanks. Anamika I tried using a dataset without missing values and changed my model statement accordingly but still get the same error. Follwoing in an excerpt. > anal.data <- read.table("nausea.txt", header=T, sep="\t") > list(names(anal.data)) [[1]] [1] "SITEID" "NAUSEA" > #anal.data <- read.csv("simuldat.csv", header=T) > attach(anal.data) The following object(s) are masked from anal.data ( position 3 ) : NAUSEA SITEID The following object(s) are masked from anal.data ( position 4 ) : NAUSEA SITEID > data.bugs <- list("NAUSEA"=NAUSEA,"SITEID"=SITEID,"n.samples"=n.samples,"n.sites"=n.sites,"n.params"=n.params) > bugsData(data.bugs, fileName = "nauseadata.txt") > inits.bugs <- list("alpha"=rep(0,n.sites), "tau"=1) > bugsInits(list(inits.bugs), fileName = "nauseainit.txt") > modelCheck("nausea_random.txt") # check model file model is syntactically correct > modelData("nauseadata.txt") # read data file expected variable name error pos 98349 *MODEL* model { for (i in 1:n.samples) {NAUSEA[i] ~ dbin(p[i],1) logit(p[i]) <- alpha[SITEID[i]]} #for(k in 1:n.params) #{b[k]~ dnorm(0.0,tau)} for (j in 1:n.sites) {alpha[j]~dnorm(0.0,1.0E-10)} tau ~ dgamma(0.001,0.001) } *Dataset:* SITEID NAUSEA 1 0 1 1 1 1 1 0 1 1 1 1 1 0 1 1 1 1 1 1 1 0 1 1 1 0 1 1 1 1 1 1 1 0 1 0 1 0 1 1 -Anamika [[alternative HTML version deleted]] ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.