stata中用固定效应模型回归有虚拟变量时为什么就omitted了 有虚拟变量的stata模型回归命令
stata\u4e2d\u53ea\u6709\u4e00\u4e2a\u865a\u62df\u53d8\u91cf\u4e3a\u4ec0\u4e48\u56de\u5f52\u4e0d\u4e86\u662f\u56e0\u4e3a\u5171\u7ebf\u6027\u95ee\u9898
\u7ed3\u679c\u7684\u524d\u4e24\u884c\u8868\u793a\u6a21\u578b\u7684\u7c7b\u522b,LZ\u91c7\u7528\u7684\u4e3arandomeffect\u968f\u673a\u6a21\u578b,\u622a\u9762\u53d8\u91cf\uff1aprovince,\u6837\u672c\u6570\u76ee310.\u7fa4\u7ec4\u6570\u76ee31,\u4e5f\u5c31\u662f\u6bcf\u7ec410\u4e2a\u89c2\u6d4b\u503c.3-5\u884c\u8868\u793a\u6a21\u578b\u7684\u62df\u5408\u4f18\u5ea6,\u5206\u522b\u4e3awithin,between,overall,\u7ec4\u5185,\u7ec4\u95f4,\u603b\u4f53\u4e09\u4e2a\u5c42\u6b21.6-7\u884c\u8868\u793a\u9488\u5bf9\u53c2\u6570\u8054\u5408\u68c0\u9a8c\u7684waldchi2\u68c0\u9a8c\u548cPvalue,p=0.000\u8868\u793a\u53c2\u6570\u6574\u4f53\u4e0a\u7070\u5e38\u663e\u8457.8-10\u884c\u8868\u793a\u89e3\u91ca\u53d8\u91cf\u7684\u4f30\u8ba1\u6743\u91cd,\u622a\u8ddd,\u6807\u51c6\u5dee,Z\u7edf\u8ba1\u91cf,P\u503c\u53ca95%\u7f6e\u4fe1\u533a\u95f4.\u8fd9\u5757\u513f\u8ddf\u622a\u9762\u56de\u5f52\u7684\u4ea7\u51fa\u7ed3\u679c\u662f\u4e00\u6837\u7684,\u5173\u4e8e\u4f60\u7684\u89e3\u91ca\u53d8\u91cfbase\u7684\u6743\u91cd\u89e3\u91ca\u662f,\u5728\u5176\u4ed6\u591a\u6709\u6761\u4ef6\u90fd\u4e0d\u53d8\u7684\u60c5\u51b5\u4e0b,base\u6bcf\u589e\u52a0\u4e00\u5355\u4f4d,city\u4f1a\u589e\u52a00.0179\u5355\u4f4d,P\u503c0.000,\u7070\u5e38\u663e\u8457.\u6700\u540e\u4e09\u884c\u5206\u522b\u662f\u968f\u673a\u6548\u5e94\u6a21\u578b\u4e2d\u4e2a\u4f53\u6548\u5e94\u548c\u968f\u673a\u5e72\u6270\u9879\u7684\u65b9\u5dee\u4f30\u8ba1\u503c,\u5206\u522b\u4e3asigma_u,sigma_e.\u4ee5\u4e0a\u4e24\u8005\u4e4b\u95f4\u7684\u5173\u7cfbrho.\u9700\u8981\u6ce8\u610f\u7684\u662f\u4f60\u7684\u6a21\u578b\u62df\u5408\u5ea6\u4e0d\u9ad8,R\u65b9\u53ea\u670926%,\u5f53\u7136\u8fd9\u8981\u770b\u5177\u4f53\u662f\u54ea\u65b9\u9762\u7684\u7814\u7a76\u4ee5\u53ca\u540c\u65b9\u5411\u5176\u4ed6\u5b66\u8005\u7684\u62df\u5408\u7ed3\u679c,\u5982\u679c\u5927\u5bb6\u90fd\u572820\u591a,\u90a3\u5c31OK.
stata中用固定效应模型回归有虚拟变量时为什么就omitted了？不要使用 xtreg 命令做 FE 回归，因为它会自动忽略 dummies.
推荐的做法是，手动生成两个 dummies. 然后，用这两个 dummies 各自乘以 size 得到两个 interaction terms.
最后，把这两个 dummies，两个 interaction terms，以及 size 放到 reg 命令中回归。
Life is so easy then.
如果加地区虚拟变量是为了控制随地区变化的因素，那么被omitted了完全没关系，地区虚拟变量能控制的因素已经被固定效应控制了。因为它也是个随时间不变的因素，而所有这样的因素都会被固定效应控制。
但如果你关心虚拟变量本身的系数，那你就不能用固定效应模型。一个合适的选择是用随机效应，但随机效应的假设更严格，你需要用hausman检验对比它跟固定效应的结果是否存在显著差异，如果存在，那么你就需要在随机效应里添加更多随时间不变的因素。
加交互项还是跟原来一样，需要添加两个虚拟变量各自的交互项。加交互项后原来的虚拟变量一般都是依然需要添加的，不过如果是固定效应模型的话，加不加就无所谓了，反正会omitted
参见 zhuyuhao.com/doc/posts/ 16和17页。
应为在stata中，i.year 这种生成变量的方式只对与单一变量有效，而且在回归方程之中不能够有运算符号。 你可以试一下使用stata自带的自动生成交叉变量的命令，Interaction expansion，或者是使用 Data > Create or change data > Other variable...
绛旓細鍔犲叆璋冭妭鏁堝簲锛氬湪鏁版嵁涓鍔犱竴涓皟鑺傚彉閲忥紝骞朵护鍏朵笌鑷彉閲忕殑鍏崇郴鏄皟鑺傛晥搴旂殑娴嬮噺锛屽苟浠ゅ叾涓庡洜鍙橀噺鐨勫叧绯讳篃鏄皟鑺傛晥搴旂殑娴嬮噺銆傚湪STATA涓杩涜妯″瀷璁剧疆锛氭偍鍙互浣跨敤xtreg鍛戒护鏉ユā鎷熷弻鍚鍥哄畾鏁堝簲妯″瀷锛屽叾涓彲浠ラ氳繃fe閫夐」鍜宺e閫夐」鏉ヨ缃浐瀹氭晥搴斿拰闅忔満鏁堝簲銆傚悓鏍凤紝鍙互浣跨敤涓粙鍙橀噺鍜岃皟鑺傚彉閲忎綔涓鸿嚜鍙橀噺锛屾潵杩涜妯″瀷...
绛旓細闈㈡澘妯″瀷寮曞叆鍥哄畾鏃堕棿鏁堝簲stata鎿嶄綔鏂规硶锛歺i: xtreg y x1 x2 x3 i.year,fe 鍙屽悜鍥哄畾鏁堝簲锛屾棦鍙互鎺у埗骞村害鏁堝簲锛屽張鍙互鐢ㄥ浐瀹氭晥搴娑堥櫎閮ㄥ垎鍐呯敓鎬 xi: xtreg y x1 x2 x3 i.year LSDV娉 灏辨槸铏氭嫙鍙橀噺鏈灏忎簩涔鍥炲綊 鍙﹀锛屽缓璁敤鑱氱被绋冲仴鏍囧噯宸,杩欐槸瑙ｅ喅寮傛柟宸殑鑹嵂 xi: xtreg y x1 ...
绛旓細1銆丠ausman妫楠岋細鍦ㄦ墽琛屽浐瀹氭晥搴旀ā鍨嬶紙FE) 鍜岄殢鏈烘晥搴旀ā鍨嬶紙RE) 涔嬪墠锛屽彲浠ヤ娇鐢╤ausman鍛戒护鏉ヨ繘琛屾楠屻傝妫楠岀殑闆跺亣璁炬槸闅忔満鏁堝簲妯″瀷鏄竴鑷翠笖鏈夋晥鐨勶紝鍗充笉瀛樺湪鍐呯敓鎬ч棶棰樸傚鏋減鍊煎皬浜0.05锛屽垯鎷掔粷闆跺亣璁撅紝琛ㄧず瀛樺湪鍐呯敓鎬ч棶棰橈紝闇瑕浣跨敤鍥哄畾鏁堝簲妯″瀷銆2銆丏WH妫楠岋細鎵ц鍥炲綊鍚庯紝鍏堝瓨鍌ㄩ殢鏈烘晥搴斾及璁￠噺鍜...
绛旓細鍒朵綔鏂规硶1銆妯″瀷鏄惁鏈夌悊璁烘敮鎸侊紝鐜版湁鏂囩尞鏄惁鍋氳繃绫讳技鐮旂┒2銆侀潰鏉挎暟鎹ā鍨嬮夋嫨鍥哄畾鏁堝簲锛屽浐瀹氭晥搴旀湁鍗曞洜绱犲拰鍙屽洜绱犱箣鍒嗭紝璇曡瘯閲囩敤鍙屽悜鍥哄畾鏁堝簲锛3銆佹ā鍨嬪彲鑳藉瓨鍦ㄥ唴鐢熸ч棶棰橈紝瀵绘壘宸ュ叿鍙橀噺锛屽埄鐢ㄤ袱闃舵鏈灏忎簩涔樻垨GMM浼拌璇曡瘯锛堝寘鎷竴姝MM,浜屾GMM,杩唬GMM,杩炵画GMM锛4銆佷笉琛岀殑璇濇崲鎺у埗鍙橀噺鍚э紝5銆侀噸鏂...
绛旓細鏄鍥哄畾鏁堝簲锛宖e锛宖ixed effect锛屼絾鏄懡浠よ搴旇鏄 xtreg x y z, fe vec(cluster province)
绛旓細鍐嶇敤xtreg鍛戒护杩涜鍥哄畾鏁堝簲闈㈡澘鏁版嵁鍥炲綊锛屽悗鍔爁閫夐」銆傚緱鍒扮粨鏋滃悗锛岀敤vif鍛戒护妫楠屾柟宸啫鑳鍥犲瓙銆傚湪鍋氬洖褰掑垎鏋愭椂鍙戠幇涓涓棶棰橈紝鍥犲彉閲弝鏈夌己澶卞硷紝濡傛灉涓嶇敤dropy==.鍛戒护锛屽洖褰掑悗vif妫楠屽皬浜10锛屽鏋滈噰鐢╠ropy==.鍛戒护锛屽洖褰掑悗vif妫楠屽间竴涓嬪瓙韫﹀埌浜27锛屽湪杩欎袱绉嶅鐞嗘柟寮鐨勫洖褰涓渶鍚庢牱鏈暟鏄浉鍚岀殑锛屽洖褰...
绛旓細stata瀵归潰鏉挎暟鎹仛闇嶆柉鏇兼楠岋紝缁撴灉鏄剧ず鐢ㄥ浐瀹氭晥搴旀ā鍨锛屼絾瑙ｉ噴鍙橀噺涓嶆樉钁楋紝鐢∣LS鍋氭樉钁楋紝瑙ｅ喅锛氳豹鏂浖妫楠岀粨鏋滅殑P鍊煎皬浜0.01锛屽嵆鎷掔粷鍘熷亣璁撅紝琛ㄦ槑搴旇閲囩敤鍥哄畾鏁堝簲銆傝豹鏂浖妫楠岀殑缁撴灉鏄憡璇夊浐瀹氭晥搴斿拰闅忔満鏁堝簲鍦ㄧ郴鏁颁及璁′笂鍑虹幇浜嗘樉钁楀樊寮傦紝鍥犳鍥哄畾鏁堝簲姣旈殢鏈烘晥搴斿ソ浣嗘槸涓嶆槸璇撮殢鏈烘晥搴斿氨涓嶈兘鐢ㄦ湁浜涙椂鍊欎负浜嗗仛...
绛旓細鍏堢敤xtset璁惧畾闈㈡澘鏁版嵁 鐒跺悗鐢▁treg锛宖e鎿嶄綔灏卞彲浠ュ仛闈㈡澘鏁版嵁鍥哄畾鏁堝簲鍟 闈㈡澘鏁版嵁鍥炲綊鍒嗘瀽鎴戝緢鐔熸倝鐨
绛旓細1銆侀鍏堬紝浣跨敤Stata鐨import鍛戒护灏嗘暟鎹鍏tata銆2銆佸叾娆★紝浣跨敤Stata鐨剅eshape鍛戒护灏嗘暟鎹浆鎹㈡垚闈㈡澘鏁版嵁鏍煎紡锛屽嵆灏嗕笉鍚屾椂闂寸偣鐨勬暟鎹爢鍙犳垚涓鍒楋紝姣忎釜涓綋瀵瑰簲涓琛屻3銆佹渶鍚庯紝闈㈡澘鏁版嵁鍥哄畾鏁堝簲妯″瀷鐨勭粨鏋滃寘鎷浐瀹氭晥搴旂郴鏁板拰鍏朵粬鍥炲綊绯绘暟锛屽浐瀹氭晥搴旂郴鏁板弽鏄犱釜浣撻棿鐨勫樊寮傦紝鍏朵粬鍥炲綊绯绘暟鍙嶆槧鑷彉閲忓鍥犲彉閲忕殑褰卞搷銆
绛旓細- 闆跺亣璁炬槸闅忔満鏁堝簲妯″瀷鏄竴鑷翠笖鏈夋晥鐨勶紝鍗充笉瀛樺湪鍐呯敓鎬ч棶棰樸- 濡傛灉p鍊煎皬浜0.05锛屾嫆缁濋浂鍋囪锛岃〃鏄庡瓨鍦ㄥ唴鐢熸ч棶棰橈紝搴浣跨敤鍥哄畾鏁堝簲妯″瀷銆2. DWH妫楠岋細- 鎵ц鍥炲綊鍒嗘瀽鍚庯紝瀛樺偍闅忔満鏁堝簲浼拌閲忓拰鍥哄畾鏁堝簲浼拌閲忥紝浣跨敤`estimates store`鍛戒护銆- 鐢熸垚闅忔満鏁堝簲浼拌閲忓拰鍥哄畾鏁堝簲浼拌閲忕殑棰勬祴鍊硷紝鍒嗗埆浣跨敤`...
