如上图片所示。不知什么原因，用这些参数的满条件分布得到的模拟结果只有mu是正确的，sigma0和sigma不对，我想是否在sigma0和sigma中推导的满条件分布不对，但我又推导了一次还是这样，求大神解答。我附上我的代码：

我用y1模拟数据作为实际数据
除了mu正确，其他两个都不对。