Example 30.2 Simple Correspondence Analysis of U.S. Population

In this example, PROC CORRESP reads an existing contingency table with supplementary observations and performs a simple correspondence analysis. The data are populations of the 50 U.S. states, grouped into regions, for each of the census years from 1920 to 1970 (U.S. Bureau of the Census; 1979). Alaska and Hawaii are treated as supplementary regions, because they were not states during this entire period and are not physically connected to the other 48 states. Consequently, it is reasonable to expect that population changes in these two states operate differently from population changes in the other states. The correspondence analysis is performed giving the supplementary points negative weight, and then the coordinates for the supplementary points are computed in the solution defined by the other points.

The initial DATA step reads the table, provides labels for the years, flags the supplementary rows with negative weights, and specifies absolute weights of 1000 for all observations since the data were originally reported in units of 1000 people.

In the PROC CORRESP statement, PRINT=PERCENT and the display options display the table of cell percentages (OBSERVED), cell contributions to the total chi-square scaled to sum to 100 (CELLCHI2), row profile rows that sum to 100 (RP), and column profile columns that sum to 100 (CP). The SHORT option specifies that the correspondence analysis summary statistics, contributions to inertia, and squared cosines should not be displayed. The option OUTC=COOR creates the output coordinate data set. Since the data are already in table form, a VAR statement is used to read the table. Row labels are specified with the ID statement, and column labels come from the variable labels. The WEIGHT statement flags the supplementary observations and restores the table values to populations.

The following statements produce Output 30.2.1:

   title 'United States Population, 1920-1970';
   
   data USPop;
   
      * Regions:
      * New England     - ME, NH, VT, MA, RI, CT.
      * Great Lakes     - OH, IN, IL, MI, WI.
      * South Atlantic  - DE, MD, DC, VA, WV, NC, SC, GA, FL.
      * Mountain        - MT, ID, WY, CO, NM, AZ, UT, NV.
      * Pacific         - WA, OR, CA.
      *
      * Note: Multiply data values by 1000 to get populations.;
   
      input Region $14. y1920 y1930 y1940 y1950 y1960 y1970;
   
      label y1920 = '1920'    y1930 = '1930'    y1940 = '1940'
            y1950 = '1950'    y1960 = '1960'    y1970 = '1970';
   
      if region = 'Hawaii' or region = 'Alaska'
         then w = -1000;       /* Flag Supplementary Observations */
         else w =  1000;
   
      datalines;
   New England        7401  8166  8437  9314 10509 11842
   NY, NJ, PA        22261 26261 27539 30146 34168 37199
   Great Lakes       21476 25297 26626 30399 36225 40252
   Midwest           12544 13297 13517 14061 15394 16319
   South Atlantic    13990 15794 17823 21182 25972 30671
   KY, TN, AL, MS     8893  9887 10778 11447 12050 12803
   AR, LA, OK, TX    10242 12177 13065 14538 16951 19321
   Mountain           3336  3702  4150  5075  6855  8282
   Pacific            5567  8195  9733 14486 20339 25454
   Alaska               55    59    73   129   226   300
   Hawaii              256   368   423   500   633   769
   ;

   ods graphics on;
   
   * Perform Simple Correspondence Analysis;
   proc corresp data=uspop print=percent observed cellchi2 rp cp
        short outc=Coor plot(flip);
      var y1920 -- y1970;
      id Region;
      weight w;
   run;
   
   ods graphics off;

The contingency table shows that the population of all regions increased over this time period. The row profiles show that population increased at a different rate for the different regions. There was a small increase in population in the Midwest, for example, but the population more than quadrupled in the Pacific region over the same period. The column profiles show that in 1920, the U.S. population was concentrated in the NY, NJ, PA, Great Lakes, Midwest, and South Atlantic regions. With time, the population shifted more to the South Atlantic, Mountain, and Pacific regions. This is also clear from the correspondence analysis. The inertia and chi-square decomposition table shows that there are five nontrivial dimensions in the table, but the association between the rows and columns is almost entirely one-dimensional.