[R-br] regressão qauntilica

geovane barbosa geovanecb em yahoo.com.br
Sexta Setembro 25 17:14:53 BRT 2015


 Olá pessoal tudo bem.
Estou trabalhando com a regressão quantilica. Abaixo segue um exemplo. Porém gostaria que a aproximação fosse por splines. Como posso fazer isso?



set.seed(1)
n <- 500x <- floor(runif(n, 18, 80))y <- 5*x - 0.05*x^2  + rnorm(n, mean=0, sd=10*sqrt(x))plot(x, y)taus <- 1:9/10
library(quantreg)modelo <- rq(y ~ x , tau=taus)
valoresPreditos <- predict(modelo)
dadosGraficos <- data.frame(valoresPre = as.numeric(valoresPreditos),                             tau = as.factor(rep(taus, each=n)),                            x = rep(x, times=length(taus)))
library(ggplot2)g <- ggplot(unique(dadosGraficos), aes(x=x, y=valoresPre, group=tau)) + theme_bw()g + geom_line(aes(linetype=tau))


Abaixo em anexo segue uma imagem do que realmente gostaria de visualizar. Como poderia fazer isso no R?
Obrigado 
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20150925/8745b357/attachment.html>
-------------- Próxima Parte ----------
Um anexo não-texto foi limpo...
Nome: ggplot.png
Tipo: image/png
Tamanho: 45236 bytes
Descrição: não disponível
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20150925/8745b357/attachment.png>
-------------- Próxima Parte ----------
Um anexo não-texto foi limpo...
Nome: quantregbin3.png
Tipo: image/png
Tamanho: 22338 bytes
Descrição: não disponível
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20150925/8745b357/attachment-0001.png>


Mais detalhes sobre a lista de discussão R-br