I'm using WinBugs 1.1.1 98 and having trouble specifying a linear mixed
model. There are 5 vendors, with 5 types of used cars, and up to 6 of
each of those types. I'm modelling the asking price of each car. I get
through the model check, and the data load, but it fails on the compile
with "variable not defined or not data for n" with respect to no line in
particular. Any ideas?
model {
for (vend in 1:5) {
for (type in 1:5) {
for (rep in 1:6) {
n<-30*(vend-1)+6*(type-1)+rep
z[vend,type,1]<-z1[n]
z[vend,type,2]<-z2[n]
z[vend,type,3]<-z3[n]
z[vend,type,4]<-z4[n]
z[vend,type,5]<-z5[n]
x[vend,type,1]<-x1[n]
x[vend,type,2]<-x2[n]
x[vend,type,3]<-x3[n]
x[vend,type,4]<-x4[n]
x[vend,type,5]<-x5[n]
x[vend,type,6]<-x6[n]
avg.price[n]<-inprod(x[vend,type,], beta)+inprod(z[vend,type,], mu)
price[n]~dnorm(avg.price[n], tau)
}
}
}
tau~dgamma(1, 100)
sigma<-1/sqrt(tau)
omega ~ dwish(R[,], 10)
beta~dmnorm(b[], one[])
mu~dmnorm(zero[], omega)
}
z1[] z2[] z3[] z4[] z5[] x1[] x2[] x3[] x4[]
x5[] x6[] price[]
1 0 0 0 0 1 0 0
0 0 11 30.495
1 0 0 0 0 1 0 0
0 0 9 22.250
1 0 0 0 0 1 0 0
0 0 7 19.995
1 0 0 0 0 1 0 0
0 0 NA NA
1 0 0 0 0 1 0 0
0 0 NA NA
1 0 0 0 0 1 0 0
0 0 NA NA
1 0 0 0 0 0 1 0
0 0 9 17.995
1 0 0 0 0 0 1 0
0 0 8 17.495
1 0 0 0 0 0 1 0
0 0 NA NA
1 0 0 0 0 0 1 0
0 0 NA NA
1 0 0 0 0 0 1 0
0 0 NA NA
1 0 0 0 0 0 1 0
0 0 NA NA
1 0 0 0 0 0 0 1
0 0 6 15.695
1 0 0 0 0 0 0 1
0 0 NA NA
1 0 0 0 0 0 0 1
0 0 NA NA
1 0 0 0 0 0 0 1
0 0 NA NA
1 0 0 0 0 0 0 1
0 0 NA NA
1 0 0 0 0 0 0 1
0 0 NA NA
1 0 0 0 0 0 0 0
1 0 10 16.295
1 0 0 0 0 0 0 0
1 0 NA NA
1 0 0 0 0 0 0 0
1 0 NA NA
1 0 0 0 0 0 0 0
1 0 NA NA
1 0 0 0 0 0 0 0
1 0 NA NA
1 0 0 0 0 0 0 0
1 0 NA NA
1 0 0 0 0 0 0 0
0 1 10 10.995
1 0 0 0 0 0 0 0
0 1 10 10.995
1 0 0 0 0 0 0 0
0 1 9 10.495
1 0 0 0 0 0 0 0
0 1 7 8.995
1 0 0 0 0 0 0 0
0 1 NA NA
1 0 0 0 0 0 0 0
0 1 NA NA
0 1 0 0 0 1 0 0
0 0 9 23.995
0 1 0 0 0 1 0 0
0 0 NA NA
0 1 0 0 0 1 0 0
0 0 NA NA
0 1 0 0 0 1 0 0
0 0 NA NA
0 1 0 0 0 1 0 0
0 0 NA NA
0 1 0 0 0 1 0 0
0 0 NA NA
0 1 0 0 0 0 1 0
0 0 9 21.995
0 1 0 0 0 0 1 0
0 0 11 21.995
0 1 0 0 0 0 1 0
0 0 NA NA
0 1 0 0 0 0 1 0
0 0 NA NA
0 1 0 0 0 0 1 0
0 0 NA NA
0 1 0 0 0 0 1 0
0 0 NA NA
0 1 0 0 0 0 0 1
0 0 6 15.995
0 1 0 0 0 0 0 1
0 0 NA NA
0 1 0 0 0 0 0 1
0 0 NA NA
0 1 0 0 0 0 0 1
0 0 NA NA
0 1 0 0 0 0 0 1
0 0 NA NA
0 1 0 0 0 0 0 1
0 0 NA NA
0 1 0 0 0 0 0 0
1 0 11 17.995
0 1 0 0 0 0 0 0
1 0 9 16.495
0 1 0 0 0 0 0 0
1 0 8 15.995
0 1 0 0 0 0 0 0
1 0 9 13.995
0 1 0 0 0 0 0 0
1 0 NA NA
0 1 0 0 0 0 0 0
1 0 NA NA
0 1 0 0 0 0 0 0
0 1 NA NA
0 1 0 0 0 0 0 0
0 1 NA NA
0 1 0 0 0 0 0 0
0 1 NA NA
0 1 0 0 0 0 0 0
0 1 NA NA
0 1 0 0 0 0 0 0
0 1 NA NA
0 1 0 0 0 0 0 0
0 1 NA NA
0 0 1 0 0 1 0 0
0 0 9 18.495
0 0 1 0 0 1 0 0
0 0 NA NA
0 0 1 0 0 1 0 0
0 0 NA NA
0 0 1 0 0 1 0 0
0 0 NA NA
0 0 1 0 0 1 0 0
0 0 NA NA
0 0 1 0 0 1 0 0
0 0 NA NA
0 0 1 0 0 0 1 0
0 0 10 19.695
0 0 1 0 0 0 1 0
0 0 NA NA
0 0 1 0 0 0 1 0
0 0 NA NA
0 0 1 0 0 0 1 0
0 0 NA NA
0 0 1 0 0 0 1 0
0 0 NA NA
0 0 1 0 0 0 1 0
0 0 NA NA
0 0 1 0 0 0 0 1
0 0 6 16.595
0 0 1 0 0 0 0 1
0 0 6 15.995
0 0 1 0 0 0 0 1
0 0 NA NA
0 0 1 0 0 0 0 1
0 0 NA NA
0 0 1 0 0 0 0 1
0 0 NA NA
0 0 1 0 0 0 0 1
0 0 NA NA
0 0 1 0 0 0 0 0
1 0 11 16.995
0 0 1 0 0 0 0 0
1 0 9 15.595
0 0 1 0 0 0 0 0
1 0 7 11.995
0 0 1 0 0 0 0 0
1 0 8 13.195
0 0 1 0 0 0 0 0
1 0 NA NA
0 0 1 0 0 0 0 0
1 0 NA NA
0 0 1 0 0 0 0 0
0 1 9 10.295
0 0 1 0 0 0 0 0
0 1 NA NA
0 0 1 0 0 0 0 0
0 1 NA NA
0 0 1 0 0 0 0 0
0 1 NA NA
0 0 1 0 0 0 0 0
0 1 NA NA
0 0 1 0 0 0 0 0
0 1 NA NA
0 0 0 1 0 1 0 0
0 0 7 18.950
0 0 0 1 0 1 0 0
0 0 NA NA
0 0 0 1 0 1 0 0
0 0 NA NA
0 0 0 1 0 1 0 0
0 0 NA NA
0 0 0 1 0 1 0 0
0 0 NA NA
0 0 0 1 0 1 0 0
0 0 NA NA
0 0 0 1 0 0 1 0
0 0 9 21.000
0 0 0 1 0 0 1 0
0 0 7 17.950
0 0 0 1 0 0 1 0
0 0 NA NA
0 0 0 1 0 0 1 0
0 0 NA NA
0 0 0 1 0 0 1 0
0 0 NA NA
0 0 0 1 0 0 1 0
0 0 NA NA
0 0 0 1 0 0 0 1
0 0 NA NA
0 0 0 1 0 0 0 1
0 0 NA NA
0 0 0 1 0 0 0 1
0 0 NA NA
0 0 0 1 0 0 0 1
0 0 NA NA
0 0 0 1 0 0 0 1
0 0 NA NA
0 0 0 1 0 0 0 1
0 0 NA NA
0 0 0 1 0 0 0 0
1 0 8 15.450
0 0 0 1 0 0 0 0
1 0 10 14.750
0 0 0 1 0 0 0 0
1 0 8 15.750
0 0 0 1 0 0 0 0
1 0 7 12.250
0 0 0 1 0 0 0 0
1 0 NA NA
0 0 0 1 0 0 0 0
1 0 NA NA
0 0 0 1 0 0 0 0
0 1 5 6.450
0 0 0 1 0 0 0 0
0 1 8 10.950
0 0 0 1 0 0 0 0
0 1 8 9.750
0 0 0 1 0 0 0 0
0 1 4 6.950
0 0 0 1 0 0 0 0
0 1 9 9.950
0 0 0 1 0 0 0 0
0 1 9 8.750
0 0 0 0 1 1 0 0
0 0 10 20.950
0 0 0 0 1 1 0 0
0 0 9 21.500
0 0 0 0 1 1 0 0
0 0 NA NA
0 0 0 0 1 1 0 0
0 0 NA NA
0 0 0 0 1 1 0 0
0 0 NA NA
0 0 0 0 1 1 0 0
0 0 NA NA
0 0 0 0 1 0 1 0
0 0 8 16.850
0 0 0 0 1 0 1 0
0 0 8 16.750
0 0 0 0 1 0 1 0
0 0 8 19.850
0 0 0 0 1 0 1 0
0 0 NA NA
0 0 0 0 1 0 1 0
0 0 NA NA
0 0 0 0 1 0 1 0
0 0 NA NA
0 0 0 0 1 0 0 1
0 0 1 9.950
0 0 0 0 1 0 0 1
0 0 NA NA
0 0 0 0 1 0 0 1
0 0 NA NA
0 0 0 0 1 0 0 1
0 0 NA NA
0 0 0 0 1 0 0 1
0 0 NA NA
0 0 0 0 1 0 0 1
0 0 NA NA
0 0 0 0 1 0 0 0
1 0 5 8.995
0 0 0 0 1 0 0 0
1 0 8 17.500
0 0 0 0 1 0 0 0
1 0 NA NA
0 0 0 0 1 0 0 0
1 0 NA NA
0 0 0 0 1 0 0 0
1 0 NA NA
0 0 0 0 1 0 0 0
1 0 NA NA
0 0 0 0 1 0 0 0
0 1 10 10.750
0 0 0 0 1 0 0 0
0 1 4 4.950
0 0 0 0 1 0 0 0
0 1 10 9.250
0 0 0 0 1 0 0 0
0 1 8 9.250
0 0 0 0 1 0 0 0
0 1 9 8.995
0 0 0 0 1 0 0 0
0 1 NA NA
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|