5

Meus dados são trazidos de uma base e eu simplesmente aplico:

ggplot (data = df_postgres, aes (x = date_time, y = duracao)) + geom_point() + stat_smooth(method = "loess") 

E obtenho a seguinte suavização:

Suavização obtida

Existe a possibilidade de obter os pontos máximos e mínimos da curva em azul?

1
  • Isso mesmo, pessoal. Das duas formas eu consegui obter o resultado que esperava. Obrigada! :)
    – Samara
    15/09/2015 às 15:14

2 Respostas 2

2

O ggplot utiliza a função loess do base R. Você pode rodar o ajuste fora do comando gráfico e obter os valores usando predict:

set.seed(4)
dados <- data.frame(date_time = 1:20, duracao = rnorm(20))
dados.loess <- loess(duracao ~ date_time, data = dados)

summary(predict(dados.loess))
#    Min.  1st Qu.   Median     Mean  3rd Qu.     Max. 
#-0.06347  0.15680  0.34110  0.35310  0.56030  0.74970

Confirmando o resultado, podemos adicionar as linhas no gráfico:

library(ggplot2)
ggplot(data = dados, aes (x = date_time, y = duracao)) +
  geom_point() + stat_smooth(method = "loess") +
  geom_hline(aes(yintercept = min(predict(dados.loess)))) +
  geom_hline(aes(yintercept = max(predict(dados.loess))))

inserir a descrição da imagem aqui

6
  • Você sabe me dizer se tem como achar, também, para cada mudança da curva, o valor máximo e mínimo?
    – Samara
    16/09/2015 às 11:45
  • Você diz as mudanças de sinal da inclinação? Se sim, você vai precisar da derivada dos dados, que pode ser obtida com diff. Isso provavelmente seria mais apropriado em uma nova pergunta, mas dê uma olhada neste código. Use o da resposta e mude apenas onde tem 20 para 1000. y <- predict(dados.loess); plot(y, type = "l"); m <- abs(diff(y, 1)); match(sort(m), m)[1:10]; abline(v=c(6953, 2621)).
    – Molx
    16/09/2015 às 13:41
  • Isso, cada mudança no sinal de inclinação, gostaria também de obter os pontos máximos e os pontos mínimos. Mas, nesse caso ele pega os máximos e mínimos dos pontos e não da curva.. Queria da curva que é gerada pelo LOESS também.
    – Samara
    16/09/2015 às 14:32
  • @Samara Sugiro que você faça uma nova pergunta para esse problema e inclua uma exemplo dos seus dados para facilitar o trabalho de quem for responder.
    – Molx
    16/09/2015 às 14:47
  • Está bem, obrigada.
    – Samara
    16/09/2015 às 15:46
2

Quando você usa a função stat_smooth do ggplot2, ele realiza os seguintes passos:

Vou usar o banco de dados mtcars, mas você pode facilmente substituir pelo seu próprio.

modelo <- loess(mpg ~ hp, data = mtcars) # ajusta o modelo  loess
mtcars$pred <- predict(modelo) # calcula as predicoes

ggplot (data = mtcars, aes (x = hp, y = mpg)) + geom_point() + 
 geom_line(aes(y = pred), colour = "blue") # plota as predicoes usando linha

Veja que este gráfico é igual ao gráfico feito usando:

ggplot (data = mtcars, aes (x = hp, y = mpg)) + geom_point() + 
 stat_smooth(method = "loess")

Portanto, para pegar o máximo valor da curva basta filtrar no seu data.frame:

mtcars[mtcars$pred == max(mtcars$pred), ]

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 .