[PYTHON] Besiege die Wahrscheinlichkeitsdichtefunktion der Normalverteilung

Zusammenfassung

f(x)=\frac{1}{\sqrt{2\piσ^2}}e^{-\frac{(x-μ)^2}{2σ^2}}

Die Normalverteilung ist eine der typischen Verteilungen und eine Verteilung, die im Bereich der Statistik häufig vorkommt. Die Wahrscheinlichkeitsdichtefunktion der Normalverteilung ist jedoch wie die ↑ Formel sehr kompliziert und hat viele Anfänger begraben. Ich erinnere mich auch an Kopfschmerzen und Schwindel, als ich diese Formel zum ersten Mal sah, als ich anfing, Statistik zu studieren.

In diesem Artikel werden der Koeffiziententeil ($ \ frac {1} {\ sqrt {2 \ piσ ^ 2}} ) und der Exponententeil ( - \ frac {(x-μ) ^ 2} {2σ ^ 2) der Formel ↑ } $) und beschreiben Sie, warum es so geformt wurde.

Indexteil

Viele Ereignisse auf der Welt haben die höchste Wahrscheinlichkeit, einen Durchschnittswert anzunehmen, und die Wahrscheinlichkeit, diesen Wert anzunehmen, nimmt mit zunehmender Entfernung vom Durchschnittswert ab.

Ein Ausdruck, der dies einfach ausdrückt, ist

f(x)=e^{-x^2}

ist. In einem Diagramm ausgedrückt

import numpy as np
import matplotlib.pyplot as plt

def normal_dist(x, ave = 0, disp=1):
    return np.exp((-(x-ave)**2)/(2*disp**2))

x = np.linspace(-3, 3)
y = normal_dist(x)

plt.plot(x,y)
plt.grid(axis="both")
plt.show()

image.png Es sieht aus wie ein Berg.

Diese Formel ist in ihrer Vielseitigkeit minderwertig, da nur ein Diagramm gezeichnet werden kann. Deshalb, ① Bewegen Sie sich links und rechts vom Diagramm ② Ändern Sie die Diagrammbreite Fügen Sie die Funktion von hinzu.

Bewegen Sie sich links und rechts vom Diagramm

Durch Ändern des $ x $ -Teils in $ x-μ $ können Sie die Position von $ x $ verschieben, die den Maximalwert annimmt.

f(x)=e^{-(x-μ)^2}

Es wird eine solche Formel. Lassen Sie uns den Wert von $ μ $ ändern und die Änderung im Diagramm sehen.

color = ["b", "g", "r", "c", "m"]
for i, col in enumerate(color):
    y = normal_dist(x, i)
    plt.plot(x, y, color=col)
    
plt.show()

image.png Der Graph verschob sich nach rechts, als $ μ $ zunahm.

Ändern Sie die Diagrammbreite

Sie können die Breite ändern, indem Sie den Exponententeil mit $ \ frac {1} {2σ ^ 2} $ multiplizieren.

f(x)=e^{-\frac{x^2}{2σ^2}}
x = np.linspace(-10, 10)
color = ["b", "g", "r", "c", "m"]
for i, col in enumerate(color):
    y = normal_dist(x, 0, i+1)
    plt.plot(x, y, color=col)
    
plt.show()

image.png

Es ist mir gelungen, die Breite zu ändern. Das Quadrat von σ soll sowohl positiv als auch negativ sein, und die Multiplikation von $ 2 $ soll die Berechnung erleichtern.

Faktor Teil

Da es sich um eine Wahrscheinlichkeitsdichtefunktion handelt, muss die Gesamtfläche $ 1 $ betragen. Multiplizieren Sie daher die Funktion mit einem geeigneten Wert. Geben Sie einen geeigneten Koeffizienten als $ c $ und ein

\int_{-\infty}^{\infty}ce^{-\frac{(x-μ)^2}{2σ^2}}dx=1

Um $ c $ zu finden.

Auf den ersten Blick mag es wie eine komplizierte Berechnung erscheinen,

\int_{-\infty}^{\infty}e^{-ax^2}dx=\sqrt{\frac{\pi}{a}}

Verwenden Sie diese Gaußsche Integralformel, wenn $ a = \ frac {1} {2σ ^ 2} $

c=\frac{1}{\sqrt{2\piσ^2}}

Es wird leicht zu berechnen. Die Lösung ist jetzt der Koeffizient der Formel oben.

Immerhin war diese Formel eine Formel, die die Wahrscheinlichkeitsdichtefunktion definiert, indem $ f (x) = e ^ {-x ^ 2} $ vielseitig gemacht und die Koeffizienten angepasst werden. Ich bin glücklich.

Referenzierte Site

https://to-kei.net/distribution/normal-distribution/density-function-derivation/ Semantisches Verständnis der Dichtefunktion der Normalverteilung

https://mathtrain.jp/gauss Zwei Beweise der Gaußschen Integralformel

https://bellcurve.jp/statistics/course/7797.html 14-1. Normalverteilung

Recommended Posts

Besiege die Wahrscheinlichkeitsdichtefunktion der Normalverteilung
Versuchen Sie, die stochastische Massenfunktion der Binomialverteilung in Python zu transkribieren
Überprüfung der Normalverteilung
[Python] Hinweis: Selbst erstellte Funktion zum Ermitteln des Bereichs der Normalverteilung
Schritte zur Berechnung der Wahrscheinlichkeit einer Normalverteilung
Tweet die Niederschlagswahrscheinlichkeit als Teil der Funktion des Bots
Überprüfen Sie die atrophische Natur der Wahrscheinlichkeitsverteilung in Python
PRML-Diagrammzeichnung Aufgabe 1.4 Nichtlineare Transformation der Wahrscheinlichkeitsdichtefunktion
Ist die Niederschlagswahrscheinlichkeit korrekt?
Testen Sie die Eignung der Verteilung
Leiten Sie den Quadrantenbereich der Standardnormalverteilung von Anfang an sorgfältig ab
Vorsichtsmaßnahmen bei der Überlagerung der Wahrscheinlichkeitsdichtefunktion und des Histogramms in matplotlib
[Python3] Schreiben Sie das Codeobjekt der Funktion neu
Über die Argumente der Setup-Funktion von PyCaret
Über die Normalgleichung der linearen Regression
Text Mining: Wahrscheinlichkeitsdichteverteilung auf Hypersphäre und Textclustering mit KMeans
Ich habe die Pivot-Table-Funktion von Pandas ausprobiert
Holen Sie sich den Aufrufer einer Funktion in Python
Passen Sie die Verteilung jeder Gruppe in Python an
Konzept des Bayes'schen Denkens (2) ... Bayes'sche Schätzung und Wahrscheinlichkeitsverteilung
Korrigieren Sie die Argumente der in map verwendeten Funktion
Berechnen Sie die Wahrscheinlichkeit von Ausreißern auf den Box-Whiskern
[CodeIQ] Ich habe die Wahrscheinlichkeitsverteilung von Würfeln geschrieben (aus dem CodeIQ-Mathematikkurs für maschinelles Lernen [Wahrscheinlichkeitsverteilung]).
Bivariate Normalverteilung
Finden Sie die kumulative Verteilungsfunktion durch Sortieren (Python-Version)
Zeichnen und verstehen Sie die multivariate Normalverteilung in Python
#Eine Funktion, die den Zeichencode einer Zeichenfolge zurückgibt
Zeichnen auf Jupyter mit der Plot-Funktion von Pandas
Ich habe das Argument class_weight von Chainers Funktion softmax_cross_entropy untersucht.
Python: Diagramm der zweidimensionalen Datenverteilung (Schätzung der Kerneldichte)
Wie man das Dokument der magischen Funktion (Linienmagie) trifft
Überprüfen Sie den Linux-Verteilungstyp und die Version
So ermitteln Sie die durchschnittliche Informationsmenge (Entropie) der ursprünglichen Wahrscheinlichkeitsverteilung aus der Stichprobe