在我之前的问题中: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
提前致谢!!!
答案 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
。