首先,需要加载 rq 包和数据集。下面以 mtcars 数据为例:

library(rq)
data(mtcars)

然后,用 rq 函数进行分位数回归拟合。下面以中位数为例,即 0.5 分位数:

fit <- rq(mpg ~ wt, data = mtcars, tau = 0.5)

其中,mpg 为因变量,wt 为自变量,tau 为分位数。

接下来,可以用 summary 函数查看拟合结果:

summary(fit)

其中,Estimate 为回归系数,Std. Error 为标准误差,t-value 为 t 值,Pr(>|t|) 为 p 值。

最后,可以用 rsq 函数计算拟合度 R 方:

rsq(fit)

这样就可以得到分位数回归拟合的 R 方值了。

R语言分位数回归 (rq) 拟合与R方计算

原文地址: https://www.cveoy.top/t/topic/osMl 著作权归作者所有。请勿转载和采集!

免费AI点我,无需注册和登录