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 non-positive 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 回答