チョコボール統計

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

第83回 チョコボール計測

今日の計測報告です。
近所のスーパーに行ったら、一個50円弱で安売りされていたので、何個かまとめて買いました。
今日はそのうちの2個を開封します。

計測結果

date best_before weight box_weight number factory shop angel net_weight mean_weight
2018-08-05 2019-02-01 34.272 4.802 17 小山工場 スーパー(さいたま市 なし 29.470 1.734
2018-08-05 2019-02-01 34.194 4.81 16 小山工場 スーパー(さいたま市 なし 29.384 1.837

今日もエンゼルさんはいませんでした。 最近全然出ないですね。。。

これで前回エンゼルを見てから、43連続で外してます。

基礎集計

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

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

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

エンゼル出現確率の予測

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

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

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

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

現在の期待値は4.96%です。
期待値5%切っちゃった。

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

広告

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

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

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

第ex13回 チョコボール計測

番外編の計測、13回目。
近所のスーパーで新作を発見したので購入しました。 「ラムネ」味です。

パッケージによると、ラムネをラムネ味のチョコレートで包んだらしいです。

実際食べてみましたが、 中のラムネは球状に加工されたものではなく、あのラムネそのままでした。
パッケージ画像の通り。

ラムネって味の主張が強いですね!

計測結果

では計測結果です。

date best_before weight box_weight number factory shop angel net_weight mean_weight
2018-08-04 2019-05-01 29.448 4.864 17 小山工場 スーパー(さいたま市 なし 24.584 1.446
2018-08-04 2019-05-01 29.319 4.888 17 小山工場 スーパー(さいたま市 なし 24.431 1.437
2018-08-04 2019-05-01 28.660 4.921 16 小山工場 スーパー(さいたま市 なし 23.739 1.484

エンゼルは今回も現れませんでした。
これで41連続で外してます。。。

基礎集計

「ラムネ」味の集計です

項目
計測データ数 3
銀のエンゼル出現数 0
金のエンゼル出現数 0
最小 中央値 最大値 平均
正味重量 23.739 24.431 24.584 24.251
個数 16.000 17.000 17.000 16.667

重量分布

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

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

エンゼルの出現確率推定

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

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

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

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

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

広告

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

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

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

第82回 チョコボール計測 + おまけ

今日の計測報告です。
同僚からもらいました。

計測結果

date best_before weight box_weight number factory shop angel net_weight mean_weight
2018-08-03 2019-01-01 34.199 4.877 17 小山工場 ドラックストア(横浜市 なし 29.322 1.725
2018-08-03 2019-01-01 34.287 4.852 16 小山工場 ドラックストア(横浜市 なし 29.435 1.840

エンゼルさんはいませんでした!

これで38連続でエンゼルさん外してます。 だいぶついてないですね。

基礎集計

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

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

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

エンゼル出現確率の予測

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

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

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

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

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

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

おまけ:次の箱でエンゼルが出てくる確率は?

最近しばらくエンゼル見てないようなので、次買ったチョコボールでエンゼルさんが出てくる確率を予想してみましょう。

アプローチ

手法は、この記事の通りです。

chocolate-ball.hatenablog.com

予想結果

前回エンゼルが出た時から連続で38回外しています。
そこで、次(前回の当たりから39箱目)にエンゼルが出てくる確率を計算してみます。 つまり、1つエンゼルが出てくるために38回連続で外す確率を負の二項分布を使って計算するということです。 エンゼルの出現確率は、簡単のために5%とします。

計算結果は、86%! ということで、次こそ出るか!?

。。。

というのは、間違いでしょう。
純粋に38個だけを選んで全て外れていたならば上記の予測は正しいのですが、 実際は174個買っていてエンゼルが8個出ています(8個のエンゼルを得るまでに166回の失敗をしている)。 これを考慮する必要があります。

では、166回の失敗で9個のエンゼルが出る、 つまり、次にエンゼルが出る確率を計算してみましょう。

結果は、51%ということで、 次当たるかはほぼ半々ということですね。

ちなみに、9個のエンゼルが当たる確率が90%を超えるためには、 249回の失敗、つまり、あと83個買う必要がありそうです。

広告

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

第81回 チョコボール計測

今日の計測報告です。
近所のコンビニで3つ手に入れてきました。

ピーナツ味もダブルチャンスで50倍チョコボールが当たるキャンペーンの商品が出まわり始めたみたいです。 (箱を展開後に写真撮ったので、一部見苦しいですが許してください。。。)

計測結果

date best_before weight box_weight number factory shop angel net_weight mean_weight
2018-08-01 2019-02-01 34.125 4.87 16 小山工場 コンビニ(さいたま市 なし 29.255 1.828
2018-08-01 2019-02-01 33.9 4.888 17 小山工場 コンビニ(さいたま市 なし 29.012 1.707
2018-08-01 2019-02-01 34.346 4.88 17 小山工場 コンビニ(さいたま市 なし 29.466 1.733

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

基礎集計

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

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

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

エンゼル出現確率の予測

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

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

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

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

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

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

金のエンゼルの予測分布に関しては、最近の変動がなくなってきました。 1%の確率であっても、172個連続で外す確率はそこそこあるので、 もっともっとデータを集めないと分布に変化が出て来ないでしょう。
レアな事象のデータ収集はお金も時間もかかりますね。

広告

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

第80回 チョコボール計測

今日の計測報告です。
同僚から4箱もらったので開封します。

久々のピーナッツ!

計測結果

date best_before weight box_weight number factory shop angel net_weight mean_weight
2018-07-28 2019-01-01 34.159 4.878 17 小山工場 ドラックストア(横浜市 なし 29.281 1.722
2018-07-28 2019-01-01 34.568 4.912 18 小山工場 ドラックストア(横浜市 なし 29.656 1.648
2018-07-28 2019-01-01 34.178 4.918 17 小山工場 ドラックストア(横浜市 なし 29.260 1.721
2018-07-28 2019-01-01 35.065 4.904 17 小山工場 ドラックストア(横浜市 なし 30.161 1.774

最近エンゼルさん見ないっすね〜

基礎集計

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

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

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

エンゼル出現確率の予測

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

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

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

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

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

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

広告

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

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

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

第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箱