第132回 チョコボール計測(ホワイトミルクピーナッツ)
本日の計測結果を報告します。 今日は新味の「ホワイトミルクピーナッツ」味を計測します。
計測結果
date | best_before | weight | box_weight | number | factory | shop | angel | net_weight | mean_weight |
---|---|---|---|---|---|---|---|---|---|
2019-02-25 | 2019-09-01 | 33.504 | 4.836 | 16 | 小山工場 | スーパー(さいたま市) | なし | 28.668 | 1.792 |
2019-02-25 | 2019-09-01 | 33.684 | 4.832 | 17 | 小山工場 | スーパー(さいたま市) | なし | 28.852 | 1.697 |
2019-02-25 | 2019-09-01 | 35.001 | 4.879 | 17 | 小山工場 | スーパー(さいたま市) | なし | 30.122 | 1.772 |
今日はエンゼルさん出ませんでした。さいきん見てない気がしますね。
基礎集計
ホワイトミルクピーナッツ味の集計結果です。
項目 | 数 |
---|---|
計測データ数 | 7 |
銀のエンゼル出現数 | 0 |
金のエンゼル出現数 | 0 |
最小 | 中央値 | 最大値 | 平均 | |
---|---|---|---|---|
正味重量 | 28.337 | 28.798 | 30.122 | 28.926 |
個数 | 15.000 | 16.000 | 17.000 | 16.000 |
この図は正味の重量のヒストグラムです。 赤い縦線が仕様(28g)を表しています。 青い太線で正規分布と仮定した最尤推定量をプロットしています。
エンゼル出現確率の予測
通常版のエンゼルの予測を行っていきます。 これまでの通常版パッケージの開封結果は次の通りです。
項目 | 数 |
---|---|
計測データ数 | 362 |
銀のエンゼル出現数 | 17 |
金のエンゼル出現数 | 0 |
この結果を使ってベイズ推定によるエンゼルの出現確率推定を行います。
はじめに銀のエンゼルの出現確率の推定です。
90%信用区間(上下それぞれ5%)は、下側が3.19%上側が7.07%という結果です。
現在の期待値は4.93%です。
次に金のエンゼルの出現確率の推定です。
90%信用区間(上下それぞれ5%)は、下側が0.01%、上側が0.78%という予測になっています。
今のところの金のエンゼルの期待値は0.26%です。
広告
Amazonの欲しいものリスト作ってみました。
チョコボールのカンパ募集中です。
チョコボールをカンパする
- 出版社/メーカー: 森永製菓
- 発売日: 2016/03/01
- メディア: 食品&飲料
- この商品を含むブログを見る
SIMERST 改良版 携帯タイプはかり ポケットデジタルスケール(秤) 0.001g-100g精密 業務用(プロ用) デジタルスケール 電子天秤
- 出版社/メーカー: Simerst
- メディア:
- この商品を含むブログを見る
第131回 チョコボール計測(ピーナッツ)
本日の計測結果を報告します。 今日は「ピーナッツ」味を計測します。
計測結果
date | best_before | weight | box_weight | number | factory | shop | angel | net_weight | mean_weight |
---|---|---|---|---|---|---|---|---|---|
2019-02-21 | 2019-08-01 | 34.579 | 4.697 | 16 | 小山工場 | コンビニ(港区) | なし | 29.882 | 1.868 |
2019-02-21 | 2019-08-01 | 33.979 | 4.699 | 17 | 小山工場 | コンビニ(港区) | なし | 29.280 | 1.722 |
今日はエンゼルさん出ませんでした。
基礎集計
ピーナッツ味の集計結果です。
項目 | 数 |
---|---|
計測データ数 | 351 |
銀のエンゼル出現数 | 12 |
金のエンゼル出現数 | 1 |
最小 | 中央値 | 最大値 | 平均 | |
---|---|---|---|---|
正味重量 | 28.594 | 29.387 | 32.232 | 29.481 |
個数 | 14.000 | 17.000 | 20.000 | 16.624 |
この図は正味の重量のヒストグラムです。 赤い縦線が仕様(28g)を表しています。 青い太線で正規分布と仮定した最尤推定量をプロットしています。
エンゼル出現確率の予測
通常版のエンゼルの予測を行っていきます。 これまでの通常版パッケージの開封結果は次の通りです。
項目 | 数 |
---|---|
計測データ数 | 359 |
銀のエンゼル出現数 | 17 |
金のエンゼル出現数 | 0 |
この結果を使ってベイズ推定によるエンゼルの出現確率推定を行います。
はじめに銀のエンゼルの出現確率の推定です。
90%信用区間(上下それぞれ5%)は、下側が3.11%上側が6.99%という結果です。
現在の期待値は4.91%です。
次に金のエンゼルの出現確率の推定です。
90%信用区間(上下それぞれ5%)は、下側が0.01%、上側が0.92%という予測になっています。
今のところの金のエンゼルの期待値は0.29%です。
広告
Amazonの欲しいものリスト作ってみました。
チョコボールのカンパ募集中です。
チョコボールをカンパする
第130回 チョコボール計測(ピーナッツ、ホワイトミルク)
本日の計測結果を報告します。 今日は新味の「ホワイトミルク ピーナッツ」味です。
通常のピーナッツと同様ですが、コーティングされているチョコがホワイトミルク味というもののようです。 仕様上の内容量は通常のピーナッツ味と同じで28gとのことです。
中身はこんなかんじです。 通常のチョコボールの色が白くなったものですね。 ただ、今回が特別なのかわかりませんが、形が均一で綺麗な気がします。
計測結果
date | best_before | weight | box_weight | number | factory | shop | angel | net_weight | mean_weight |
---|---|---|---|---|---|---|---|---|---|
2019-02-19 | 2019-09-01 | 33.498 | 4.845 | 16 | 小山工場 | コンビニ(さいたま市) | なし | 28.653 | 1.791 |
2019-02-19 | 2019-09-01 | 33.889 | 4.837 | 16 | 小山工場 | コンビニ(さいたま市) | なし | 29.052 | 1.816 |
2019-02-19 | 2019-09-01 | 33.646 | 4.848 | 15 | 小山工場 | コンビニ(さいたま市) | なし | 28.798 | 1.920 |
2019-02-19 | 2019-09-01 | 33.202 | 4.865 | 15 | 小山工場 | コンビニ(さいたま市) | なし | 28.337 | 1.889 |
今日はエンゼルさん出ませんでした。
基礎集計
ホワイトミルクピーナッツ味の集計結果です。
項目 | 数 |
---|---|
計測データ数 | 4 |
銀のエンゼル出現数 | 0 |
金のエンゼル出現数 | 0 |
最小 | 中央値 | 最大値 | 平均 | |
---|---|---|---|---|
正味重量 | 28.337 | 28.726 | 29.052 | 28.710 |
個数 | 15.000 | 15.500 | 16.000 | 15.500 |
この図は正味の重量のヒストグラムです。 赤い縦線が仕様(28g)を表しています。 青い太線で正規分布と仮定した最尤推定量をプロットしています。
まだ4箱なので、傾向ははっきりとは見えてきませんが、通常のピーナッツよりも仕様に対してのマージンが少ないような気がしますね。 通常ピーナッツとホワイトミルクピーナッツの重量分布をプロットしてみます。 なお、仕様上の重量がどちらも同じなので、右図と左図は同じになります。
今後の調査に期待です。
エンゼル出現確率の予測
通常版のエンゼルの予測を行っていきます。 これまでの通常版パッケージの開封結果は次の通りです。
項目 | 数 |
---|---|
計測データ数 | 357 |
銀のエンゼル出現数 | 17 |
金のエンゼル出現数 | 0 |
この結果を使ってベイズ推定によるエンゼルの出現確率推定を行います。
はじめに銀のエンゼルの出現確率の推定です。
90%信用区間(上下それぞれ5%)は、下側が3.25%上側が6.99%という結果です。
現在の期待値は5.02%です。
次に金のエンゼルの出現確率の推定です。
90%信用区間(上下それぞれ5%)は、下側が0.02%、上側が0.80%という予測になっています。
今のところの金のエンゼルの期待値は0.28%です。
広告
Amazonの欲しいものリスト作ってみました。
チョコボールのカンパ募集中です。
チョコボールをカンパする
第129回 チョコボール計測(チョコバナナ)
本日の計測結果を報告します。 今日はチョコバナナ味を3箱計測します。
計測結果
date | best_before | weight | box_weight | number | factory | shop | angel | net_weight | mean_weight |
---|---|---|---|---|---|---|---|---|---|
2019-02-18 | 2019-10-01 | 28.095 | 4.824 | 15 | 小山工場 | コンビニ(さいたま市) | なし | 23.271 | 1.551 |
2019-02-18 | 2019-10-01 | 27.001 | 4.831 | 14 | 小山工場 | コンビニ(さいたま市) | なし | 22.170 | 1.584 |
2019-02-18 | 2019-10-01 | 29.159 | 4.821 | 16 | 小山工場 | コンビニ(さいたま市) | なし | 24.338 | 1.521 |
今日はエンゼルさん出ませんでした。
基礎集計
これまでのチョコバナナ味の集計結果です。
項目 | 数 |
---|---|
計測データ数 | 23 |
銀のエンゼル出現数 | 1 |
金のエンゼル出現数 | 0 |
最小 | 中央値 | 最大値 | 平均 | |
---|---|---|---|---|
正味重量 | 22.107 | 22.531 | 24.338 | 22.748 |
個数 | 14.000 | 15.000 | 16.000 | 14.783 |
この図は正味の重量のヒストグラムです。 赤い縦線が仕様(22g)を表しています。 青い太線で正規分布と仮定した最尤推定量をプロットしています。
仕様に近い量が多いですが、2g以上多いケースも散見されますね。
工業製品なのにこんなにばらつくんですね。
時期によって傾向が変わったりするかも知れないので、この辺りは追加調査をご期待下さい。
エンゼル出現確率の予測
通常版のエンゼルの予測を行っていきます。 これまでの通常版パッケージの開封結果は次の通りです。
項目 | 数 |
---|---|
計測データ数 | 353 |
銀のエンゼル出現数 | 17 |
金のエンゼル出現数 | 0 |
この結果を使ってベイズ推定によるエンゼルの出現確率推定を行います。
はじめに銀のエンゼルの出現確率の推定です。
90%信用区間(上下それぞれ5%)は、下側が3.27%上側が7.37%という結果です。
現在の期待値は5.15%です。
次に金のエンゼルの出現確率の推定です。
90%信用区間(上下それぞれ5%)は、下側が0.01%、上側が0.80%という予測になっています。
今のところの金のエンゼルの期待値は0.25%です。
広告
Amazonの欲しいものリスト作ってみました。
チョコボールのカンパ募集中です。
チョコボールをカンパする
第128回 チョコボール計測(いちご)
本日の計測結果を報告します。 今日はいちご味を3箱計測します。
計測結果
date | best_before | weight | box_weight | number | factory | shop | angel | net_weight | mean_weight |
---|---|---|---|---|---|---|---|---|---|
2019-02-14 | 2019-10-01 | 30.877 | 4.818 | 15 | 小山工場 | スーパー(さいたま市) | なし | 26.059 | 1.737 |
2019-02-14 | 2019-10-01 | 30.869 | 4.819 | 15 | 小山工場 | スーパー(さいたま市) | なし | 26.050 | 1.737 |
2019-02-14 | 2019-10-01 | 30.932 | 4.820 | 15 | 小山工場 | スーパー(さいたま市) | なし | 26.112 | 1.741 |
今回はエンゼルさんは出ませんでした。
基礎集計
これまでのいちご味の集計結果です。
項目 | 数 |
---|---|
計測データ数 | 37 |
銀のエンゼル出現数 | 1 |
金のエンゼル出現数 | 0 |
最小 | 中央値 | 最大値 | 平均 | |
---|---|---|---|---|
正味重量 | 25.862 | 26.624 | 27.654 | 26.671 |
個数 | 15.000 | 15.000 | 16.000 | 15.405 |
この図は正味の重量のヒストグラムです。 赤い縦線が仕様(25g)を表しています。 青い太線で正規分布と仮定した最尤推定量をプロットしています。
エンゼル出現確率の予測
通常版のエンゼルの予測を行っていきます。 これまでの通常版パッケージの開封結果は次の通りです。
項目 | 数 |
---|---|
計測データ数 | 350 |
銀のエンゼル出現数 | 17 |
金のエンゼル出現数 | 0 |
この結果を使ってベイズ推定によるエンゼルの出現確率推定を行います。
はじめに銀のエンゼルの出現確率の推定です。
90%信用区間(上下それぞれ5%)は、下側が3.27%上側が7.31%という結果です。
現在の期待値は5.14%です。
次に金のエンゼルの出現確率の推定です。
90%信用区間(上下それぞれ5%)は、下側が0.02%、上側が0.90%という予測になっています。
今のところの金のエンゼルの期待値は0.29%です。
広告
Amazonの欲しいものリスト作ってみました。
チョコボールのカンパ募集中です。
チョコボールをカンパする
第127回 チョコボール計測(ピーナツ)
本日の計測結果を報告します。 今日はピーナツ味を2箱計測します。
計測結果
date | best_before | weight | box_weight | number | factory | shop | angel | net_weight | mean_weight |
---|---|---|---|---|---|---|---|---|---|
2019-02-13 | 2019-08-01 | 34.065 | 4.777 | 14 | 小山工場 | コンビニ(さいたま市) | なし | 29.288 | 2.092 |
2019-02-13 | 2019-08-01 | 34.937 | 4.755 | 16 | 小山工場 | コンビニ(さいたま市) | なし | 30.182 | 1.886 |
今回はエンゼルさんは出ませんでした。
基礎集計
これまでのピーナツ味の集計結果です。
項目 | 数 |
---|---|
計測データ数 | 349 |
銀のエンゼル出現数 | 12 |
金のエンゼル出現数 | 1 |
最小 | 中央値 | 最大値 | 平均 | |
---|---|---|---|---|
正味重量 | 28.594 | 29.387 | 32.232 | 29.481 |
個数 | 14.000 | 17.000 | 20.000 | 16.625 |
この図は正味の重量のヒストグラムです。 赤い縦線が仕様(28g)を表しています。 青い太線で正規分布と仮定した最尤推定量をプロットしています。
エンゼル出現確率の予測
通常版のエンゼルの予測を行っていきます。 これまでの通常版パッケージの開封結果は次の通りです。
項目 | 数 |
---|---|
計測データ数 | 347 |
銀のエンゼル出現数 | 17 |
金のエンゼル出現数 | 0 |
この結果を使ってベイズ推定によるエンゼルの出現確率推定を行います。
はじめに銀のエンゼルの出現確率の推定です。
90%信用区間(上下それぞれ5%)は、下側が3.36%上側が7.06%という結果です。
現在の期待値は5.17%です。
次に金のエンゼルの出現確率の推定です。
90%信用区間(上下それぞれ5%)は、下側が0.01%、上側が0.83%という予測になっています。
今のところの金のエンゼルの期待値は0.27%です。
広告
Amazonの欲しいものリスト作ってみました。
チョコボールのカンパ募集中です。
チョコボールをカンパする
- 出版社/メーカー: 森永製菓
- 発売日: 2016/03/01
- メディア: 食品&飲料
- この商品を含むブログを見る
SIMERST 改良版 携帯タイプはかり ポケットデジタルスケール(秤) 0.001g-100g精密 業務用(プロ用) デジタルスケール 電子天秤
- 出版社/メーカー: Simerst
- メディア:
- この商品を含むブログを見る
チョコボールに学ぶ実践的ベイズ統計モデリング入門 #1
【概要】
【目次】
はじめに
2019/02/09に下記のセミナーを実施しました。 私は、2時間枠でのセミナーというのは初体験だったもので、 最後ちょっと駆け足になってしまったかなと反省しております。
この記事では、セミナーのフォローアップを目的として、 資料や資料の中で言い足りなかったことなどを補足します。
資料
当日使用した資料は下記Slideshareにアップしています。 なお、公開用に一部修正を加えています。
ハンズオンで使用したコードは下記のGithubリポジトリに上がっております。 MITライセンスの範囲で自由に利用していただいて構いません。
統計モデリングの目的
統計モデリングとは何か?
「統計モデリング」とは?という疑問に対して、明確な定義は無いように見受けます。 Googleで検索するといくつかヒットしますが、 私自身は「統計学の手法を使ってシステムを理解するための枠組み」と捉えています。
「物事を理解する」というのには色々やり方があると思いますが、 特に「活用*1」をするために、 統計学の手法(確率分布)を利用するのだと考えています。 純粋に楽しいというのはあるにしろ、 活用が出来なければビジネスとしては利用価値が無いと考えます。 そのため、過度に複雑なモデルを作ることにはあまり意味がなく、 できる限りシンプルなモデルの方が利用効率が高いと思っています。
具体的な活用の仕方という点については、当然ながら問題によってケースバイケースです。 モデリングの結果から、異常値の定義のようなものを直接導いても良いですし、 ニューラルネットワークを使うという判断をする場合もあると思います。
「AI」と統計モデリング
上記では、統計モデリングの結果は分類や予測などの問題に活用すると述べてきました。 このような問題は「AI」を使って解決できると言われていますが、 「AI」と統計モデリングの関係について考えてみます。
結論から言うと、「AI」と言う技術は存在せず、統計モデリングだって「AI」系の技術の一つです。
データの構造を推定すると言う点で、統計モデリングはこのような問題を考える上で最も基礎となる考え方です。
ニューラルネットワークなど注目されている技術は、
データ構造を明示的に推定する部分を飛ばして、
前節で言うところの「活用」を直接行なっているものと考えています。
この部分については、私がやっている別のブログの記事でも言及しておりますので、
よければこちらも読んでいただけたらと思います。
セミナー資料で識別境界の例を載せていますが、 こちらはirisデータセットのうち2種(versicolor, virginica)を分類するために NN*2を使って学習した境界線です。 この図の計算も上記のGithubリポジトリに上がっていますので、 よければ参照ください*3。
ここまでの参考資料等
ここまでを理解するための資料として、 統計学入門が統計学についてかっちりと書かれているのでおすすめです。 これがちょっと難しい場合には、道具としてのベイズ統計など易しめの本を一読すると理解が進むと思います。
モデリングや活用についての具体的なイメージについては、その問題、数理モデルが解決しますがとても イメージしやすくて良いと思いました。
ハンズオンでも利用する「ベイズ統計」について、 成り立ちや利用シーンについて、異端の統計学 ベイズを一読されることをおすすめします。 この書籍は、理論書ではないので数理的な話はありませんが、活用されてきた場面を知ることで、 どのような問題に適しているのかがイメージできると思います。
最尤推定ハンズオン
今回のセミナーでは、モデリングの具体的な問題としてチョコボールの銀のエンゼルの出現率(含有率)を推定する問題を扱います。
問題設定
銀のエンゼルの出現はベルヌーイ試行*4であると仮定しています。 ベルヌーイ試行を仮定することで、エンゼルの出現を二項分布でモデル化できます。 当たり外れのような2値の問題でも、ベルヌーイ試行を仮定できない場合はあると思いますが、 そのような場合にはモデルを工夫する必要があります。 操作できるならデータは取得段階でランダム化した方が良いですし、どこかでランダム性を仮定することが必要になります。
最尤推定とは
最尤推定とは、推定対象のパラメータを変数に持った「尤度関数」を定義し、 尤度関数を最大化するを求める方法です。 尤度関数とは、モデル(二項分布などの確率分布)に具体的なデータがどのくらい当てはまるかを定める関数です。
今回は単純な二項分布をモデルにしているので、最尤推定量は解析的に計算ができて、標本平均となります。 モデルが複雑で最尤推定量を解析的に計算できない場合には、 勾配法などの探索的なアルゴリズムを利用してパラメータを推定します。
ここまでの参考資料等
ここまでを理解するための資料としては、前章で述べた以下の2冊が参考になります。
また、Googleなどで検索してもたくさん出てきます。
ベイズ推定ハンズオン
前章で、銀のエンゼルの出現は二項分布でモデル化することにしました。 ここでは、二項分布のパラメータ(エンゼルの出現確率)をベイズ推定することを試してみます。
ベイズ推定とは
最尤推定では、推定値(エンゼルの出現確率)は一つの値として得られました(結局、標本平均でした)。 一方ベイズ推定では、推定結果は事後分布と言う関数(分布)で得られます。
推定結果が分布になると理解が難しいと感じるかもしれませんが、 データが10個しかない場合と1,000個ある場合で、 推定結果の信頼性はどっちが上かと問われたら当然後者であると理解できるかと思います *5。 結果が分布で得られるとは、結果の信頼性がどのくらいあるのかを併せて推定していることになります。 こうしてみると、分布で結果が得られるとは、人間の直感に適合するように感じられるかと思います。
ベイズ推定の具体的な計算
ベイズ推定の計算は、最尤推定のように単純に平均を出すようには行きません。 解析的に解く方法と数値計算的に解く方法の二つがあります。
解析的に解くためには、尤度と事前分布の組み合わせを工夫する必要があります。 エンゼルの含有率推定では、二項分布をモデルに選択したので、ベータ分布を事前分布にすることで、 事後分布はベータ分布になります(当然、ベータ分布のパラメータはデータによって事前分布のパラメータとは異なる)。 このような事前分布を共役事前分布と呼んでおり、尤度として選ぶ関数(確率分布)毎に共役事前分布が決まっています。
尤度と共役事前分布のペアを利用すれば解析的に事後分布が計算できるのですが、 事前分布を複雑にしたり*6、共役ではない事前分布を選びたい場合があると思います。 そのような場合に、MCMC(マルコフチェーンモンテカルロ)と呼ばれる手法を利用することで、数値計算的に事後分布を推定できます。
事後分布の活用
MCMCの結果は、事後分布からサンプルしたデータの集合になります。 そのため、サンプル系列のヒストグラムを見れば事後分布の形が想像できますし、 平均は事後分布の期待値になります。 信用区間を計算するのも、2.5%点と97.5%点(パーセンタイル)を求めれば95%信用区間を得ることができます。
統計モデリングの「活用」としては、期待値や最頻値を使って点推定と同様の活用ができます。 分布の広がりを考慮して信用区間の幅を活用することもできます。 エンゼルの含有率では、90%信用区間は3%〜8%と推定されました。
この結果を利用して、5個のエンゼルを得るまでに必要な外れチョコボールの個数を予測したところ以下となりました。
この結果から、 エンゼル5個を50%の確率で得るまでに平均的には77個買えば良いが、 リスクを考慮すると130個程度買えば良いであろうと言うことが読み取れます。
ここまでの参考資料等
ハンズオンではPythonを利用し、Pythonの確率プログラミング言語であるpymc3を利用しました。 pymc3を利用した参考書としては、Pythonによるベイズ統計モデリングが参考になるかと思います。 こちらは理論的な解説も載っているので、まずこちらの書籍を利用して触りながらベイズ統計を理解すると良いかもしれません。
ベイズ統計の実践的参考書としてはPythonで体験するベイズ推論もあり、 こちらは事例が豊富で参考になるシーンは多いかと思います。 しかし、こちらの書籍ではpymc2を利用しているため、コードの記述方法を多少読み替える必要があります。
終わりに
実は、当初はモデリング対象を2つ用意していたのですが、 「どう考えても2時間では収まらない」という助言をもらったので、 今回はエンゼルの含有率の推定のみを扱いました。 実際、2時間のセミナーではこのボリュームでちょうど良かったです。
せっかく作ったし、また、二項分布のパラメータ推定の問題だけでは現実問題をカバーしきれないので、 第2回も開催します。 第1回に来れなかった方も是非参加いただけたらと思います。
参考文献
統計学入門
有名な統計学の入門書。かっちり書いてあるので、統計学を学ぶにはすごく良い。- 作者: 東京大学教養学部統計学教室
- 出版社/メーカー: 東京大学出版会
- 発売日: 1991/07/09
- メディア: 単行本
- 購入: 158人 クリック: 3,604回
- この商品を含むブログ (79件) を見る
道具としてのベイズ統計
とりあえず使ってみてから理論的なことを勉強しようという場合にはちょうど良いかも。 初学者向け。- 作者: 涌井良幸
- 出版社/メーカー: 日本実業出版社
- 発売日: 2009/11/19
- メディア: 単行本
- 購入: 10人 クリック: 71回
- この商品を含むブログ (18件) を見る
その問題、数理モデルが解決します
モデリングの考え方が物語として(事例ベースで)書かれている。 とても理解しやすいので、モデリングとは何かを知りたい場合におすすめ。- 作者: 浜田宏
- 出版社/メーカー: ベレ出版
- 発売日: 2018/12/12
- メディア: 単行本
- この商品を含むブログを見る
- どのような場面でベイズ統計が使われてきたのか、歴史的な背景などを記した読み物。 ベイズを扱うか扱わないかに関わらず、是非一読してほしいオススメの書籍。 読めばベイズ推定を利用したくなると思います。
- 作者: シャロン・バーチュマグレイン,Sharon Bertsch McGrayne,冨永星
- 出版社/メーカー: 草思社
- 発売日: 2013/10/23
- メディア: 単行本
- この商品を含むブログ (28件) を見る
- pymc3を使った入門書。 理論面も多少書かれており、こちらで手を動かしながら学習するのが良いと思う。
Pythonによるベイズ統計モデリング: PyMCでのデータ分析実践ガイド
- 作者: オズワルドマーティン,Osvaldo Martin,金子武久
- 出版社/メーカー: 共立出版
- 発売日: 2018/06/22
- メディア: 単行本
- この商品を含むブログを見る
- pymcを利用した事例ベースの参考書。 事例が豊富で参考になるシーンは出てくるかもしれないが、階層モデルについては記載がない。 また、ptmc2なので、pymc3と記法が少し違うので注意。
Pythonで体験するベイズ推論 PyMCによるMCMC入門
- 作者: キャメロンデビッドソン=ピロン,玉木徹
- 出版社/メーカー: 森北出版
- 発売日: 2017/04/06
- メディア: 単行本(ソフトカバー)
- この商品を含むブログ (1件) を見る
広告
Amazonの欲しいものリスト作ってみました。
チョコボールのカンパ募集中です。
チョコボールをカンパする
- 出版社/メーカー: 森永製菓
- 発売日: 2016/03/01
- メディア: 食品&飲料
- この商品を含むブログを見る
SIMERST 改良版 携帯タイプはかり ポケットデジタルスケール(秤) 0.001g-100g精密 業務用(プロ用) デジタルスケール 電子天秤
- 出版社/メーカー: Simerst
- メディア:
- この商品を含むブログを見る