在随机游走中加入停止功能

时间:2014-10-16 22:36:24

标签: r random random-walk

在我之前的问题中:How do I put arena limits on a random walk?社区帮助在集合竞技场中创建随机漫步功能。此功能用于模拟在区域内移动的鱼,但现在我需要决定何时在满足某个条件时停止。

我认为这会很简单 {{if(z>P)break}}放在循环函数之前。我想要理解的是“如果满足这个条件就停止,否则一直走到最大步数。

相反,它导致我的随机游走变得确定(我总是得到相同的路径,它在step.max之前永远不会停止)。

主要问题:如果z> P?

,如何判断随机游走是否停止?

供参考:

step.max<-125
step.prob<-function(n.times=step.max){
draw=sample(0:100,1,replace=T)
CS<-sample(draw,size=1,replace=TRUE)
CS.max<-100
step.num<-15
SP<-((CS/CS.max)*(1-(step.num/step.max))+(step.num/step.max))*100
if(SP>P){stop('Settled at step number',P)}else{SP
     }
}
z<-step.prob(1) #renaming the above function to be easier to reference later
P<-80 #preset cutoff point for value z, ranges from 0-100
walkE <- function(n.times=125,
               xlim=c(524058,542800),
               ylim=c(2799758,2818500),
               start=c(525000,2810000),
               stepsize=c(4000,4000)) {
    plot(c(0,0),type="n",xlim=xlim,ylim=ylim,
           xlab="Easting",ylab="Northing") 
    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,-steps[c(1,5,6)],0)
    points(x,y,pch=16,col="red",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="blue")
        x <- newx
        y <- newy
if(z>P){stop(points(newx,newy,col="green",cex=1))} 
                        #this is where I want it to stop if z>P
         else
if(z<P){points(newx,newy,pch=1,col="blue",cex=1)}
    else
if(step.max){points(newx,newy,pch=16,col="green",cex=1)}
set.seed(101)}
}
walkE(step.max) #run above random walk function walkE looped for the step.max number

提前致谢!!!

1 个答案:

答案 0 :(得分:0)

这非常简单,可以通过在用户定义的stop(...)函数中插入step.prob函数来实现。

step.prob<-function(n.times=step.max, p){
  draw=sample(0:100,1,replace=T)
  CS<-sample(draw,size=1,replace=TRUE)
  CS.max<-100
  CS.max
  step.num<-15
  SP<-((CS/CS.max)*(1-(step.num/step.max))+(step.num/step.max))*100
  if(SP > p) {
    stop('Your random walk exceeded ', p)
  } else {
  SP
  }
}

如果不这样做,请查看break命令。

因此,当随机游走值是>号码:

step.prob(p=300000)
# Error in step.prob(p = 3) : Your random walk exceeded 3

如果你想将函数返回的值设置为p,你可以在SP <- p命令之前添加stop