	model
	{
		beta0 ~ dnorm(0.0, 0.001)
		beta1 ~ dnorm(0.0, 0.001)
		for (i in 1:N) {
			Y[i, 1:M] ~ dmnorm(mu[], Omega[ , ]) 
		}                                  
		for(j in 1:M) { 
			mu[j] <- beta0 + beta1* age[j]
		}
		Omega[1 : M , 1 : M]  ~ dwish(R[ , ], 4)
		Sigma[1 : M , 1 : M] <- inverse(Omega[ , ])

	}
