機械学習基礎理論独習

誤りがあればご指摘いただけると幸いです。数式が整うまで少し時間かかります。リンクフリーです。

勉強ログです。リンクフリーです
目次へ戻る

最尤推定 - 正規分布

尤度

同一の正規分布から標本{\bf x}=\{x_1,\ldots,x_n\}が互いに独立に観察されたとします(独立同分布)。

この時の尤度関数は

\begin{eqnarray}
f({\bf x}|\mu,\sigma^2)&=&\prod_{i=1}^nf(x_i|\mu,\sigma^2)\\
&=&\prod_{i=1}^n\frac{1}{\sqrt{2\pi}\sigma}\exp\left(-\frac{1}{2\sigma^2}(x_i-\mu)^2\right)\tag{1}\\
\end{eqnarray}

となります。
(2)に対数をとります。

\begin{eqnarray}
\ln f({\bf x}|\mu,\sigma^2) &=& \ln \prod_{i=1}^n\frac{1}{\sqrt{2\pi}\sigma}\exp\left(-\frac{1}{2\sigma^2}(x_i-\mu)^2\right)\\
&=&\ln\left(\frac{1}{\sqrt{2\pi}}\right)+\ln\left(\frac{1}{\sigma}\right)+\sum_{i=1}^n-\frac{1}{2\sigma^2}(x_i-\mu)^2\\
&=&-\frac{n}{2}\ln2\pi-\frac{n}{2}\ln\sigma^2-\frac{1}{2\sigma^2}\sum_{i=1}^n(x_i-\mu)^2\tag{2}
\end{eqnarray}

最尤推定

(2)\mu微分して=0とおきます。

\begin{eqnarray}
&&\frac{\rm d}{{\rm d}{\mu}}\ln f({\bf x}|\mu,\sigma^2)=0\\
&&\Leftrightarrow \frac{\rm d}{{\rm d}{\mu}}\left(-\frac{1}{2\sigma^2}\sum_{i=1}^n(x_i-\mu)^2\right)=0\\
&&\Leftrightarrow -\frac{1}{2\sigma^2}\frac{\rm d}{{\rm d}{\mu}}\sum_{i=1}^n(x_i-\mu)^2=0\\
&&\Leftrightarrow \sum_{i=1}^n(x_i-\mu)=0\\
&&\Leftrightarrow \mu=\frac{1}{2}\sum_{i=1}^nx_i\tag{3}
\end{eqnarray}

(3)は標本平均と一致します。

(2)\sigma^2微分して=0とおきます。

\begin{eqnarray}
&&\frac{\rm d}{{\rm d}{\sigma^2}}\ln f({\bf x}|\mu,\sigma^2)=0\\
&&\Leftrightarrow \frac{\rm d}{{\rm d}{\sigma^2}}\left(-\frac{n}{2}\ln\sigma^2-\frac{1}{2\sigma^2}\sum_{i=1}^n(x_i-\mu)^2\right)=0\\
&&\Leftrightarrow \frac{\rm d}{{\rm d}{\sigma^2}}\left(-\frac{n}{2}\ln\sigma^2\right)-\frac{\rm d}{{\rm d}{\sigma^2}}\left(\frac{1}{2\sigma^2}\sum_{i=1}^n(x_i-\mu)^2\right)=0\\
&&\Leftrightarrow -\frac{n}{2\sigma^2}+\frac{1}{2\sigma^4}\sum_{i=1}^2(x_i-\mu)^2=0\\
&&\Leftrightarrow \sigma^2=\frac{1}{n}\sum_{i=1}^n(x_i-\mu)^2\tag{4}
\end{eqnarray}

(4)は標本分散と一致します。

対数尤度関数のグラフを示します。左が3次元のグラフで、右が等高線のグラフです。
{\bf x}=\{162,172,178,154,173,174,166,166,166,164\}とします。
標本平均は167.5、標本分散は43.45です。

f:id:olj611:20210415224447p:plain:w800

偉人の名言

f:id:olj611:20210416000457p:plain:w300
君にできる最大の冒険は、君の夢に生きる事だ。
オプラ・ウィンフリー

参考文献

基礎からのベイズ統計学

動画

目次へ戻る