Puede usar la biblioteca de visualización de R ggplot2 para trazar un modelo de regresión lineal ajustado usando la siguiente sintaxis básica:
ggplot (datos, aes (x, y)) + geom_point () + geom_smooth (método = ' lm ')
El siguiente ejemplo muestra cómo utilizar esta sintaxis en la práctica.
Ejemplo: trazar una línea de regresión lineal en ggplot2
Supongamos que ajustamos un modelo de regresión lineal simple al siguiente conjunto de datos:
#create dataset
data <- data.frame (y = c (6, 7, 7, 9, 12, 13, 13, 15, 16, 19, 22, 23, 23, 25, 26),
x = c (1, 2, 2, 3, 4, 4, 5, 6, 6, 8, 9, 9, 11, 12, 12))
#ajustar el modelo de regresión lineal al conjunto de datos y ver el modelo de resumen del
modelo <- lm (y ~ x, data = data)
resumen (modelo)
Llamada:
lm (fórmula = y ~ x, datos = datos)
Derechos residuales de autor:
Mín. 1T Mediana 3T Máx.
-1,4444 -0,8013 -0,2426 0,5978 2,2363
Coeficientes:
Estimar Std. Valor t de error Pr (> | t |)
(Intercepción) 4.20041 0.56730 7.404 5.16e-06 ***
x 1.84036 0.07857 23.423 5.13e-12 ***
---
Signif. códigos: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0,1 pulg. 1
Error estándar residual: 1.091 en 13 grados de libertad
R cuadrado múltiple: 0,9769, R cuadrado ajustado: 0,9751
Estadístico F: 548,7 en 1 y 13 DF, valor de p: 5,13e-12
El siguiente código muestra cómo visualizar el modelo de regresión lineal ajustado:
biblioteca (ggplot2) #create plot para visualizar el modelo de regresión lineal ajustado ggplot (datos, aes (x, y)) + geom_point () + geom_smooth (método = ' lm ')
De forma predeterminada, ggplot2 agrega líneas de error estándar al gráfico. Puede deshabilitarlos utilizando el argumento se = FALSE de la siguiente manera:
biblioteca (ggplot2) #crear gráfico de regresión sin líneas de error estándar ggplot (datos, aes (x, y)) + geom_point () + geom_smooth (método = ' lm ', se = FALSE )
Por último, podemos personalizar algunos aspectos del gráfico para hacerlo más atractivo visualmente:
biblioteca (ggplot2) #create un gráfico de regresión con un estilo personalizado ggplot (datos, aes (x, y)) + geom_point () + geom_smooth (método = ' lm ', se = FALSE , color = ' turquoise4 ') + theme_minimal () + laboratorios (x = ' Valores X ', y = ' Valores Y ', título = ' Gráfica de regresión lineal ') + tema (plot.title = element_text (hjust = 0.5, size = 20, face = ' bold '))
Consulte esta publicación para obtener una guía completa de los mejores temas de ggplot2.
Recursos adicionales
Introducción a la regresión lineal múltiple en R
Cómo trazar un intervalo de confianza en R
- https://r-project.org
- https://www.python.org/
- https://www.stata.com/
