相关随机效应面板分位数回归的问题:rqpd

时间:2020-06-25 07:01:41

标签: panel quantile-regression

我有一项研究的108名参与者的Panel数据,历时7年。回归方程式使用“工资”作为因变量,并使用了受教育年限,年龄,经验以及行业状况(工作)变量,该变量是未观察到的变量,但会影响因变量。一年中在培训上花费的天数是另一个重要变量。没有类别添加到模型。所有这些都是连续变量,包括行业条件。

rqpd软件包可与以下R模型的固定效果模型配合使用:

fe.form <-薪水〜教育程度+年龄+经验+行业+培训| countid |

fixed.fit <-rqpd(fe.form,panel(lambda = 5),data = quantile7)

个人ID的取值范围为1001到1008,年份的时间ID的取值范围为1到7。联合声明是通过将个人ID声明为因子来定义的。

对于相关随机效应模型,规范为

cre.form <-工资〜教育程度+年龄+经验+行业+培训|县|行业+培训

random.fit <-rqpd(cre.form,panel(method =“ cre”,cre =“ ad”),data = quantile7,subset = quantile7 $ seqno%in%1:7)

此处,countid是个人ID和年份ID的串联。因此数字将变为10011、10012、10013,...,10017、10021、10022,...,10027、10031 ...这就是如何将rqpd中出生体重影响示例的momid3传递到模型中。

这是一个平衡的数据,因此使用cre =“ ad”。变量“ seqno”标识1到7年。

R程序因致命错误而崩溃或引发以下错误

切换(错误,调用提取时,存储空间不足(工作空间)...

数据太小还是我传递数据的方式有问题(也许未正确定义countid)。感谢您解决此错误的帮助。

0 个答案:

没有答案