多変量正規分布の条件つき分布① ~ブロック行列の逆行列と行列式~
ベイズ統計学を勉強していると、多変量正規分布の条件つき分布の使用は避けては通れません。 つまり、$\boldsymbol{x}=(\boldsymbol{x_1}, \boldsymbol{x_2})$が多変量正規分布$ … 続きを読む
Pyroで非負値行列分解を実装してみる
前回の記事(『Numpyroで主成分分析を実装してみる』)では、確率的プログラミング言語(PPL)の一種であるNumpyroを使って、MCMCによる確率的主成分分析(PPCA)の実装方法を紹介しました。
今回は、PPCAと同じ次元削減の手法のひとつである非負値行列分解(NMF)を、PPLの一種であり、Numpyroの親戚であるPyroを使って実装してみます。
Numpyroで主成分分析を実装してみる
みなさんは普段、どのようにベイズモデリングを実装されていますか？ アカデミックの方や理論・プログラミングに明るい人は「自分でMHアルゴリズムなどMCMCを一から実装している」という選択肢をとられているかもしれません。 が … 続きを読む
理系でもビジネストランスレーターになれるのか？
ビジネストランスレーターとは？ みなさん「ビジネストランスレーター」という言葉を聞いたことありますか？私は恥ずかしながら最近知ったのですが、その定義を一言で表すと『ビジネスとデータサイエンスの橋渡し』を担うポジションのこ … 続きを読む
統計調査士の合格体験記
先日、統計調査士試験を受験し、無事合格することができました!! そこで忘れないうちに、私がどのような勉強を行ってきたかをこのブログに記録しておこうと思います。 統計調査士とは？ 統計検定の一種であり、国勢調査や家計調査な … 続きを読む
DeepARで複数系列の予測をする
業務で時系列データの予測モデルを作る時、私は基本的にProphetを最初に選択します。なかなかの精度のモデルが、手早く簡単に作れるからです。 しかし、時系列データが数百系列あるような場合、いかにProphetといえどもモ … 続きを読む
AWSクラウドプラクティショナーを受験してきた
数日前に受験してきまして、何とか合格できたみたいなので、忘れないうちにブログにいろいろ書き残しておこうと思います。(ホントはこのブログでは技術関連の投稿だけやるつもりだったんですが、それだけだと中々ヘビーなので…) AW … 続きを読む
Pyroで混合ベルヌーイによるクラスタリングをやってみる
混合分布モデルでは、離散値の場合は混合ポアソンモデル、連続値の場合は混合ガウスモデルが使われることが多いですが、この2つ以外にも適当な確率分布を使った「混合○○モデル」を定義することはできます。 今回は、データが0 or … 続きを読む
HDIの解説とRでの計算
HDIとは? ベイズ統計では、MCMCで発生させたサンプルを用いて、各変数の取り得る値の区間を計算することがあります。 例を以下のシンプルな単回帰で考えてみます。 \begin{gather}y_i = \alpha + … 続きを読む