use ART.dta, replace
// Fix missing values
recode avgemtrans (-1 -3 = .)
recode emhum (-1 = .)
// Calculate r-squared
qui regress avgemtrans emhum
local r2 = round(e(r2), .01) // Round coeff.
// Create first panel
twoway (scatter avgemtrans emhum, mlab(countrycode) mlabpos(0) msymbol(i)) ///
(lfit avgemtrans emhum) ///
, xlabel(, format(%6.1f)) ylabel(, format(%6.1f)) ///
ytitle("Avg. no. of fresh non-donor embryos" "per fresh embryo transfer cycle") ///
xtitle("Avg. level disagreement that" "embryo is a human being") ///
legend(order(2 "Linear fit R{char 178} = `r2'") ring(0) pos(8)) ///
name(figure2a, replace)
// Fix missing values
recode afford clonehelpinf (-1 = .)
// Listwise deletion to control y-axis range
preserve
keep if !missing(afford, clonehelpinf)
// Calculate r-squared
regress afford clonehelpinf
local r2 = round(e(r2), .01) // Round coeff.
// Create second panel
twoway (scatter afford clonehelpinf, mlab(countrycode) mlabpos(0) msymbol(i)) ///
(lfit afford clonehelpinf) ///
, xlabel(, format(%6.1f)) ylabel(0 (5) 25, format(%6.1f)) ///
ytitle("Affordability (net cost of a fresh ART cycle" "as % of annual disposable income)") ///
xtitle("Avg. level disagreement with" "cloning to help infertile couples") ///
legend(order(2 "Linear fit R{char 178} = `r2'") ring(0) pos(8)) ///
name(figure2b, replace)
restore
// Create final Figure
graph combine figure2a figure2b, row(1)
