高校数学の質問スレPart407at MATH高校数学の質問スレPart407 - 暇つぶし2ch■コピペモード□スレを通常表示□オプションモード□このスレッドのURL■項目テキスト400:132人目の素数さん 20/09/02 04:55:24.53 KVQ8oPlE.net >>382-383 頭悪過ぎw 401:132人目の素数さん 20/09/02 05:49:23.58 62Rt+r6G.net >>378 今、手元にPCないのでスマホで100回シミュレーションして平均値を出してみた。 sim <- function(n=900,p=0.02,m=100){ i = 0 # trial number j = 0 # failure k = 0# success while(i<=n){ a=ifelse(j==(m-1),1,rbinom(1,1,p)) # ceiling or not j=ifelse(a==1,0,j+1) # success or failure k=k+a # total success i=i+1 } return(k) } mean(replicate(1e2,sim())) 実行結果 > mean(replicate(1e2,sim())) [1] 21.26 > 402:132人目の素数さん 20/09/02 06:10:57.22 62Rt+r6G.net >>386 求めるのは16回以上当たる確率だった 1万回のシミュレーション sim <- function(n=900,p=0.02,m=100){ i = 0 # 試行回数 j = 0 # 現在の外れの連続数 k = 0 # 当たった総数 while(i<=n){ a=ifelse(j==(m-1),1,rbinom(1,1,p)) # 天井なら1でそれ以外はpで当たる j=ifelse(a==1,0,j+1) # 当たりか否かでjをリセット k=k+a # 現在の当たりの総数 i=i+1 # 試行回数 } return(k) # 当たりの数を返す } K=replicate(1e4,sim()) mean(K>=16) # 16回以上の割合 > mean(K>=16) # 16回以上の割合 [1] 0.9376 次ページ最新レス表示レスジャンプ類似スレ一覧スレッドの検索話題のニュースおまかせリストオプションしおりを挟むスレッドに書込スレッドの一覧暇つぶし2ch