|
728 | 728 | "cell_type": "markdown", |
729 | 729 | "metadata": { |
730 | 730 | "heading_collapsed": true, |
731 | | - "hidden": true |
| 731 | + "hidden": true, |
| 732 | + "jp-MarkdownHeadingCollapsed": true |
732 | 733 | }, |
733 | 734 | "source": [ |
734 | 735 | "### 説明" |
|
877 | 878 | { |
878 | 879 | "cell_type": "markdown", |
879 | 880 | "metadata": { |
880 | | - "hidden": true |
| 881 | + "hidden": true, |
| 882 | + "jp-MarkdownHeadingCollapsed": true |
881 | 883 | }, |
882 | 884 | "source": [ |
883 | 885 | "### 推定" |
|
893 | 895 | "\n", |
894 | 896 | "(注意)\n", |
895 | 897 | "\n", |
896 | | - "時間ダミー変数の代わりに`TimeEffects`を使わないように。入れることができますが,そのような仕様になっていません。" |
| 898 | + "時間ダミー変数の代わりに`TimeEffects`を使わないように。\n", |
| 899 | + "入れてもエラーは出ませんが,`TimeEffects`は無視されて,時間ダミーがない結果と同じになります。" |
897 | 900 | ] |
898 | 901 | }, |
899 | 902 | { |
|
906 | 909 | "source": [ |
907 | 910 | "formula_re = 'lwage ~ 1 + married + union + expersq \\\n", |
908 | 911 | " + exper + educ + black + hisp \\\n", |
909 | | - " +d81+d82+d83+d84+d85+d86+d87'" |
| 912 | + " + d81 + d82 + d83 +d84 + d85 + d86 + d87'" |
910 | 913 | ] |
911 | 914 | }, |
912 | 915 | { |
|
1155 | 1158 | "source": [ |
1156 | 1159 | "formula_cre = 'lwage ~ 1 + married + union + expersq \\\n", |
1157 | 1160 | " + married_mean + union_mean + expersq_mean \\\n", |
1158 | | - " +d81+d82+d83+d84+d85+d86+d87'\n", |
| 1161 | + " + d81 + d82 + d83 + d84 + d85 + d86 + d87'\n", |
1159 | 1162 | "\n", |
1160 | 1163 | "result_cre = RandomEffects.from_formula(formula_cre, data=wagepan).fit()\n", |
1161 | 1164 | "\n", |
|
1199 | 1202 | "source": [ |
1200 | 1203 | "FEモデルとREモデルのどちらが適しているかを調べることができるHausman検定というものがある。CREモデルを使うことにより,同様の検定が簡単に行える。(式6)を考えよう。\n", |
1201 | 1204 | "* $\\gamma=0$の場合,REモデルの(式5)と同じになり,REモデルが妥当ということになる。\n", |
1202 | | - "* $\\gamma\\neq 0$の場合,$a_i$と$\\bar{x}_{i}$は相関することになり,これは$\\text{Cov}\\left(a_i,{x}_{it}\\right)\\neq 0$を意味し,REモデルが妥当な推定方法となる。\n", |
| 1205 | + "* $\\gamma\\neq 0$の場合,$a_i$と$\\bar{x}_{i}$は相関することになり,これは$\\text{Cov}\\left(a_i,{x}_{it}\\right)\\neq 0$を意味し,FEモデルが妥当な推定方法となる。\n", |
1203 | 1206 | "\n", |
1204 | 1207 | "この考えを利用して,次の帰無仮説と対立仮説のもとで$\\gamma$の優位性を調べる。\n", |
1205 | 1208 | "* $\\text{H}_0:\\;\\text{Cov}\\left(a_i,x_{it}\\right)=0$\n", |
|
1293 | 1296 | "formula_cre2 = 'lwage ~ 1 + married + union + expersq \\\n", |
1294 | 1297 | " + exper + educ + black + hisp \\\n", |
1295 | 1298 | " + married_mean + union_mean + expersq_mean \\\n", |
1296 | | - " +d81+d82+d83+d84+d85+d86+d87'\n", |
| 1299 | + " + d81 + d82 + d83 + d84 + d85 + d86 + d87'\n", |
1297 | 1300 | "\n", |
1298 | 1301 | "result_cre2 = RandomEffects.from_formula(formula_cre2, data=wagepan).fit()\n", |
1299 | 1302 | "\n", |
|
2378 | 2381 | "metadata": { |
2379 | 2382 | "celltoolbar": "Tags", |
2380 | 2383 | "kernelspec": { |
2381 | | - "display_name": "Python 3 (ipykernel)", |
| 2384 | + "display_name": "py11", |
2382 | 2385 | "language": "python", |
2383 | | - "name": "python3" |
| 2386 | + "name": "py11" |
2384 | 2387 | }, |
2385 | 2388 | "language_info": { |
2386 | 2389 | "codemirror_mode": { |
|
2392 | 2395 | "name": "python", |
2393 | 2396 | "nbconvert_exporter": "python", |
2394 | 2397 | "pygments_lexer": "ipython3", |
2395 | | - "version": "3.12.7" |
| 2398 | + "version": "3.11.9" |
2396 | 2399 | }, |
2397 | 2400 | "toc": { |
2398 | 2401 | "base_numbering": 1, |
|
0 commit comments