Comment tracer une distribution normale de journal dans R

Pour tracer la fonction de densité de probabilité pour une distribution log-normale dans R, nous pouvons utiliser les fonctions suivantes :

• dlnorm(x, meanlog = 0, sdlog = 1) pour créer la fonction de densité de probabilité.
• courbe(fonction, de = NULL, à = NULL) pour tracer la fonction de densité de probabilité.

Par exemple, le code suivant illustre comment tracer une fonction de densité de probabilité pour une distribution log-normale avec moyenne = 0 et écart type = 1 (sur une échelle logarithmique) où l’axe des X du tracé va de 0 à 10 :

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

Par défaut, meanlog = 0 et sdlog =1, ce qui signifie que nous pouvons produire exactement le même tracé sans spécifier ces paramètres dans la fonction dlnorm() :

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

Nous pouvons ajouter un titre, modifier l’étiquette de l’axe Y, augmenter la largeur de la ligne et même changer la couleur de la ligne pour rendre l’intrigue plus esthétique :

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 

Nous pouvons également ajouter plusieurs courbes au graphique pour comparer les distributions log-normales avec différents écarts types. Par exemple, le code suivant crée des tracés de distribution normale avec sdlog = 0,3, sdlog = 0,5 et 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)

Nous pouvons ajouter une légende au tracé en utilisant la fonction legend() , qui adopte la syntaxe suivante :

légende(x, y=NULL, légende, remplissage, col, bg, lty, cex)

où:

• x, y : les coordonnées x et y utilisées pour positionner la légende
• légende : le texte à mettre dans la légende
• fill : couleur de remplissage à l’intérieur de la légende
• col : la liste des couleurs à utiliser pour les lignes à l’intérieur de la légende
• bg : la couleur de fond de la légende
• lty : style de ligne
• cex : taille du texte dans la légende

Dans notre exemple nous utiliserons la syntaxe suivante pour créer une légende :

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