我试图创建一个基本的plm回归。问题是我试图同时运行数千个不同的回归(美国每个县一个回归),所以我无法自定义每个回归。 我想要运行的命令如下所示:
pm1 <- plm(y~x+unemploymentRate+rnetmig+percentageAbove65+xmas+xmas1+newyears+valentines+easter+thanksgiving+mothersday+halloween+superbowl+julyfourth,index=c("upc","week_end"),data=counties[[c]])
从圣诞到七月的变量都是假人,其他一切都不是。 以下适用于所有县:
pm2 <- plm(y~x+unemploymentRate+rnetmig+percentageAbove65,index=c("upc","week_end"),data=counties[[c]])
但是,一旦我添加了假人,在第一个plm中,我会收到此错误或类似错误:
series fips_state_code, xfips_county_code, xupc_ver_uc, xSTNAME, xCTYNAME, xxmas, xnewyears, xvalentines, xeaster, xsuperbowl, xjulyfourth are constants and have been removed
Error in eval(expr, envir, enclos) : object 'xmas' not found
我意识到这是因为每个傻瓜在某些县都有所不同。在这种情况下,圣诞节始终为0,因此导致回归崩溃。 有没有办法指示plm只包含那些非常量的回归量,而不是收到这个错误?我不能手动完成每个回归,因为有数千个回归。我也无法合理地创造&#34;如果&#34;每个虚拟变量的语句,因为我很快就会添加更多。 有谁知道会为我做这个命令吗? 非常感谢。