走者ベクトルを付け加えるやり方で進めることにした。
ここで問題となるのが、超低い確率だが3アウトまで行きつかねぇってときにうわ〜。
ということで制限設けて強制終了することで、終了しない事態、を防ぎましょう。
p<- seq(from=0,to=1,by=0.01) #アウトかヒットしか打撃パターンがないとして、pはアウトになる確率 k<- 3:15 #kは1イニングで打席に立つ打者数。 #今、15人までを上限として足し合わせると、16人以上打席に立つのはこれの余事象 data<- rep(0,length(p)) for(i in 1:length(p)){ a<- (k-1)*(k-2)/2*((1-p[i])^(k-3))*p[i]^3 #がちゃがちゃ計算。二項分布。 cumsum(a)[length(a)] -> data[i] #末尾が総足し算結果 } plot(p,data) #どや?
どや?とかいっても先生にスクリプト書いてもらったのだが。
どうがんばってもpは0.6〜0.8くらいだから、これでいこうか…
大天使「そんな評価で大丈夫か?」
イーなんとかさん「大丈夫だ、問題ない」
神は言っている。16人以上打席に立ったら排除してカウントしておけと。