我通过应用 rq 估计了之前由 OLS 估计的线性回归模型的 0.1 条件分位数。原始数据由 8760 行组成(一年的每小时数据)。
> require("quantreg")
> Qr_0.1 = rq(Reg_OLS, tau=0.1)
我现在想使用估计的模型来预测更小的时间跨度:
> Qrpredict_0.1 = predict.rq(Qr_0.1, newdata)
但它会导致以下警告消息:
Warning message:'newdata' had 192 rows but variables found have 8760 rows
我得到的结果包含 8760 行。似乎它只是使用旧数据。
这里到底是什么问题,有没有办法解决这个问题?
编辑:当我使用具有相同行数的新数据时,它工作得很好。