トップページ > 数学 > 2018年11月09日 > Qla0VxTD

書き込み順位&時間帯一覧

25 位/124 ID中時間01234567891011121314151617181920212223Total
書き込み数0010000100000000000000002


使用した名前一覧書き込んだスレッド一覧
132人目の素数さん
【R言語】統計解析フリーソフトR 第6章【GNU R】 [無断転載禁止]©2ch.net

書き込みレス一覧

【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


※このページは、『2ちゃんねる』の書き込みを基に自動生成したものです。オリジナルはリンク先の2ちゃんねるの書き込みです。
※このサイトでオリジナルの書き込みについては対応できません。
※何か問題のある場合はメールをしてください。対応します。