* s0.do		1/06
* full model
* sreq 19 quaniles
* raw scale $1000
* 2001 only

tempfile t
set matsize 400

global X age age2 blk hsp asn scl hsg nhs mh fh sg ot rural

forvalues j=2/2 {
use cinc $X year if year==`j' using d0, clear
drop year
replace cinc=cinc/1000

* centering
sum $X
tokenize $X
  while "`1'"~="" {
    egen m=mean(`1')
    replace `1'=`1'-m
    drop m
    macro shift
}
sum $X

sqreg cinc $X, reps(500) q(.05 .10 .15 .20 .25 .30 .35 .40 .45 .50 .55 .60 .65 .70 .75 .80 .85 .90 .95)
outreg using s0`j', nolabel bdec(3) coefast se replace
mstore b`j', from(e(b))
mstore v`j', from(e(V))
keep age
keep if _n<11
save s0`j', replace
}