除了对轨迹的标准检查（相关性，是否混合良好等），对于“单元测试”每个变量的采样，例如，在我的 Gibbs 采样器中，我执行以下操作：

• 创建一个小型玩具数据集。

• 在我的 Gibbs 采样器中，注释掉（我使用if True和if False声明）除一个变量之外的所有变量，其中固定变量（这次未采样的变量）设置为它们的真实值。

• 查看一个采样变量是否收敛到其真实值。

对所有变量一一进行。然后对两个变量进行采样（尝试不同的对），依此类推，直到我全部采样。然后我可以使用更复杂的玩具数据集，或者只是尝试真实数据。

（不过，我很想看看其他人是如何做到的，以及在编程方面是否有一些“好的做法”）