rqpd パッケージを使用してパネルで分位点回帰を実行すると、係数の上限/下限を取得できますか?
以下はセットアップです（私はRの初心者です）：
s = as.factor(rep(1:30,rep(63,30)))
s = as.matrix(s)
fit = rqpd(formula = y ~ x | s, panel = panel(lambda = 1, taus=c(0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 0.95), tauw=c(0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 0.95)), data = pdata, control = list(tmpmax = 500))
capture.output(summary(fit), file = "QR.txt")
コードはうまく機能し、出力は通常の係数、t 値、p 値などを提供します。しかし、95% の有意性を得るために、SE の上限と下限を把握したいと考えています。
パッケージ quantreg の rq を使用して、次の出力をキャプチャすることで上限/下限を抽出できます。
fit = rq(y~z, tau = taeus)
plot.fit = summary(fit)
capture.output(plot.fit, file = "BOUNDS.txt")
ただし、これは現在見ているパネルデータには当てはまりません。皆さん、何か指針はありますか？