tenho o seguinte conjunto de dados no R:
x <- c(0.1, 3, 4, 5, 9, 12, 13, 19, 22, 25)
y <- c(5, 12, 17, 23, 28, 39, 26, 31, 38, 40)
bd <- data.frame(x, y)
Minha pergunta é como eu faço no R para gerar um modelo de regressão que melhor se ajuste de forma que 90% dos dados fiquem abaixo da linha de regressão e que o modelo estime a origem (zero).
Parece que o modelo geométrico se ajusta melhor a este caso. Tentei utilizar a exponencial conforme a seguir, ele começa na origem mais os dados não ficam 90% abaixo da curva.
library(ggplot2)
ggplot(bd,aes(x = x, y = y)) +
geom_point() +
stat_smooth(method = 'nls', formula = 'y~a*x^b',
method.args = list(start = list(a = 1, b = 1)),
se = FALSE)
fit0 <- quantreg::rq(y ~ 0 + x, tau = 0.90)
?ggplot(bd, aes(x = log(x), y = log(y))) + geom_point() + stat_smooth(method = 'quantreg::rq', formula = 'y ~ x', method.args = list(tau = 0.9), se = FALSE)