clear
cd D:\teaching\stata\multi_reg		/* ワーキングディレクトリの指定 */

set logtype text
log using wage2.txt, replace		/* logファイル名の指定 */

import delimited wage2.csv, clear	/* データファイルの読込 */

gen lwage = ln(wage)	/* 変数"ln(wage)"を発生させ"lwage"と定義 */
gen agesq = age^2		/* 変数"age^2"を発生させ"agesq"と定義 */
gen em = educ*married	/* 変数"educ*married"を発生させ"em"と定義 */

reg lwage educ age agesq			/* 最小2乗推定:1(i) */
estimates store reg1				/* 結果を"reg1"の名でメモリに記憶 */

reg lwage educ age agesq married	/* 最小2乗推定:1(ii) */
estimates store reg2				/* 結果を"reg2"の名でメモリに記憶 */

reg lwage educ age agesq married em	/* 最小2乗推定:1(iii) */
estimates store reg3				/* 結果を"reg3"の名でメモリに記憶 */

estimates table reg1 reg2 reg3, b(%9.4f) se stats(N r2 r2_a) stfmt(%9.4f) style(noline)              /* 推定値・標準誤差・サンプルサイズ等を一覧表形式で表示 */

test (married=0)(em=0)	/* 仮説検定:2 */

log close