* s_m0.do	1/06
* matrix operation 
* 13 covariates + cons
* graphs for beta's (19 QR)
* 500 bootstrap se
* 2001 only

set matsize 400

* 13 covariate + cons
local k=14
* 19 quantiles times k
local k1=`k'*19

forvalues m=2/2 {
use s0`m', clear

qui mstore b`m'
qui mstore v`m'
mat ren b`m' b
mat ren v`m' v

* 95%ci
* dimension `k' x 1 

mat vv=vecdiag(v)
mat vv=vv'
svmat vv
mat drop vv
qui replace vv1=sqrt(vv1)
mkmat vv1 if _n<=`k1', mat(v)
drop vv1

mat b=b'
mat l=b-1.96*v
mat u=b+1.96*v

* 19 quantiles
mat q=(.05\.10\.15\.20\.25\.30\.35\.40\.45\.50\.55\.60\.65\.70\.75\.80\.85\.90\.95)

* reorganize matrix by variable

forvalues j=1/`k' {
forvalues i=1/19 {
local l=`k'*(`i'-1)+`j'
mat x`j'q`i'=q[`i',1],b[`l',1],l[`l',1],u[`l',1],v[`l',1]
}
}

forvalues j=1/`k' {
mat x`j'=x`j'q1
forvalues i=2/19 {
mat x`j'=x`j'\x`j'q`i'
}
* q b l u v
mat list x`j', format(%8.3f)
svmat x`j'

mat a1=x`j'[1...,2]
mat a2=x`j'[1...,5]
mat xx`j'=q,a1,a2
* q b v
mat list xx`j', format(%8.3f)
mat drop a1 a2 xx`j'
}

* graphs using the same scale for a concept
* age
forvalues j=1/1 {
line x`j'2 x`j'1 || line x`j'3 x`j'1 || line x`j'4 x`j'1, ylabel(0(1)3) ytitle(income $1000) xlabel(0(.1)1) xtitle(p) legend(off)
graph export x`m'`j'.ps, as(ps) logo(off) tmargin(0) lmargin(0) replace
}

* age2
forvalues j=2/2 {
line x`j'2 x`j'1 || line x`j'3 x`j'1 || line x`j'4 x`j'1, ylabel(-.03(.01)0) ytitle(income $1000) xlabel(0(.1)1) xtitle(p) legend(off)
graph export x`m'`j'.ps, as(ps) logo(off) tmargin(0) lmargin(0) replace
}

* eth
forvalues j=3/5 {
line x`j'2 x`j'1 || line x`j'3 x`j'1 || line x`j'4 x`j'1, ylabel(-20(5)25) ytitle(income $1000) xlabel(0(.1)1) xtitle(p) legend(off)
graph export x`m'`j'.ps, as(ps) logo(off) tmargin(0) lmargin(0) replace
}

* ced
forvalues j=6/8 {
line x`j'2 x`j'1 || line x`j'3 x`j'1 || line x`j'4 x`j'1, ylabel(-100(20)0) ytitle(income $1000) xlabel(0(.1)1) xtitle(p) legend(off)
graph export x`m'`j'.ps, as(ps) logo(off) tmargin(0) lmargin(0) replace
}

* htype
forvalues j=9/12 {
line x`j'2 x`j'1 || line x`j'3 x`j'1 || line x`j'4 x`j'1, ylabel(-80(10)0) ytitle(income $1000) xlabel(0(.1)1) xtitle(p) legend(off)
graph export x`m'`j'.ps, as(ps) logo(off) tmargin(0) lmargin(0) replace
}

* rural
forvalues j=13/13 {
line x`j'2 x`j'1 || line x`j'3 x`j'1 || line x`j'4 x`j'1, ylabel(-18(3)0) ytitle(income $1000) xlabel(0(.1)1) xtitle(p) legend(off)
graph export x`m'`j'.ps, as(ps) logo(off) tmargin(0) lmargin(0) replace
}

* constant (mean of all covariates)
forvalues j=14/14 {
line x`j'2 x`j'1 || line x`j'3 x`j'1 || line x`j'4 x`j'1, ylabel(10(20)110) ytitle(income $1000) xlabel(0(.1)1) xtitle(p) legend(off)
graph export x`m'`j'.ps, as(ps) logo(off) tmargin(0) lmargin(0) replace
}

drop x*
matrix drop _all
}