9

Suponha um modelo de regressão linear como o seguinte:

set.seed(1)
x <- rnorm(100)
y <- rnorm(100) + 2*x +10
modelo <- lm(y~x)

Se eu plotar y contra x, como faço para incluir a linha de regressão no gráfico?

plot(y~x)

inserir a descrição da imagem aqui

4 Respostas 4

9

Outra alternativa é utilizar o pacote ggplot2:

set.seed(1)
x <- rnorm(100)
y <- rnorm(100) + 2*x +10

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() # Informa que eu quero um gráfico de dispersão.
p

Gráfico 1

p1 <- p + geom_smooth(method=lm) # Acrescenta a linha de tendência e o intervalo de confiança de predição
p1

Gráfico 2

p2 <- p + geom_smooth(method=lm, se=FALSE) # Acrescenta a linha de tendência, sem o intervalo de predição
p2

Gráfico 3

4
  • 1
    o método tem que estar entre aspas. geom_smooth(method="lm") Commented 18/02/2014 às 15:37
  • 1
    Certeza? Aqui funciona sem as aspas.
    – Rcoster
    Commented 18/02/2014 às 16:46
  • 1
    Sim, certeza. Pode ser problemas com as versões. R Under development (unstable) (2014-02-06 r64930) Platform: x86_64-unknown-linux-gnu (64-bit) other attached packages: [1] ggplot2_0.9.3.1.99 são as partes que eu achei relevantes de sessionInfo(). Vou conferir se pode ser um erro de instalação da minha parte. Commented 18/02/2014 às 17:00
  • 1
    Aqui: R version 3.0.2 (2013-09-25); Platform: x86_64-w64-mingw32/x64 (64-bit); ggplot2_0.9.3.1
    – Rcoster
    Commented 18/02/2014 às 20:02
8

Você pode usar a função abline junto com coef para extrair os coeficientes do modelo e plotar a linha:

plot(y~x)
abline(coef(modelo))

linha de regressão

5

Uma outra maneira é usar a função curve.

curve(coef(modelo)[1]+coef(modelo)[2]*x,add=TRUE,col = "blue",lwd=2)

Uma vantagem dessa função é que ela também serve caso o modelo envolva componentes não lineares. Como no caso abaixo:

n=1000
x1 = rnorm(n)
y = 5 +  4*x1 + 2*x1^2 + rnorm(n)
modelo <- lm( y ~ x1 + I(x1^2) )
plot(y~x1)
curve(coef(modelo)[1]+coef(modelo)[2]*x+coef(modelo)[3]*x^2,add=TRUE,col = "blue",lwd=2)

Modelo não linear

2

É possível também trabalhar com gráfico interativo, caso necessário. Para tal, fica a dica do pacote highcharter.

set.seed(1)
x <- rnorm(100)
y <- rnorm(100) + 2*x +10
dados <- data.frame(x = x, y = y)

library(highcharter)
library(ggplot2)
library(dplyr)

highcharter::hchart(
  dados,   # data frame
  "point", # estilo do plot
  highcharter::hcaes(
    x = x, # columa de dados para eixo x
    y = y # columa de dados para eixo y
    ),
  regression = TRUE) %>% 
  highcharter::hc_colors("#d35400") %>% 
  highcharter::hc_add_dependency("plugins/highcharts-regression.js")

inserir a descrição da imagem aqui

Você deve fazer log-in para responder a esta pergunta.

Esta não é a resposta que você está procurando? Pesquise outras perguntas com a tag .