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)

Gráfico de distribución de chi-cuadrado en R con 5 grados de libertad

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

Gráfico de densidad de chi-cuadrado en R

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)

Gráfico de distribución de chi-cuadrado en R con los valores rellenados

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)

Distribución de chi-cuadrado con gráfico de 10 grados de libertad

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)

Distribución de chi-cuadrado con 95% de valores fuera de la distribuciónPor ú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)

Distribución chi-cuadrado media 95%

  • https://r-project.org
  • https://www.python.org/
  • https://www.stata.com/

Redactor del artículo

  • Luis Benites
    Director de Statologos.com

    Tengo una Maestría en Ciencias en Estadística Aplicada y he trabajado en algoritmos de aprendizaje automático para empresas profesionales tanto en el sector de la salud como en el comercio minorista.

    Ver todas las entradas

¿Te hemos ayudado?

Ayudanos ahora tú, dejanos un comentario de agradecimiento, nos ayuda a motivarnos y si te es viable puedes hacer una donación:

La ayuda no cuesta nada

Por otro lado te rogamos que compartas nuestro sitio con tus amigos, compañeros de clase y colegas, la educación de calidad y gratuita debe ser difundida, recuerdalo:

Deja un comentario

Para trazar una distribución normal en Python, puede usar la siguiente sintaxis: # El eje x oscila entre -3 y…
statologos comunidad-2

Compartimos información EXCLUSIVA y GRATUITA solo para suscriptores (cursos privados, programas, consejos y mucho más)

You have Successfully Subscribed!