clear
set more off
cap log close
// log using sj-graphs.log, replace
// cd D:\home\stephenj\MYProjects\cappellari\deprivation\mvprobit
version 8.1

* sj-graphs.do
* SPJ 20 August 2003
* Create figures for LC/SPJ SJ -mvprobit- article
* [reads in v7-format data; creates graphs in v8-format]

use rho21-gen, clear

twoway (connected rho21d draws if draws >= 1 & draws <= 150  , /// 
	sort msymbol(circle)  ), yscale(range(.1 .3)) yline(.266729) ///
	ylabel(none .1 .15 .2 .25 .3) ///
	xlabel(1 10(20)150) subtitle(Seed = 123456789) scheme(sj) name(a)
	
twoway (connected rho21s draws if draws >= 1 & draws <= 150 , /// 
	sort msymbol(circle) ), yscale(range(.1 .3)) yline(.266729) ///
	ylabel(none .1 .15 .2 .25 .3) ///
	xlabel(1 10(20)150) subtitle(Seed = 999) scheme(sj) name(b)

graph combine a b, xsize(5) ysize(2.6)

graph export ../../eps/mvprobit1.eps, replace

use rho21, clear

twoway (connected rho21d draws if draws >= 1 & draws <= 150 , /// 
	sort msymbol(circle) ), yscale(range(-.35 0)) yline(-.2696186) ///
	ylabel(none -.35 -.30 -.25 -.20 -.15 -.10 -.05 0) /// 
	xlabel(1 10(20)150) subtitle(Seed = 123456789) scheme(sj) name(c)

* NB	ylabel(-.35(.05)0) lead to "-1.30e-17" printed instead of "0"
* in this and next command

twoway (connected rho21s draws if draws >= 1 & draws <= 150 , /// 
	sort msymbol(circle) ), yscale(range(-.35 0)) yline(-.2696186) ///
	ylabel(none -.35 -.30 -.25 -.20 -.15 -.10 -.05 0) /// 
	xlabel(1 10(20)150) subtitle(Seed = 999) scheme(sj) name(d)

graph combine c d, xsize(5) ysize(2.6)

graph export ../../eps/mvprobit2.eps, replace

/// log close