* working directory is c:\tmp\ * dataset will be downloaded from http://www.louischauvel.org cd "C:\tmp\" use "http://www.louischauvel.org/hpiv2015.dta", clear gen yq=real(substr(time,1,4)) replace yq=yq+(real(substr(ti,7,1))-1)/4 two (li rhpi_lu y) (li rpdi_lu y) /// (li rhpi_fr y) (li rpdi_fr y) /// (li rhpi_us y) (li rpdi_us y) /// (li rhpi_de y) (li rpdi_de y) keep if yq==int(yq) *reshape long rhpi_ rpdi_ , i(yq) j(co) local ff " rhpi_au rhpi_be rhpi_ca rhpi_ch rhpi_de rhpi_dk rhpi_es rhpi_fi rhpi_fr rhpi_uk rhpi_ir rhpi_it rhpi_jp rhpi_kr rhpi_lu rhpi_ne rhpi_no rhpi_nz rhpi_se rhpi_us rhpi_za rhpi_hr rhpi_il rpdi_au rpdi_be rpdi_ca rpdi_ch rpdi_de rpdi_dk rpdi_es rpdi_fi rpdi_fr rpdi_uk rpdi_ir rpdi_it rpdi_jp rpdi_kr rpdi_lu rpdi_ne rpdi_no rpdi_nz rpdi_se rpdi_us rpdi_za rpdi_hr rpdi_il " foreach cc in `ff' { su `cc' if y==2000 replace `cc' = `cc'/r(mean) *100 } rename y year local ff "de fr uk au ca us" local ii=0 foreach cc in `ff' { local ii=`ii'+1 two (sca rhpi_`cc' y) (li rhpi_`cc' y) (li rpdi_`cc' y), title("`cc'") saving(g`ii', replace) legend(off) } graph combine g1.gph g2.gph g3.gph g4.gph g5.gph g6.gph , ycomm imargin(zero)