* check that institutionalization of mental health patients can't explain findings
use country numISO year oecd prisrate mental using paneldata.dta, clear
replace prisrate=L.prisrate if mi(prisrate) // ICPS 9 covers 2010, WHO covers 2011
keep if year==2011 & !mi(mental)
gen institutionalized=mental+prisrate
sum mentalh institutionalized prisrate if oecd
gsort -institut
list country instit mental prisrate if oecd // US tops the list, 4 s.d. above mean
gsort -mentalh
list country mentalh if oecd // US ranked 19/30, 1/2 s.d. below mean; DK, SWE even lower
