チョコボール統計

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

第102回 チョコボール計測

ちょっと時間空いてしまいましたが、 この間に以下のイベントに本ブログ関連の作品を出展してました。

ma2018.we-are-ma.jp

https://festa2018.we-are-ma.jp/festa2018.we-are-ma.jp

イベントについては、別途書くかもしれません。

ということで本日の計測です。 今日は6箱計測します。

計測結果

date best_before weight box_weight number factory shop angel net_weight mean_weight
2018-12-05 2019-05-01 33.998 4.840 18 小山工場 ドラックストア(横浜市 なし 29.158 1.620
2018-12-05 2019-05-01 34.180 4.913 17 小山工場 ドラックストア(横浜市 なし 29.267 1.722
2018-12-05 2019-05-01 34.001 4.904 17 小山工場 ドラックストア(横浜市 なし 29.097 1.712
2018-12-05 2019-05-01 34.071 4.865 16 小山工場 コンビニ(東京都) なし 29.206 1.825
2018-12-05 2019-07-01 34.139 4.871 17 小山工場 コンビニ(東京都) 29.268 1.722
2018-12-05 2019-05-01 34.040 4.857 16 小山工場 コンビニ(東京都) なし 29.183 1.824

銀のエンゼルさん出てくれました〜

基礎集計

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

項目
計測データ数 297
銀のエンゼル出現数 8
金のエンゼル出現数 1
最小 中央値 最大値 平均
正味重量 28.750 29.389 32.232 29.490
個数 14.000 16.000 20.000 16.559

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

エンゼル出現確率の予測

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

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

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

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

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

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

金のエンゼルの出現確率の予測分布がだいぶ収束してきて予測値が動かなくなってきました。
もう少しMCMCのサンプル数を増やして表現力を上げたほうが良いかも。

広告

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