- 【R言語】統計解析フリーソフトR 第6章【GNU R】 [無断転載禁止]©2ch.net
559 :132人目の素数さん[sage]:2018/11/09(金) 02:01:36.95 ID:Qla0VxTD - >>558
ご指摘ありがとうございました。 修正しました。 rle1 <- function(N=100,n=5,p=0.5){ r=rle(rbinom(N,1,p)) max(r$len[which(r$val==1)]) } 結果 > system.time(mean(replicate(1e4,max(rle1()>=5)))) user system elapsed 4.430 0.000 4.546 > system.time(mean(replicate(1e4,seqn()))) user system elapsed 4.490 0.010 4.609 > system.time(mean(replicate(1e4,mhs(rbinom(100,1,0.5))>=5))) user system elapsed 7.440 0.000 7.656 > system.time(searchseq()) user system elapsed 1.950 0.000 2.066 >
| - 【R言語】統計解析フリーソフトR 第6章【GNU R】 [無断転載禁止]©2ch.net
560 :132人目の素数さん[sage]:2018/11/09(金) 07:25:00.97 ID:Qla0VxTD - 無理矢理1行にして実行
system.time(mean(replicate(1e4,any(diff(cumsum(rbinom(100,1,0.5)),5)==5)))) user system elapsed 1.820 0.000 1.886 > > system.time(mean(replicate(1e4,with(rle(rbinom(100,1,0.5)), max(lengths[wh <e(1e4,with(rle(rbinom(100,1,0.5)), max(lengths[whi ch(values==1)])>=5)))) user system elapsed 4.370 0.010 4.478
|
|