So zeichnen sie eine logarithmische normalverteilung in r auf


Um die Wahrscheinlichkeitsdichtefunktion für eine logarithmische Normalverteilung in R darzustellen, können wir die folgenden Funktionen verwenden:

  • dlnorm(x, meanlog = 0, sdlog = 1) , um die Wahrscheinlichkeitsdichtefunktion zu erstellen.
  • Kurve(Funktion, von = NULL, bis = NULL), um die Wahrscheinlichkeitsdichtefunktion darzustellen.

Der folgende Code veranschaulicht beispielsweise, wie eine Wahrscheinlichkeitsdichtefunktion für eine logarithmische Normalverteilung mit Mittelwert = 0 und Standardabweichung = 1 (auf einer logarithmischen Skala) dargestellt wird, wobei die x-Achse des Diagramms von 0 bis 10 reicht:

 curve(dlnorm(x, meanlog=0, sdlog=1), from=0, to=10)

Zeichnen Sie ein Normalverteilungsdiagramm in R

Standardmäßig sind „meanlog = 0“ und „sdlog =1“, was bedeutet, dass wir genau denselben Plot erstellen können, ohne diese Parameter in der Funktion dlnorm() anzugeben:

 curve(dlnorm(x), from=0, to=10)

Wir können einen Titel hinzufügen, die Beschriftung der Y-Achse ändern, die Zeilenbreite erhöhen und sogar die Zeilenfarbe ändern, um den Plot ästhetisch ansprechender zu gestalten:

 curve(dlnorm(x), from=0, to=10, 
    main = 'Log Normal Distribution', #add title
    ylab = 'Density', #change y-axis label
    lwd = 2, #increase line width to 2
    col = 'steelblue') #change line color to steelblue 

Protokollieren Sie das Normalverteilungsdiagramm mit dem Titel in R

Wir können dem Diagramm auch mehrere Kurven hinzufügen, um Lognormalverteilungen mit unterschiedlichen Standardabweichungen zu vergleichen. Der folgende Code erstellt beispielsweise Normalverteilungsdiagramme mit sdlog = 0,3, sdlog = 0,5 und sdlog = 1:

 curve(dlnorm(x, meanlog=0, sdlog=.3), from=0, to=10, col='blue')
curve(dlnorm(x, meanlog=0, sdlog=.5), from=0, to=10, col='red', add=TRUE)
curve(dlnorm(x, meanlog=0, sdlog=1), from=0, to=10, col='purple', add=TRUE)

Mit der Funktion legend() können wir dem Plot eine Legende hinzufügen, die die folgende Syntax hat:

legend(x, y=NULL, legend, fill, col, bg, lty, cex)

Gold:

  • x, y: die x- und y-Koordinaten, die zur Positionierung der Legende verwendet werden
  • Legende: Der Text, der in die Legende eingefügt werden soll
  • Füllung: Füllfarbe innerhalb der Legende
  • col: die Liste der Farben, die für die Linien innerhalb der Legende verwendet werden sollen
  • bg: die Hintergrundfarbe der Legende
  • lty: Linienstil
  • cex: Größe des Textes in der Legende

In unserem Beispiel verwenden wir die folgende Syntax, um eine Legende zu erstellen:

 #create density plots
curve(dlnorm(x, meanlog=0, sdlog=.3), from=0, to=10, col='blue')
curve(dlnorm(x, meanlog=0, sdlog=.5), from=0, to=10, col='red', add=TRUE)
curve(dlnorm(x, meanlog=0, sdlog=1), from=0, to=10, col='purple', add=TRUE)

#add legend
legend(6, 1.2, legend=c("sdlog=.3", "sdlog=.5", "sdlog=1"),
       col=c("blue", "red", "purple"), lty=1, cex=1.2) 

Mehrere logarithmische Normaldichtefunktionen in einem Diagramm in R

Einen Kommentar hinzufügen

Deine E-Mail-Adresse wird nicht veröffentlicht. Erforderliche Felder sind mit * markiert