在R上绘制随机游动图

我运行了一个Gibbs取样器，并以$X_1$和$X_2$获得了一个样本。我试图重现这样一个情节：

---

如何在R上重建漫游部分？

我能想到的重建漫游部分最简单的方法如下：

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()

旁白：我做了一次简单的散步

结果是：

---

欢迎来到SO，汤姆！如果这样的问题（特别是在R中）是可复制和自包含的，那么它们会做得更好。我的意思是包括尝试的代码（请明确说明非基本包）和样本代表性数据（可能通过

dput（head（x））

或以编程方式构建数据（例如，

data.frame（…）

），可能在

set.seed（1）

之后随机）。参考文献：，和。