在WinBUGS中，有model如下：
model
{for(i in 1:N)
  {
    y[i] ~ dnorm(mean.y[i], prec.y)
    mean.y[i] <- beta3 + beta2*m[i] + beta1*x[i]
  }
  beta2 ~ dnorm(0, 1.0E-6)
  beta3 ~ dnorm(0, 1.0E-6)
  beta1 ~ dnorm(0, 1.0E-6)   
  prec.y ~ dgamma(0.009, 0.001)
}
如果我想求多元线性回归方程Y=beta3 + beta2*m + beta1*x的回归系数beta2和beta1的协方差，应该如何写model？谢谢！请各位专家指教！