不balance的multivarite regression怎么做MCMC?

大家好，我有一个关于手写R codes做MCMC的问题。

假设说我有两个不同的y（y1,y2)
假设我有这两个y共同的x1,x2,x3

目标是做  y=Xb+e,即y1=b0+b1x1+b2x2+b3x3,y2=b0+b1x1+b2x2+b3x3，但y1和y2同时还有covariance（假设是multivariate normal）


我知道用R（不用任何bugs或者jags）如何写codes，但只限于y1和y2用的是同一组的X，如果y1和y2的X不同（假设有n人是回归y1的，m人是回归y2的）怎么写codes的呢？还没有怎么接触这种multivariate的东西， 但是有个程序现在需要用到。

请各位大神们指点迷津！谢谢了！！！！！！