作业问题:
考虑一维伊辛模型。
让。为 -1 或 +1
设计一个 gibbs 采样算法以近似从目标分布生成样本。
我的尝试:
随机选择值(-1 或 1)来填充向量。所以也许。所以这是。
所以现在我们需要继续进行第一次迭代。我们必须分别为绘制 40 个不同的 x 。所以...
从绘制
从绘制
从绘制
ETC..
所以让我绊倒的部分是我们如何从条件分布中实际提取。如何发挥作用?也许一次平局的例子可以解决问题。
作业问题:
考虑一维伊辛模型。
让。为 -1 或 +1
设计一个 gibbs 采样算法以近似从目标分布生成样本。
我的尝试:
随机选择值(-1 或 1)来填充向量。所以也许。所以这是。
所以现在我们需要继续进行第一次迭代。我们必须分别为绘制 40 个不同的 x 。所以...
从绘制
从绘制
从绘制
ETC..
所以让我绊倒的部分是我们如何从条件分布中实际提取。如何发挥作用?也许一次平局的例子可以解决问题。
先看这个案例。删除不依赖于的术语,我们有。
x_1 <- sample(c(-1, 1), 1, prob = c(exp(-x_2), exp(x_2)) / (2*cosh(x_2)))
将其推广到(注意差异;请参阅下面 Ilmari 的评论)。
您可以使用 Ising 的分析结果来检查您的模拟吗?