パッケージrq()内の関数を使用して分位点回帰を実行しようとしています。quantreg私は自分のデータをセットアップし、モデルなどを作成しました (別のサイトの指示に従います)。下記参照; REとWTPCは連続変数でArmoredあり、以下で定義されます。
# Convert Treatment to 1/0 Binary factor where Armored=1
Armored<-as.numeric(Treatment=="A")
# Define variables
Y<-WTPC
X<-cbind(RE,Armored)
# Quantile Regression @ 25th
QR25<-rq(Y~X,tau=0.25)
summary(QR25)
summary()ただし、rq()モデルを呼び出すと、これが得られます
> summary(QR25)
Call: rq(formula = Y ~ X, tau = 0.25)
tau: [1] 0.25
Coefficients:
coefficients lower bd upper bd
(Intercept) 11.09456 9.34170 18.22058
XRE -1.81530 -3.46350 -0.78062
XArmored -3.68480 -14.29227 -0.74389
これは素晴らしい情報ですが、通常の線形回帰のように標準誤差、p 値などが必要です。ドキュメントと他の例を確認しましたが、これがこれを行う方法のようです。モデル ( QR25) を作成してから呼び出すsummary(QR25)と、エラーと p 値を含む回帰出力が得られます。
これは分位点回帰の最初の試みであるため、以前のステップが欠落しているか、何か他のものを指定する必要がある可能性があります。何かご意見は？これを RStudio バージョン 0.98.501 で実行しています。