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)
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
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)