第147回 チョコボール計測(チョコバナナ)
本日の計測結果を報告します。 Amazonで箱買いしたチョコバナナ味を4箱計測します。
Amazon箱買いのチョコバナナ味、これで最後でした。 チョコバナナ味はもう販売終了したので、たぶんこれで最後でしょう。
計測結果
date | best_before | weight | box_weight | number | factory | shop | angel | net_weight | mean_weight |
---|---|---|---|---|---|---|---|---|---|
2019-04-18 | 2019-11-01 | 28.683 | 4.860 | 16 | 小山工場 | amazon | なし | 23.823 | 1.489 |
2019-04-18 | 2019-11-01 | 28.117 | 4.842 | 15 | 小山工場 | amazon | なし | 23.275 | 1.552 |
2019-04-18 | 2019-11-01 | 27.769 | 4.851 | 15 | 小山工場 | amazon | なし | 22.918 | 1.528 |
2019-04-18 | 2019-11-01 | 28.079 | 4.867 | 15 | 小山工場 | amazon | なし | 23.212 | 1.547 |
今回はエンゼルさんは現れませんでした。
基礎集計
チョコバナナ味の集計結果です。
項目 | 数 |
---|---|
計測データ数 | 43 |
銀のエンゼル出現数 | 2 |
金のエンゼル出現数 | 0 |
最小 | 中央値 | 最大値 | 平均 | |
---|---|---|---|---|
正味重量 | 22.107 | 22.915 | 24.403 | 22.929 |
個数 | 14.000 | 15.000 | 16.000 | 15.023 |
この図は正味の重量のヒストグラムです。 赤い縦線が仕様(22g)を表しています。 青い太線で正規分布と仮定した最尤推定量をプロットしています。
チョコバナナ味のここまでの計測結果の分析は、後ほど記事としてまとめたいと思います。
エンゼル出現確率の予測
通常版のエンゼルの予測を行っていきます。 これまでの通常版パッケージの開封結果は次の通りです。
項目 | 数 |
---|---|
計測データ数 | 417 |
銀のエンゼル出現数 | 19 |
金のエンゼル出現数 | 0 |
この結果を使ってベイズ推定によるエンゼルの出現確率推定を行います。
はじめに銀のエンゼルの出現確率の推定です。
90%信用区間(上下それぞれ5%)は、下側が3.20%上側が6.66%という結果です。
現在の期待値は4.77%です。
次に金のエンゼルの出現確率の推定です。
90%信用区間(上下それぞれ5%)は、下側が0.01%、上側が0.73%という予測になっています。
今のところの金のエンゼルの期待値は0.24%です。
広告
Amazonの欲しいものリスト作ってみました。
チョコボールのカンパ募集中です。
チョコボールをカンパする
- 出版社/メーカー: 森永製菓
- 発売日: 2016/03/01
- メディア: 食品&飲料
- この商品を含むブログを見る
SIMERST 改良版 携帯タイプはかり ポケットデジタルスケール(秤) 0.001g-100g精密 業務用(プロ用) デジタルスケール 電子天秤
- 出版社/メーカー: Simerst
- メディア:
- この商品を含むブログを見る
第146回 チョコボール計測(チョコバナナ)
本日の計測結果を報告します。 今日もAmazonで箱買いしたチョコバナナ味です。
チョコバナナ味、次回で終了かな(残り4箱)。
計測結果
date | best_before | weight | box_weight | number | factory | shop | angel | net_weight | mean_weight |
---|---|---|---|---|---|---|---|---|---|
2019-04-17 | 2019-11-01 | 27.899 | 4.831 | 15 | 小山工場 | amazon | なし | 23.068 | 1.538 |
2019-04-17 | 2019-11-01 | 27.723 | 4.864 | 16 | 小山工場 | amazon | なし | 22.859 | 1.429 |
2019-04-17 | 2019-11-01 | 27.527 | 4.845 | 15 | 小山工場 | amazon | なし | 22.682 | 1.512 |
2019-04-17 | 2019-11-01 | 29.232 | 4.829 | 16 | 小山工場 | amazon | なし | 24.403 | 1.525 |
今回はエンゼルさんは現れませんでした。
基礎集計
チョコバナナ味の集計結果です。
項目 | 数 |
---|---|
計測データ数 | 39 |
銀のエンゼル出現数 | 2 |
金のエンゼル出現数 | 0 |
最小 | 中央値 | 最大値 | 平均 | |
---|---|---|---|---|
正味重量 | 22.107 | 22.872 | 24.403 | 22.891 |
個数 | 14.000 | 15.000 | 16.000 | 15.000 |
この図は正味の重量のヒストグラムです。 赤い縦線が仕様(22g)を表しています。 青い太線で正規分布と仮定した最尤推定量をプロットしています。
エンゼル出現確率の予測
通常版のエンゼルの予測を行っていきます。 これまでの通常版パッケージの開封結果は次の通りです。
項目 | 数 |
---|---|
計測データ数 | 413 |
銀のエンゼル出現数 | 19 |
金のエンゼル出現数 | 0 |
この結果を使ってベイズ推定によるエンゼルの出現確率推定を行います。
はじめに銀のエンゼルの出現確率の推定です。
90%信用区間(上下それぞれ5%)は、下側が3.22%上側が6.65%という結果です。
現在の期待値は4.82%です。
次に金のエンゼルの出現確率の推定です。
90%信用区間(上下それぞれ5%)は、下側が0.01%、上側が0.72%という予測になっています。
今のところの金のエンゼルの期待値は0.24%です。
広告
Amazonの欲しいものリスト作ってみました。
チョコボールのカンパ募集中です。
チョコボールをカンパする
第145回 チョコボール計測(チョコバナナ)
本日の計測結果を報告します。 ようやくメインのUbuntuマシンのセッティング(し直し)がほぼほぼ終わったので、平常運転で計測できそう。
計測結果
date | best_before | weight | box_weight | number | factory | shop | angel | net_weight | mean_weight |
---|---|---|---|---|---|---|---|---|---|
2019-04-08 | 2019-11-01 | 27.804 | 4.875 | 15 | 小山工場 | amazon | なし | 22.929 | 1.529 |
2019-04-08 | 2019-11-01 | 27.835 | 4.855 | 15 | 小山工場 | amazon | なし | 22.980 | 1.532 |
2019-04-08 | 2019-11-01 | 28.407 | 4.815 | 16 | 小山工場 | amazon | 銀 | 23.592 | 1.474 |
2019-04-08 | 2019-11-01 | 27.777 | 4.835 | 16 | 小山工場 | amazon | なし | 22.942 | 1.434 |
久しぶりに銀のエンゼルさん出てくれました〜!
基礎集計
チョコバナナ味の集計結果です。
項目 | 数 |
---|---|
計測データ数 | 35 |
銀のエンゼル出現数 | 2 |
金のエンゼル出現数 | 0 |
最小 | 中央値 | 最大値 | 平均 | |
---|---|---|---|---|
正味重量 | 22.107 | 22.872 | 24.338 | 22.849 |
個数 | 14.000 | 15.000 | 16.000 | 14.943 |
この図は正味の重量のヒストグラムです。 赤い縦線が仕様(22g)を表しています。 青い太線で正規分布と仮定した最尤推定量をプロットしています。
エンゼル出現確率の予測
通常版のエンゼルの予測を行っていきます。 これまでの通常版パッケージの開封結果は次の通りです。
項目 | 数 |
---|---|
計測データ数 | 409 |
銀のエンゼル出現数 | 19 |
金のエンゼル出現数 | 0 |
この結果を使ってベイズ推定によるエンゼルの出現確率推定を行います。
はじめに銀のエンゼルの出現確率の推定です。
90%信用区間(上下それぞれ5%)は、下側が3.26%上側が6.72%という結果です。
現在の期待値は4.87%です。
次に金のエンゼルの出現確率の推定です。
90%信用区間(上下それぞれ5%)は、下側が0.01%、上側が0.72%という予測になっています。
今のところの金のエンゼルの期待値は0.25%です。
広告
Amazonの欲しいものリスト作ってみました。
チョコボールのカンパ募集中です。
チョコボールをカンパする
- 出版社/メーカー: 森永製菓
- 発売日: 2016/03/01
- メディア: 食品&飲料
- この商品を含むブログを見る
SIMERST 改良版 携帯タイプはかり ポケットデジタルスケール(秤) 0.001g-100g精密 業務用(プロ用) デジタルスケール 電子天秤
- 出版社/メーカー: Simerst
- メディア:
- この商品を含むブログを見る
第144回 チョコボール計測(チョコバナナ)
本日の計測結果を報告します。 今週、マシンの調子が悪く、ログインできない状態だったので、2回分の計測をまとめます。
計測結果
date | best_before | weight | box_weight | number | factory | shop | angel | net_weight | mean_weight |
---|---|---|---|---|---|---|---|---|---|
2019-04-04 | 2019-11-01 | 28.002 | 4.869 | 15 | 小山工場 | amazon | なし | 23.133 | 1.542 |
2019-04-04 | 2019-11-01 | 27.763 | 4.857 | 15 | 小山工場 | amazon | なし | 22.906 | 1.527 |
2019-04-04 | 2019-11-01 | 28.043 | 4.872 | 15 | 小山工場 | amazon | なし | 23.171 | 1.545 |
2019-04-04 | 2019-11-01 | 27.914 | 4.848 | 15 | 小山工場 | amazon | なし | 23.066 | 1.538 |
エンゼルさん現れず。
基礎集計
チョコバナナ味の集計結果です。
項目 | 数 |
---|---|
計測データ数 | 31 |
銀のエンゼル出現数 | 1 |
金のエンゼル出現数 | 0 |
最小 | 中央値 | 最大値 | 平均 | |
---|---|---|---|---|
正味重量 | 22.107 | 22.814 | 24.338 | 22.815 |
個数 | 14.000 | 15.000 | 16.000 | 14.871 |
この図は正味の重量のヒストグラムです。 赤い縦線が仕様(22g)を表しています。 青い太線で正規分布と仮定した最尤推定量をプロットしています。
計測ミスなんじゃないかってくらいバラつきが大きいですね。
エンゼル出現確率の予測
通常版のエンゼルの予測を行っていきます。 これまでの通常版パッケージの開封結果は次の通りです。
項目 | 数 |
---|---|
計測データ数 | 405 |
銀のエンゼル出現数 | 18 |
金のエンゼル出現数 | 0 |
この結果を使ってベイズ推定によるエンゼルの出現確率推定を行います。
はじめに銀のエンゼルの出現確率の推定です。
90%信用区間(上下それぞれ5%)は、下側が3.05%上側が6.43%という結果です。
現在の期待値は4.67%です。
次に金のエンゼルの出現確率の推定です。
90%信用区間(上下それぞれ5%)は、下側が0.01%、上側が0.75%という予測になっています。
今のところの金のエンゼルの期待値は0.24%です。
広告
Amazonの欲しいものリスト作ってみました。
チョコボールのカンパ募集中です。
チョコボールをカンパする
- 出版社/メーカー: 森永製菓
- 発売日: 2016/03/01
- メディア: 食品&飲料
- この商品を含むブログを見る
SIMERST 改良版 携帯タイプはかり ポケットデジタルスケール(秤) 0.001g-100g精密 業務用(プロ用) デジタルスケール 電子天秤
- 出版社/メーカー: Simerst
- メディア:
- この商品を含むブログを見る
Chainer Meetup #09 でLT:チョコボール数の傾向を分析
【概要】
- 3/30に開催のChainer Meetupに参加してきました
- ただ参加するだけではつまらないので、LTしてきました
- 本記事ではLTの概要と資料をまとめます
【目次】
はじめに
今回で2回目の参加となりましたChainer Meetupですが、今回もとても興味深い話が聞けて良かったです。 Chainerの中の人の開発情報が聴けるというのはもちろん良かったのですが、 nVidiaの中の人の講演が今回もとても興味深かったです。 前回もパフォーマンスチューニングについてとても興味深く聞いていましたが、今回は自分でも取り組めそうな話で希望が持てました。 資料はconnpassのイベントページに公開されているので、興味ある方はご覧になられたらと思います。
本記事ですが、 前回に引き続き今回もLTをさせてもらったので、 どのような話をしたのか概要をまとめておくものです。
資料
当日使用した資料は以下です。
LT概要
今回のLTでは、チョコボールの個数を分析した話をしました。
チョコボール(ピーナッツ味)は、内容量が28gと定められています*1。 内容量(重量)が定められているのですが、 ピーナッツは自然のものであるため大きさや重さを厳密に一定にすることはできません*2。 そのため、一箱に入っているチョコボールの個数は箱ごとに異なります。
一方、せっかくチョコボールを買うなら美味しいものを食べたいのですが、 小さかったり形がいびつなチョコボールは味もいまいちな印象です(個人的な感想であり品質を保証するものではありません)。
そこで、チョコボールの粒が大きく揃っているものを購入するにはどのようにすれば良いのかについて、 個数の傾向から考察してみたという話です。
データと分析アプローチ
データ収集効率化の取り組み
(詳しくは上記の資料を参照ください)
個数の傾向を分析するには、各商品に個数がいくつ入っていたのかを計測し、 データを利用できる形に整理することが必要です。 しかし、一つ一つ数えるのはなかなかの手間です*3。
そこで、Faster R-CNNを使ってチョコボールを検出するWebサービスを開発しました。 これで、写真を撮るだけで個数が自動で計測できます! (ここまでは、前回のLTでお話ししたところです)
しかしながら実際のところ、 このシステムは計測の場面ではまだ実運用できておらず、 手で数えています*4。 ChainerのイベントでChainerについて一切触れないのはマズかろうということで(大人の事情)。
実験に利用したデータの概要
ということで手で数えて収集したデータですが、2017年11月から2019年3月まででピーナッツ味373箱分のデータを利用します。
このデータを賞味期限が同じチョコボールをひとまとめにして、個数の分布(箱ひげ図)を確認してみてみます。 なお、チョコボールには製造年月日の記載が無く、賞味期限が月単位で表示されています。 そのため、賞味期限が同じであれば、同じ時期に製造されたものであると仮定して、データを分析します。
上記の図からわかるように、2018年の8~10月の個数が少ないように見えますね。 また、2019年の8~10月もデータ数が少し少ないですが、同様に個数が少ないように見える気がします。
そのため今回は、賞味期限が8~10月の個数が他の期間と比べて有意に数が少ないのかを確認してみます。 この期間の数が少ないことが統計的に確認できたならば、 賞味期限が8~10月のチョコボールを狙って購入することで、粒の大きいチョコボールを手に入れられる可能性が高まる と期待できます。
他にも分析する軸は考えられますが、今回は時期による個数の変化にだけ注目します。
モデル
個数の傾向を分析するに当たって、下記の統計モデルを設定しました*5。
一箱に入っているチョコボールの個数yは、ポアソン分布に従う確率変数と考えます。 ポアソン分布のパラメータは各月の平均個数を示すパラメータであり、 ガンマ分布を事前分布とします。 ガンマ分布のパラメータは任意に設定するのではなく、 それぞれ階層事前分布を設定してデータから推定させることにします。 ガンマ分布のパラメータは、 平均と標準偏差を示すパラメータを使って以下のように算出されます。
今回は、8~10月の期間のデータから推定した事前分布()と、 その他の期間のデータから推定した事前分布()を比較し、 統計的に差があるのかを確認します。
結果
上記のモデルを使い、パラメータを推定しました。
パラメータは、MCMCアルゴリズム(NUTS)を利用して数値的に推定しました。
推定した二つのガンマ分布は以下のようなものです。
95%信用区間で分布の平均に差があるのかを確認するために、 平均のパラメータの事後分布の差の分布を以下に示します。
ということで、二つの分布に差がないことを示す0の位置が90%信用区間から外れているため、 二つの分布の平均には差があるであろうということがわかりました。
その差は、約1.3個。
以上の結果、賞味期限が8~10月のチョコボールは個数が少ないということがわかりました。 資料p.22にポアソン分布の各月のパラメータの事後分布を載せていますが、 2018年だけでなく、2019年の8~10月も同様に個数は少ないと言えそうです。 そのため、賞味期限が8~10月のチョコボールを狙って買えば、粒の大きいチョコボールを手に入れられそうです!
終わりに
ということでチョコボールの個数の傾向を分析し、 特定の月の個数に差がありそうだということまではわかりました。
しかし、今回は8~10月を予め切り分けてパラメータを推定したので、
8~10月だけが少ないのかはわからないです。
データから変化点を合わせて推定させることが必要と考えています*6。
(実は、変化点を推定する実験を最初にやってみたのですが、推定結果が安定せずでして、プログラム上に問題があるのか、モデルがまずいのか検討中の状態です。。。)
また、チョコボールの大きさが大きくて粒が揃っているというのは、 個数だけで特徴付けられるのかについて、はっきりとした確証はありません。 そのため、大きさのばらつきについても今後調査していきたいと考えています。
参考文献
データ解析のための統計モデリング入門
一般化線形モデルを軸に階層ベイズモデルまで丁寧に記載されている。 Rとstanを利用した実装例もあるが、実装部分は無視しても理解できる。データ解析のための統計モデリング入門――一般化線形モデル・階層ベイズモデル・MCMC (確率と情報の科学)
- 作者: 久保拓弥
- 出版社/メーカー: 岩波書店
- 発売日: 2012/05/19
- メディア: 単行本
- 購入: 16人 クリック: 163回
- この商品を含むブログ (29件) を見る
- pymc3を使った入門書。 理論面も多少書かれており、こちらで手を動かしながら学習するのが良いと思う。
Pythonによるベイズ統計モデリング: PyMCでのデータ分析実践ガイド
- 作者: オズワルドマーティン,Osvaldo Martin,金子武久
- 出版社/メーカー: 共立出版
- 発売日: 2018/06/22
- メディア: 単行本
- この商品を含むブログを見る
- pymcを利用した事例ベースの参考書。 事例が豊富で参考になるシーンは出てくるかもしれないが、階層モデルについては記載がない。 また、pymc2なので、pymc3と記法が少し違うので注意。
Pythonで体験するベイズ推論 PyMCによるMCMC入門
- 作者: キャメロンデビッドソン=ピロン,玉木徹
- 出版社/メーカー: 森北出版
- 発売日: 2017/04/06
- メディア: 単行本(ソフトカバー)
- この商品を含むブログ (1件) を見る
広告
Amazonの欲しいものリスト作ってみました。
チョコボールのカンパ募集中です。
チョコボールをカンパする
- 出版社/メーカー: 森永製菓
- 発売日: 2016/03/01
- メディア: 食品&飲料
- この商品を含むブログを見る
SIMERST 改良版 携帯タイプはかり ポケットデジタルスケール(秤) 0.001g-100g精密 業務用(プロ用) デジタルスケール 電子天秤
- 出版社/メーカー: Simerst
- メディア:
- この商品を含むブログを見る
*1:設定値は28gですが、実際の重量傾向についてはこちらの記事を参照ください
*2:製造時のばらつきがあるため、どの味でも厳密に大きさを一定にすることはできません
*3:資料p.5の画像から16個だってすぐにはわからないですよね
*4:正直、手で数えた方が速い。僕の作りが悪いので改良していきます
*5:このモデルは階層ベイズモデルと呼ばれています。データ解析のための統計モデリング入門やPythonによるベイズ統計モデリングが参考になると思います。
*6:変化点検知については、Pythonで体験するベイズ推論に記載があります。
第143回 チョコボール計測(チョコバナナ)
本日の計測結果を報告します。 今回からチョコバナナ味の計測が続きます。
チョコバナナ味の重量傾向を調査したかったので、在庫がなくなる前にAmazonで購入しておきました。
今日は2箱計測します。
計測結果
date | best_before | weight | box_weight | number | factory | shop | angel | net_weight | mean_weight |
---|---|---|---|---|---|---|---|---|---|
2019-03-27 | 2019-11-01 | 27.829 | 4.859 | 15 | 小山工場 | amazon | なし | 22.970 | 1.531 |
2019-03-27 | 2019-11-01 | 27.665 | 4.854 | 15 | 小山工場 | amazon | なし | 22.811 | 1.521 |
今日もエンゼルさん現れずでした。
基礎集計
チョコバナナ味の集計結果です。
項目 | 数 |
---|---|
計測データ数 | 25 |
銀のエンゼル出現数 | 1 |
金のエンゼル出現数 | 0 |
最小 | 中央値 | 最大値 | 平均 | |
---|---|---|---|---|
正味重量 | 22.107 | 22.599 | 24.338 | 22.759 |
個数 | 14.000 | 15.000 | 16.000 | 14.800 |
この図は正味の重量のヒストグラムです。 赤い縦線が仕様(22g)を表しています。 青い太線で正規分布と仮定した最尤推定量をプロットしています。
エンゼル出現確率の予測
通常版のエンゼルの予測を行っていきます。 これまでの通常版パッケージの開封結果は次の通りです。
項目 | 数 |
---|---|
計測データ数 | 399 |
銀のエンゼル出現数 | 18 |
金のエンゼル出現数 | 0 |
この結果を使ってベイズ推定によるエンゼルの出現確率推定を行います。
はじめに銀のエンゼルの出現確率の推定です。
90%信用区間(上下それぞれ5%)は、下側が3.08%上側が6.54%という結果です。
現在の期待値は4.73%です。
次に金のエンゼルの出現確率の推定です。
90%信用区間(上下それぞれ5%)は、下側が0.01%、上側が0.74%という予測になっています。
今のところの金のエンゼルの期待値は0.24%です。
広告
Amazonの欲しいものリスト作ってみました。
チョコボールのカンパ募集中です。
チョコボールをカンパする
第142回 チョコボール計測(ピーナッツ)
本日の計測結果を報告します。 今回も通常版のピーナッツを4箱計測します。
↓これの残りでした(やっと今年のバレンタインが終った)
うちの妻から、バレンタインだってさ pic.twitter.com/E1ClZOuCee
— yoichi_t (@yoichi_t) 2019年2月14日
計測結果
date | best_before | weight | box_weight | number | factory | shop | angel | net_weight | mean_weight |
---|---|---|---|---|---|---|---|---|---|
2019-03-26 | 2019-08-01 | 34.040 | 4.732 | 16 | 小山工場 | スーパー(さいたま市) | なし | 29.308 | 1.832 |
2019-03-26 | 2019-08-01 | 34.223 | 4.707 | 15 | 小山工場 | スーパー(さいたま市) | なし | 29.516 | 1.968 |
2019-03-26 | 2019-08-01 | 35.153 | 4.760 | 15 | 小山工場 | スーパー(さいたま市) | なし | 30.393 | 2.026 |
2019-03-26 | 2019-08-01 | 34.214 | 4.742 | 16 | 小山工場 | スーパー(さいたま市) | なし | 29.472 | 1.842 |
今日もエンゼルさん現れずでした。
基礎集計
通常版のピーナッツ味の集計結果です。
項目 | 数 |
---|---|
計測データ数 | 379 |
銀のエンゼル出現数 | 13 |
金のエンゼル出現数 | 1 |
最小 | 中央値 | 最大値 | 平均 | |
---|---|---|---|---|
正味重量 | 28.594 | 29.382 | 32.232 | 29.476 |
個数 | 14.000 | 17.000 | 20.000 | 16.538 |
この図は正味の重量のヒストグラムです。 赤い縦線が仕様(28g)を表しています。 青い太線で正規分布と仮定した最尤推定量をプロットしています。
エンゼル出現確率の予測
通常版のエンゼルの予測を行っていきます。 これまでの通常版パッケージの開封結果は次の通りです。
項目 | 数 |
---|---|
計測データ数 | 397 |
銀のエンゼル出現数 | 18 |
金のエンゼル出現数 | 0 |
この結果を使ってベイズ推定によるエンゼルの出現確率推定を行います。
はじめに銀のエンゼルの出現確率の推定です。
90%信用区間(上下それぞれ5%)は、下側が3.14%上側が6.74%という結果です。
現在の期待値は4.80%です。
次に金のエンゼルの出現確率の推定です。
90%信用区間(上下それぞれ5%)は、下側が0.01%、上側が0.80%という予測になっています。
今のところの金のエンゼルの期待値は0.26%です。
広告
Amazonの欲しいものリスト作ってみました。
チョコボールのカンパ募集中です。
チョコボールをカンパする
- 出版社/メーカー: 森永製菓
- 発売日: 2016/03/01
- メディア: 食品&飲料
- この商品を含むブログを見る
SIMERST 改良版 携帯タイプはかり ポケットデジタルスケール(秤) 0.001g-100g精密 業務用(プロ用) デジタルスケール 電子天秤
- 出版社/メーカー: Simerst
- メディア:
- この商品を含むブログを見る