model{ for(i in 1:n){ FEV[i] ~ dnorm(mu[i],tau) mu[i] <- beta[1] + beta[2]*Age[i] + beta[3]*Age[i]*Age[i] + beta[4]*Smoke[i] + beta[5]*Age[i]*Smoke[i] + beta[6]*Age[i]*Age[i]*Smoke[i] CPOinv[i] <- sqrt(2*3.14159/tau)*exp(0.5*tau*pow(FEV[i]-mu[i],2)) } beta[1] ~ dnorm(0,0.001) beta[2] ~ dnorm(0,0.001) beta[3] ~ dnorm(0,0.001) beta[4] ~ dnorm(0,0.001) beta[5] ~ dnorm(0,0.001) beta[6] ~ dnorm(0,0.001) tau ~ dgamma(0.001,0.001) }