Pyroで非負値行列分解を実装してみる
前回の記事(『Numpyroで主成分分析を実装してみる』)では、確率的プログラミング言語(PPL)の一種であるNumpyroを使って、MCMCによる確率的主成分分析(PPCA)の実装方法を紹介しました。
今回は、PPCAと同じ次元削減の手法のひとつである非負値行列分解(NMF)を、PPLの一種であり、Numpyroの親戚であるPyroを使って実装してみます。
前回の記事(『Numpyroで主成分分析を実装してみる』)では、確率的プログラミング言語(PPL)の一種であるNumpyroを使って、MCMCによる確率的主成分分析(PPCA)の実装方法を紹介しました。
今回は、PPCAと同じ次元削減の手法のひとつである非負値行列分解(NMF)を、PPLの一種であり、Numpyroの親戚であるPyroを使って実装してみます。
みなさんは普段、どのようにベイズモデリングを実装されていますか? アカデミックの方や理論・プログラミングに明るい人は「自分でMHアルゴリズムなどMCMCを一から実装している」という選択肢をとられているかもしれません。 が … 続きを読む
先日、統計調査士試験を受験し、無事合格することができました!! そこで忘れないうちに、私がどのような勉強を行ってきたかをこのブログに記録しておこうと思います。 統計調査士とは? 統計検定の一種であり、国勢調査や家計調査な … 続きを読む
混合分布モデルでは、離散値の場合は混合ポアソンモデル、連続値の場合は混合ガウスモデルが使われることが多いですが、この2つ以外にも適当な確率分布を使った「混合○○モデル」を定義することはできます。 今回は、データが0 or … 続きを読む
HDIとは? ベイズ統計では、MCMCで発生させたサンプルを用いて、各変数の取り得る値の区間を計算することがあります。 例を以下のシンプルな単回帰で考えてみます。 \begin{gather}y_i = \alpha + … 続きを読む