人工置顶

LB <- 0
UB <- 10
k <- 0.25
ARL0 <- 200 #step1
repeat {
    RL <- NULL
    for(i in 1:10000) {
        z = 0
        count <- 0 #step2
        h <- (LB + UB) / 2
        while(z < h) {
            x <- rnorm(1)
            z <- max(z + x - k, 0)
            count <- count + 1
        } #this loop is step3
        RL[i] <- count
    } #this loop is step4
    ARL <- mean(RL)
    if(abs(ARL0 - ARL) <= 1) break
    if(ARL < ARL0 - 1) {
        LB <- h
    } else if(ARL > ARL0 + 1) {
        UB <- h
    }
}#this loop is step5
> h
[1] 5.576172
> ARL
[1] 199.9728
感觉跑的好慢啊

ARL <- mean(RL)
    if(abs(ARL0 - ARL) <= 1) break
    if(ARL < ARL0 - 1) {
        LB <- h
    } else if(ARL > ARL0 + 1) {
        UB <- h
    }
}#this loop is step5
大神，我想问一下，您写的这个step5中我没看明白，题中更新LB=h后，不是还有个重复步骤2-5，直到h+=h的这个，这个是怎么循环的呢
我复制到我得Rstudio中，没有跑出来结果

就是repeat啊，蒙特卡洛速度太慢，你可以少模拟几次看看

太感谢了，我之前模拟的次数太多，电脑跑不出来，谢谢~现在出结果了