如何用R函数进行mcmc抽样估计联合分布参数值

求各位大神相助：      已知几个参数的联合先验分布和似然函数，推导得到这些参数的联合后验分布形式（但是只是成比例的，并不是真实的分布形式），然后我希望求得其中一个参数在某个区间的概率？要用MCMC抽样吗？具体要怎么实现？各位亲们帮帮忙！已知xi和yi服从二元正态分布：


概率密度函数：

似然函数为：

其中

参数的联合先验分布为：

参数的联合后验分布为：
联合分布是积分形式的，就我如果参数的联合分布得到delta的分布，或者是在某个区间的概率？
R中有现成的用于贝叶斯估计的函数吗？
请各位大神帮帮忙，在下感激不尽！

图片显示好像有问题

求各位大神相助：      已知几个参数的联合先验分布和似然函数，推导得到这些参数的联合后验分布形式（但是只是成比例的，并不是真实的分布形式），然后我希望求得其中一个参数在某个区间的概率？要用MCMC抽样吗？具体要怎么实现？各位亲们帮帮忙！已知xi和yi服从二元正态分布：


概率密度函数：

似然函数为：

其中

参数的联合先验分布为：

参数的联合后验分布为：
联合分布是积分形式的，就我如果参数的联合分布得到delta的分布，或者是在某个区间的概率？

R中有现成的用于贝叶斯估计的函数吗？