Para crear una gráfica de densidad para una distribución de Chi-cuadrado en R, podemos usar las siguientes funciones:
- dchisq () para crear la función de densidad de probabilidad
- curva () para trazar la función de densidad de probabilidad
Todo lo que necesitamos hacer para crear la gráfica es especificar los grados de libertad para dchisq () junto con los puntos hacia y desde la curva () .
Por ejemplo, el siguiente código ilustra cómo crear una gráfica de densidad para una distribución de chi-cuadrado con 10 grados de libertad donde el eje x de la gráfica varía de 0 a 40:
curva (dchisq (x, gl = 10), desde = 0, hasta = 40)
Modificación de la gráfica de densidad
También podemos modificar la gráfica de densidad agregando un título, cambiando la etiqueta del eje y, aumentando el ancho de la línea y modificando el color de la línea:
curva (dchisq (x, gl = 10), desde = 0, hasta = 40,
main = 'Chi-Square Distribution (df = 10)', #add title
ylab = 'Density', #cambiar la etiqueta del eje y
lwd = 2, #increase line width to 2
col = 'steelblue') #cambiar el color de la línea azul acero
Relleno de la gráfica de densidad
Además de crear la gráfica de densidad, podemos completar parte de la gráfica usando la función polygon () basada en un valor inicial y final.
El siguiente código ilustra cómo completar la parte de la gráfica de densidad para los valores de x que van de 10 a 40:
#creación de la curva de densidad
(dchisq (x, df = 10), from = 0, to = 40,
main = 'Chi-Square Distribution (df = 10)',
ylab = 'Density',
lwd = 2)
#crear vector de valores x
x_vector <- seq (10, 40)
#crear vector de valores de densidad chi-cuadrado
p_vector <- dchisq (x_vector, df = 10)
# rellene parte de la gráfica de densidad de 0 a 40
polígono (c (x_vector, rev (x_vector)), c (p_vector, rep (0, length (p_vector))),
col = ajustecolor ('rojo', alfa = 0.3), borde = NA)
El siguiente código ilustra cómo completar la parte de la gráfica de densidad para los valores de x que van de 0 a 10:
#creación de la curva de densidad
(dchisq (x, df = 10), from = 0, to = 40,
main = 'Chi-Square Distribution (df = 10)',
ylab = 'Density',
lwd = 2)
#crear vector de valores x
x_vector <- seq ( 0, 10 )
#crear vector de valores de densidad chi-cuadrado
p_vector <- dchisq (x_vector, df = 10)
# rellene una parte de la gráfica de densidad de 0 a 10
polígono (c (x_vector, rev (x_vector)), c (p_vector, rep (0, length (p_vector))),
col = ajustecolor ('rojo', alfa = 0.3), borde = NA)
El siguiente código ilustra cómo completar la parte de la gráfica de densidad para los valores x que se encuentran fuera del 95% medio de la distribución:
#crear curva de densidad de curva
(dchisq (x, df = 10), de = 0, a = 40,
main = 'Distribución de chi-cuadrado (df = 10)',
ylab = 'Densidad',
lwd = 2)
# encontrar valores superiores e inferiores para el 95% medio de la distribución
inferior 95 <- qchisq (.025, 10)
superior 95 <- qchisq (.975, 10)
#create vector de valores x
x_lower95 <- seq (0, lower95)
#create vector de valores de densidad chi-cuadrado
p_lower95 <- dchisq (x_lower95, df = 10)
# rellene una parte del gráfico de densidad de 0 a un valor inferior del 95% del
polígono (c (x_lower95, rev (x_lower95)), c (p_lower95, rep (0, length (p_lower95))),
col = ajustecolor ('rojo', alfa = 0.3), borde = NA)
#create vector de valores x
x_upper95 <- seq (upper95, 40)
#crear vector de valores de densidad chi-cuadrado
p_upper95 <- dchisq (x_upper95, df = 10)
# rellene la parte de la gráfica de densidad para el valor superior del 95% hasta el final del polígono de la gráfica
(c (x_upper95, rev (x_upper95)), c (p_upper95, rep (0, length (p_upper95))),
col = ajustecolor ('rojo', alfa = 0.3), borde = NA)
Por último, el siguiente código ilustra cómo completar la parte de la gráfica de densidad para los valores x que se encuentran dentro del 95% medio de la distribución:
#crear curva de densidad de curva
(dchisq (x, df = 10), de = 0, a = 40,
main = 'Distribución de chi-cuadrado (df = 10)',
ylab = 'Densidad',
lwd = 2)
# encontrar valores superiores e inferiores para el 95% medio de la distribución
inferior 95 <- qchisq (.025, 10)
superior 95 <- qchisq (.975, 10)
#crear vector de valores x
x_vector <- seq (lower95, upper95)
#crear vector de valores de densidad chi-cuadrado
p_vector <- dchisq (x_vector, df = 10)
# rellenar el polígono de la parcela de densidad
(c (x_vector, rev (x_vector)), c (p_vector, rep (0, length (p_vector))),
col = ajustecolor ('rojo', alfa = 0.3), borde = NA)
- https://r-project.org
- https://www.python.org/
- https://www.stata.com/
