【求助】这段代码运行不出来

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

这段二分法的代码只要迭代9步，但是运行之后一直busy，不知道哪里出了问题

blci<-function(r,n,alpha=.05,eps=1e-3){
  if(r<=0) pl<-0 else{
    plu<-r/n
    pll<-0
    pl<-(pll+plu)/2
    u<-1-pbinom(r-1,n,pl)
    i<-1
    while((plu-pll)/plu>eps){
      i<-i+1
      if(u>alpha){
        plu<-pl
        pl<-(pll+plu)/2
        u<-1-pbinom(r-1,n,pl)
      }
        else{
        pll<pl
        pl<-(pll+plu)/2
        u<-1-pbinom(r-1,n,pl)
      }
    }
  }
    c(lcl=pl,neval=i)
}

blci(5,15,alpha=.05,eps=1e-2)

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：BUS 二分法 不知道 二分法

就是普通的二分法代码，但就是求解不出来

正在看。

in the else{} part of the while loop,

pll<pl  

should be

pll<- pl