quantreg~2000 の観測値のデータベースがあり、パッケージを使用して 95 パーセンタイルで分位点回帰を行いました。
さらに分析を行うために、95 パーセンタイル回帰の勾配と切片の計算に実際に使用された観測値を特定したいと考えました。それを行う方法はありますか？
これは私がこれまでに使用したコードですquantreg：
datos<-quantreg.example
library(quantreg)
rq(y ~ x, tau=0.95, data=datos, method="br", model = TRUE)
ここにデータファイルがあります: http://www.filedropper.com/quantregexample
わかりました、私はあなたの質問が求めていると思うものに対する簡単な解決策を考え出そうと思いました:
set.seed(123)
library(dplyr) #data transformation
library(quantreg) #quantile regression
#make dummy data
df <- data.frame(x = sample(1:10, 200, replace = T))
df$y <- df$x + rnorm(200)
#fit quantile regression
my_q <- rq(y~x, data = df, tau = 0.95)
#use dplyr to get 95% quantile at each x
df_q <- df %>% group_by(x) %>% summarise(yq = quantile(y, probs = .95))
#quick viz with red points being 95% quantiles
with(df, plot(x,y))
legend('topleft',legend = '95% Conditional Quantiles',col = 'red',pch = 19, bty = 'n')
with(df_q, points(x, yq, col = 'red', pch = 19))
abline(reg = my_q)
お役に立てれば。