我估计了先前由OLS通过应用rq估计的线性回归模型的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行。它似乎只是使用旧数据。
这里到底出了什么问题,是否有办法解决这个问题?
编辑:当我使用具有相同行数的新数据时,它完全正常。