/* immorigin.sas */
/* the notation "region : $ 14." in the input
     statement indicates that the character 
     variable "region" has values up to 14 
     characters long --
     $ indicates character values -- this will
        truncate at 8 characters
     : $ 14. indicates up to 14 characters  */
data immorigin;
   input region : $ 14. yrs31_40 yrs61_70 yrs91_00;   
cards;
Europe 347566 1123492 1359737
Asia 16595 427692 2795672
NorthAmerica 130871 886891 2441448
Caribbean 15502 470213 978787
CentralAmerica 5861 101330 526915
SouthAmerica 7803 257940 539656
Africa 1750 28954 354939
Oceania 2483 25122 55845
;
proc print data=immorigin;
title 'immigrant region of last residence data';
proc sgplot data=immorigin;
  title 'immigrant region of last residence 1931-1940';
  hbar region / response=yrs31_40
                stat=percent
                datalabel;
proc sgplot data=immorigin;
  title 'immigrant region of last residence 1961-1970';
  hbar region / response=yrs61_70
                stat=percent
                datalabel;
proc sgplot data=immorigin;
  title 'immigrant region of last residence 1991-2000'; 
  hbar region / response=yrs91_00
                stat=percent
                datalabel;
run;                
title;
data decade1;
    set immorigin;
    decade='yrs31_40';
    number=yrs31_40;
    keep region number decade;
data decade2;
    set immorigin;
    decade='yrs61_70';
    number=yrs61_70;
    keep region number decade;
data decade3;
    set immorigin;
    decade='yrs91_00';
    number=yrs91_00;
    keep region number decade;
data immorigin2;
    set decade1 decade2 decade3;
proc print data=immorigin2;
title 'immigrant region of last residence data';
proc sgpanel data=immorigin2 pctlevel=cell;
title 'Immigrant region of last residence by decade (bar graphs)';
  panelby decade / columns=1;
  hbar region / response=number
                stat=percent
                datalabel;
run;
title; 
proc sgplot data=immorigin2 pctlevel=group;
title 'Immigrant region of last residence by decade (segmented bar graphs)';
  hbar decade / response=number stat=percent
     group=region;
run;
title;
proc export data=immorigin
   outfile='/folders/myfolders/CBSASfiles/immorigin.jmp'
   dbms=jmp replace;
proc export data=immorigin2
   outfile='/folders/myfolders/CBSASfiles/immorigin2.jmp'
   dbms=jmp replace;