Exercises 3.4.1: Use the read.table function to read in the syphilis. Evaluate structure of data frame. Do not attach the data frame (yet). Create a 3-dimensional array using both the table or xtabs function. Now attach the data frame using the attach function. Create the same 3-dimensional array using both the table or xtabs function.

std<-read.csv("https://raw.githubusercontent.com/taragonmd/data/master/syphilis89c.csv")

Without attach Using Table

table(std$Race, std$Age, std$Sex)
## , ,  = Female
## 
##        
##         <=14  >55 15-19 20-24 25-29 30-34 35-44 45-54
##   Black  165   92  2257  4503  3590  2628  1505   392
##   Other   11   15   158   307   283   167   149    40
##   White   14   24   253   475   433   316   243    55
## 
## , ,  = Male
## 
##        
##         <=14  >55 15-19 20-24 25-29 30-34 35-44 45-54
##   Black   31  823  1412  4059  4121  4453  3858  1619
##   Other    7  108   210   654   633   520   492   202
##   White    2  216    88   407   550   564   654   323

Without attach Using xtabs function

z <- xtabs(~ Age + Race + Sex, data = std)
z
## , , Sex = Female
## 
##        Race
## Age     Black Other White
##   <=14    165    11    14
##   >55      92    15    24
##   15-19  2257   158   253
##   20-24  4503   307   475
##   25-29  3590   283   433
##   30-34  2628   167   316
##   35-44  1505   149   243
##   45-54   392    40    55
## 
## , , Sex = Male
## 
##        Race
## Age     Black Other White
##   <=14     31     7     2
##   >55     823   108   216
##   15-19  1412   210    88
##   20-24  4059   654   407
##   25-29  4121   633   550
##   30-34  4453   520   564
##   35-44  3858   492   654
##   45-54  1619   202   323

With attach

attach(std)
table(Race, Age, Sex)
## , , Sex = Female
## 
##        Age
## Race    <=14  >55 15-19 20-24 25-29 30-34 35-44 45-54
##   Black  165   92  2257  4503  3590  2628  1505   392
##   Other   11   15   158   307   283   167   149    40
##   White   14   24   253   475   433   316   243    55
## 
## , , Sex = Male
## 
##        Age
## Race    <=14  >55 15-19 20-24 25-29 30-34 35-44 45-54
##   Black   31  823  1412  4059  4121  4453  3858  1619
##   Other    7  108   210   654   633   520   492   202
##   White    2  216    88   407   550   564   654   323
xtabs(~ Race + Age + Sex)
## , , Sex = Female
## 
##        Age
## Race    <=14  >55 15-19 20-24 25-29 30-34 35-44 45-54
##   Black  165   92  2257  4503  3590  2628  1505   392
##   Other   11   15   158   307   283   167   149    40
##   White   14   24   253   475   433   316   243    55
## 
## , , Sex = Male
## 
##        Age
## Race    <=14  >55 15-19 20-24 25-29 30-34 35-44 45-54
##   Black   31  823  1412  4059  4121  4453  3858  1619
##   Other    7  108   210   654   633   520   492   202
##   White    2  216    88   407   550   564   654   323
detach(std)

3.4.2: Use the apply function to get marginal totals for the syphilis 3-dimensional array.

z <- xtabs(~ Age + Race + Sex, data = std)
#### 2-D tables
z.ar <- apply(z, c(1, 2), sum); z.ar
##        Race
## Age     Black Other White
##   <=14    196    18    16
##   >55     915   123   240
##   15-19  3669   368   341
##   20-24  8562   961   882
##   25-29  7711   916   983
##   30-34  7081   687   880
##   35-44  5363   641   897
##   45-54  2011   242   378
z.as <- apply(z, c(1, 3), sum); z.as
##        Sex
## Age     Female Male
##   <=14     190   40
##   >55      131 1147
##   15-19   2668 1710
##   20-24   5285 5120
##   25-29   4306 5304
##   30-34   3111 5537
##   35-44   1897 5004
##   45-54    487 2144
z.rs <- apply(z, c(2, 3), sum); z.rs
##        Sex
## Race    Female  Male
##   Black  15132 20376
##   Other   1130  2826
##   White   1813  2804
#### 1-D tables
z.a <- apply(z, 1, sum); z.a
##  <=14   >55 15-19 20-24 25-29 30-34 35-44 45-54 
##   230  1278  4378 10405  9610  8648  6901  2631
z.r <- apply(z, 2, sum); z.r
## Black Other White 
## 35508  3956  4617
z.s <- apply(z, 3, sum); z.s
## Female   Male 
##  18075  26006

3.4.3:Use the sweep and apply functions to get marginal and joint distributions for a 3-D array.

z <- xtabs(~ Age + Race + Sex, data = std)
rowt <- apply(z, c(1, 3), sum) # row distrib
rowd <- sweep(z, c(1, 3), rowt, "/"); rowd
## , , Sex = Female
## 
##        Race
## Age          Black      Other      White
##   <=14  0.86842105 0.05789474 0.07368421
##   >55   0.70229008 0.11450382 0.18320611
##   15-19 0.84595202 0.05922039 0.09482759
##   20-24 0.85203406 0.05808893 0.08987701
##   25-29 0.83372039 0.06572225 0.10055736
##   30-34 0.84474446 0.05368049 0.10157506
##   35-44 0.79335793 0.07854507 0.12809700
##   45-54 0.80492813 0.08213552 0.11293634
## 
## , , Sex = Male
## 
##        Race
## Age          Black      Other      White
##   <=14  0.77500000 0.17500000 0.05000000
##   >55   0.71752398 0.09415867 0.18831735
##   15-19 0.82573099 0.12280702 0.05146199
##   20-24 0.79277344 0.12773438 0.07949219
##   25-29 0.77696078 0.11934389 0.10369532
##   30-34 0.80422612 0.09391367 0.10186021
##   35-44 0.77098321 0.09832134 0.13069544
##   45-54 0.75513060 0.09421642 0.15065299
apply(rowd, c(1, 3), sum) # confirm
##        Sex
## Age     Female Male
##   <=14       1    1
##   >55        1    1
##   15-19      1    1
##   20-24      1    1
##   25-29      1    1
##   30-34      1    1
##   35-44      1    1
##   45-54      1    1
colt <- apply(z, c(2, 3), sum) #  col distrib
cold <- sweep(z, c(2, 3), colt, "/"); cold
## , , Sex = Female
## 
##        Race
## Age            Black        Other        White
##   <=14  0.0109040444 0.0097345133 0.0077220077
##   >55   0.0060798308 0.0132743363 0.0132377275
##   15-19 0.1491541105 0.1398230088 0.1395477110
##   20-24 0.2975812847 0.2716814159 0.2619966906
##   25-29 0.2372455723 0.2504424779 0.2388306674
##   30-34 0.1736716891 0.1477876106 0.1742967457
##   35-44 0.0994581020 0.1318584071 0.1340319912
##   45-54 0.0259053661 0.0353982301 0.0303364589
## 
## , , Sex = Male
## 
##        Race
## Age            Black        Other        White
##   <=14  0.0015213977 0.0024769993 0.0007132668
##   >55   0.0403906557 0.0382165605 0.0770328103
##   15-19 0.0692972124 0.0743099788 0.0313837375
##   20-24 0.1992049470 0.2314225053 0.1451497860
##   25-29 0.2022477424 0.2239915074 0.1961483595
##   30-34 0.2185414213 0.1840056617 0.2011412268
##   35-44 0.1893404005 0.1740976645 0.2332382311
##   45-54 0.0794562230 0.0714791224 0.1151925820
apply(cold, c(2, 3), sum) # confirm
##        Sex
## Race    Female Male
##   Black      1    1
##   Other      1    1
##   White      1    1
jtd <- apply(z, 3, sum) # joint distrib
jtd <- sweep(z, 3, jtd, "/"); jtd
## , , Sex = Female
## 
##        Race
## Age            Black        Other        White
##   <=14  9.128631e-03 6.085754e-04 7.745505e-04
##   >55   5.089903e-03 8.298755e-04 1.327801e-03
##   15-19 1.248686e-01 8.741355e-03 1.399723e-02
##   20-24 2.491286e-01 1.698479e-02 2.627939e-02
##   25-29 1.986169e-01 1.565698e-02 2.395574e-02
##   30-34 1.453942e-01 9.239281e-03 1.748271e-02
##   35-44 8.326418e-02 8.243430e-03 1.344398e-02
##   45-54 2.168741e-02 2.213001e-03 3.042877e-03
## 
## , , Sex = Male
## 
##        Race
## Age            Black        Other        White
##   <=14  1.192033e-03 2.691687e-04 7.690533e-05
##   >55   3.164654e-02 4.152888e-03 8.305776e-03
##   15-19 5.429516e-02 8.075060e-03 3.383834e-03
##   20-24 1.560794e-01 2.514804e-02 1.565023e-02
##   25-29 1.584634e-01 2.434054e-02 2.114897e-02
##   30-34 1.712297e-01 1.999539e-02 2.168730e-02
##   35-44 1.483504e-01 1.891871e-02 2.514804e-02
##   45-54 6.225486e-02 7.767438e-03 1.242021e-02
apply(jtd, 3, sum) # confirm
## Female   Male 
##      1      1
distr <- list(rowd, cold, jtd); distr 
## [[1]]
## , , Sex = Female
## 
##        Race
## Age          Black      Other      White
##   <=14  0.86842105 0.05789474 0.07368421
##   >55   0.70229008 0.11450382 0.18320611
##   15-19 0.84595202 0.05922039 0.09482759
##   20-24 0.85203406 0.05808893 0.08987701
##   25-29 0.83372039 0.06572225 0.10055736
##   30-34 0.84474446 0.05368049 0.10157506
##   35-44 0.79335793 0.07854507 0.12809700
##   45-54 0.80492813 0.08213552 0.11293634
## 
## , , Sex = Male
## 
##        Race
## Age          Black      Other      White
##   <=14  0.77500000 0.17500000 0.05000000
##   >55   0.71752398 0.09415867 0.18831735
##   15-19 0.82573099 0.12280702 0.05146199
##   20-24 0.79277344 0.12773438 0.07949219
##   25-29 0.77696078 0.11934389 0.10369532
##   30-34 0.80422612 0.09391367 0.10186021
##   35-44 0.77098321 0.09832134 0.13069544
##   45-54 0.75513060 0.09421642 0.15065299
## 
## 
## [[2]]
## , , Sex = Female
## 
##        Race
## Age            Black        Other        White
##   <=14  0.0109040444 0.0097345133 0.0077220077
##   >55   0.0060798308 0.0132743363 0.0132377275
##   15-19 0.1491541105 0.1398230088 0.1395477110
##   20-24 0.2975812847 0.2716814159 0.2619966906
##   25-29 0.2372455723 0.2504424779 0.2388306674
##   30-34 0.1736716891 0.1477876106 0.1742967457
##   35-44 0.0994581020 0.1318584071 0.1340319912
##   45-54 0.0259053661 0.0353982301 0.0303364589
## 
## , , Sex = Male
## 
##        Race
## Age            Black        Other        White
##   <=14  0.0015213977 0.0024769993 0.0007132668
##   >55   0.0403906557 0.0382165605 0.0770328103
##   15-19 0.0692972124 0.0743099788 0.0313837375
##   20-24 0.1992049470 0.2314225053 0.1451497860
##   25-29 0.2022477424 0.2239915074 0.1961483595
##   30-34 0.2185414213 0.1840056617 0.2011412268
##   35-44 0.1893404005 0.1740976645 0.2332382311
##   45-54 0.0794562230 0.0714791224 0.1151925820
## 
## 
## [[3]]
## , , Sex = Female
## 
##        Race
## Age            Black        Other        White
##   <=14  9.128631e-03 6.085754e-04 7.745505e-04
##   >55   5.089903e-03 8.298755e-04 1.327801e-03
##   15-19 1.248686e-01 8.741355e-03 1.399723e-02
##   20-24 2.491286e-01 1.698479e-02 2.627939e-02
##   25-29 1.986169e-01 1.565698e-02 2.395574e-02
##   30-34 1.453942e-01 9.239281e-03 1.748271e-02
##   35-44 8.326418e-02 8.243430e-03 1.344398e-02
##   45-54 2.168741e-02 2.213001e-03 3.042877e-03
## 
## , , Sex = Male
## 
##        Race
## Age            Black        Other        White
##   <=14  1.192033e-03 2.691687e-04 7.690533e-05
##   >55   3.164654e-02 4.152888e-03 8.305776e-03
##   15-19 5.429516e-02 8.075060e-03 3.383834e-03
##   20-24 1.560794e-01 2.514804e-02 1.565023e-02
##   25-29 1.584634e-01 2.434054e-02 2.114897e-02
##   30-34 1.712297e-01 1.999539e-02 2.168730e-02
##   35-44 1.483504e-01 1.891871e-02 2.514804e-02
##   45-54 6.225486e-02 7.767438e-03 1.242021e-02

3.4.4: Review and read in the group-level, tabular data set of primary and secondary syphilis cases in the United States in 1989. Use the rep function on the data frame fields to recreate the individual-level data frame with over 40,000 observations.

stdb <- read.csv("https://raw.githubusercontent.com/taragonmd/data/master/syphilis89b.csv",
as.is = c(FALSE, FALSE, TRUE, FALSE))
str(stdb)
## 'data.frame':    48 obs. of  4 variables:
##  $ Sex : Factor w/ 2 levels "Female","Male": 2 1 2 1 2 1 2 1 2 1 ...
##  $ Race: Factor w/ 3 levels "Black","Other",..: 3 3 1 1 2 2 3 3 1 1 ...
##  $ Age : chr  "<=14" "<=14" "<=14" "<=14" ...
##  $ Freq: int  2 14 31 165 7 11 88 253 1412 2257 ...
stdb
##       Sex  Race   Age Freq
## 1    Male White  <=14    2
## 2  Female White  <=14   14
## 3    Male Black  <=14   31
## 4  Female Black  <=14  165
## 5    Male Other  <=14    7
## 6  Female Other  <=14   11
## 7    Male White 15-19   88
## 8  Female White 15-19  253
## 9    Male Black 15-19 1412
## 10 Female Black 15-19 2257
## 11   Male Other 15-19  210
## 12 Female Other 15-19  158
## 13   Male White 20-24  407
## 14 Female White 20-24  475
## 15   Male Black 20-24 4059
## 16 Female Black 20-24 4503
## 17   Male Other 20-24  654
## 18 Female Other 20-24  307
## 19   Male White 25-29  550
## 20 Female White 25-29  433
## 21   Male Black 25-29 4121
## 22 Female Black 25-29 3590
## 23   Male Other 25-29  633
## 24 Female Other 25-29  283
## 25   Male White 30-34  564
## 26 Female White 30-34  316
## 27   Male Black 30-34 4453
## 28 Female Black 30-34 2628
## 29   Male Other 30-34  520
## 30 Female Other 30-34  167
## 31   Male White 35-44  654
## 32 Female White 35-44  243
## 33   Male Black 35-44 3858
## 34 Female Black 35-44 1505
## 35   Male Other 35-44  492
## 36 Female Other 35-44  149
## 37   Male White 45-54  323
## 38 Female White 45-54   55
## 39   Male Black 45-54 1619
## 40 Female Black 45-54  392
## 41   Male Other 45-54  202
## 42 Female Other 45-54   40
## 43   Male White   >55  216
## 44 Female White   >55   24
## 45   Male Black   >55  823
## 46 Female Black   >55   92
## 47   Male Other   >55  108
## 48 Female Other   >55   15
expand.rows <- rep(1:nrow(stdb),stdb$Freq)
stdb.df <- stdb[expand.rows,]
agelab<-c("<=14","15-19","20-24","25-29","30-34","35-44","45-54",">55")
stdb.df$Age <- factor(stdb.df$Age, levels = agelab)
str(stdb.df)
## 'data.frame':    44081 obs. of  4 variables:
##  $ Sex : Factor w/ 2 levels "Female","Male": 2 2 1 1 1 1 1 1 1 1 ...
##  $ Race: Factor w/ 3 levels "Black","Other",..: 3 3 3 3 3 3 3 3 3 3 ...
##  $ Age : Factor w/ 8 levels "<=14","15-19",..: 1 1 1 1 1 1 1 1 1 1 ...
##  $ Freq: int  2 2 14 14 14 14 14 14 14 14 ...
table(stdb.df$Age)
## 
##  <=14 15-19 20-24 25-29 30-34 35-44 45-54   >55 
##   230  4378 10405  9610  8648  6901  2631  1278

3.4.5:Working with population estimates can be challenging because of the amount of data manipulation. Study the 2000 population estimates for California Counties: [~/data/calpop/CalCounties2000.txt]. Now, study and implement this R code. For each expression or group of expressions, explain in words what the R code is doing. Be sure to display intermediate objects to understand each step.

#### 1 Read county names
cty <- scan("https://raw.githubusercontent.com/taragonmd/data/master/calpop/calcounty.txt", what="")
#### 2 Read county population estimates
calpop = read.csv("https://raw.githubusercontent.com/taragonmd/data/master/calpop/CalCounties2000.csv", header = T)
#### 2 Replace county number with county name
for(i in 1:length(cty)){
calpop$County[calpop$County==i] <- cty[i]
}
#### 3 Discretize age into categories
calpop$Agecat <- cut(calpop$Age, c(0,20,45,65,100),
include.lowest = TRUE, right = FALSE)
#### 4 Create combined API category
calpop$AsianPI <- calpop$Asian + calpop$Pac.Isl
#### 5 Shorten selected ethnic labels
names(calpop)[c(6, 9, 10)] = c("Latino", "AfrAmer", "AmerInd")
#### 6 Index Bay Area Counties
baindex <- calpop$County=="Alameda" | calpop$County=="San Francisco"
bapop <- calpop[baindex,]
bapop
##             County Year Sex Age White Latino Asian Pac.Isl AfrAmer AmerInd
## 1          Alameda 2000   F   0  2751   2910  1921      63    1346      32
## 2          Alameda 2000   F   1  2673   2837  1820      70    1343      37
## 3          Alameda 2000   F   2  2654   2770  1930      64    1442      37
## 4          Alameda 2000   F   3  2646   2798  1998      64    1455      39
## 5          Alameda 2000   F   4  2780   2762  2029      85    1558      38
## 6          Alameda 2000   F   5  2826   2738  1998      93    1632      45
## 7          Alameda 2000   F   6  2871   2872  1984      78    1657      53
## 8          Alameda 2000   F   7  2942   2837  1993      94    1732      51
## 9          Alameda 2000   F   8  3087   2756  1957      85    1817      43
## 10         Alameda 2000   F   9  3070   2603  2031      83    1819      44
## 11         Alameda 2000   F  10  3076   2465  1967      89    1820      39
## 12         Alameda 2000   F  11  3056   2315  1945      77    1802      40
## 13         Alameda 2000   F  12  3167   2236  1897      71    1762      47
## 14         Alameda 2000   F  13  3077   2168  1862      78    1570      40
## 15         Alameda 2000   F  14  2979   2094  1736      82    1573      37
## 16         Alameda 2000   F  15  3027   2023  1855      93    1524      47
## 17         Alameda 2000   F  16  2972   1987  1869      79    1543      51
## 18         Alameda 2000   F  17  2880   2045  1896      69    1552      33
## 19         Alameda 2000   F  18  2787   2210  2056      93    1442      45
## 20         Alameda 2000   F  19  2839   2312  2307      92    1376      41
## 21         Alameda 2000   F  20  2870   2371  2522      86    1469      49
## 22         Alameda 2000   F  21  2869   2417  2526      92    1462      51
## 23         Alameda 2000   F  22  3079   2581  2399      86    1445      45
## 24         Alameda 2000   F  23  3260   2629  2432      83    1490      50
## 25         Alameda 2000   F  24  3430   2675  2573      77    1540      42
## 26         Alameda 2000   F  25  3694   2771  2731      73    1597      48
## 27         Alameda 2000   F  26  3701   2748  2790      74    1636      48
## 28         Alameda 2000   F  27  3874   2861  2966      88    1696      48
## 29         Alameda 2000   F  28  4037   2837  3022      76    1773      64
## 30         Alameda 2000   F  29  4449   2734  3125      87    1915      66
## 31         Alameda 2000   F  30  4978   2680  3324      86    2097      59
## 32         Alameda 2000   F  31  4792   2540  3251      85    1847      51
## 33         Alameda 2000   F  32  4521   2493  3062     101    1667      49
## 34         Alameda 2000   F  33  4462   2383  2887      82    1719      47
## 35         Alameda 2000   F  34  4612   2341  2916      73    1877      60
## 36         Alameda 2000   F  35  5100   2362  3077      90    2110      67
## 37         Alameda 2000   F  36  5252   2295  3195      92    2030      75
## 38         Alameda 2000   F  37  5183   2172  3145      91    1909      71
## 39         Alameda 2000   F  38  5330   2154  2923      90    1922      54
## 40         Alameda 2000   F  39  5484   1961  2815      95    1924      61
## 41         Alameda 2000   F  40  5586   1934  2847      99    2044      68
## 42         Alameda 2000   F  41  5453   1790  2645      92    1872      60
## 43         Alameda 2000   F  42  5460   1697  2607      83    1885      51
## 44         Alameda 2000   F  43  5456   1569  2538      70    1754      54
## 45         Alameda 2000   F  44  5394   1558  2484      70    1757      63
## 46         Alameda 2000   F  45  5445   1532  2485      62    1835      60
## 47         Alameda 2000   F  46  5328   1423  2337      68    1676      59
## 48         Alameda 2000   F  47  5310   1306  2266      64    1765      55
## 49         Alameda 2000   F  48  5204   1197  2190      51    1679      60
## 50         Alameda 2000   F  49  5179   1181  2174      54    1629      55
## 51         Alameda 2000   F  50  5274   1179  2246      69    1728      46
## 52         Alameda 2000   F  51  5099   1126  2085      58    1597      64
## 53         Alameda 2000   F  52  5206   1089  1870      43    1508      53
## 54         Alameda 2000   F  53  5387   1001  1721      36    1433      52
## 55         Alameda 2000   F  54  4387    884  1479      35    1186      46
## 56         Alameda 2000   F  55  4132    841  1363      30    1164      40
## 57         Alameda 2000   F  56  3921    759  1284      28    1059      34
## 58         Alameda 2000   F  57  3776    755  1202      25     987      26
## 59         Alameda 2000   F  58  3308    668  1166      20     973      30
## 60         Alameda 2000   F  59  2915    659  1144      25     918      36
## 61         Alameda 2000   F  60  2775    694  1103      29     884      32
## 62         Alameda 2000   F  61  2658    603  1065      25     783      22
## 63         Alameda 2000   F  62  2464    573  1082      26     751      15
## 64         Alameda 2000   F  63  2274    558  1099      24     725      19
## 65         Alameda 2000   F  64  2163    506  1079      24     668      18
## 66         Alameda 2000   F  65  2219    505  1106      24     735      14
## 67         Alameda 2000   F  66  2075    445  1065      26     650       9
## 68         Alameda 2000   F  67  2016    462   986      18     637      17
## 69         Alameda 2000   F  68  1994    468   958      21     577      15
## 70         Alameda 2000   F  69  1984    459   949      15     589      12
## 71         Alameda 2000   F  70  2154    465   943      12     648      10
## 72         Alameda 2000   F  71  2074    396   888      13     593      11
## 73         Alameda 2000   F  72  2046    392   841      18     588      12
## 74         Alameda 2000   F  73  2160    397   815      15     613      13
## 75         Alameda 2000   F  74  2307    368   824      12     638      12
## 76         Alameda 2000   F  75  2301    336   767      11     628      13
## 77         Alameda 2000   F  76  2272    317   715      13     585      15
## 78         Alameda 2000   F  77  2209    302   615       8     541      10
## 79         Alameda 2000   F  78  2246    276   539       4     521       9
## 80         Alameda 2000   F  79  2127    248   533       5     490       7
## 81         Alameda 2000   F  80  1980    216   455       6     501       9
## 82         Alameda 2000   F  81  1800    201   393       7     441       5
## 83         Alameda 2000   F  82  1679    178   347       9     414       5
## 84         Alameda 2000   F  83  1542    158   307       4     358       6
## 85         Alameda 2000   F  84  1445    143   281       3     300       8
## 86         Alameda 2000   F  85  1355    127   262       2     303       7
## 87         Alameda 2000   F  86  1294    108   234       2     263       5
## 88         Alameda 2000   F  87  1113    110   188       3     213       8
## 89         Alameda 2000   F  88   990     98   145       2     182       5
## 90         Alameda 2000   F  89   861     68   130       2     149       3
## 91         Alameda 2000   F  90   717     65   102       3     121       2
## 92         Alameda 2000   F  91   580     52    78       2     102       2
## 93         Alameda 2000   F  92   503     39    56       1      91       0
## 94         Alameda 2000   F  93   438     43    46       0      76       1
## 95         Alameda 2000   F  94   360     24    36       0      56       1
## 96         Alameda 2000   F  95   254     21    34       0      50       0
## 97         Alameda 2000   F  96   208     12    24       0      46       0
## 98         Alameda 2000   F  97   146     13    18       0      39       0
## 99         Alameda 2000   F  98    88     15    21       0      25       0
## 100        Alameda 2000   F  99    87     10    14       1      22       0
## 101        Alameda 2000   F 100   138     11    27       0      49       0
## 102        Alameda 2000   M   0  2905   3025  2033      57    1344      30
## 103        Alameda 2000   M   1  2788   2983  1985      76    1346      37
## 104        Alameda 2000   M   2  2793   2913  2007      82    1452      30
## 105        Alameda 2000   M   3  2841   2896  2132      77    1539      36
## 106        Alameda 2000   M   4  2914   2877  2141      75    1620      32
## 107        Alameda 2000   M   5  3082   2971  2178      95    1632      44
## 108        Alameda 2000   M   6  3068   2904  2214      87    1720      38
## 109        Alameda 2000   M   7  3113   2934  2128      82    1785      44
## 110        Alameda 2000   M   8  3228   2801  2177      88    1800      55
## 111        Alameda 2000   M   9  3282   2748  2108      81    1874      48
## 112        Alameda 2000   M  10  3287   2669  2052      76    1968      48
## 113        Alameda 2000   M  11  3214   2445  2110      89    1825      45
## 114        Alameda 2000   M  12  3212   2398  2002      90    1706      49
## 115        Alameda 2000   M  13  3174   2223  1876      86    1603      47
## 116        Alameda 2000   M  14  3211   2168  1882      95    1561      40
## 117        Alameda 2000   M  15  3120   2194  1963      79    1568      40
## 118        Alameda 2000   M  16  3044   2244  1970      89    1553      34
## 119        Alameda 2000   M  17  3125   2365  1947      76    1557      50
## 120        Alameda 2000   M  18  3014   2559  2038      84    1488      42
## 121        Alameda 2000   M  19  3031   2780  2322      93    1336      39
## 122        Alameda 2000   M  20  3163   3058  2539      91    1347      37
## 123        Alameda 2000   M  21  3075   3086  2454      84    1278      49
## 124        Alameda 2000   M  22  3179   3161  2233      81    1178      53
## 125        Alameda 2000   M  23  3261   3120  2266      70    1216      49
## 126        Alameda 2000   M  24  3458   3257  2297      73    1271      36
## 127        Alameda 2000   M  25  3752   3349  2494      70    1312      38
## 128        Alameda 2000   M  26  3715   3352  2566      75    1256      39
## 129        Alameda 2000   M  27  3910   3374  2811      82    1317      47
## 130        Alameda 2000   M  28  4100   3351  2897      74    1340      42
## 131        Alameda 2000   M  29  4555   3481  2988      82    1491      43
## 132        Alameda 2000   M  30  5165   3412  3078      84    1731      56
## 133        Alameda 2000   M  31  4873   3032  2995      85    1579      59
## 134        Alameda 2000   M  32  4627   2944  2827      75    1441      70
## 135        Alameda 2000   M  33  4641   2817  2701      78    1452      41
## 136        Alameda 2000   M  34  4823   2802  2749      84    1558      46
## 137        Alameda 2000   M  35  5355   2815  3011      86    1762      74
## 138        Alameda 2000   M  36  5464   2561  3048     104    1624      58
## 139        Alameda 2000   M  37  5444   2384  3001      95    1617      48
## 140        Alameda 2000   M  38  5550   2306  2722      90    1666      48
## 141        Alameda 2000   M  39  5580   2104  2611      84    1651      65
## 142        Alameda 2000   M  40  5856   2175  2779      98    1823      55
## 143        Alameda 2000   M  41  5537   1984  2652      90    1613      57
## 144        Alameda 2000   M  42  5635   1843  2509      84    1529      56
## 145        Alameda 2000   M  43  5492   1777  2412      69    1464      53
## 146        Alameda 2000   M  44  5535   1688  2351      74    1460      44
## 147        Alameda 2000   M  45  5520   1630  2381      72    1591      53
## 148        Alameda 2000   M  46  5316   1416  2226      63    1491      46
## 149        Alameda 2000   M  47  5327   1276  2044      56    1403      49
## 150        Alameda 2000   M  48  5207   1235  1993      55    1300      47
## 151        Alameda 2000   M  49  4939   1140  1929      58    1289      44
## 152        Alameda 2000   M  50  5217   1163  2032      56    1453      38
## 153        Alameda 2000   M  51  5066   1036  1870      59    1332      39
## 154        Alameda 2000   M  52  5127    998  1744      53    1237      42
## 155        Alameda 2000   M  53  5282    914  1525      49    1146      52
## 156        Alameda 2000   M  54  4377    799  1346      40     975      37
## 157        Alameda 2000   M  55  4066    797  1256      36     937      28
## 158        Alameda 2000   M  56  3884    720  1154      41     911      38
## 159        Alameda 2000   M  57  3892    657  1054      30     866      27
## 160        Alameda 2000   M  58  3334    625   949      21     770      32
## 161        Alameda 2000   M  59  2890    558   980      25     751      32
## 162        Alameda 2000   M  60  2754    554  1023      23     754      24
## 163        Alameda 2000   M  61  2573    521   926      26     697      19
## 164        Alameda 2000   M  62  2403    457   904      25     606      18
## 165        Alameda 2000   M  63  2152    457   863      15     615      11
## 166        Alameda 2000   M  64  2006    449   823      16     613      11
## 167        Alameda 2000   M  65  1985    426   854      22     588      14
## 168        Alameda 2000   M  66  1851    385   829      18     499      12
## 169        Alameda 2000   M  67  1813    366   782      17     490       4
## 170        Alameda 2000   M  68  1786    362   754      20     471       5
## 171        Alameda 2000   M  69  1735    361   751      15     460       7
## 172        Alameda 2000   M  70  1749    333   733      12     432       7
## 173        Alameda 2000   M  71  1672    333   672      17     390       9
## 174        Alameda 2000   M  72  1648    311   634      10     405      10
## 175        Alameda 2000   M  73  1655    278   619       7     411      11
## 176        Alameda 2000   M  74  1618    252   561       4     435       7
## 177        Alameda 2000   M  75  1599    268   578       8     388       7
## 178        Alameda 2000   M  76  1577    254   528       7     346       8
## 179        Alameda 2000   M  77  1498    200   476       4     344       9
## 180        Alameda 2000   M  78  1526    191   408       3     346       5
## 181        Alameda 2000   M  79  1384    194   361       3     314       7
## 182        Alameda 2000   M  80  1256    147   368       3     312       5
## 183        Alameda 2000   M  81  1118    118   274       2     227       3
## 184        Alameda 2000   M  82  1002     98   232       2     209       4
## 185        Alameda 2000   M  83   884     82   218       5     162       3
## 186        Alameda 2000   M  84   810     75   174       3     161       2
## 187        Alameda 2000   M  85   777     62   139       4     142       3
## 188        Alameda 2000   M  86   659     58   138       2     108       2
## 189        Alameda 2000   M  87   513     54   119       1      87       2
## 190        Alameda 2000   M  88   411     52   105       1      77       0
## 191        Alameda 2000   M  89   342     36    88       2      67       3
## 192        Alameda 2000   M  90   280     29    64       1      49       3
## 193        Alameda 2000   M  91   231     29    43       0      33       2
## 194        Alameda 2000   M  92   165     23    48       0      20       1
## 195        Alameda 2000   M  93   125     10    36       1      20       0
## 196        Alameda 2000   M  94   107     10    22       0      24       0
## 197        Alameda 2000   M  95    73      7    13       0      17       0
## 198        Alameda 2000   M  96    49      3    10       1      14       1
## 199        Alameda 2000   M  97    32      3    10       0      12       0
## 200        Alameda 2000   M  98    20      3     7       0       9       0
## 201        Alameda 2000   M  99    17      2     4       0       6       0
## 202        Alameda 2000   M 100    27      1     8       0      12       0
## 7475 San Francisco 2000   F   0  1248    782   949      26     279       9
## 7476 San Francisco 2000   F   1  1012    693   936      22     285       8
## 7477 San Francisco 2000   F   2   784    667   936      22     266       9
## 7478 San Francisco 2000   F   3   731    694   950      29     283      11
## 7479 San Francisco 2000   F   4   660    656   982      35     307       7
## 7480 San Francisco 2000   F   5   622    669  1097      30     320       8
## 7481 San Francisco 2000   F   6   582    659  1102      35     380       9
## 7482 San Francisco 2000   F   7   602    675  1148      25     387       6
## 7483 San Francisco 2000   F   8   616    677  1168      22     395      12
## 7484 San Francisco 2000   F   9   608    671  1152      38     386       6
## 7485 San Francisco 2000   F  10   651    654  1133      36     398       7
## 7486 San Francisco 2000   F  11   639    592  1205      31     361       8
## 7487 San Francisco 2000   F  12   623    616  1184      25     391      11
## 7488 San Francisco 2000   F  13   651    612  1130      26     358      11
## 7489 San Francisco 2000   F  14   626    572  1148      33     348       5
## 7490 San Francisco 2000   F  15   633    588  1213      28     368       5
## 7491 San Francisco 2000   F  16   621    600  1260      30     353       6
## 7492 San Francisco 2000   F  17   617    631  1297      37     367       9
## 7493 San Francisco 2000   F  18   822    738  1254      42     370      12
## 7494 San Francisco 2000   F  19  1007    805  1404      45     384      14
## 7495 San Francisco 2000   F  20  1157    857  1573      39     361      14
## 7496 San Francisco 2000   F  21  1336    865  1619      40     350      15
## 7497 San Francisco 2000   F  22  2114    943  1779      40     342      17
## 7498 San Francisco 2000   F  23  3278   1032  2006      31     380      16
## 7499 San Francisco 2000   F  24  3993   1125  2195      31     379      22
## 7500 San Francisco 2000   F  25  4590   1202  2323      40     376      20
## 7501 San Francisco 2000   F  26  4747   1185  2394      40     362      24
## 7502 San Francisco 2000   F  27  4872   1170  2458      41     357      20
## 7503 San Francisco 2000   F  28  5259   1166  2418      37     409      26
## 7504 San Francisco 2000   F  29  5721   1198  2415      39     430      28
## 7505 San Francisco 2000   F  30  6054   1230  2448      42     484      24
## 7506 San Francisco 2000   F  31  5140   1067  2223      40     422      25
## 7507 San Francisco 2000   F  32  4370   1017  2071      41     392      25
## 7508 San Francisco 2000   F  33  3934    930  1951      30     376      20
## 7509 San Francisco 2000   F  34  3622    890  1973      35     394      20
## 7510 San Francisco 2000   F  35  3577    891  2100      34     428      18
## 7511 San Francisco 2000   F  36  3136    858  2107      31     422      18
## 7512 San Francisco 2000   F  37  2739    782  2113      33     441      23
## 7513 San Francisco 2000   F  38  2560    758  1962      26     444      18
## 7514 San Francisco 2000   F  39  2457    711  1877      32     429      23
## 7515 San Francisco 2000   F  40  2524    748  1973      32     489      26
## 7516 San Francisco 2000   F  41  2230    761  1901      31     440      29
## 7517 San Francisco 2000   F  42  2112    713  1909      20     486      18
## 7518 San Francisco 2000   F  43  2123    676  1908      19     456      19
## 7519 San Francisco 2000   F  44  2121    683  1943      16     435      18
## 7520 San Francisco 2000   F  45  2219    709  1971      23     485      22
## 7521 San Francisco 2000   F  46  2111    619  1893      17     459      21
## 7522 San Francisco 2000   F  47  2223    589  1874      15     433      14
## 7523 San Francisco 2000   F  48  2222    553  1866      18     405      12
## 7524 San Francisco 2000   F  49  2188    528  1924      19     388      18
## 7525 San Francisco 2000   F  50  2327    574  2206      25     407      20
## 7526 San Francisco 2000   F  51  2212    530  2091      15     392      21
## 7527 San Francisco 2000   F  52  2284    472  1800      17     396      21
## 7528 San Francisco 2000   F  53  2341    441  1600      20     357      17
## 7529 San Francisco 2000   F  54  1914    441  1408      19     313      16
## 7530 San Francisco 2000   F  55  1793    420  1322      14     291      17
## 7531 San Francisco 2000   F  56  1695    381  1187      11     319       9
## 7532 San Francisco 2000   F  57  1619    379  1127      10     316      10
## 7533 San Francisco 2000   F  58  1413    385  1183      13     273       9
## 7534 San Francisco 2000   F  59  1277    365  1220      15     283       8
## 7535 San Francisco 2000   F  60  1265    373  1305      10     317       9
## 7536 San Francisco 2000   F  61  1213    368  1249      11     261       6
## 7537 San Francisco 2000   F  62  1164    353  1320      11     263       8
## 7538 San Francisco 2000   F  63  1098    356  1293      15     259      12
## 7539 San Francisco 2000   F  64  1039    348  1329      12     273      12
## 7540 San Francisco 2000   F  65   979    348  1325      10     259       6
## 7541 San Francisco 2000   F  66   932    339  1362       4     234       6
## 7542 San Francisco 2000   F  67   937    309  1289       5     245       6
## 7543 San Francisco 2000   F  68   958    298  1213       9     224       7
## 7544 San Francisco 2000   F  69  1016    335  1317       9     231       7
## 7545 San Francisco 2000   F  70  1087    320  1396      10     232       4
## 7546 San Francisco 2000   F  71  1081    303  1278       8     231       5
## 7547 San Francisco 2000   F  72  1129    281  1260       8     243       9
## 7548 San Francisco 2000   F  73  1138    254  1277       8     254       5
## 7549 San Francisco 2000   F  74  1208    260  1262       5     237       4
## 7550 San Francisco 2000   F  75  1257    245  1135       3     258       4
## 7551 San Francisco 2000   F  76  1261    223   996       4     249       4
## 7552 San Francisco 2000   F  77  1265    242   913       6     240       8
## 7553 San Francisco 2000   F  78  1198    232   855       3     239       5
## 7554 San Francisco 2000   F  79  1168    193   773       3     199       3
## 7555 San Francisco 2000   F  80  1091    211   718       2     196       2
## 7556 San Francisco 2000   F  81   995    156   628       2     171       2
## 7557 San Francisco 2000   F  82   925    145   540       5     159       2
## 7558 San Francisco 2000   F  83   882    139   487       3     157       3
## 7559 San Francisco 2000   F  84   839    126   454       3     156       4
## 7560 San Francisco 2000   F  85   870    137   395       3     136       4
## 7561 San Francisco 2000   F  86   815    120   343       1     119       4
## 7562 San Francisco 2000   F  87   748    118   305       1      88       6
## 7563 San Francisco 2000   F  88   644     85   271       2      77       4
## 7564 San Francisco 2000   F  89   548     65   221       3      69       1
## 7565 San Francisco 2000   F  90   475     58   177       1      50       1
## 7566 San Francisco 2000   F  91   408     48   153       2      48       1
## 7567 San Francisco 2000   F  92   310     39   122       4      44       1
## 7568 San Francisco 2000   F  93   278     36   105       1      34       2
## 7569 San Francisco 2000   F  94   215     31    85       0      25       0
## 7570 San Francisco 2000   F  95   151     23    65       0      22       1
## 7571 San Francisco 2000   F  96   106     14    52       0       9       0
## 7572 San Francisco 2000   F  97    86     10    45       0       8       0
## 7573 San Francisco 2000   F  98    67      9    31       0       6       0
## 7574 San Francisco 2000   F  99    50     10    20       0       7       0
## 7575 San Francisco 2000   F 100    98     11    47       1      16       0
## 7576 San Francisco 2000   M   0  1245    850   981      29     269      12
## 7577 San Francisco 2000   M   1  1043    749   921      27     279       8
## 7578 San Francisco 2000   M   2   876    723   941      24     314       4
## 7579 San Francisco 2000   M   3   790    721  1054      31     327       7
## 7580 San Francisco 2000   M   4   732    696  1080      31     335       8
## 7581 San Francisco 2000   M   5   708    713  1141      35     336       8
## 7582 San Francisco 2000   M   6   683    732  1147      38     332       9
## 7583 San Francisco 2000   M   7   627    713  1123      40     330      11
## 7584 San Francisco 2000   M   8   661    688  1172      34     389       8
## 7585 San Francisco 2000   M   9   633    696  1211      29     382      10
## 7586 San Francisco 2000   M  10   669    692  1263      39     409       9
## 7587 San Francisco 2000   M  11   671    644  1276      41     383       8
## 7588 San Francisco 2000   M  12   648    651  1263      37     356       8
## 7589 San Francisco 2000   M  13   616    591  1245      29     385       8
## 7590 San Francisco 2000   M  14   654    585  1270      33     388       5
## 7591 San Francisco 2000   M  15   653    624  1303      34     342       7
## 7592 San Francisco 2000   M  16   653    670  1343      45     346       7
## 7593 San Francisco 2000   M  17   686    778  1365      37     360       7
## 7594 San Francisco 2000   M  18   765    902  1357      32     346      10
## 7595 San Francisco 2000   M  19   868   1062  1411      29     351      11
## 7596 San Francisco 2000   M  20  1076   1190  1546      26     348      16
## 7597 San Francisco 2000   M  21  1247   1189  1529      21     328      19
## 7598 San Francisco 2000   M  22  1749   1308  1653      25     302      19
## 7599 San Francisco 2000   M  23  2783   1351  1829      41     322      21
## 7600 San Francisco 2000   M  24  3630   1378  2020      35     342      23
## 7601 San Francisco 2000   M  25  4593   1575  2155      39     351      30
## 7602 San Francisco 2000   M  26  4913   1529  2190      28     345      37
## 7603 San Francisco 2000   M  27  5117   1554  2251      25     351      34
## 7604 San Francisco 2000   M  28  5729   1602  2252      31     403      36
## 7605 San Francisco 2000   M  29  6456   1613  2210      32     428      38
## 7606 San Francisco 2000   M  30  7205   1767  2294      38     501      37
## 7607 San Francisco 2000   M  31  6476   1533  2102      29     483      39
## 7608 San Francisco 2000   M  32  5829   1412  1879      31     446      35
## 7609 San Francisco 2000   M  33  5315   1285  1763      31     439      33
## 7610 San Francisco 2000   M  34  5116   1369  1810      30     443      38
## 7611 San Francisco 2000   M  35  5331   1332  1940      37     539      40
## 7612 San Francisco 2000   M  36  4780   1244  1877      35     499      35
## 7613 San Francisco 2000   M  37  4246   1056  1934      36     500      34
## 7614 San Francisco 2000   M  38  4045    994  1813      38     502      35
## 7615 San Francisco 2000   M  39  3703    990  1762      36     535      35
## 7616 San Francisco 2000   M  40  3803   1005  1850      41     586      37
## 7617 San Francisco 2000   M  41  3307    914  1740      26     549      32
## 7618 San Francisco 2000   M  42  3229    858  1732      25     570      22
## 7619 San Francisco 2000   M  43  3156    825  1734      24     522      33
## 7620 San Francisco 2000   M  44  2964    732  1737      17     477      24
## 7621 San Francisco 2000   M  45  3117    756  1814      19     522      24
## 7622 San Francisco 2000   M  46  2914    689  1696      18     495      28
## 7623 San Francisco 2000   M  47  2830    657  1643      19     474      21
## 7624 San Francisco 2000   M  48  2846    617  1546      18     495      27
## 7625 San Francisco 2000   M  49  2840    603  1639      15     486      21
## 7626 San Francisco 2000   M  50  2974    590  1916      22     536      25
## 7627 San Francisco 2000   M  51  2702    490  1811      14     474      19
## 7628 San Francisco 2000   M  52  2697    481  1569      14     415      20
## 7629 San Francisco 2000   M  53  2732    450  1328      13     391      24
## 7630 San Francisco 2000   M  54  2334    402  1118      10     315      17
## 7631 San Francisco 2000   M  55  2199    370  1109       6     323      13
## 7632 San Francisco 2000   M  56  1999    353  1007      10     294      17
## 7633 San Francisco 2000   M  57  1849    316   937      15     311      16
## 7634 San Francisco 2000   M  58  1608    286   972       9     257      19
## 7635 San Francisco 2000   M  59  1496    307   957      10     268      11
## 7636 San Francisco 2000   M  60  1478    302  1046      10     293      16
## 7637 San Francisco 2000   M  61  1392    275  1023       6     262      13
## 7638 San Francisco 2000   M  62  1325    260   998       3     257       8
## 7639 San Francisco 2000   M  63  1234    230   997       9     253      10
## 7640 San Francisco 2000   M  64  1128    240  1027      11     231       5
## 7641 San Francisco 2000   M  65  1133    256  1104      10     242       8
## 7642 San Francisco 2000   M  66  1035    237  1139      10     209       5
## 7643 San Francisco 2000   M  67   988    205  1077      10     205       5
## 7644 San Francisco 2000   M  68   967    196   975       4     175       6
## 7645 San Francisco 2000   M  69  1044    226  1005       7     178       5
## 7646 San Francisco 2000   M  70  1086    228  1038       8     192       6
## 7647 San Francisco 2000   M  71  1015    207   955      10     167       3
## 7648 San Francisco 2000   M  72   992    194   913       9     182       2
## 7649 San Francisco 2000   M  73   972    182   900       6     180       4
## 7650 San Francisco 2000   M  74   935    177   876       6     163       4
## 7651 San Francisco 2000   M  75  1000    164   912       6     137       2
## 7652 San Francisco 2000   M  76   962    149   810       5     140       1
## 7653 San Francisco 2000   M  77   836    133   821       5     138       2
## 7654 San Francisco 2000   M  78   825    124   736       3     138       2
## 7655 San Francisco 2000   M  79   794     99   579       3     127       1
## 7656 San Francisco 2000   M  80   741     96   538       1     130       2
## 7657 San Francisco 2000   M  81   648     81   418       2     111       3
## 7658 San Francisco 2000   M  82   578     83   335       3      96       2
## 7659 San Francisco 2000   M  83   524     45   305       3      64       1
## 7660 San Francisco 2000   M  84   518     45   295       1      53       2
## 7661 San Francisco 2000   M  85   462     50   248       1      45       1
## 7662 San Francisco 2000   M  86   387     48   203       2      42       2
## 7663 San Francisco 2000   M  87   346     35   166       3      39       1
## 7664 San Francisco 2000   M  88   293     30   148       1      40       1
## 7665 San Francisco 2000   M  89   238     30   141       0      32       2
## 7666 San Francisco 2000   M  90   187     29   110       0      15       2
## 7667 San Francisco 2000   M  91   145     15    87       0      15       0
## 7668 San Francisco 2000   M  92   117     15    70       0      16       0
## 7669 San Francisco 2000   M  93    86     10    45       0      16       0
## 7670 San Francisco 2000   M  94    60      8    40       1      12       1
## 7671 San Francisco 2000   M  95    49      7    23       0       9       0
## 7672 San Francisco 2000   M  96    31      9    11       0       9       0
## 7673 San Francisco 2000   M  97    22      5    12       1       3       0
## 7674 San Francisco 2000   M  98    13      2     8       1       0       0
## 7675 San Francisco 2000   M  99    15      0     8       0       2       0
## 7676 San Francisco 2000   M 100    28      8    17       0       7       0
##      Multirace   Agecat AsianPI
## 1          711   [0,20)    1984
## 2          574   [0,20)    1890
## 3          579   [0,20)    1994
## 4          588   [0,20)    2062
## 5          575   [0,20)    2114
## 6          576   [0,20)    2091
## 7          548   [0,20)    2062
## 8          514   [0,20)    2087
## 9          545   [0,20)    2042
## 10         537   [0,20)    2114
## 11         524   [0,20)    2056
## 12         489   [0,20)    2022
## 13         475   [0,20)    1968
## 14         463   [0,20)    1940
## 15         433   [0,20)    1818
## 16         431   [0,20)    1948
## 17         402   [0,20)    1948
## 18         384   [0,20)    1965
## 19         388   [0,20)    2149
## 20         384   [0,20)    2399
## 21         361  [20,45)    2608
## 22         357  [20,45)    2618
## 23         350  [20,45)    2485
## 24         356  [20,45)    2515
## 25         356  [20,45)    2650
## 26         328  [20,45)    2804
## 27         318  [20,45)    2864
## 28         318  [20,45)    3054
## 29         320  [20,45)    3098
## 30         333  [20,45)    3212
## 31         338  [20,45)    3410
## 32         308  [20,45)    3336
## 33         293  [20,45)    3163
## 34         281  [20,45)    2969
## 35         280  [20,45)    2989
## 36         313  [20,45)    3167
## 37         288  [20,45)    3287
## 38         272  [20,45)    3236
## 39         277  [20,45)    3013
## 40         295  [20,45)    2910
## 41         285  [20,45)    2946
## 42         273  [20,45)    2737
## 43         265  [20,45)    2690
## 44         254  [20,45)    2608
## 45         239  [20,45)    2554
## 46         233  [45,65)    2547
## 47         219  [45,65)    2405
## 48         227  [45,65)    2330
## 49         199  [45,65)    2241
## 50         182  [45,65)    2228
## 51         199  [45,65)    2315
## 52         178  [45,65)    2143
## 53         172  [45,65)    1913
## 54         166  [45,65)    1757
## 55         146  [45,65)    1514
## 56         133  [45,65)    1393
## 57         127  [45,65)    1312
## 58         124  [45,65)    1227
## 59         106  [45,65)    1186
## 60          91  [45,65)    1169
## 61          96  [45,65)    1132
## 62          86  [45,65)    1090
## 63          87  [45,65)    1108
## 64          79  [45,65)    1123
## 65          72  [45,65)    1103
## 66          71 [65,100]    1130
## 67          69 [65,100]    1091
## 68          65 [65,100]    1004
## 69          62 [65,100]     979
## 70          62 [65,100]     964
## 71          63 [65,100]     955
## 72          54 [65,100]     901
## 73          49 [65,100]     859
## 74          48 [65,100]     830
## 75          51 [65,100]     836
## 76          43 [65,100]     778
## 77          40 [65,100]     728
## 78          42 [65,100]     623
## 79          41 [65,100]     543
## 80          48 [65,100]     538
## 81          42 [65,100]     461
## 82          28 [65,100]     400
## 83          23 [65,100]     356
## 84          17 [65,100]     311
## 85          16 [65,100]     284
## 86          12 [65,100]     264
## 87          13 [65,100]     236
## 88           9 [65,100]     191
## 89           9 [65,100]     147
## 90           7 [65,100]     132
## 91           8 [65,100]     105
## 92           4 [65,100]      80
## 93           3 [65,100]      57
## 94           3 [65,100]      46
## 95           3 [65,100]      36
## 96           3 [65,100]      34
## 97           2 [65,100]      24
## 98           1 [65,100]      18
## 99           1 [65,100]      21
## 100          1 [65,100]      15
## 101          1 [65,100]      27
## 102        755   [0,20)    2090
## 103        606   [0,20)    2061
## 104        602   [0,20)    2089
## 105        585   [0,20)    2209
## 106        576   [0,20)    2216
## 107        596   [0,20)    2273
## 108        550   [0,20)    2301
## 109        532   [0,20)    2210
## 110        528   [0,20)    2265
## 111        550   [0,20)    2189
## 112        524   [0,20)    2128
## 113        478   [0,20)    2199
## 114        470   [0,20)    2092
## 115        453   [0,20)    1962
## 116        453   [0,20)    1977
## 117        387   [0,20)    2042
## 118        379   [0,20)    2059
## 119        362   [0,20)    2023
## 120        364   [0,20)    2122
## 121        352   [0,20)    2415
## 122        312  [20,45)    2630
## 123        305  [20,45)    2538
## 124        292  [20,45)    2314
## 125        257  [20,45)    2336
## 126        284  [20,45)    2370
## 127        301  [20,45)    2564
## 128        297  [20,45)    2641
## 129        286  [20,45)    2893
## 130        293  [20,45)    2971
## 131        309  [20,45)    3070
## 132        308  [20,45)    3162
## 133        283  [20,45)    3080
## 134        257  [20,45)    2902
## 135        256  [20,45)    2779
## 136        267  [20,45)    2833
## 137        271  [20,45)    3097
## 138        262  [20,45)    3152
## 139        253  [20,45)    3096
## 140        253  [20,45)    2812
## 141        257  [20,45)    2695
## 142        279  [20,45)    2877
## 143        238  [20,45)    2742
## 144        243  [20,45)    2593
## 145        225  [20,45)    2481
## 146        207  [20,45)    2425
## 147        220  [45,65)    2453
## 148        204  [45,65)    2289
## 149        174  [45,65)    2100
## 150        164  [45,65)    2048
## 151        159  [45,65)    1987
## 152        162  [45,65)    2088
## 153        151  [45,65)    1929
## 154        147  [45,65)    1797
## 155        147  [45,65)    1574
## 156        120  [45,65)    1386
## 157        126  [45,65)    1292
## 158        113  [45,65)    1195
## 159        101  [45,65)    1084
## 160         98  [45,65)     970
## 161         77  [45,65)    1005
## 162         83  [45,65)    1046
## 163         72  [45,65)     952
## 164         64  [45,65)     929
## 165         66  [45,65)     878
## 166         58  [45,65)     839
## 167         51 [65,100]     876
## 168         49 [65,100]     847
## 169         53 [65,100]     799
## 170         46 [65,100]     774
## 171         45 [65,100]     766
## 172         35 [65,100]     745
## 173         39 [65,100]     689
## 174         34 [65,100]     644
## 175         35 [65,100]     626
## 176         35 [65,100]     565
## 177         32 [65,100]     586
## 178         25 [65,100]     535
## 179         26 [65,100]     480
## 180         29 [65,100]     411
## 181         30 [65,100]     364
## 182         28 [65,100]     371
## 183         24 [65,100]     276
## 184         15 [65,100]     234
## 185          8 [65,100]     223
## 186         11 [65,100]     177
## 187         10 [65,100]     143
## 188         10 [65,100]     140
## 189          7 [65,100]     120
## 190          5 [65,100]     106
## 191          6 [65,100]      90
## 192          6 [65,100]      65
## 193          3 [65,100]      43
## 194          3 [65,100]      48
## 195          3 [65,100]      37
## 196          2 [65,100]      22
## 197          1 [65,100]      13
## 198          2 [65,100]      11
## 199          2 [65,100]      10
## 200          0 [65,100]       7
## 201          0 [65,100]       4
## 202          1 [65,100]       8
## 7475       267   [0,20)     975
## 7476       202   [0,20)     958
## 7477       174   [0,20)     958
## 7478       183   [0,20)     979
## 7479       156   [0,20)    1017
## 7480       143   [0,20)    1127
## 7481       132   [0,20)    1137
## 7482       145   [0,20)    1173
## 7483       153   [0,20)    1190
## 7484       139   [0,20)    1190
## 7485       131   [0,20)    1169
## 7486       124   [0,20)    1236
## 7487       132   [0,20)    1209
## 7488       126   [0,20)    1156
## 7489       124   [0,20)    1181
## 7490       127   [0,20)    1241
## 7491       115   [0,20)    1290
## 7492       119   [0,20)    1334
## 7493       122   [0,20)    1296
## 7494       126   [0,20)    1449
## 7495       117  [20,45)    1612
## 7496       114  [20,45)    1659
## 7497       139  [20,45)    1819
## 7498       177  [20,45)    2037
## 7499       181  [20,45)    2226
## 7500       199  [20,45)    2363
## 7501       204  [20,45)    2434
## 7502       202  [20,45)    2499
## 7503       212  [20,45)    2455
## 7504       221  [20,45)    2454
## 7505       208  [20,45)    2490
## 7506       179  [20,45)    2263
## 7507       163  [20,45)    2112
## 7508       147  [20,45)    1981
## 7509       136  [20,45)    2008
## 7510       146  [20,45)    2134
## 7511       134  [20,45)    2138
## 7512       117  [20,45)    2146
## 7513       103  [20,45)    1988
## 7514        95  [20,45)    1909
## 7515       104  [20,45)    2005
## 7516        94  [20,45)    1932
## 7517        95  [20,45)    1929
## 7518        87  [20,45)    1927
## 7519        82  [20,45)    1959
## 7520        83  [45,65)    1994
## 7521        78  [45,65)    1910
## 7522        78  [45,65)    1889
## 7523        76  [45,65)    1884
## 7524        66  [45,65)    1943
## 7525        79  [45,65)    2231
## 7526        74  [45,65)    2106
## 7527        69  [45,65)    1817
## 7528        66  [45,65)    1620
## 7529        58  [45,65)    1427
## 7530        55  [45,65)    1336
## 7531        55  [45,65)    1198
## 7532        47  [45,65)    1137
## 7533        35  [45,65)    1196
## 7534        35  [45,65)    1235
## 7535        41  [45,65)    1315
## 7536        38  [45,65)    1260
## 7537        34  [45,65)    1331
## 7538        40  [45,65)    1308
## 7539        37  [45,65)    1341
## 7540        35 [65,100]    1335
## 7541        26 [65,100]    1366
## 7542        23 [65,100]    1294
## 7543        28 [65,100]    1222
## 7544        26 [65,100]    1326
## 7545        26 [65,100]    1406
## 7546        27 [65,100]    1286
## 7547        33 [65,100]    1268
## 7548        29 [65,100]    1285
## 7549        28 [65,100]    1267
## 7550        33 [65,100]    1138
## 7551        25 [65,100]    1000
## 7552        23 [65,100]     919
## 7553        24 [65,100]     858
## 7554        14 [65,100]     776
## 7555        17 [65,100]     720
## 7556        18 [65,100]     630
## 7557        17 [65,100]     545
## 7558        12 [65,100]     490
## 7559        11 [65,100]     457
## 7560        11 [65,100]     398
## 7561        13 [65,100]     344
## 7562        11 [65,100]     306
## 7563         9 [65,100]     273
## 7564         5 [65,100]     224
## 7565         4 [65,100]     178
## 7566         5 [65,100]     155
## 7567         5 [65,100]     126
## 7568         5 [65,100]     106
## 7569         3 [65,100]      85
## 7570         3 [65,100]      65
## 7571         2 [65,100]      52
## 7572         1 [65,100]      45
## 7573         1 [65,100]      31
## 7574         0 [65,100]      20
## 7575         1 [65,100]      48
## 7576       278   [0,20)    1010
## 7577       213   [0,20)     948
## 7578       187   [0,20)     965
## 7579       182   [0,20)    1085
## 7580       159   [0,20)    1111
## 7581       133   [0,20)    1176
## 7582       134   [0,20)    1185
## 7583       138   [0,20)    1163
## 7584       132   [0,20)    1206
## 7585       126   [0,20)    1240
## 7586       129   [0,20)    1302
## 7587       133   [0,20)    1317
## 7588       121   [0,20)    1300
## 7589       109   [0,20)    1274
## 7590       119   [0,20)    1303
## 7591       123   [0,20)    1337
## 7592       108   [0,20)    1388
## 7593       105   [0,20)    1402
## 7594       110   [0,20)    1389
## 7595       112   [0,20)    1440
## 7596       105  [20,45)    1572
## 7597       103  [20,45)    1550
## 7598       131  [20,45)    1678
## 7599       145  [20,45)    1870
## 7600       162  [20,45)    2055
## 7601       171  [20,45)    2194
## 7602       172  [20,45)    2218
## 7603       186  [20,45)    2276
## 7604       196  [20,45)    2283
## 7605       198  [20,45)    2242
## 7606       213  [20,45)    2332
## 7607       199  [20,45)    2131
## 7608       170  [20,45)    1910
## 7609       154  [20,45)    1794
## 7610       144  [20,45)    1840
## 7611       160  [20,45)    1977
## 7612       165  [20,45)    1912
## 7613       141  [20,45)    1970
## 7614       126  [20,45)    1851
## 7615       132  [20,45)    1798
## 7616       140  [20,45)    1891
## 7617       125  [20,45)    1766
## 7618       119  [20,45)    1757
## 7619       111  [20,45)    1758
## 7620        98  [20,45)    1754
## 7621       106  [45,65)    1833
## 7622       101  [45,65)    1714
## 7623        96  [45,65)    1662
## 7624        80  [45,65)    1564
## 7625        84  [45,65)    1654
## 7626       100  [45,65)    1938
## 7627        78  [45,65)    1825
## 7628        62  [45,65)    1583
## 7629        65  [45,65)    1341
## 7630        57  [45,65)    1128
## 7631        51  [45,65)    1115
## 7632        43  [45,65)    1017
## 7633        42  [45,65)     952
## 7634        43  [45,65)     981
## 7635        38  [45,65)     967
## 7636        40  [45,65)    1056
## 7637        32  [45,65)    1029
## 7638        33  [45,65)    1001
## 7639        35  [45,65)    1006
## 7640        34  [45,65)    1038
## 7641        29 [65,100]    1114
## 7642        19 [65,100]    1149
## 7643        23 [65,100]    1087
## 7644        20 [65,100]     979
## 7645        22 [65,100]    1012
## 7646        30 [65,100]    1046
## 7647        24 [65,100]     965
## 7648        28 [65,100]     922
## 7649        29 [65,100]     906
## 7650        23 [65,100]     882
## 7651        20 [65,100]     918
## 7652        20 [65,100]     815
## 7653        17 [65,100]     826
## 7654        18 [65,100]     739
## 7655        16 [65,100]     582
## 7656        15 [65,100]     539
## 7657        15 [65,100]     420
## 7658        12 [65,100]     338
## 7659        10 [65,100]     308
## 7660         9 [65,100]     296
## 7661         9 [65,100]     249
## 7662         7 [65,100]     205
## 7663         9 [65,100]     169
## 7664         5 [65,100]     149
## 7665         4 [65,100]     141
## 7666         3 [65,100]     110
## 7667         2 [65,100]      87
## 7668         2 [65,100]      70
## 7669         3 [65,100]      45
## 7670         3 [65,100]      41
## 7671         2 [65,100]      23
## 7672         1 [65,100]      11
## 7673         0 [65,100]      13
## 7674         0 [65,100]       9
## 7675         0 [65,100]       8
## 7676         1 [65,100]      17
#### 7 Labels for later use
agelabs <- names(table(bapop$Agecat))
sexlabs <- c("Female", "Male")
racen <- c("White", "AfrAmer", "AsianPI", "Latino", "Multirace",
"AmerInd")
ctylabs <- names(table(bapop$County))
#### 8 Aggregate
bapop2 <- aggregate(bapop[,racen], list(Agecat = bapop$Agecat,
Sex = bapop$Sex, County = bapop$County), sum)
bapop2
##      Agecat Sex        County  White AfrAmer AsianPI Latino Multirace
## 1    [0,20)   F       Alameda  58160   31765   40653  49738     10120
## 2   [20,45)   F       Alameda 112326   44437   72923  58553      7658
## 3   [45,65)   F       Alameda  82205   24948   33236  18534      2922
## 4  [65,100]   F       Alameda  49762   12834   16004   7548      1014
## 5    [0,20)   M       Alameda  61446   32277   42922  53097     10102
## 6   [20,45)   M       Alameda 115745   36976   69053  69233      6795
## 7   [45,65)   M       Alameda  81332   20737   29841  17402      2506
## 8  [65,100]   M       Alameda  33994    8087   11855   5416       711
## 9    [0,20)   F San Francisco  14355    6986   23265  13251      2940
## 10  [20,45)   F San Francisco  85766   10284   52479  23458      3656
## 11  [45,65)   F San Francisco  35617    6890   31478   9184      1144
## 12 [65,100]   F San Francisco  27215    5172   23044   5773       554
## 13   [0,20)   M San Francisco  14881    6959   24541  14480      2851
## 14  [20,45)   M San Francisco 105798   11111   48379  31605      3766
## 15  [45,65)   M San Francisco  43694    7352   26404   8674      1220
## 16 [65,100]   M San Francisco  20072    3329   17190   3428       450
##    AmerInd
## 1      839
## 2     1401
## 3      822
## 4      246
## 5      828
## 6     1263
## 7      687
## 8      156
## 9      173
## 10     526
## 11     282
## 12     121
## 13     165
## 14     782
## 15     354
## 16      76
#### 9 Temp matrix of counts
tmp <- as.matrix(cbind(bapop2[1:4,racen], bapop2[5:8,racen],
bapop2[9:12,racen], bapop2[13:16,racen]))
#### 10 Convert into final array
bapop3 <- array(tmp, c(4, 6, 2, 2))
dimnames(bapop3) <- list(agelabs, racen, sexlabs, ctylabs)
bapop3
## , , Female, Alameda
## 
##           White AfrAmer AsianPI Latino Multirace AmerInd
## [0,20)    58160   31765   40653  49738     10120     839
## [20,45)  112326   44437   72923  58553      7658    1401
## [45,65)   82205   24948   33236  18534      2922     822
## [65,100]  49762   12834   16004   7548      1014     246
## 
## , , Male, Alameda
## 
##           White AfrAmer AsianPI Latino Multirace AmerInd
## [0,20)    61446   32277   42922  53097     10102     828
## [20,45)  115745   36976   69053  69233      6795    1263
## [45,65)   81332   20737   29841  17402      2506     687
## [65,100]  33994    8087   11855   5416       711     156
## 
## , , Female, San Francisco
## 
##          White AfrAmer AsianPI Latino Multirace AmerInd
## [0,20)   14355    6986   23265  13251      2940     173
## [20,45)  85766   10284   52479  23458      3656     526
## [45,65)  35617    6890   31478   9184      1144     282
## [65,100] 27215    5172   23044   5773       554     121
## 
## , , Male, San Francisco
## 
##           White AfrAmer AsianPI Latino Multirace AmerInd
## [0,20)    14881    6959   24541  14480      2851     165
## [20,45)  105798   11111   48379  31605      3766     782
## [45,65)   43694    7352   26404   8674      1220     354
## [65,100]  20072    3329   17190   3428       450      76

The R Code separates, categorize and displays the data by county, year, sex, age and race