午前中の日誌ではMicrosoft Excelで理論計算しましたが、今回は1万回くじを引くケースを1万通りシミュレーションしました。
結果、1〜4等の獲得メダル合計は次の通りとなりました。
最小:1137
最大:1745
平均:1418.9
標準偏差:77.6
尚、特等は1万通りでは誤差が大きくなる恐れがあるので計算に入れていません。
このデータを元にnorm.inv関数とnorm.dist関数から逆算すると次の通りとなります。
99%確定値:1238.5
99.99%確定値:1130.5
99.9999%確定値:1050.3
1000個確定確率:99.9999967%
1100個確定確率:99.9980%
1200個確定確率:99.76%
ただし、1万通りのデータなので、確率予測値に0.01%台の誤差はあると考えられます。
この結果からは、誤差を考慮しても、1万回(約12万円)で99.9%以上確定、8千回強で99%確定と言えそうです。
午前中の日誌の計算結果と大差ないことから、計算に大きな間違いはなかったと推定できます。