Em outra pergunta viu-se como colocar a linha de regressão em um gráfico. Mas, como colocar a equação da regressão no gráfico? Por exemplo:
Ou
Em outra pergunta viu-se como colocar a linha de regressão em um gráfico. Mas, como colocar a equação da regressão no gráfico? Por exemplo:
Ou
Segue uma possível solução, usando geom_text
:
set.seed(1)
x <- rnorm(100)
y <- rnorm(100) + 2*x +10
modelo <- lm(y~x)
coeficientes <- modelo$coefficients
texto <- sprintf('y = %.2f + %.2fx, r² = %.2f', coeficientes[1], coeficientes[2], summary(modelo)$r.squared)
require(ggplot2)
dados <- data.frame(x=x, y=y) # O ggplot2 exige que os dados estejam em um data.frame
p <- ggplot(dados, aes(x=x, y=y)) + # Informa os dados a serem utilizadps
geom_point() +
geom_text(aes(x=min(x), y=max(y), label=texto), hjust=0, vjust=1)
p
O hjust=0
serve para informar que o x indicado é o do limite esquerdo do texto (padrão é o centro, .5) e vjust=1
serve para informar que o y indicado é o limite superior do texto (mesmo padrao do hjust).
Também é possível automatizar esses output com a função stat_poly_eq
do pacote ggpmisc.
library(ggplot2)
library(ggpmisc)
set.seed(1)
x <- rnorm(100)
y <- rnorm(100) + 2*x +10
dados <- data.frame(x = x, y = y)
formula <- y ~ x
ggplot(dados, aes(x = x, y = y)) +
geom_point() +
geom_smooth(method = "lm", se = T, color = "black", formula = formula) +
ggpmisc::stat_poly_eq(aes(label = paste(..eq.label.., ..adj.rr.label..,
sep = "~~~~")),
formula = formula, parse = TRUE)
ggplot(dados, aes(x = x, y = y)) +
geom_point() +
ggpmisc::stat_poly_eq(aes(label = paste(..eq.label.., ..adj.rr.label..,
sep = "~~~~")),
formula = formula, parse = TRUE)
Para ficar registrada mais uma opção: a função ggpubr::stat_regline_equation
. A sintaxe é a mesma da ggpmisc::stat_poly_eq
:
library(ggplot2)
library(ggpubr)
ggplot(dados, aes(x, y)) +
geom_point() +
geom_smooth(method = "lm") +
stat_regline_equation(aes(label = paste(..eq.label.., ..rr.label.., sep = "*plain(\",\")~~")))
Se a fórmula não for especificada, assume-se y ~ x.