<< coding stata | Main | 15° F & Real Feel -1° F >>

2005年11月17日

panel by R

Rでpanelができないかということだけれど,昨晩,Rのhelpをいじっていて,方法発見:

例えば,
state_base <- paste("state", 1:51, sep="")
formula <- as.formula(paste("y ~ ", paste(state_base, collapse= "+")))
とすると,formulaに,state FEが入ります(もちろん,このままだと,rank conditionで落ちるので,どれか一つダミー落とす)。2段階でちょっと面倒だけれど,これで,state trendとかも入れることができる。stateだと,51しかないからまだいいけれど,company FEとかで300とか500とかあった場合には,大いにコスト削減になるでしょう。
Reconometricsに加筆しておきます。

ソスソスソスeソスソス hatsuru : 2005年11月17日 08:16

ソスgソスソスソスbソスNソスoソスbソスN

ソスソスソスフエソスソスソスgソスソスソス[ソスフトソスソスソスbソスNソスoソスbソスNURL:
http://www.law.tohoku.ac.jp/~hatsuru/cgi-local/mt/mt-tb.cgi/494