確率について

（ ソース ）

機械学習、特にニューラルネットワークの仕組みを他の人に伝えなければならないこともあります。 通常、勾配降下と線形回帰から始め、徐々に多層パーセプトロン、自動エンコーダー、畳み込みネットワークに移行します。 全員が頭をうなずきながらうなずきますが、ある時点で、抜け目のない人は必ず尋ねます：

そして、なぜ線形回帰の変数が独立していることがそれほど重要なのでしょうか？

または

また、通常の完全に接続されたネットワークではなく、画像に畳み込みネットワークが使用されるのはなぜですか？

「ああ、簡単です」と答えたいです。 -「変数が依存している場合、変数間の条件付き確率分布をモデル化する必要があるため」または「小さなローカルエリアでピクセルの共同分布を学習する方がはるかに簡単であるため」 しかし、ここに問題があります。生徒はまだ確率分布と確率変数について何も知らないので、他の方法で出て、より複雑に説明しますが、概念と用語は少なくなります。 また、バッチに関する正規化モデルまたは生成モデルについて説明するように求められた場合はどうすればよいかわかりません。

ですから、自分自身や他人を苦しめるのではなく、確率論の基本概念を覚えておいてください。

ランダム変数

年齢、身長、性別、子供の数が示されている人々のプロファイルがあると想像してください。

このようなテーブルの各行はオブジェクトです。 各セルは、このオブジェクトを特徴付ける変数の値です。 たとえば、最初の人は32歳で、2人目の人は身長180 cmです。 しかし、すべてのオブジェクトに対して変数を一度に記述したい場合、つまり 列全体を取りますか？ この場合、1つの特定の値ではなく、一度に複数の値があり、それぞれに固有の発生頻度があります。 可能な値のリスト+対応する確率は、 ランダム変数 （ランダム変数、rv）と呼ばれます。

離散および連続確率変数

これを頭に留めておくために、もう一度繰り返します。ランダム変数は、その値の確率分布によって完全に決定されます。 ランダム変数には、離散型と連続型の2つの主なタイプがあります。

離散変数は、明確に分離可能な値のセットを取ることができます。 通常、私はそれらを次のように描写します（確率質量関数、pmf）：

Pkg.add("Plots") using Plots plotly() plot(["0","1"], [0.3, 0.7], linetype=:bar, legend=false) 

テキストでは、これは通常次のように書かれています（g-性別）：

つまり サンプルからランダムな人が女性になる確率（ $$$g = 0$ ）は0.3で、男（ $$$g = 1$ ）-0.7、これは女性の30％と男性の70％がサンプルに含まれていたという事実に相当します。

離散変数には、1人あたりの子供の数、テキスト内の単語の出現頻度、映画が視聴された回数などが含まれます。 ちなみに、有限数のクラスへの分類の結果も、離散確率変数です。

連続変数は、特定の間隔で任意の値を取ることができます。 たとえば、身長が175cmであると記録した場合でも、つまり 1センチメートルに丸めます。実際には175.8231 cmになります。 連続変数は通常、確率密度関数（pdf）曲線を使用して描かれます。

 Pkg.add("Distributions") using Distributions xs = 140:0.1:200 ys = [pdf(Normal(172, 10), x) for x in xs] plot(xs, ys; xlabel="h", ylabel="p(h)", legend=false, show=true) 

確率密度のグラフは注意が必要です：各列の高さがそのような値を得る確率を直接示す離散変数の確率質量のグラフとは異なり、確率密度は特定のポイント周辺の相対的な確率量を示します。 この場合、確率自体は区間についてのみ計算できます。 たとえば、この例では、サンプルから無作為に抽出された人の身長が160〜170 cmになる確率は約0.3です。

 d = Normal(172, 10) prob = cdf(d, 170) - cdf(d, 160) 

質問：ある時点での確率密度は1よりも大きくなりますか？ もちろん、答えはイエスです。主なことは、グラフの下の総面積（または数学的に言えば、確率密度積分）が1に等しいことです。

連続変数のもう1つの難点は、その確率密度が常にうまく表現できるとは限らないことです。 離散変数の場合、値のテーブル->確率がありました。 連続の場合、これらは一般に無限の数の意味を持つため、機能しません。 したがって、彼らは通常、よく研究されたパラメトリック分布によってデータセットを近似しようとします。 たとえば、上のグラフはいわゆる例です。 正規分布。 その確率密度は次の式で与えられます。

どこで $$$\ mu$ （mat。期待、平均）および $$$\シグマ^ 2$ （分散）-分布パラメーター。 つまり 2つの数値しかないため、分布を完全に記述し、任意の点での確率密度または2つの値間の合計確率を計算できます。 残念ながら、どのデータセットからも遠く離れて、それを美しく表現できる分布があります。 これに対処するには多くの方法があります（少なくとも正規分布の混合を取ります）が、これは完全に異なるトピックです。

継続的分布の他の例：人の年齢、画像内のピクセルの強度、サーバーからの応答時間など。

共同分布、周辺分布、条件付き分布

通常、オブジェクトのプロパティは一度に1つではなく、他のプロパティと組み合わせて考慮されます。ここでは、いくつかの変数の共同分布の概念が表示されます。 2つの離散変数の場合、テーブルの形式でそれを表すことができます（g-性別、c-子の数）：

この分布によると、データセットで2人の子供を持つ女性に会う確率は次のとおりです。 $$$p（g = 0、c = 2）= 0.1$ と子供のない男- $$$p（g = 1、c = 0）= 0.2$ 。

たとえば、身長と年齢などの2つの連続変数の場合、再び分析分布関数を定義する必要があります $$$p（h、a）$ それを概算し、
たとえば、 多次元法線 。 これをテーブルに書くことはできませんが、描くことができます：

 d = MvNormal([172.0, 30.0], [10 0; 0 5]) xs = 160:0.1:180 ys = 22:1:38 zs = [pdf(d, [x, y]) for x in xs, y in ys] surface(zs) 

共同分布がある場合、残りの変数を単純に合計（離散の場合）または積分（連続の場合）することで、各変数の分布を個別に見つけることができます。

これは、テーブルの各行または列の合計として表され、結果をテーブルのフィールドに入れることができます。

だから再び $$$p（g = 0）= 0.3$ そして $$$p（g = 1）= 0.7$ 。 マージンプロセスは、結果の分布自体に名前を与えます-限界確率。

しかし、変数の1つの値がすでにわかっている場合はどうでしょうか？ たとえば、目の前に男性がいて、その子供の数の確率分布を取得したいことがわかりますか？ 共同確率表もここで役立ちます：前に男がいることを確実に知っているので、 $$$g = 1$ 、他のすべてのオプションを考慮から破棄し、1行のみを考慮することができます。

確率は何らかの方法で1つに合計する必要があるため、結果の値を正規化する必要があります。

既知の値を持つ別の変数の分布は、条件付き確率と呼ばれます。

チェーンルール

そして、これらのすべての確率は、チェーンルールと呼ばれる1つの簡単な公式によって結び付けられています（チェーンルール、差別化においてチェーンルールと混同しないでください）。

この式は対称であるため、これも実行できます。

ルールの解釈は非常に簡単です：if $$$p（x）$ -私が赤信号に行く確率、そして $$$p（y | x）$ -赤信号に目を向ける人がヒットする確率、赤信号に移動してヒットする共同確率は、これら2つのイベントの確率の積に正確に等しくなります。 しかし、一般的に、緑になります。

従属変数と独立変数

すでに述べたように、共同分布表がある場合、システムに関するすべてを知っています。変数の限界確率を計算したり、ある変数を別の既知の変数で条件付きで分布したりできます。 残念ながら、実際には、このようなテーブルをコンパイルする（または連続分布のパラメーターを計算する）ことはできません。 たとえば、1000個の単語の出現の共同分布を計算する場合は、次の表が必要です。

107150860718626732094842504906000181056140481170553360744375038837035105112493612
249319837881569585812759467291755314682518714528569231404359845775746985748039345
677748242309854210746050623711418779541821530464749835819412673987675591655439460
77062914571196477686542167660429831652624386837205668069376

（1e301を少し超える）セル。 比較のために、観測可能な宇宙の原子数は約1e81です。 おそらく、追加のメモリバーを購入するだけでは十分ではありません。

ただし、詳細が1つあります。すべての変数が互いに依存しているわけではありません。 明日雨が降る確率は、私が道路を横断して赤信号になるかどうかにほとんど依存しません。 独立変数の場合、一方から他方への条件付き分布は単純に周辺分布です。

正直に言うと、1000語の結合確率は次のように記述されます。

しかし、単語が互いに独立していると「単純に」仮定すると、式は次のようになります。

そして、確率を保つために $$$p（w_i）$ 1000ワードの場合、1000セルのみのテーブルが必要です。これはまったく問題ありません。

では、すべての変数を独立と見なさないのはなぜですか 残念ながら、大量の情報が失われます。 喉の痛みと発熱という2つの変数に応じて、患者がインフルエンザにかかる確率を計算したいとします。 それとは別に、喉の痛みは、病気と患者が大声で歌っているだけのことを示します。 温度が別に上昇している場合は、病気と、走りから戻ったばかりの事実の両方を示します。 しかし、体温とのどの痛みを同時に観察する場合、これは患者の病気休暇を処方する重大な理由です。

対数

文献では、確率だけでなくその対数が使用されることがよくあります。 なんで？ すべてがかなり平凡です：

1. 対数は単調に増加する関数です。 のために $$$p（x_1）$ そして $$$p（x_2）$ もし $$$p（x_1）> p（x_2）$ それから $$$\ log p（x_1）> \ log p（x_2）$ 。
2. 積の対数は、対数の合計に等しくなります。 $$$\ log（p（x_1）p（x_2））= \ log p（x_1）+ \ log p（x_2）$ 。

単語を含む例では、任意の単語に出会う確率 $$$p（w_i）$ 、原則として、統一よりもはるかに少ない。 計算精度が限られているコンピューターで多くの小さな確率を掛けようとすると、どうなるのでしょうか？ ええ、私たちの確率はすぐにゼロに丸められます。 ただし、多数の個別の対数を追加すると 、計算の精度の限界を超えることは事実上不可能になります。

関数としての条件付き確率

これらのすべての例の後、特定の値が発生する回数をカウントすることで条件付き確率が常に計算されるという印象がある場合、私はこのエラーを払拭することを急いでいます：一般的な場合、条件付き確率は別のランダム変数の関数です：

どこで $$$\イプシロン$ -これはノイズです。 ノイズの種類-これは別のトピックでもありますが、ここでは取り上げませんが、機能については $$$f（x）$ もっと詳しく見てみましょう。 上記の離散変数の例では、関数として単純な出現回数を使用しました。 これは、それ自体で多くの場合、たとえば、テキストまたはユーザーの振る舞いに対する単純なベイジアン分類器でうまく機能します。 もう少し複雑なモデルは線形回帰です：

ここでも、変数は $$$x_i$ 互いに独立しているが分布 $$$p（y | \ mathbf {x}）$ パラメーターが線形関数を使用して既にモデル化されている $$$\ mathbf {\ theta}$ 見つける必要があります。

多層パーセプトロンも機能ですが、すべての入力変数の影響を一度に受ける中間層のおかげで、MLPを使用すると、個々の変数だけでなく、入力の組み合わせに対する出力変数の依存性をシミュレートすることができます（喉と温度の例を思い出してください）。

畳み込みネットワークは、フィルターサイズでカバーされるローカルエリアのピクセル分布で動作します。 リカレントネットワークは、前のデータと入力データからの次の状態の条件付き分布と、現在の状態からの出力変数をモデル化します。 まあ、一般的に、あなたはアイデアを得る。

ベイズの定理と連続変数の乗算

ネットワークルールを覚えていますか？

左側を削除すると、単純で明白な同等性が得られます。

そして、今転送する場合 $$$p（x）$ 右に、有名なベイズの公式を取得します。

式は非常に打ち解かれているため、各部分には独自の名前が付いています。

• $$$p（y）$ 事前配布（事前）と呼ばれます。 これは、特定のオブジェクト（たとえば、期限内にローンを支払った人の総数）を見る前からわかっていることです。
• $$$p（x | y）$ 尤度と呼ばれます。 これは、そのようなオブジェクトを見る確率です（変数によって記述されます $$$x$ ）出力変数のこの値で $$$y$ 。 たとえば、融資を行った人に2人の子供がいる可能性。
• $$$p（x）= \ int p（x、y）dy$ -限界妥当性、一般的にそのようなオブジェクトを見る確率。 みんな同じです $$$y$ 、したがって、ほとんどの場合考慮されず、ベイズの式の分子によって単純に最大化されます。
• $$$p（y | x）$ -事後分布（事後）。 これは変数の確率分布です $$$y$ オブジェクトを見た後。 たとえば、2人の子供を持つ人が期限内にローンを返済する可能性。

ベイジアン統計はとてつもなく興味深いものですが、これからは入りません。 私が触れたい唯一の質問は、連続変数の2つの分布の乗算です。これは、たとえば、ベイズ式の分子、および実際には連続変数の2番目の式すべてにあります。

2つの分布があるとしましょう $$$p_1（y）$ そして $$$p_2（y）$ ：

 d1 = Normal(175, 5) d2 = Normal(168, 5) space = 150:0.1:200 y1 = [pdf(d1, y) for y in space] y2 = [pdf(d2, y) for y in space] plot(space, y1, label="p_1(y)") plot!(space, y2, label="p_2(y)") 

そして、彼らの製品を入手したいのです。

各ポイントでの両方の分布の確率密度を知っているので、正直なところ、一般的な場合、各ポイントで密度を掛ける必要があります。 しかし、私たちがうまく行けば、 $$$p_1（y）$ そして $$$p_2（y）$ たとえば、2つの数値による正規分布（期待値と分散）のパラメーターを指定し、それらの製品の各ポイントでの確率を​​考慮する必要がありますか？

幸いなことに、多くの既知の分布の積は、簡単に計算可能なパラメーターを持つ別の既知の分布を与えます。 ここでのキーワードは、 共役事前です。

計算方法に関係なく、2つの正規分布の積により、もう1つの正規分布が得られます（正規化されていません）。

 #        # ,      plot(space, y1 .* y2, label="p_1(y)p_2(y)") 

さて、比較のために、3つの正規分布の混合の分布：

 plot(space, [pdf(Normal(130, 5), x) for x in space] .+ [pdf(Normal(150, 20), x) for x in space] .+ [pdf(Normal(190, 3), x) for x in space]) 

ご質問

これはチュートリアルであり、おそらく誰かがここに書かれた内容を思い出したいと思うでしょうから、ここに資料を修正するためのいくつかの質問があります。

人間の成長をパラメータを持つ正規分布のランダム変数とする $$$\ mu = 172$ そして $$$\シグマ^ 2 = 10$ 。 身長178cmの人と出会う確率はどのくらいですか？

させる $$$x$ -ローンの借り手である子供の数（3つの可能な値）、 $$$y$ -人がローンを提供したかどうかのサイン（2つの可能な値）。 ベイズの公式を使用して、子供が1人いる特定の顧客がローンを提供するかどうかを予測します。 先験的分布と事後分布は、尤度と限界尤度と同様に、いくつの値を取ることができますか？

2つの分布間の距離を最適化するニューラルネットワーク $$$q（x）$ そして $$$p（x）$ 多くの場合、最適化の目標としてクロスエントロピーまたはKullback-Leibler発散距離を使用します。 後者は次のように定義されます：

$$$\ int q（x）（。）dx$ -それは仲間です。 待っている $$$q（x）$ 、そしてなぜ主な部分で- $$$\ log \ frac {q（x）} {p（x）}$ -2つの関数の密度の差だけでなく、除算が使用されます $$$q（x）-p（x）$ ？