clear
set obs 1
generate var1 = 66.1 in 1
set obs 2
replace var1 = 67.3 in 2
set obs 3
replace var1 = 68.6 in 3
set obs 4
replace var1 = 55.1 in 4
set obs 5
replace var1 = 65.8 in 5
set obs 6
replace var1 = 52.4 in 6
generate var2 = 70 in 1
replace var2 = 73.7 in 2
replace var2 = 63.4 in 3
replace var2 = 85.1 in 4
replace var2 = 66.6 in 5
replace var2 = 35.4 in 6
rename var1 gs
rename var2 gm
gen name=""
replace name="LAT" if _n==1
replace name="NYT" if _n==2
replace name="USAT" if _n==3
replace name="WSJ" if _n==4
replace name="WP" if _n==5
replace name="WT" if _n==6

twoway (scatter gm gs, xtitle("Gentzkow-Shapiro Slant Quotient estimates") ylabel(30(10)90) ymtick(35(10)85) xlabel(50(5)70) text(90 65.5 "Data from Groseclose (2011: 175)",size(*.8)) xmtick(52.5(5)67.5) ytitle("Groseclose-Milyo Slant Quotient estimates") mlabel(name) legend(order(2 3) label(2 "Excludes Wall Street Journal") size(*.8) label(3 "Excludes Washington Times"))) (lfit gm gs if name!="WSJ",lcolor(black) lpattern(solid)) (lfit gm gs if name!="WT", scheme(s1mono) lcolor(gs10) lpattern(dash))

graph export "gmgs.eps",replace

pwcorr gm gs 
pwcorr gm gs if name!="WSJ"
pwcorr gm gs if name!="WT"