上次论坛上有个大神帮忙根据第一个图中的内容编写了程序,我类比将第二张图的公式加进去,进行新的蒙特卡洛模拟,可是50次都出不来结果,求教这个程序是不是进入死循环了?
LB <- 0
UB <- 10
k=0.05
a=1
ARL0 <- 200
repeat {
RL <- NULL
for(i in 1:100) {
e = 0
count <- 0
h <- (LB + UB) / 2
while(e < h ) {
x <- rnorm(1)
e <- max(k*x+(1-k)*e, 0)
count <- count + 1
}
RL
<- count
}
ARL <- mean(RL)
if(abs(ARL0 - ARL) <= 1) break
if(ARL < ARL0 - 1) {
LB <- h
} else if(ARL > ARL0 + 1) {
UB <- h
}
}