在 R 上绘制随机游走

Question

我已经运行了一个 Gibbs 采样器并获得了 $X_1$ 和 $X_2$ 的样本。我正在尝试重新创建一个这样的情节：

如何在 R 上重新创建步行部分？

Answer 1

我能想到的最直接的重建步行部分的方法是这样的：

x = c(0)

for(i in 2:100){
    x[i] = x[i-1] + floor(runif(1, -5, 5))
}

pdf(file = "My Plot.pdf", width = 4, height = 4)

plot(x, type="l")

dev.off()

Obs.: 我做了一个简单的步行

结果是这个：

在 R 上绘制随机游走

Plotting a random walk on R

r

stat