基本的なplm回帰を作成しようとしています。問題は、何千もの異なる回帰 (米国の郡ごとに 1 つ) を一度に実行しようとしているため、それぞれをカスタマイズできないことです。実行したいコマンドは次のようになります。
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]])
xmas から julyfourth までの変数はすべてダミーであり、それ以外はすべてダミーではありません。以下は、すべての郡で問題なく実行されます。
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
これは、それぞれのダミーが一部の郡でのみ異なるためであることに気付きました。この場合、xmas は常に 0 だったため、回帰がクラッシュしました。このエラーを受け取るのではなく、定数ではないリグレッサーのみを含めるように plm に指示する方法はありますか? 何千もの回帰があるため、手動で回帰を行うことはできません。また、ダミー変数ごとに「if」ステートメントを合理的に作成することもできません。おそらく、すぐにさらに多くの変数を追加するからです。誰かが私のためにこれを行うコマンドを知っていますか? よろしくお願いします。