如果计算出一个分布在不同损失函数下的贝叶斯估计，怎么利用R软件编程MCMC模拟，计算估计误差(bias)，和其他的估计方法想比较？