Cómo trazar fácilmente una distribución de chi-cuadrado en R

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/