	model
	{
		for (i in 1:N) {
			for (j in 1:M) {
	# Survival times bounded below by censoring times:
				t[i,j] ~ dweib(r,mu[i,j]) I(t.cen[i,j],);
				log(mu[i,j]) <- alpha + beta.age*age[i,j] 
						+ beta.sex *sex[i]  
						+ beta.dis[disease[i]] + b[i];
			}
	# Random effects:
			b[i] ~ dnorm(0.0, tau)   
		}
	# Priors:
		alpha ~ dnorm(0.0, 0.0001);
		beta.age ~ dnorm(0.0, 0.0001);
		beta.sex ~ dnorm(0.0, 0.0001);
	#	beta.dis[1] <- 0;  # corner-point constraint
		for(k in 2 : 4) {
			beta.dis[k] ~ dnorm(0.0, 0.0001);
			}
		tau ~ dgamma(1.0E-3, 1.0E-3);
		r ~ dgamma(1.0, 1.0E-3); 
		sigma <- 1/sqrt(tau); # s.d. of random effects
	}
