2

このquantregパッケージを使用して、R で分位点回帰を計算しています。次のコマンドを使用して QR を呼び出します。

quantGsReg15 <- rq(gsRMSD ~ kMeanGrp + medianDurationMS + flightHours + 
                   flightHoursType + landings30days + privPilot + afterKSS,
                   tau = .15 , data = accumulatedNM , method="fn")
summary(quantGsReg15)

戻り値を見ると、変数 flightHours と flightHoursType の係数はゼロ/NULL です。

tau: [1] 0.15

Coefficients:
                 coefficients   lower bd       upper bd      
(Intercept)        1.080000e-03  -1.381000e-02   1.618000e-02
kMeanGrp2          1.510000e-03 -1.797693e+308   1.510000e-03
kMeanGrp3          6.170000e-03  -2.630000e-03  1.797693e+308
medianDurationMS  -1.000000e-05  -4.000000e-05   3.000000e-05
flightHours        0.000000e+00   0.000000e+00   0.000000e+00
flightHoursType    0.000000e+00   0.000000e+00   1.000000e-05
landings30days     3.000000e-05  -6.600000e-04   3.600000e-04
privPilotyes       4.440000e-03 -1.797693e+308   1.380000e-02
afterKSS           3.000000e-05  -1.130000e-03   2.250000e-03

追加するse = "nid"と、警告メッセージも表示されます。

警告メッセージ: summary.rq(quantGsReg15, se = "nid") : 7 非正の fis

以前に計算した OLS 回帰は、そうではない可能性が高いことを示しています。さまざまな QR メソッド ( method = "fn") を試しましたが、成功しませんでした。

0.5 分位数に対して同じコマンドを呼び出すと、同様の出力が返されます。ただし、別の警告メッセージが表示されます。

警告メッセージ: rq.fit.br(x, y, tau = tau, ci = TRUE, ...) : 解が一意でない可能性があります

私は何を間違っていますか？私は今30分間グーグルで検索しましたが、成功しませんでした。

4

1 に答える 1

0

PS: さらにテストを行い、従属変数に 10000 を掛けて、同じ回帰を計算したところ、出力は正常に見えました。コンマの後ろの数桁を計算しているため、何らかの丸めが行われている可能性がありますか?

Call: rq(formula = gsRMSD10k ~ kMeanGrp + medianDurationMS + flightHours + 
    flightHoursType + landings30days + privPilot + afterKSS, 
    tau = 0.15, data = accumulatedNM, method = "fn")

tau: [1] 0.15

Coefficients:
                 Value     Std. Error t value   Pr(>|t|) 
(Intercept)       10.80968 101.20390    0.10681   0.91563
kMeanGrp2         15.09255  24.98586    0.60404   0.55021
kMeanGrp3         61.70173  28.31814    2.17888   0.03706
medianDurationMS  -0.06609   0.15616   -0.42323   0.67505
flightHours        0.00948   0.00183    5.19156   0.00001
flightHoursType    0.02038   0.00595    3.42483   0.00175
landings30days     0.28104   1.91373    0.14685   0.88420
privPilotyes      44.41189  28.79167    1.54253   0.13309
afterKSS           0.31738   6.38128    0.04974   0.96065
Warning message:
In summary.rq(quantGsReg15, se = "nid") : 5 non-positive fis
于 2014-08-02T09:45:12.437 に答える