如何在随机游走中包含计步器?

时间:2014-10-29 22:15:11

标签: r loops random-walk

在我的随机游走循环中,我需要合并一个使用步骤编号作为变量的函数。

如何在随机游走过程中让R产生一个步数(循环计数器可能是一个更好的术语)?

我需要在此代码的末尾,我目前有虚拟代码step.num

作为参考我的随机游走代码是:

    walkW <- function(n.times=125,
               xlim=c(524058,542800),
               ylim=c(2799758,2818500),
               start=c(542800,2815550),
               stepsize=c(4000,4000)) {

    plot(c(0,0),type="n",xlim=xlim,ylim=ylim,
           xlab="Easting",ylab="Northing",col="white",col.lab="white")#arena 
    x <- start[1]
    y <- start[2]
    steps <- 1/c(1,2,4,8,12,16)
    steps.y <- c(steps,-steps,0)
    steps.x <- c(steps[c(1,5,6)],-steps,0)
    points(x,y,pch=16,col="green",cex=1)

for (i in 1:n.times) {
        repeat {
           xi <- stepsize[1]*sample(steps.x,1)
           yi <- stepsize[2]*sample(steps.y,1)
           newx <- x+xi
           newy <- y+yi
           if (newx>xlim[1] && newx<xlim[2] &&
               newy>ylim[1] && newy<ylim[2]) break
        }
        lines(c(x,newx),c(y,newy),col="white")
         x <- newx
           y <- newy
}
  step.num<-(your magical step counter function) #I need the step counter here
 }

set.seed(101)
walkW(s)

1 个答案:

答案 0 :(得分:3)

在&#34;破坏&#34;的for循环之后,索引变量仍然存在,所以你的神奇代码就是:

 step.num <- i

虽然for循环在技术上是功能,但有时它们似乎有不同的范围规则,这就是这样的时间。您的函数会返回&#39; step.num&#39;价值,退出walkW后,您将无法访问“我”,但如果for循环在.GlobalEnv中完成,那么索引&#34; ; I&#34;当for循环通过break结束时,仍然不会被重置或空转。