R量子模型不再现分位数:为什么?
我正在使用R量子模型不再现分位数:为什么?,r,quantile,quantreg,R,Quantile,Quantreg,我正在使用quantreg包来预测分位数及其置信区间。我不明白为什么预测的分位数与使用quantile()直接从数据计算的分位数不同 库(tidyverse) 图书馆(量子力学) 数据您应该意识到,预测只有10个项目的数据集的第90个百分位数实际上完全是基于两个最高值。您应该查看分位数的帮助页面,在那里可以找到该术语的多个定义 当我运行此命令时,我看到: yqr1<- predict(qr1, newdata=tibble(data=c(1)) ) yqr1 1 8.
quantreg
包来预测分位数及其置信区间。我不明白为什么预测的分位数与使用quantile()
直接从数据计算的分位数不同
库(tidyverse)
图书馆(量子力学)
数据您应该意识到,预测只有10个项目的数据集的第90个百分位数实际上完全是基于两个最高值。您应该查看分位数的帮助页面,在那里可以找到该术语的多个定义
当我运行此命令时,我看到:
yqr1<- predict(qr1, newdata=tibble(data=c(1)) )
yqr1
1
8.525812
杰出的我忘了quantile()
有type
参数。当我设置type=1
时,我从两种方法中得到了相同的答案,这对于我的目的来说已经足够了。
yqr1<- predict(qr1, newdata=tibble(data=c(1)) )
yqr1
1
8.525812
data
# A tibble: 10 x 1
data
<dbl>
1 8.52581158
2 7.73959380
3 4.53000680
4 0.03431813
5 2.13842058
6 5.60713159
7 6.17525537
8 8.76262959
9 5.30750304
10 4.61817190
quantile(data$data, .9)
# 90%
#8.549493
?quantile