Gráfico de la función de densidad de la distribución normal en R

En R, vamos a crear el gráfico de la función de densidad de la distribución normal usando dos funciones: plot y curve.

Función plot

Creamos un vector con las coordenadas x, en los argumentos de seq (secuencia) especificamos los límites y el número de valores intermedios o el incremento. En las coordenadas y, la función de densidad dnorm, por defecto una distribución normal estándar, N(0, 1), con una media de 0 y una varianza de 1.

x <- seq(-3, 3, 0.001) # Incremento de la secuencia: 6001 valores
# O bien
x <- seq(-3, 3, length = 1000) # 1000 valores intermedios

Usamos la función plot:

plot(x,          # Coordenadas x
     dnorm(x),   # Coordenadas y
     type = "l", # Tipo de gráfico líneas
     las = 1)    # Alineación horizontal de etiquetas

Función curve

No es necesario crear un vector con las coordenadas x. Solamente especificar el nombre de la función que evaluará: dnorm. Por defecto, curve evalua 101 valores (n = 101).

curve(dnorm(x), xlim = c(-3, 3), las = 1)

Más ejemplos

Dos ejemplos con distribuciones distintas de la distribución normal estándar. En este caso, sí es necesario suministrar los argumentos: media y desviación típica (raíz cuadrada de la varianza). Añadimos el parámetro xaxp para forzar el número de marcas del eje x entre las dos coordenadas de los extremos.

N(-2, 0.5)

x <- seq(-5, 5, 0.001)
plot(x, dnorm(x, -2, sqrt(.5)), type = "l", las = 1, xaxp =c(-5, 5, 10))

N(0, 0.2)

curve(dnorm(x, 0, sqrt(.2)), -5, 5, las = 1, xaxp = c(-5, 5, 10))

Referencias:

?dnorm
?plot
?curve
?par