チョコボール統計

チョコボールの秘密を統計解析で明らかにしていく。おもちゃのカンヅメ欲しい。

第ex12回 チョコボール計測

今日も番外編です。12回目。
同僚から「パチパチ」味をもらったので開封します。

計測結果

では計測結果です。

date best_before weight box_weight number factory shop angel net_weight mean_weight
2018-07-25 2019-02-01 32.421 4.902 16 小山工場 ドラックストア(横浜市 なし 27.519 1.720
2018-07-25 2019-02-01 32.393 4.862 16 小山工場 ドラックストア(横浜市 なし 27.531 1.721
2018-07-25 2019-02-01 31.578 4.841 16 小山工場 ドラックストア(横浜市 なし 26.737 1.671
2018-07-25 2019-02-01 32.139 4.886 16 小山工場 ドラックストア(横浜市 なし 27.253 1.703

今回もエンゼルさん現れませんでした。

基礎集計

「パチパチ」味の集計です

項目
計測データ数 9
銀のエンゼル出現数 0
金のエンゼル出現数 0
最小 中央値 最大値 平均
正味重量 26.737 27.519 28.819 27.579
個数 16.000 16.000 17.000 16.111

重量分布

f:id:hippy-hikky:20180726002803p:plain

この図は正味の重量のヒストグラムです。 赤い縦線が仕様を表しています。 青い太線で正規分布と仮定した最尤推定量をプロットしています。

エンゼルの出現確率推定

通常版のエンゼルの予測を行っていきます。 これまでの通常版パッケージの開封結果は次の通りです

項目
計測データ数 165
銀のエンゼル出現数 8
金のエンゼル出現数 0

この結果を使ってベイズ推定によるエンゼルの出現確率推定を行います。

はじめに銀のエンゼルの出現確率の推定です。
90%信用区間(上下それぞれ5%)は、下側が2.93%上側が8.63%という結果です。
f:id:hippy-hikky:20180726002951p:plain

現在の期待値は5.36%です。

次に金のエンゼルの出現確率の推定です。 90%信用区間(上下それぞれ5%)は、下側が0.03%、上側が1.84%という予測になっています。
f:id:hippy-hikky:20180726003022p:plain

広告

Amazonの欲しいものリスト作ってみました。 チョコボールのカンパ募集中です。
チョコボールをカンパする

第ex11回 チョコボール計測

お久しぶりです。 最近の酷暑によりチョコボールを部屋に置いておくと溶けちゃうんですよね。 なので、最近あんまり買ってなくて更新が滞ってます。。。(言い訳)

今日は番外編の11回目です。
「パチパチ」味という新作が出ていたので試してみました。

この新作が出たタイミングからなのかわかりませんが、 エンゼル無しのクチバシを集めて抽選で当たるキャンペーンが始まりました。

体積比50倍のチョコボールが当たるらしいです。
当然応募しますよ!

ちなみに、何が“パチパチ”かというと、駄菓子の「わたパチ」の中にあるような弾けるキャンディ(?)がチョコの中に入っています。 食べるとパチパチします。
妻曰く、「一回食べたらもういらないかな〜」とのこと。

計測結果

では計測結果です。

date best_before weight box_weight number factory shop angel net_weight mean_weight
2018-07-19 2019-02-01 32.366 4.852 16 小山工場 コンビニ(さいたま市 なし 27.514 1.720
2018-07-19 2019-02-01 32.325 4.845 16 小山工場 コンビニ(さいたま市 なし 27.480 1.718
2018-07-19 2019-02-01 32.592 4.858 16 小山工場 コンビニ(さいたま市 なし 27.734 1.733
2018-07-19 2019-02-01 32.478 4.852 16 小山工場 コンビニ(さいたま市 なし 27.626 1.727
2018-07-19 2019-02-01 33.691 4.872 17 小山工場 コンビニ(さいたま市 なし 28.819 1.695

ということで、エンゼルは現れませんでした。

基礎集計

「パチパチ」味の集計です

項目
計測データ数 5
銀のエンゼル出現数 0
金のエンゼル出現数 0
最小 中央値 最大値 平均
正味重量 27.480 27.626 28.819 27.835
個数 16.000 16.000 17.000 16.200

重量分布

f:id:hippy-hikky:20180722183630p:plain

この図は正味の重量のヒストグラムです。 赤い縦線が仕様を表しています。 青い太線で正規分布と仮定した最尤推定量をプロットしています。

パチパチ味もその他の味と同様に、 仕様上の重量よりも1gほど重めにコントロールされて出荷されているようですね。

エンゼルの出現確率推定

通常版のエンゼルの予測を行っていきます。 これまでの通常版パッケージの開封結果は次の通りです

項目
計測データ数 161
銀のエンゼル出現数 8
金のエンゼル出現数 0

この結果を使ってベイズ推定によるエンゼルの出現確率推定を行います。

はじめに銀のエンゼルの出現確率の推定です。
90%信用区間(上下それぞれ5%)は、下側が2.76%上側が8.70%という結果です。
f:id:hippy-hikky:20180722184054p:plain

現在の期待値は5.44%です。

次に金のエンゼルの出現確率の推定です。 90%信用区間(上下それぞれ5%)は、下側が0.02%、上側が1.86%という予測になっています。
f:id:hippy-hikky:20180722184129p:plain

広告

Amazonの欲しいものリスト作ってみました。 チョコボールのカンパ募集中です。
チョコボールをカンパする

森永製菓 チョコボール<ピーナッツ> 28g×20箱

森永製菓 チョコボール<ピーナッツ> 28g×20箱

で、結局、チョコボール何個買ったらエンゼルが出るの?

概要

  • 当ブログではベイズ推定によってエンゼルの出現確率の分布を推定しています
  • 出現確率が予測できたら、次は、いくら投資すればおもちゃのカンヅメがもらえるのか知りたいですよね?
  • ということで、チョコボールを何個買ったらエンゼルが当たるのかを計算してみます
  • 結果は、85個買えば銀のエンゼルが5個当たる確率が50%を超えます。85箱買うのだ!

【トップに戻る】

目次

はじめに

当ブログでは、日々チョコボールの金/銀エンゼルの出現確率を予測しています。 最近の調査では、銀のエンゼルについては、予測分布がだいぶ収束してきた感じがします。 金のエンゼルはキャンペーンの期間外ではまだ出てない。。。

しかし、エンゼルの出現確率の分布が予測できただけでは不十分ですよね?
おもちゃのカンヅメを手に入れることが目的の我々にとっては、 エンゼルの出現確率分布が予測できたなら次は、 おもちゃのカンヅメが得られるまでに投資するべき金額の見積もりを立てる必要があります。 じゃないと予算申請できません(どこに?)

というわけで今回は、 エンゼル出現確率の予測分布に基づいて、 チョコボールを何個買ったらおもちゃのカンヅメが手に入るのかを計算してみます。

【トップに戻る】

先行調査との比較

2018/01/28の記事で、金のエンゼルが78箱目で得られた確率について計算していました。
chocolate-ball.hatenablog.com

この時の調査では、幾何分布のパラメータを一般の先行研究のデータを参考にして、 点推定したものを利用していました。
現在、エンゼルの予測分布がかなり収束してきているので、予測分布を活用して 何個買えばエンゼルが当たるのかを計算してみます。

点推定、分布推定については、 LTしたときの資料や下記参考資料、 webなどにたくさん解説があるので、そちらを参考にしてください。

【トップに戻る】

当たりをk回引くまでにx回の試行が必要となる確率

この問題は、チョコボールだけでなく、ガチャや宝くじのように確率的な事象に対して いくら投資すれば当たりが期待できるのか見積もりを立てるために応用できます。
クジ全般は確率的な事象なので、以下の計算で赤字になるようなら手を出さないのが鉄則です。 もしも確率以外の現象を仮定するなら、つまり、そのクジは不正に操作されていると言えるので、 もっと手を出しちゃいけませんね。

基本

以前の記事にも書いたように、 成功(確率p)か失敗(確率1-p)という二値の事象で、 毎回の確率pは変化しないという試行のことをベルヌーイ試行といいます。 チョコボールのエンゼルなどクジは基本的にこの仮定を前提にしています。

ベルヌーイ試行に従う事象で成功が起こるまでの回数の分布は幾何分布(下式)に従うことが知られています。
{ \displaystyle
p(n) = (1-p)^{n-1}p
}

金のエンゼルは一つで良いので幾何分布を考えれば良いのですが、 銀のエンゼルは5つ当たらなければおもちゃのカンヅメと交換できません。 そこで、「負の二項分布(下式)」を利用します。
{ \displaystyle
f(x) = {}_{k+x-1}C_{x}p^{k}(1-p)^{x}
}

ここで、pは成功確率、kは成功数、xはk回成功するまでの失敗回数です。 幾何分布と比較すると、k回成功するのでp^{k}となっていますね。 この式から、k=1であれば幾何分布と同様の式になることがわかります。 ただし負の二項分布では、k回成功するまでの失敗回数の分布であることがちょっと違うので注意が必要です。

負の二項分布について詳しくは「参考文献1」などを参照ください。 成功確率p=10%、成功数5回の場合の負の二項分布の例を図1に示します。

f:id:hippy-hikky:20180708000055p:plain
図1. 成功確率p=10%で5回の成功が得られるまでの失敗数の分布。左が確率密度、右が累積確率。

図1右の累積確率をみることで、 確率が10%であるなら、41回までの失敗回数で5回の成功が得られる確率が50%を超えるということがわかりました。

このように、確率が一点でわかっている場合には、負の二項分布から直接失敗回数の見積もりを立てることができます。

確率が分布で与えられる場合

成功確率を点推定で予測する場合には、上記の通り負の二項分布を直接利用すれば良いです。 しかし、ここでは予測の曖昧さも含めて成功確率が分布として得られていた場合にどのように計算すれば良いのかを考えます。

成功確率\thetaでk回の成功が得られるまでの失敗回数xの分布をf(x | \theta, k)とします (この分布は上記の通り負の二項分布で表現されます)。
次に、データDを用いて予測した成功確率\thetaの事後分布をp(\theta | D)とします (これはいつも当ブログでベイズ推定により求めているエンゼル確率の予測分布などです)。

成功確率が一点で与えられない場合には、成功確率の確率を用いて期待値を計算すれば良いです。 つまり、これら2つの分布の積をとり\thetaに対して周辺化(積分)します。

{ \displaystyle
f(x | D) = \int p(\theta | D) f(x | \theta, k) d\theta
}

積分の簡略化

上記の式を計算すれば、試行回数の見積もりが立てられることがわかったのですが、 実際に計算しようとすると大きな2つの問題があります。

  • p(\theta | D)の式ってなに?
  • 積分の計算が超大変(かも)

一つ目の問題は、\thetaの予測方法に依存しています。 ベイズ推定で事後分布がよく知られた分布になるように共役事前分布を利用すれば数式で分布を表現できます。 しかし、MCMCなどで予測した場合には、数式で表現できません。

二つ目の問題は、p(\theta | D)がたとえよく知られた分布であったとしても、 負の二項分布との積の積分を解析的に解くのはなかなか面倒です。

これらの問題は、実は、確率の事後分布p(\theta | D)MCMC法によって サンプルの集合で表現することで、簡単に解決することができます。

MCMCで得られたサンプルは、事後分布p(\theta | D)を離散値の集合で表現します。 そこで、前節に書いた周辺化(積分)は、MCMCサンプル毎の平均で近似されるのです。

{ \displaystyle
f(x | D) \simeq \frac{1}{N}\sum_{\theta \in \Theta} f(x | \theta, k)
}
ここで、\Theta\thetaのサンプル集合、Nはサンプルの数を表します。

これで複雑な積分を計算する手間が削減されました。 たとえ共役事前分布を使って解析的にベイズ推定しても、その後の計算が複雑になってしまう場合はよくあると思います。 なので、近似ではありますが、 MCMCのようなサンプルベースの手法で分布の表現を持っておくと便利な場合があります。

【トップに戻る】

エンゼル出現までの購入個数の予測

ここまでで、 エンゼルが出るまでに買わなければならないチョコボールの個数の見積もり方について解説してきました。 実際に現時点で得られたデータを使って、いくら買えば良いのか見積もって見ましょう。

データ概要

2018年7月8日現在、金のエンゼル二倍キャンペーンを除いたデータ数は以下の通りです。

項目
計測データ数 156
銀のエンゼル出現数 8
金のエンゼル出現数 0

最尤推定量(サンプル平均)では、 銀のエンゼルの出現確率は、5.13%、 金のエンゼルは0.00%となっています。

このデータをベイズ推定で予測すると、事後分布は以下のようになります。

f:id:hippy-hikky:20180708151627p:plain
図2. ベイズ推定によるエンゼル出現確率の事後分布。左が銀のエンゼルの出現確率。右が金のエンゼルの出現確率。

銀のエンゼルの場合

上記の銀のエンゼル出現確率の事後分布を利用して、 銀のエンゼルが5つ当たるまでに必要なチョコボールの購入数を計算します。 計算に利用したコードは下に載せておきますので、必要な方はご参照ください。 5章に書いたことをそのままコードにしているだけですが。

結果は以下のとおりでした。

f:id:hippy-hikky:20180708152300p:plain
図3. 銀のエンゼルの予測分布に基づき5個当たるまでに必要な失敗数の分布。左が確率密度。右が累積確率。

図3右の赤の縦線は、累積確率が50%を超えるラインを示していて、80個でした。 つまり、85箱チョコボールを買えば、銀のエンゼルが5つ当たる確率が50%を超えるということです。 (失敗数80個+銀のエンゼル5箱=85箱)

一応、最尤推定量を使って計算した結果も図4に載せます。

f:id:hippy-hikky:20180708152332p:plain
図4. 銀のエンゼルの点推定結果に基づき5個当たるまでに必要な失敗数の分布。左が確率密度。右が累積確率。

図4右の累積確率が50%を超えるラインは86箱であり、ベイズ推定の結果と比較すると少し弱気な予測となりました。 おそらく、ベイズ推定による事後分布が少し右に偏った分布になっているからと考えています。

金のエンゼルの場合

続いて、金のエンゼルの場合にも同様の計算をしてみましょう。 はたして、いくら買えば金は当たるのでしょうか??

f:id:hippy-hikky:20180708153243p:plain
図5. 金のエンゼルの予測分布に基づき金のエンゼルが当たるまでに必要な失敗数の分布。左が確率密度。右が累積確率。

図5右の累積確率が50%を超えるラインは147箱でした。 が、結局金のエンゼルはまだ一つも出ておらず、これはだいぶ甘い予測になっていると考えられます。

以下に、以前の記事で引用させて頂いた 他の方の調査を載せた3つの記事の結果を使ってみます。 上記の記事をまとめると、金のエンゼルの出現確率の最尤推定量は0.14%となります。

f:id:hippy-hikky:20180708153705p:plain
図4. 金のエンゼルの点推定結果に基づき金のエンゼルが当たるまでに必要な失敗数の分布。左が確率密度。右が累積確率。図3とはx軸のレンジが異なっていることに注意。

累積確率が50%超えるまでには、494箱必要との見積もりになります。

金のエンゼルについては、 そもそも一つも金のエンゼルが出ていない状態ですので予測分布が収束していないものと思われます。 予測の事前分布に一様分布を用いているのも問題で、 先行研究で0.14%という予測が出ているので、事前分布をもう少し工夫するべきかなと思っています。

実験コード

【トップに戻る】

まとめ

ということで、予測分布を利用して、何個買えばエンゼルが現れるのかを見積もってみたところ、 50%の確率で5枚の銀のエンゼルを獲得するためには85箱必要との結果が出ました。 チョコボールのコンビニ価格は80円なので、6,800円の投資でおもちゃのカンヅメが手に入りそうです!

ここで冷静に、メルカリでおもちゃのカンヅメがいくらで出ているのか調べて見ました。
だいたい1,000円〜2,000円ですね(笑)

どこに価値を置くかが重要なので、メルカリでいくらで出品されているかなんて関係ないんです!

【トップに戻る】

おわりに

ということで、予測分布を利用して、何個買えばエンゼルが現れるのかを見積もってみました。 これで根拠を持って予算申請できますね!

また、 MCMCを使って事後確率を予測していたことで、 積分が簡単になるメリットが享受できましたね。 どうしても厳密解にこだわらないのであれば、 MCMCのようなサンプルベースの手法が応用しやすくて良いのかなと思っています。 サンプル数を増やせば厳密解に近づいていきますし。

【トップに戻る】

参考文献

  1. 統計学入門

    統計学入門 (基礎統計学?)

    統計学入門 (基礎統計学?)

    有名な統計学の入門書。かっちり書いてあるので、統計学を学ぶにはすごく良い。

  2. 自然科学の統計学

    自然科学の統計学 (基礎統計学)

    自然科学の統計学 (基礎統計学)

    ベイズ推定はもちろん、回帰など統計の醍醐味がいろいろ載っている。 上記の入門書を読んでから入ると良いかも。

  3. データ解析のための統計モデリング入門

    通称「みどり本」。 MCMCについての理論的な部分はこのくらいは読んでおけば良いかも。

  4. Pythonで体験するベイズ推論

    Pythonで体験するベイズ推論 PyMCによるMCMC入門

    Pythonで体験するベイズ推論 PyMCによるMCMC入門

    pythonのMCMCパッケージであるpymcの使い方を学ぶにはちょうどよい。 ただし、プログラムとして書くためのノウハウが主なので、理論的な部分はほとんど載っていない。 理論は上記の参考書などを読むと良い。

【トップに戻る】

第79回 チョコボール計測

今日の計測報告です。
同僚から4箱もらったうちの3,4箱目を開封します。

計測結果

date best_before weight box_weight number factory shop angel net_weight mean_weight
2018-07-04 2019-01-01 34.913 4.937 17 小山工場 ドラックストア(横浜市 なし 29.976 1.763
2018-07-04 2019-01-01 34.303 4.94 15 小山工場 ドラックストア(横浜市 なし 29.363 1.958

エンゼルは今日も現れぬ。。。

基礎集計

この集計はピーナツ味のチョコボールの集計結果です。

項目
計測データ数 232
銀のエンゼル出現数 6
金のエンゼル出現数 1
最小 中央値 最大値 平均
正味重量 28.913 29.396 32.232 29.501
個数 14.000 16.000 20.000 16.534

f:id:hippy-hikky:20180704233455p:plain この図は正味の重量のヒストグラムです。 赤い縦線が仕様(28g)を表しています。 青い太線で正規分布と仮定した最尤推定量をプロットしています。

エンゼル出現確率の予測

通常版のエンゼルの予測を行っていきます。 これまでの通常版パッケージの開封結果は次の通りです。

項目
計測データ数 156
銀のエンゼル出現数 8
金のエンゼル出現数 0

この結果を使ってベイズ推定によるエンゼルの出現確率推定を行います。

はじめに銀のエンゼルの出現確率の推定です。
90%信用区間(上下それぞれ5%)は、下側が2.97%上側が8.98%という結果です。
f:id:hippy-hikky:20180704233537p:plain

現在の期待値は5.68%です。

次に金のエンゼルの出現確率の推定です。 90%信用区間(上下それぞれ5%)は、下側が0.03%、上側が1.78%という予測になっています。
f:id:hippy-hikky:20180704233611p:plain

広告

Amazonの欲しいものリスト作ってみました。 チョコボールのカンパ募集中です。
チョコボールをカンパする

森永製菓 チョコボール<ピーナッツ> 28g×20箱

森永製菓 チョコボール<ピーナッツ> 28g×20箱

第78回 チョコボール計測

今日の計測報告です。
同僚から4箱もらったので、そのうち2箱を開封します。

計測結果

date best_before weight box_weight number factory shop angel net_weight mean_weight
2018-07-02 2019-01-01 34.568 4.937 15 小山工場 ドラックストア(横浜市 なし 29.631 1.975
2018-07-02 2019-01-01 34.321 4.913 17 小山工場 ドラックストア(横浜市 なし 29.408 1.730

今日もエンゼルは現れませんでした。

ここ最近暑いので、チョコボールを置いておくと溶けちゃうので危険です。

基礎集計

この集計はピーナツ味のチョコボールの集計結果です。

項目
計測データ数 230
銀のエンゼル出現数 6
金のエンゼル出現数 1
最小 中央値 最大値 平均
正味重量 28.913 29.396 32.232 29.500
個数 14.000 16.000 20.000 16.539

f:id:hippy-hikky:20180702235644p:plain この図は正味の重量のヒストグラムです。 赤い縦線が仕様(28g)を表しています。 青い太線で正規分布と仮定した最尤推定量をプロットしています。

エンゼル出現確率の予測

通常版のエンゼルの予測を行っていきます。 これまでの通常版パッケージの開封結果は次の通りです。

項目
計測データ数 154
銀のエンゼル出現数 8
金のエンゼル出現数 0

この結果を使ってベイズ推定によるエンゼルの出現確率推定を行います。

はじめに銀のエンゼルの出現確率の推定です。
90%信用区間(上下それぞれ5%)は、下側が3.14%、上側が9.30%という結果です。
f:id:hippy-hikky:20180702235715p:plain

現在の期待値は5.86%です。

次に金のエンゼルの出現確率の推定です。 90%信用区間(上下それぞれ5%)は、下側が0.04%、上側が1.91%という予測になっています。
f:id:hippy-hikky:20180702235748p:plain

広告

Amazonの欲しいものリスト作ってみました。 チョコボールのカンパ募集中です。
チョコボールをカンパする

森永製菓 チョコボール<ピーナッツ> 28g×20箱

森永製菓 チョコボール<ピーナッツ> 28g×20箱

第77回 チョコボール計測

今日の計測報告です。
先日同僚からもらったうちの最後の2箱を開封します。

計測結果

date best_before weight box_weight number factory shop angel net_weight mean_weight
2018-06-29 2018-11-01 34.338 4.951 16 小山工場 ドラッグストア(上尾市 なし 29.387 1.837
2018-06-29 2018-11-01 34.367 4.928 16 小山工場 ドラッグストア(上尾市 なし 29.439 1.840

今回はエンゼルさん出ませんでした!

基礎集計

この集計はピーナツ味のチョコボールの集計結果です。

項目
計測データ数 228
銀のエンゼル出現数 6
金のエンゼル出現数 1
最小 中央値 最大値 平均
正味重量 28.913 29.395 32.232 29.500
個数 14.000 16.000 20.000 16.544

f:id:hippy-hikky:20180630000840p:plain この図は正味の重量のヒストグラムです。 赤い縦線が仕様(28g)を表しています。 青い太線で正規分布と仮定した最尤推定量をプロットしています。

エンゼル出現確率の予測

通常版のエンゼルの予測を行っていきます。 これまでの通常版パッケージの開封結果は次の通りです。

項目
計測データ数 152
銀のエンゼル出現数 8
金のエンゼル出現数 0

この結果を使ってベイズ推定によるエンゼルの出現確率推定を行います。

はじめに銀のエンゼルの出現確率の推定です。
90%信用区間(上下それぞれ5%)は、下側が3.08%、上側が9.07%という結果です。
f:id:hippy-hikky:20180630000914p:plain

現在の期待値は5.72%です。

次に金のエンゼルの出現確率の推定です。 90%信用区間(上下それぞれ5%)は、下側が0.04%、上側が1.88%という予測になっています。
f:id:hippy-hikky:20180630000939p:plain

広告

Amazonの欲しいものリスト作ってみました。 チョコボールのカンパ募集中です。
チョコボールをカンパする

森永製菓 チョコボール<ピーナッツ> 28g×20箱

森永製菓 チョコボール<ピーナッツ> 28g×20箱

第76回 チョコボール計測

今日の計測報告です。
前回同僚からもらったうちの3箱を開封します。

これの他に残り2箱あるので、それは後日開封します。

計測結果

date best_before weight box_weight number factory shop angel net_weight mean_weight
2018-06-28 2018-11-01 34.423 4.858 17 小山工場 ドラッグストア(上尾市 なし 29.565 1.739
2018-06-28 2018-11-01 34.402 4.925 16 小山工場 ドラッグストア(上尾市 なし 29.477 1.842
2018-06-28 2018-11-01 35.128 5.001 16 小山工場 ドラッグストア(上尾市 なし 30.127 1.883

エンゼルさん出ませんでした〜

基礎集計

この集計はピーナツ味のチョコボールの集計結果です。

項目
計測データ数 226
銀のエンゼル出現数 6
金のエンゼル出現数 1
最小 中央値 最大値 平均
正味重量 28.913 29.395 32.232 29.501
個数 14.000 16.000 20.000 16.549

f:id:hippy-hikky:20180628155605p:plain この図は正味の重量のヒストグラムです。 赤い縦線が仕様(28g)を表しています。 青い太線で正規分布と仮定した最尤推定量をプロットしています。

エンゼル出現確率の予測

通常版のエンゼルの予測を行っていきます。 これまでの通常版パッケージの開封結果は次の通りです。

項目
計測データ数 150
銀のエンゼル出現数 8
金のエンゼル出現数 0

この結果を使ってベイズ推定によるエンゼルの出現確率推定を行います。

はじめに銀のエンゼルの出現確率の推定です。
90%信用区間(上下それぞれ5%)は、下側が3.28%、上側が9.41%という結果です。
f:id:hippy-hikky:20180628155635p:plain

現在の期待値は5.95%です。

次に金のエンゼルの出現確率の推定です。 90%信用区間(上下それぞれ5%)は、下側が0.03%、上側が1.88%という予測になっています。
f:id:hippy-hikky:20180628155708p:plain

広告

Amazonの欲しいものリスト作ってみました。 チョコボールのカンパ募集中です。
チョコボールをカンパする

森永製菓 チョコボール<ピーナッツ> 28g×20箱

森永製菓 チョコボール<ピーナッツ> 28g×20箱