求大神看看我这个R语言程序

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

根据这个图片中的过程，写出编程过程，我是R语言初学者，好多都不怎么会，就自己编了一个过程，想求教大神们，过程有什么要改进的么？步骤5我也不会循环。如果我编的过程太惨目忍睹，能不能帮忙编一编

[size=0.83em]半小时前 上传 下载附件 [size=0.83em](71.5 KB) 这是我的过程 i<-1:100;j<-1:10000; n<-length(i);m<-length(j); z<-runif(n,0,1) LB<-0 UB<-10 h<-(LB:UB) k<-0.25 ARL0<-200 z[1]=0 count=0 h=(LB+UB)/2; for(z in 1:n) {    for(RL in 1:m)    {      if(z>=h)      {        RL[j]=count;      }      if(z<h)      {        x<-rnorm(1,0,1)        z<-max(z+x-k,0)        count=count+1;z=z+1;        while(z>h)        {          RL[j]=count;        }      }      ARL=0      ARL=ARL+RL[j]    } } ARL if(abs(ARL-200)<=1) {   h+=h }

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：R语言 R语言初学者 size 初学者

人工置顶

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
感觉跑的好慢啊

jinkelazzz 发表于 2017-12-28 17:55
LB

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中，没有跑出来结果

一朵小墨墨 发表于 2017-12-28 21:07
ARL

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

jinkelazzz 发表于 2017-12-29 09:57
就是repeat啊，蒙特卡洛速度太慢，你可以少模拟几次看看

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