10/27 プログラミングセミナー

走者ベクトルを付け加えるやり方で進めることにした。
ここで問題となるのが、超低い確率だが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人以上打席に立ったら排除してカウントしておけと。