|
刚才我在用rqpd回归时也遇到了这个问题。
rqpd(formula, panel = panel(), data = parent.frame(),
na.action, subset, contrasts = NULL, control = NULL, ...)
tmpmax是control里面的一个参数,解释为upper bound of the working array in Cholesky factorization;computed automatically inside the routine if missing
这些内容是在动态面板分位数回归的帮助文件里找到的,这个问题的解决方法是在
http://blog.sina.com.cn/s/blog_65b2919d0102vhtq.html
这个网址内找到的
然后我把代码改为
fe.fit <- rqpd(fe.form, p, data=paneldata, control=list(tmpmax=1000))
就可以了
另外lambda就是λ,是Koenker的面板惩罚分位数回归公式中的调节参数,用来调节个体效应。具体请见论文“面板数据分位数回归模型求解及应用研究”,我是初学者,这些问题理解的不透彻。
|