- 【R言語】統計解析フリーソフトＲ 第6章【GNU R】 [無断転載禁止]©2ch.net
742 ：１３２人目の素数さん[sage]：2019/12/14(土) 03:40:27.74 ID:L1XaepqW - 分からない問題スレから、
>>
1回3.6%で激レアが出るガチャを10回回した確率って
36%なのでしょうか？
それとも0.964*0.964*0.964(略 0.964を10回電卓にかけた数なのでしょうか？
教えてください。
<<
百万回シミュレーション
p=3.6/100
N=10
sim <- function() any(rbinom(N,1,p)==1)
mean(replicate(1e6,sim()))
> mean(replicate(1e6,sim()))
[1] 0.306628
| - 【R言語】統計解析フリーソフトＲ 第6章【GNU R】 [無断転載禁止]©2ch.net
743 ：１３２人目の素数さん[sage]：2019/12/14(土) 03:51:08.05 ID:L1XaepqW - シミュレーション その２
gacha=c(rep(1,36),rep(0,1000-36))
sim2 <-function() any(sample(gacha,10,replace=TRUE)==1)
mean(replicate(1e6,sim2()))
> mean(replicate(1e6,sim2()))
[1] 0.306904
何故か0.36にならない、どうしてだろ？
| - 【R言語】統計解析フリーソフトＲ 第6章【GNU R】 [無断転載禁止]©2ch.net
744 ：１３２人目の素数さん[sage]：2019/12/14(土) 04:10:22.87 ID:L1XaepqW - シミュレーションその３ （処理速度の関係で１０万回の平均）
gacha=c(rep(1,36),rep(0,1000-36))
sim3 <- function() any(replicate(10,sample(gacha,1))==1)
mean(replicate(1e5,sim3()))
> mean(replicate(1e5,sim3()))
[1] 0.30691
これも0.307弱だな。
何が悪いんだろ？ 俺の頭かな？？
| - 【R言語】統計解析フリーソフトＲ 第6章【GNU R】 [無断転載禁止]©2ch.net
746 ：１３２人目の素数さん[sage]：2019/12/14(土) 08:18:13.50 ID:L1XaepqW - >>745
ありがとう、
やっぱ、
これだったｗ
何が悪いんだろ？ 俺の頭かな？？
| - 【R言語】統計解析フリーソフトＲ 第6章【GNU R】 [無断転載禁止]©2ch.net
747 ：１３２人目の素数さん[sage]：2019/12/14(土) 08:19:41.64 ID:L1XaepqW - >>746
シミュレーションにバグはなかったｗ
|
|