Python【Box-Muller法】標準正規分布にしたがう乱数生成 ボックス=ミュラー法(Box-Muller's method)は、一様分布にしたがう確率変数から、標準正規分布にしたがう確率変数を生成する手法です。 コンピュータで標準で生成できる乱数から、正規分布にしたがう乱数(正規乱数)を作るために使用...2020.11.16Python確率・統計
Python指数分布の性質と、最尤推定・パラメータ変換・乱数生成まで 指数分布を定義し、平均や分散といった基本的な性質を示す。その後、最尤推定によりパラメータの値を求める、パラメータを変換してデータをスケーリングする、逆関数法により乱数を生成するという、指数関数の実用テクニックを網羅的に解説する。また、文末にはPython言語による計算例を示し、上記の実用テクニックを一通り体験できるようにした。2020.11.15Python確率・統計
Matplotlibヒストグラムと確率分布の同時プロット―フィッティング精度の検証用に 多数のデータ点が得られたとき のようなフィッティングを行えば、それらのデータ点が従う分布を推定することができる。この場合、データのヒストグラムはその分布に類似した形状をとるはずであり、そのことを視覚的に確かめたい(または、その推定が正しいこ...2020.11.14Matplotlib確率・統計
Python逆関数法とは?その導出と乱数生成のための実装 この記事では、サンプリング法の1つである逆関数法を導出する。逆関数法では、乱数を生成したい累積分布関数の逆関数にもとづき、標準一様乱数を任意の乱数に変換する。PRML(パターン認識と機械学習)の「第11章サンプリング法 11.1.1 標準的な分布」の前半の解説としても有用である。2020.11.13Python数学
機械学習・AIPRML(パターン認識と機械学習)演習問題【難問】解答集 C.M.ビショップ著、PRML(パターン認識と機械学習―ベイズ理論による統計的予測、)の演習問題において、【難問】に指定されている問題の解答集。解答を当サイトに投稿し次第、またはどこかで見つけ次第、リンクを増やしていく。2020.11.12機械学習・AI確率・統計
機械学習・AI【PRML10章 演習10.39解答】EP法による事後分布の更新 PRML10章 演習10.39の解答・解説。EP法によって更新した事後分布の平均と分散を導出する。それに際し、パラメータの期待値を計算する。その後近似因子のパラメータを求め、モデルエビデンスが近似される値を導く。2020.11.12機械学習・AI確率・統計
機械学習・AIPRML10章 演習10.39解答(その4:(10.223),(10.224)の導出) 問題、記号の意味、利用した式については を参照。 目標 モデルエビデンスが $$p(\mathcal{D})\simeq\left(\frac{\nu^{\mathrm {new}}}{b}\right)^{D/2}\left\exp{\l...2020.11.12機械学習・AI確率・統計
機械学習・AIPRML10章 演習10.39解答(その3:(10.220)-(10.222)の導出) 問題、記号の意味、利用した式については を参照。 目標 改良された近似因子 \(\tilde{f}_{n}(\boldsymbol{\theta})\) のパラメータが $$\nu_{n}^{-1}=(\nu^{\mathrm{new}})...2020.11.12機械学習・AI確率・統計
機械学習・AIPRML10章 演習10.39解答(その2:(10.218)の導出) 問題、記号の意味、利用した式については を参照。 目標 \(q^{\mathrm{new}}(\boldsymbol{\theta})\) の平均と分散を \(q^{\backslash n}(\boldsymbol{\theta})f_{...2020.11.12機械学習・AI確率・統計
機械学習・AIPRML10章 演習10.39解答(その1:(10.217),(10.219)の導出) 問題、記号の意味、利用した式については を参照。 目標 \(q^{\mathrm{new}}(\boldsymbol{\theta})\) の平均と分散を \(q^{\backslash n}(\boldsymbol{\theta})f_{...2020.11.11機械学習・AI確率・統計