Random graphs (65): Dot plots with confidence intervals

preserve
// Keep those countries for which three waves are available
keep if inlist(cntry, "AT", "AU", "BG", "CZ", "DEE", "DEW", "ES", "GB", "HU") | ///
        inlist(cntry, "IE", "IL", "JP", "NL", "NO", "PH", "PL", "RU", "SE")   | ///
  inlist(cntry, "SI", "US")

// Define temporary objects
tempname foo
tempname foox
postfile `foo' str3 cntry year perc perc_ll perc_ul using `foox', replace

levelsof cntry, local(levels1)
levelsof year, local(levels2)

// Calculate proportions by country and year
foreach year of local levels2 {
 foreach country of local levels1 {

    capture proportion separate if cntry == "`country'" & year == `year'
    *matrix list r(table)
    matrix fcoefs  = r(table)
    local fperc    = fcoefs[1,2] * 100
    local fperc_ll = fcoefs[5,2] * 100
    local fperc_ul = fcoefs[6,2] * 100

    di "`country'"  _skip(2) `year' _skip(2)  `fperc_ll' _skip(2) `fperc' _skip(2) `fperc_ul'

    if "`fperc'" != "" {    // Make sure that the loop doesn't break if empty
    post `foo' ("`country'") (`year') (`fperc') (`fperc_ll') (`fperc_ul')
    }
 }
}
postclose `foo'

// Use posted data set
use `foox', clear

// Generate country name variable using -kountry-
kountry cntry, from(iso2c)
rename NAMES_STD country
replace country = "Germany (West)" if country == "dew"
replace country = "Germany (East)" if country == "dee"

// Plot average change over time
twoway (scatter perc year, connect(direct) msymbol(o)) ///
       (rcap perc_ll perc_ul year) ///
      , by(country, ///
           note("{it:Source:} ISSP 'Family and Changing Gender Roles II{c 150}IV. {it:Note:} Error bars denote 95% confidence intervals.", span) ///
           legend(off)) ///
        xlabel(1994 2002 2012) xtitle("") ///
        ytitle("% of couples keeping incomes separate") ///
        ylabel(0 10 20 30 40 50, gstyle(minor)) yscale(r(0))
        // yscale(r(0)) adds missing gridline according to
        // http://www.stata.com/statalist/archive/2013-04/msg00904.html    
restore

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