小编use*_*124的帖子

Xn可以取值-1或1,概率为0.5.并且Sn = Sn-1 + Xn如何计算在Sn = X1 + X2 + ::: + Xn给出的时刻n观察到的部分和.我想在这里模拟随机游走.我做了以下但我不确定它是对的:

rw <- function(n){
    x=numeric(n)
    xdir=c(TRUE, FALSE)
    step=c(1,-1)
    for (i in 2:n)
    if (sample(xdir,1)) {
        x[i]=x[i-1]+sample(step,1)
    } else {
        x[i]=x[i-1]
    }
    list(x=x)
}

请帮忙!