连老师:
您好!
::高级部分::
* 计量分析与Stata应用
* ==========================
* 第八讲 STATA编程进阶
* ==========================
* 8.3 输入项
* 第二部分
第544-565行
cap program drop mysepreg3
program define mysepreg3
version 9.2
syntax varlist [if] [in]
gettoken y x: varlist
tokenize "`x'"
local i = 1
while "``i''" != ""{
dis _n _n in y "Regression for " in w "``i''"
reg `y' ``i'' `if' `in' /*!!注意此处的变化*/
est store r_``i++''
}
dis _n _n in w "Summary of the Results: "
est table r_*, stat(r2) star(0.01 0.05 0.1)
end
sysuse auto, clear
mysepreg3 price wei len turn foreign if price<10000
mysepreg3 price wei len in 20/60
第564行mysepreg3 price wei len turn foreign if price<10000
Stata结果显示
------------------------------------------------------------------------------
Variable | r_weight r_length r_turn r_foreign
-------------+----------------------------------------------------------------
weight | .45870576*
length | 18.249028**
turn | 18.013473
foreign | 894.92273**
_cons | 3846.8138*** 1789.4139 4454.7079** 4878.9773***
-------------+----------------------------------------------------------------
r2 | .04943905 .07727094 .00283984 .08756693
------------------------------------------------------------------------------
legend: * p<.1; ** p<.05; *** p<.01
这个好理解。
第565行 mysepreg3 price wei len in 20/60
Stata结果显示
------------------------------------------------------------------------------
Variable | r_weight r_length r_turn r_foreign
-------------+----------------------------------------------------------------
weight | 1.7157901***
length | 50.831656***
turn | 18.013473
foreign | 894.92273**
_cons | 659.07758 -3749.7836 4454.7079** 4878.9773***
-------------+----------------------------------------------------------------
r2 | .26003544 .18006048 .00283984 .08756693
------------------------------------------------------------------------------
legend: * p<.1; ** p<.05; *** p<.01
问题:
根据第565行 mysepreg3 price wei len in 20/60的命令,以及mysepreg3(自编程序),应该是price分别对weight和length作回归,两次回归两个结果,为什么却有四个结果(包括 r_turn 和 r_foreign)?
谢谢连老师!


雷达卡





京公网安备 11010802022788号







