|
do文件内容如下:请各位,问题出在哪里?
set seed 10101
quietly set obs 500
gen x2=rnormal(0,5)
gen x3=rnormal(0,5)
gen e=rnormal(0,5)
gen u=sqrt(exp(-1+0.2*x2))*e
gen y=1+1*x2+1*x3+u
sum
reg y x2 x3
predict uhat,resid
gen absu=abs(uhat)
quietly twoway (scatter absu x2)(lowess absu x2,bw(0.4) lw(thick)),scale(1.2) xscale(titleg(*5)) yscale(titleg(*5))plotr(style(none))
graph rename gls1
quietly twoway (scatter absu x3)(lowess absu x3,bw(0.4) lw(thick)),scale(1.2) xscale(titleg(*5)) yscale(titleg(*5))plotr(style(none))
graph rename gls2
graph combine gls1 gls2
|