problem in gen. initial in winbugs

#1
dear all members
i have a small problem in the code below. in the last step gen initial and after press on this key i cannot get on result because the program still generate these values and i cannot get on any resutls after then i must close the program because it is not work and restart it.
many thanks

thanoon





model{
#Group one
for(i in 1:N1){
#measurement equation model
for(j in 1:p){y1[i,j]~dnorm(mu1[i,j],psi1[j])I(low[z1[i,j]+1],high[z1[i,j]+1])}
mu1[i,1]<-xi1[i,1]
mu1[i,2]<-lam1[1]*xi1[i,1]
mu1[i,3]<-lam1[2]*xi1[i,1]
mu1[i,4]<-lam1[3]*xi1[i,1]
mu1[i,5]<-xi1[i,2]
mu1[i,6]<-lam1[4]*xi1[i,2]
mu1[i,7]<-xi1[i,3]
mu1[i,8]<-lam1[5]*xi1[i,3]
mu1[i,9]<-eta1
mu1[i,10]<-lam1[6]*eta1
#structural equation model
xi1[i,1:4]~dmnorm(u1[1:4],phi1[1:4,1:4])
eta1~dnorm(nu1,psd1)
nu1<-gam1[1]*xi1[i,1]+gam1[2]*xi1[i,2]+gam1[3]*xi1[i,3]+gam1[4]*xi1[i,1]*xi1[i,2]
}# end of i

#Group two
for(i in 1:N2){
#measurement equation model
for(j in 1:p){y2[i,j]~dnorm(mu2[i,j],psi2[j])I(low[z2[i,j]+1],high[z2[i,j]+1])}
mu2[i,1]<-xi2[i,1]
mu2[i,2]<-lam2[1]*xi2[i,1]
mu2[i,3]<-lam2[2]*xi2[i,1]
mu2[i,4]<-lam2[3]*xi2[i,1]
mu2[i,5]<-xi2[i,2]
mu2[i,6]<-lam2[4]*xi2[i,2]
mu2[i,7]<-xi2[i,3]
mu2[i,8]<-lam2[5]*xi2[i,3]
mu2[i,9]<-eta2
mu2[i,10]<-lam2[6]*eta2

#structural equation model
xi2[i,1:4]~dmnorm(u2[1:4],phi2[1:4,1:4])
eta2~dnorm(nu2,psd2)
nu2<-gam2[1]*xi2[i,1]+gam2[2]*xi2[i,2]+gam2[3]*xi2[i,3]+gam2[4]*xi2[i,1]*xi2[i,2]
}# end of i
for(i in 1:4){zero4<-0.0}

#priors on loadings and coefficients
for(i in 1:10){mu.y1~dnorm(0.0,4.0)}
for(i in 1:6){lam1~dnorm(0.8,4.0)}
for(i in 1:4){gam1~dnorm(0.6,4.0)}
for(i in 1:10){mu.y2~dnorm(0.0,4.0)}
for(i in 1:6){lam2~dnorm(0.8,4.0)}
for(i in 1:4){gam2~dnorm(0.6,4.0)}
#priors on precisions
for(j in 1:p){
psi1[j]~dgamma(10,8)
sgm1[j]<-1/psi1[j]
}
psd1~dgamma(10,8)
sgd1<-1/psd1
phi1[1:4,1:4]~dwish(R[1:4,1:4], 30)
phx1[1:4,1:4]<-inverse(phi1[1:4,1:4])
for(j in 1:p){
psi2[j]~dgamma(10,8)
sgm2[j]<-1/psi2[j]
}
psd2~dgamma(10,8)
sgd2<-1/psd2
phi2[1:4,1:4]~dwish(R[1:4,1:4], 30)
phx2[1:4,1:4]<-inverse(phi2[1:4,1:4])
} #end of model