求教怎么改进这个程序

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

上次论坛上有个大神帮忙根据第一个图中的内容编写了程序，我类比将第二张图的公式加进去，进行新的蒙特卡洛模拟，可是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
  }
}

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：蒙特卡洛模拟 蒙特卡洛 蒙特卡

人工置顶~

有人么

##RL <- count
RL[i] <- count

jinkelazzz 发表于 2018-1-3 13:30
##RL

您好，我改了后还是出不来结果，是不是计算量太大了，还是程序有哪不对的呢