Достижение аналогичных результатов с использованием Quantreg, как это было сделано с лессом - PullRequest
0 голосов
/ 03 июня 2019

Я пытаюсь рассчитать процентили для массы жира и мышечной массы.Набор данных имеет ограниченный размер выборки.Использование функции LOESS обеспечивает отличное представление данных.Однако LOESS не предоставляет квантили, такие как 10-й, 25-й, 50-й, 75-й или 90-й процентиль.Тем не менее, как я могу достичь с Quantreg Tau = 0,5 аналогичный результат по сравнению с LOESS?

data(cars)
lmodel <- loess(dist ~ speed, data = cars, span = 0.9 )
plot(cars)
lines( x=4:25 , y=predict(lmodel, newdata= data.frame(speed=4:25)) ,col="Blue")

enter image description here

Как я мог изменить код с помощью Quantregполучить аналогичные результаты / выход по сравнению с потерями?Можно ли изменить зигзагообразную линию на плавную линию, похожую на линию LOESS?

library(quantreg)
data(cars)
xx <- seq(min(cars$speed), max(cars$speed), length.out = nrow(cars))
taus <- c( 0.5)
rqmodel <- rq(dist ~ speed, tau = taus, data = cars)
rq_yy <- predict(rqmodel)
plot(cars)
matlines(xx, rq_yy, col = "black",lwd=3)

enter image description here

...