使用quantreg获得与黄土一样的相似结果

时间:2019-06-03 18:55:42

标签: r percentile loess quantreg quantile-regression

我正在尝试计算脂肪质量和瘦肉质量的百分位数。数据集的样本量有限。使用LOESS功能可以很好地表示数据。但是,LOESS不提供第10、25、50、75或90%的分位数。

也许quantreg可能有用。但是,用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

与LOESS相比,如何使用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

0 个答案:

没有答案