Welcome to OStack Knowledge Sharing Community for programmer and developer-Open, Learning and Share
Welcome To Ask or Share your Answers For Others

Categories

0 votes
245 views
in Technique[技术] by (71.8m points)

stat - Plotting a random walk on R

I've run a Gibbs sampler and obtained a sample for $X_1$ and $X_2$. I'm trying to recreate a plot like this one:

enter image description here

How do I recreate the walk part on R?

question from:https://stackoverflow.com/questions/65887658/plotting-a-random-walk-on-r

与恶龙缠斗过久,自身亦成为恶龙;凝视深渊过久,深渊将回以凝视…
Welcome To Ask or Share your Answers For Others

1 Answer

0 votes
by (71.8m points)

The most straightforward way I can think of to recreate the walk part would be like this:

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.: I made up a simple walk

The result was this one:

enter image description here


与恶龙缠斗过久,自身亦成为恶龙;凝视深渊过久,深渊将回以凝视…
Welcome to OStack Knowledge Sharing Community for programmer and developer-Open, Learning and Share
Click Here to Ask a Question

...