So zeichnen Sie einfach eine Chi-Quadrat-Verteilung in R

Um ein Dichtediagramm für eine Chi-Quadrat-Verteilung in R zu erstellen, können wir die folgenden Funktionen verwenden:

dchisq(), um die Wahrscheinlichkeitsdichtefunktion zu erstellen
curve() zur Darstellung der Wahrscheinlichkeitsdichtefunktion

Alles, was wir tun müssen, um den Plot zu erstellen, ist die Freiheitsgrade für dchisq() angeben zusammen mit dem zu und von den Punkten für curve().

Der folgende Code veranschaulicht beispielsweise, wie ein Dichtediagramm für eine Chi-Quadrat-Verteilung mit 10 Freiheitsgraden erstellt wird, wobei die x-Achse des Diagramms zwischen 0 und 40 liegt:

curve(dchisq(x, df = 10), from = 0, to = 40)

Sie können das Dichtediagramm auch ändern, indem Sie einen Titel hinzufügen, die Beschriftung der y-Achse ändern, die Linienbreite erhöhen und die Linienfarbe ändern:

curve(dchisq(x, df = 10), from = 0, to = 40,
      main = 'Chi-Square Distribution (df = 10)', #Titel hinzufügen
      ylab = 'Density', #y-Achsenlabel ändern
      lwd = 2, #Linenstärke auf 2 erhöhen
      col = 'steelblue') #Linienfarbe zu steelblue ändern

Ausfüllen des Dichtediagramms

Zusätzlich zum Erstellen des Dichtediagramms können wir einen Teil des Diagramms mithilfe der Funktion polygon() basierend auf einem Start- und Endwert ausfüllen.

Der folgende Code zeigt, wie der Teil des Dichtediagramms für die x-Werte zwischen 10 und 40 ausgefüllt wird:

#Dichtekurve erstellen
curve(dchisq(x, df = 10), from = 0, to = 40,
      main = 'Chi-Square Distribution (df = 10)',
      ylab = 'Density',
      lwd = 2)

#Vektor von x-Werten erstellen
x_vector <- seq(10, 40)
#Erstellen Sie einen Vektor mit Chi-Quadrat-Dichtewerten
p_vector <- dchisq(x_vector, df = 10)
# Füllen Sie einen Teil des Dichtediagramms von 0 bis 40 aus
polygon(c(x_vector, rev(x_vector)), c(p_vector, rep(0, length(p_vector))),
        col = adjustcolor('red', alpha=0.3), border = NA)

Der folgende Code zeigt, wie der Teil des Dichtediagramms für die x-Werte im Bereich von 0 bis 10 ausgefüllt wird:

#Dichtekurve erstellen
curve(dchisq(x, df = 10), from = 0, to = 40,
      main = 'Chi-Square Distribution (df = 10)',
      ylab = 'Density',
      lwd = 2) 

#Vektor von x-Werten erstellen
x_vector <- seq(0, 10)

#Erstellen Sie einen Vektor mit Chi-Quadrat-Dichtewerten
p_vector <- dchisq(x_vector, df = 10)
# Füllen Sie einen Teil des Dichtediagramms von 0 bis 10 aus
polygon(c(x_vector, rev(x_vector)), c(p_vector, rep(0, length(p_vector))),
        col = adjustcolor('red', alpha=0.3), border = NA)

Der folgende Code zeigt, wie der Teil des Dichtediagramms für die x-Werte ausgefüllt wird, die außerhalb der mittleren 95% der Verteilung liegen:

#Dichtekurve erstellen
curve(dchisq(x, df = 10), from = 0, to = 40,
        main = 'Chi-Square Distribution (df = 10)',
        ylab = 'Density',
        lwd = 2)

#Finden Sie obere und untere Werte für mittlere 95% der Verteilung
lower95 <- qchisq(.025, 10)
upper95 <- qchisq(.975, 10)

#Vektor von x-Werten erstellen
x_lower95 <- seq(0, lower95)
#Erstellen Sie einen Vektor mit Chi-Quadrat-Dichtewerten
p_lower95 <- dchisq(x_lower95, df = 10)
# Füllen Sie einen Teil des Dichtediagramms von 0 bis zu einem niedrigeren 95%-Wert aus
polygon(c(x_lower95, rev(x_lower95)), c(p_lower95, rep(0, length(p_lower95))),
        col = adjustcolor('red', alpha=0.3), border = NA)

# Vektor von x-Werten erstellen
x_upper95 <- seq(upper95, 40)
#Erstellen Sie einen Vektor mit Chi-Quadrat-Dichtewerten
p_upper95 <- dchisq(x_upper95, df = 10)
# Füllen Sie einen Teil des Dichtediagramms für den oberen 95%-Wert bis zum Ende des Diagramms aus
polygon(c(x_upper95, rev(x_upper95)), c(p_upper95, rep(0, length(p_upper95))),
        col = adjustcolor('red', alpha=0.3), border = NA)

Der folgende Code zeigt schließlich, wie der Teil des Dichtediagramms für die x-Werte ausgefüllt wird, die innerhalb der mittleren 95% der Verteilung liegen:

#Dichtekurve erstellen
curve(dchisq(x, df = 10), from = 0, to = 40,
        main = 'Chi-Square Distribution (df = 10)',
        ylab = 'Density',
        lwd = 2)

# Finden Sie obere und untere Werte für mittlere 95% der Verteilung
lower95 <- qchisq(.025, 10)
upper95 <- qchisq(.975, 10)

# Vektor von x-Werten erstellen
x_vector <- seq(lower95, upper95)
#Erstellen Sie einen Vektor mit Chi-Quadrat-Dichtewerten
p_vector <- dchisq(x_vector, df = 10)
#Dichtediagramm ausfüllen
polygon(c(x_vector, rev(x_vector)), c(p_vector, rep(0, length(p_vector))),
        col = adjustcolor('red', alpha=0.3), border = NA)

Ausfüllen des Dichtediagramms

Statistik: Der Weg zur Datenanalyse

So zeichnen Sie eine Gamma-Verteilung in Python (mit Beispielen)

So verwenden Sie die einheitliche Verteilung in Python