Random graphs (35): Line plots

clear
import excel "art births across countries.xlsx", sheet("Data") firstrow case(lower) clear

ren a geo // Rename country ID

// Rename variables: perc_YEAR format
foreach v of varlist b-o {
   local x : variable label `v'
   rename `v' perc`x'
}

// Reshape data set to long format
reshape long perc, i(geo) j(year)

// Cross-sectional bar graph at latest point in time
by geo: egen lastobs = max(cond(!missing(perc), year, .))  // Identify year of latest observation
gen byte last = year == lastobs                            // ID variable of latest observation
// According to http://www.stata.com/support/faqs/data-management/first-and-last-occurrences/



graph hbar perc if year >= 2007 & last == 1, over(geo, sort(perc) descending) nofill ///
      ytitle("% ART live births among all births") ylabel(0(1)6) ///
      note("{it:Note:} Data refer to 2010 or latest available (no earlier than 2007)", span) ///
      /*title("Country differences in ART birth rates")*/ ///
      name(crosssec, replace)



// Create classical twoway line plot
twoway (line perc year if geo == "DK", cmissing(no)) ///
       (line perc year if geo == "SI", cmissing(no)) ///
       (line perc year if geo == "BE", cmissing(no)) ///
       (line perc year if geo == "SE", cmissing(no)) ///
       (line perc year if geo == "DE", cmissing(no)) ///
       (line perc year if geo == "FR", cmissing(no)) ///
       (line perc year if geo == "IT", cmissing(no)) ///
       (line perc year if geo == "UK", cmissing(no)) ///
       (line perc year if geo == "US", cmissing(no)) ///
     , xlabel(1997(1)2010, ang(h) labsize(small)) xtitle("") ///
       ylabel(0(1)6) ytitle("% ART live births among all births") ///
       legend(label(1 "DK") label(2 "SI") label(3 "BE") ///
              label(4 "SE") label(5 "DE") label(6 "FR") ///
              label(7 "IT") label(8 "UK") label(9 "US") ///
              pos(2))

//cmissing(no) makes sure that gaps are not filled

 


// Use panel line plot functionality
keep if inlist(geo, "DK", "SI", "BE", "SE", "DE", "FR", "IT", "UK", "US")

encode geo, gen(country)
xtset country year, yearly
xtline perc, overlay xlabel(1997(1)2010, ang(h) labsize(small)) xtitle("") ///
      ylabel(0(1)6) ytitle("% ART live births among all births") ///

             legend(pos(2))
// Note that the gaps are much harder to spot in this plot

// Small multiples version
twoway (line perc year, by(geo, note("") /*title("Development of ART birth rate, selected countries")*/) cmissing(no)), ///
        xlabel(1998 (2) 2010, ang(h) labsize(small) alternate) xtitle("") ///
        ylabel(0(1)6) yscale(r(0)) ytitle("% ART live births among all births") ///
        name(longit, replace)
        // yscale(r(0)) adds missing gridline according to
        // http://www.stata.com/statalist/archive/2013-04/msg00904.html