Question

我正在使用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"时，我还会收到一条警告消息：

警告讯息：总结.rq（quantGsReg15，se =＆＃34; nid＆＃34;）：7个非正fis

我之前计算过的OLS回归表明很可能不是这种情况。我尝试了不同的QR方法（方法=＆＃34; fn＆＃34;）但没有成功。

为0.5分位数调用相同的命令会返回类似的输出。但是，我收到了一条不同的警告信息：

警告讯息：在rq.fit.br（x，y，tau = tau，ci = TRUE，......）：解决方案可能不是唯一的

我做错了什么？我现在用谷歌搜索30分钟但没有成功。

Answer 1

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

R使用quantreg分位数回归包返回零/ NULL系数

1 个答案: