/* to illustrate impact of China */

/* note the code in each line excludes countries with most recent marriage data from before 2000 (lstyr<2000),
   and Jamaica, with 28% marriage (lstmar<40) */

twoway (scatter gii lstmar if lstmar>40 & lstyr>1999, msymbol(o) mlc(white) mfc(black) msize(medium) mlw(.1)                     ///
        scheme(s1mono) xlab(40(10)100) legend(off)                                                                               ///
        xti("Women ages 30-34 married (%)") yti("Gender inequality") scale(.8) plotregion(style(none)) )		                 ///
       (scatter gii lstmar if lstmar>40 & lstyr>1999 & pop<1000000 & pop>48999, mlabposition(0) msymbol(i) mlabel(country) mlabs(small)  ///
        mlabp(6) mlabc(red) xlab(40(10)100) )                                                                                    ///
       (scatter gii lstmar if lstmar>40 & lstmar<75 & gii>.42 & lstyr>1999 & country != "South Africa", mlabposition(0) msymbol(i) mlabel(country) mlabs(small)       ///
        mlabp(12) mlabc(blue) xlab(40(10)100) )                                                                                  ///
       (scatter gii lstmar if gii>.71 & lstyr>1999, mlabposition(0) msymbol(i) mlabel(country) mlabs(small)                      ///
        mlabp(12) mlabc(blue) xlab(40(10)100) )                                                                                  ///
       (scatter gii lstmar if gii<.08 & lstmar>75 & lstyr>1999, mlabposition(0) msymbol(i) mlabel(country) mlabs(small)          ///
        mlabp(6) mlabc(green) xlab(40(10)100) )                                                                                  ///
       (scatter gii lstmar if lstmar<50 & lstmar>40 & lstyr>1999, mlabposition(0) msymbol(i) mlabel(country) mlabs(small)        ///
        mlabp(6) mlabc(green) xlab(40(10)100) )                                                                                  ///
       (qfit gii lstmar if lstmar>40 & lstyr>1999, mlabposition(0) xlab(40(10)100) lp(shortdash) lc(gray)) ///
	   (scatter gii lstmar if lstmar>40 & lstyr>1999 & pop>999999, msymbol(i) mlabel(country) mlabs(large)  mlabp(12) mlabc(red) xlab(40(10)100) ///
	   title("Marriage and gender inequality in 124 countries", size(medium)) note("{it:Chart by PN Cohen, from U.N. statistics (marriage most recent before 2008, gender from 2015 report)}") ///
	   		t2title("{it:Quadratic fit unweighted; countries over 49 million population in red}", size(medsmall)))

twoway (scatter gii lstmar if lstmar>40 & lstyr>1999 [w=pop], msymbol(o) mlc(black) mfc(none) msize(medium) mlw(.1)                  ///
        scheme(s1mono) xlab(40(10)100) ylab(0(.2).8) legend(off)                                                                     ///
        xti("Women ages 30-34 married (%)") yti("Gender inequality") scale(.8) plotregion(style(none)) )		                     ///
       (qfit gii lstmar if lstmar>40 & lstyr>1999 & country !="China" [w=pop], mlabposition(0) xlab(40(10)100) lp(shortdash))        ///
       (qfit gii lstmar if lstmar>40 & lstyr>1999                     [w=pop], mlabposition(0) xlab(40(10)100) lp(shortdash) lc(red) ///
        title("Marriage and gender inequality in 124 countries", size(medium))                                                       ///
		t2title("{it:Quadratic fit weighted by population: Gray line excludes China, red line includes China}", size(medsmall))      ///
		note("{it:Chart by PN Cohen, from U.N. statistics (marriage most recent before 2008, gender from 2015 report)}"))