全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
870 5
2018-01-02
上次论坛上有个大神帮忙根据第一个图中的内容编写了程序,我类比将第二张图的公式加进去,进行新的蒙特卡洛模拟,可是50次都出不来结果,求教这个程序是不是进入死循环了? 1 2


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
  }
}

二维码

扫码加我 拉你入群

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

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

全部回复
2018-1-2 16:31:07
人工置顶~
二维码

扫码加我 拉你入群

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

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

2018-1-2 17:14:13
有人么
二维码

扫码加我 拉你入群

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

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

2018-1-2 19:38:01
二维码

扫码加我 拉你入群

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

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

2018-1-3 13:30:34
##RL <- count
RL[i] <- count
二维码

扫码加我 拉你入群

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

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

2018-1-3 15:36:11
jinkelazzz 发表于 2018-1-3 13:30
##RL
您好,我改了后还是出不来结果,是不是计算量太大了,还是程序有哪不对的呢
二维码

扫码加我 拉你入群

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

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

相关推荐
栏目导航
热门文章
推荐文章

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群