mutate more efficiently

Erica · July 31, 2020, 10:33pm

Hi,

I'm using dplyr::mutate to aggregate a large dataset (census government finances).... they provide an excel sheet that shows what items should be aggregated together into groups......

The dataset is 8,000 rows x 300 columns, I want to aggregate columns together by summing so I end up with the 300 columns + 141 new aggregated columns.

I was originally doing this manually ie

agg2017 %>% mutate("Total Revenue"=(B01 + B21 + B22 + B30 + B42 + B46 + B50 + B59 + B79 + B80 + B89 + B91 + B92 + B93 + B94 +
                                    C21 + C30 + C42 + C46 + C50 + C79 + C80 + C89 + C91 + C92 + C93 + C94 +
                                    D21 + D30 + D42 + D46 + D50 + D79 + D80 + D89 + D91 + D92 + D93 + D94 +
                                    T01 + T09 + T10 +  T11 + T12+ T13 + T14 + T15 + T16 + T19 + T20 + T21 + T22 + T23 +
                                    T24 + T25 + T27 + T28 + T29 + T40 + T41 + T50 + T51 + T53 + T99 + 
                                    A01 + A03 + A09 + A10 + A12 + A16 + A18 + A21 + A36 + A44 + A45 +
                                    A50 + A56 + A59 + A60 + A61 + A80 + A81 + A87 + A89 + 
                                    U01 + U11 + U20 + U21 + U30 + U40 + U41 + U50 +U95 + U99 + 
                                    A90 + A91 + A92 + A93 + A94 +
                                    X01 + X02 + X05 + X08 +
                                    Y01 + Y02 + Y04 + Y11 + Y12 + Y51 + Y52)

but there must be a more efficient way right?

Thanks so much!!!!

nirgrahamuk · July 31, 2020, 10:38pm

yeah, a manual but much faster way would be to recognise you can sum() with comma seperated variable names, so you could directly paste the text as provided sum(B01, B21 ,... etc)

next level more automation would be to get the excel info in as a dataframe in itself, and parse the values from the cell, to save you copy and pasting.

Erica · July 31, 2020, 11:14pm

Thanks!

When I use sum in mutate, the resulting column has the same value for every row.

agg2017 %>% mutate("Total Revenue"=sum(B01, B21, B22, B30, B42, B46, B50, B59, B79, B80, B89, B91, B92, B93, B94, C21,
                                       C30, C42, C46, C50, C79, C80, C89, C91, C92, C93, C94, D21, D30, D42, D46, D50, D79,
                                       D80, D89, D91, D92, D93, D94, T01, T09, T10, T11, T12, T13, T14, T15, T16, T19, T20, 
                                       T21, T22, T23, T24, T25, T27, T28, T29, T40, T41, T50, T51, T53, T99, A01, A03, A09, A10,
                                       A12, A16, A18, A21, A36, A44, A45, A50, A59, A60, A61, A80, A81, A87, A89, U01, 
                                       U11, U20, U30, U40, U41, U50, U95, U99, A90, A91,
                                       A92, A93, A94, X01, X02, X05, X08, Y01, Y02, Y11, Y12, Y51, Y52)) -> agg2017

head(agg2017$`Total Revenue`, 50)
 [1] 4524495818 4524495818 4524495818 4524495818 4524495818 4524495818 4524495818 4524495818 4524495818 4524495818
[11] 4524495818 4524495818 4524495818 4524495818 4524495818 4524495818 4524495818 4524495818 4524495818 4524495818
[21] 4524495818 4524495818 4524495818 4524495818 4524495818 4524495818 4524495818 4524495818 4524495818 4524495818
[31] 4524495818 4524495818 4524495818 4524495818 4524495818 4524495818 4524495818 4524495818 4524495818 4524495818
[41] 4524495818 4524495818 4524495818 4524495818 4524495818 4524495818 4524495818 4524495818 4524495818 4524495818

Is my syntax wrong? Thanks!!

Hlynur · August 1, 2020, 12:18am

Hi, you're right. As of dplyr 1.0, there's a really efficient way to do what you want to do, using the rowwise approach.

First let's check that we do in fact have dplyr 1.0 installed, and then I'm just gonna make some dummy data to work on.

packageVersion("dplyr")
#> [1] '1.0.0'

library(dplyr)
#> 
#> Attaching package: 'dplyr'
#> The following objects are masked from 'package:stats':
#> 
#>     filter, lag
#> The following objects are masked from 'package:base':
#> 
#>     intersect, setdiff, setequal, union
demodata <- tibble(col_1 = rep(1, 10),
               col_2 = rep(2, 10),
               col_3 = rep(3, 10),
               col_4 = rep(4, 10))

demodata
#> # A tibble: 10 x 4
#>    col_1 col_2 col_3 col_4
#>    <dbl> <dbl> <dbl> <dbl>
#>  1     1     2     3     4
#>  2     1     2     3     4
#>  3     1     2     3     4
#>  4     1     2     3     4
#>  5     1     2     3     4
#>  6     1     2     3     4
#>  7     1     2     3     4
#>  8     1     2     3     4
#>  9     1     2     3     4
#> 10     1     2     3     4

^{Created on 2020-08-01 by the reprex package (v0.3.0)}

There are two main things to have in mind:

We use rowwise() to set up rowwise calculations
We use the colon operator : for columns that are in succession.

col_1:col_4 means do this for action for every column from col_1 to col_4, both included.

Also, as a sidenote, we use ungroup() to turn of the rowwise calculations again.

demodata %>% 
  rowwise() %>% 
  mutate(`Total revenue` = sum(c(col_1:col_4))) %>% 
  ungroup()
#> # A tibble: 10 x 5
#>    col_1 col_2 col_3 col_4 `Total revenue`
#>    <dbl> <dbl> <dbl> <dbl>           <int>
#>  1     1     2     3     4              10
#>  2     1     2     3     4              10
#>  3     1     2     3     4              10
#>  4     1     2     3     4              10
#>  5     1     2     3     4              10
#>  6     1     2     3     4              10
#>  7     1     2     3     4              10
#>  8     1     2     3     4              10
#>  9     1     2     3     4              10
#> 10     1     2     3     4              10

demodata %>% 
  rowwise() %>% 
  mutate(`Total revenue` = sum(c(col_1, col_3:col_4))) %>% 
  ungroup()
#> # A tibble: 10 x 5
#>    col_1 col_2 col_3 col_4 `Total revenue`
#>    <dbl> <dbl> <dbl> <dbl>           <dbl>
#>  1     1     2     3     4               8
#>  2     1     2     3     4               8
#>  3     1     2     3     4               8
#>  4     1     2     3     4               8
#>  5     1     2     3     4               8
#>  6     1     2     3     4               8
#>  7     1     2     3     4               8
#>  8     1     2     3     4               8
#>  9     1     2     3     4               8
#> 10     1     2     3     4               8

As you can see, you can put together the successions of columns within the c() function call selecting the columns. Like in the example above, skipping `col_2``

So what you could do for your data, assuming all the columns for Total revenue are in succession from B01 to Y52, would be.

agg2017 %>% 
  rowwise() %>% 
  mutate(`Total revenue` = sum(c(B01:Y52))) %>% 
  ungroup()

Hope that helps

Erica · August 2, 2020, 11:16pm

Hlynur:

demodata %>% 
  rowwise() %>% 
  mutate(`Total revenue` = sum(c(col_1, col_3:col_4))) %>% 
  ungroup()
#> # A tibble: 10 x 5
#>    col_1 col_2 col_3 col_4 `Total revenue`
#>    <dbl> <dbl> <dbl> <dbl>           <dbl>
#>  1     1     2     3     4               8
#>  2     1     2     3     4               8
#>  3     1     2     3     4               8
#>  4     1     2     3     4               8
#>  5     1     2     3     4               8
#>  6     1     2     3     4               8
#>  7     1     2     3     4               8
#>  8     1     2     3     4               8
#>  9     1     2     3     4               8
#> 10     1     2     3     4               8

Thank you so much for this in depth answer! When I try it on my data I get NA's.... how do I avoid this? I know I have NAs in some columns but na.rm=TRUE didn't help

I'm attaching a subset of the data

structure(list(Geo_ID = c("01000000000000", "01100100100000", 
"01100200200000", "01100300300000", "01100400400000", "01100500500000", 
"01100600600000", "01100700700000", "01100800800000", "01100900900000"
), Name = c("ALABAMA                                                         ", 
"AUTAUGA COUNTY                                                  ", 
"BALDWIN COUNTY                                                  ", 
"BARBOUR COUNTY                                                  ", 
"BIBB COUNTY                                                     ", 
"BLOUNT COUNTY                                                   ", 
"BULLOCK COUNTY                                                  ", 
"BUTLER COUNTY                                                   ", 
"CALHOUN COUNTY                                                  ", 
"CHAMBERS COUNTY                                                 "
), State_FIPS = c("01", "01", "01", "01", "01", "01", "01", "01", 
"01", "01"), County_FIPS = c("   ", "001", "003", "005", "007", 
"009", "011", "013", "015", "017"), Place_FIPS = c("     ", "99001", 
"99003", "99005", "99007", "99009", "99011", "99013", "99015", 
"99017"), Population = c("  4858979", "    55347", "   203709", 
"    26489", "    22583", "    57673", "    10696", "    20154", 
"   115620", "    34123"), Government_type = c("0", "1", "1", 
"1", "1", "1", "1", "1", "1", "1"), `19U` = c(8137124, 13276, 
95987, 1875, 4472, NA, 1472, 5032, 6893, 484), `24T` = c(38958, 
NA, NA, NA, NA, NA, NA, NA, NA, NA), `29U` = c(1410000, NA, 5348, 
NA, NA, NA, NA, NA, NA, NA), `34T` = c(98959, NA, 280, NA, NA, 
NA, 70, NA, NA, NA), `39U` = c(1399849, 1185, 10682, 633, 486, 
NA, 242, 877, 1031, 484), `44T` = c(411065, NA, 4475, NA, NA, 
NA, 645, NA, NA, NA), `49U` = c(8325192, 12090, 90653, 1242, 
3986, NA, 1230, 4155, 5862, NA), `61V` = c(54357, NA, NA, NA, 
NA, NA, NA, NA, NA, NA), `64V` = c(57220, NA, NA, NA, NA, NA, 
NA, NA, NA, NA), A01 = c(47, NA, NA, NA, NA, NA, NA, NA, NA, 
NA), A16 = c(540132, NA, NA, NA, NA, NA, NA, NA, NA, NA), A18 = c(2095136, 
NA, NA, NA, NA, NA, NA, NA, NA, NA), A21 = c(28367, NA, NA, NA, 
NA, NA, NA, NA, NA, NA), A36 = c(2103921, NA, NA, NA, NA, NA, 
NA, NA, NA, NA), A44 = c(1597, NA, 62, NA, 1030, NA, 57, NA, 
NA, NA), A50 = c(44096, NA, NA, NA, NA, NA, NA, NA, NA, NA), 
    A59 = c(5157, NA, NA, 1, NA, NA, NA, NA, NA, NA), A61 = c(1891, 
    NA, 17, NA, NA, NA, NA, NA, NA, NA), A81 = c(6535, NA, 13226, 
    820, NA, NA, 285, NA, 2530, NA), A87 = c(125872, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA), A89 = c(88027, 1143, 13948, 
    NA, 725, NA, 361, 1015, 5542, NA), A90 = c(335417, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA), B01 = c(18386, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA), B21 = c(1676402, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA), B22 = c(54109, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA), B30 = c(87706, NA, NA, NA, NA, NA, NA, NA, NA, NA
    ), B42 = c(266937, NA, NA, NA, NA, NA, NA, NA, NA, NA), B46 = c(829164, 
    NA, NA, 43, NA, NA, NA, NA, NA, NA), B50 = c(32920, NA, NA, 
    4, NA, NA, NA, NA, NA, NA), B59 = c(87747, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA), B79 = c(6402819, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA), B80 = c(14588, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA), B89 = c(370921, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA), B91 = c(15330, NA, NA, NA, NA, NA, NA, NA, NA, NA), 
    D21 = c(13176, NA, NA, NA, NA, NA, NA, NA, NA, NA), D42 = c(7349, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA), D46 = c(46548, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA), D89 = c(25145, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA), E01 = c(1033, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA), E04 = c(439585, 2737, 13535, 83, NA, 1749, NA, 
    691, NA, 1546), E05 = c(53752, NA, 23, NA, NA, 88, NA, NA, 
    NA, NA), E16 = c(628349, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA), E18 = c(3868594, NA, NA, NA, NA, NA, NA, NA, NA, NA), 
    E21 = c(339568, NA, NA, NA, NA, NA, NA, NA, NA, NA), E22 = c(97126, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA), E23 = c(308672, 3725, 
    7970, 387, NA, 739, 114, 291, NA, NA), E25 = c(196940, NA, 
    4814, 392, NA, 148, 113, 428, NA, NA), E29 = c(52777, 589, 
    9580, NA, 1786, 396, 606, 278, 12577, 457), E32 = c(404599, 
    128, 2753, 3, 183, NA, 40, 24, 385, NA), E36 = c(2109215, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA), E44 = c(566855, 3562, 
    15891, 1582, 4719, 4981, 2999, 4658, 11257, NA), E50 = c(37178, 
    NA, NA, NA, NA, NA, 171, NA, NA, NA), E52 = c(3782, NA, 91, 
    NA, NA, 30, NA, 54, NA, NA), E59 = c(226576, NA, NA, NA, 
    NA, NA, NA, 13, NA, NA), E61 = c(12900, 338, 1354, NA, NA, 
    201, 18, 11, 447, NA), E62 = c(152219, 2271, 12129, 1914, 
    2239, 2804, 504, 748, 13665, NA), E66 = c(82916, NA, 567, 
    99, NA, NA, NA, NA, NA, NA), E74 = c(6273571, NA, NA, NA, 
    NA, NA, NA, NA, NA, NA), E75 = c(206605, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA), E77 = c(52729, NA, NA, 1, NA, NA, NA, 
    NA, NA, NA), E79 = c(627836, 624, 504, NA, NA, 25, 261, 6, 
    100, NA), E85 = c(5874, NA, NA, NA, NA, NA, NA, NA, NA, NA
    ), E87 = c(76987, NA, NA, NA, NA, NA, NA, NA, NA, NA), E89 = c(525053, 
    NA, 7478, 3275, NA, NA, 251, 732, NA, NA), E90 = c(329044, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA), E91 = c(75, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA), F01 = c(3, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA), F04 = c(2540, NA, 70, NA, NA, NA, NA, NA, 
    NA, NA), F16 = c(84749, NA, NA, NA, NA, NA, NA, NA, NA, NA
    ), F18 = c(346940, NA, NA, NA, NA, NA, NA, NA, NA, NA), F21 = c(766, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA), F29 = c(529, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA), F31 = c(875, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA), F32 = c(3006, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA), F36 = c(18767, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA), F44 = c(1011648, NA, 20545, NA, NA, 323, NA, NA, NA, 
    NA), F59 = c(4373, NA, NA, NA, NA, NA, NA, NA, NA, NA), F61 = c(676, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA), F62 = c(40, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA), F87 = c(11963, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA), F89 = c(9636, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA), G01 = c(88, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA), G04 = c(14062, NA, 471, NA, NA, NA, NA, NA, NA, NA), 
    G05 = c(5384, NA, NA, NA, NA, NA, NA, NA, NA, NA), G16 = c(9975, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA), G18 = c(173344, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA), G21 = c(10973, NA, NA, NA, 
    NA, NA, NA, NA, NA, NA), G22 = c(1347, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA), G23 = c(3018, NA, 158, NA, NA, NA, NA, NA, 
    NA, NA), G25 = c(1378, NA, NA, NA, NA, NA, NA, NA, NA, NA
    ), G29 = c(4088, 275, 82, NA, NA, NA, NA, NA, NA, NA), G32 = c(7238, 
    NA, 1, NA, NA, NA, NA, NA, NA, NA), G36 = c(91706, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA), G44 = c(69568, 1219, NA, NA, 
    NA, NA, 63, 1110, NA, NA), G50 = c(338, NA, NA, NA, NA, NA, 
    6, NA, NA, NA), G52 = c(1251, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA), G59 = c(49176, NA, NA, NA, NA, NA, NA, NA, NA, NA
    ), G61 = c(4615, NA, 141, NA, NA, NA, NA, NA, NA, NA), G62 = c(14758, 
    NA, 689, NA, NA, NA, NA, 46, NA, NA), G66 = c(2324, NA, 45, 
    NA, NA, NA, NA, NA, NA, NA), G77 = c(1651, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA), G79 = c(4958, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA), G85 = c(134, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA), G87 = c(2268, NA, NA, NA, NA, NA, NA, NA, NA, NA), G89 = c(12842, 
    NA, 2257, NA, NA, NA, 164, NA, NA, NA), I89 = c(389784, 333, 
    4339, 75, 198, NA, 113, 202, 547, 19), J19 = c(749580, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA), J67 = c(41337, NA, NA, NA, 
    NA, NA, NA, NA, NA, NA), J68 = c(1233518, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA), M01 = c(21764, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA), M04 = c(35996, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA), M05 = c(2625, NA, NA, NA, NA, NA, NA, NA, NA, NA), M21 = c(178852, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA), M23 = c(75558, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA), M25 = c(1883, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA), M29 = c(14, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA), M30 = c(356061, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA), M32 = c(136784, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA), M36 = c(655587, NA, NA, NA, NA, NA, NA, NA, NA, NA), 
    M44 = c(11548, NA, NA, NA, NA, NA, NA, NA, NA, NA), `19T` = c(465299, 
    NA, 4755, NA, NA, NA, 715, NA, NA, NA), M52 = c(5403, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA), M59 = c(15191, NA, NA, NA, 
    NA, NA, NA, NA, NA, NA), M61 = c(4623, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA), M62 = c(4590, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA), M66 = c(2836, NA, NA, NA, NA, NA, NA, NA, NA, NA
    ), M79 = c(104850, NA, NA, NA, NA, NA, NA, NA, NA, NA), M89 = c(207190, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA), Q12 = c(5080209, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA), T01 = c(395930, 5254, 34590, 
    2129, 1979, 5974, 1650, 2802, 11566, 3792), T09 = c(2654883, 
    2227, 12431, 935, 1085, 1617, 760, 1014, 6365, 1753), T10 = c(217790, 
    183, 116, 123, NA, 466, NA, 26, NA, 177), T12 = c(339719, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA), T13 = c(588992, NA, 
    7850, 1655, NA, 2963, NA, NA, NA, 956), T14 = c(1226, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA), T15 = c(688756, NA, 251, 
    NA, NA, NA, NA, NA, NA, NA), T16 = c(188783, NA, 1403, 102, 
    NA, 409, NA, NA, NA, 13), T19 = c(567145, 224, 2609, NA, 
    NA, 122, NA, NA, NA, NA), T20 = c(4086, NA, NA, 8, NA, NA, 
    NA, 11, NA, NA), T22 = c(173265, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA), T23 = c(22956, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA), T24 = c(219393, NA, NA, 201, NA, NA, NA, 58, NA, NA), 
    T25 = c(33730, NA, NA, NA, NA, NA, NA, NA, NA, NA), T27 = c(14346, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA), T28 = c(60763, 60, 1079, 
    87, NA, NA, NA, 17, NA, NA), T40 = c(3624543, NA, NA, NA, 
    NA, NA, NA, NA, NA, NA), T41 = c(520113, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA), T51 = c(46150, NA, 1708, 177, NA, NA, 
    NA, 43, NA, NA), T53 = c(54515, 171, NA, NA, NA, NA, NA, 
    NA, NA, NA), U20 = c(214873, 56, 317, 0, NA, NA, 31, 27, 
    2120, NA), U30 = c(207184, NA, 40, NA, NA, 240, NA, 3, NA, 
    NA), U40 = c(6076, NA, 491, 716, NA, 200, NA, NA, NA, NA), 
    U41 = c(56287, NA, 9, 17, NA, NA, NA, NA, NA, NA), U50 = c(337349, 
    495, 209, 13, NA, 22, NA, 1, NA, NA), U99 = c(542151, 94, 
    5335, 365, 51, 1420, 248, NA, 2231, NA), W01 = c(1331356, 
    NA, 9031, NA, NA, NA, 1018, 108, NA, NA), W31 = c(330053, 
    NA, NA, NA, NA, NA, 1115, NA, NA, NA), W61 = c(12819587, 
    13213, 81785, NA, 2746, 13086, 1378, 5104, 27871, NA), X01 = c(434900, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA), X02 = c(288641, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA), X05 = c(1026823, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA), X08 = c(4159557, NA, NA, NA, 
    NA, NA, NA, NA, NA, NA), X11 = c(3206779, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA), X12 = c(107813, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA), X21 = c(1370957, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA), X30 = c(1196546, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA), X44 = c(4818607, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA), X47 = c(2033, NA, NA, NA, NA, NA, NA, NA, NA, NA), Y01 = c(220211, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA), Y02 = c(11833, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA), Y05 = c(176348, NA, NA, NA, 
    NA, NA, NA, NA, NA, NA), Y06 = c(1, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA), Y07 = c(574485, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA), Y08 = c(574485, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA), Z00 = c(4899183, 6066, 40366, 3182, NA, NA, 1849, 2548, 
    11809, 1148), Z77 = c(4764854, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA), Z78 = c(21074107, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA), M50 = c(30062, NA, NA, NA, NA, NA, NA, NA, NA, NA), 
    C46 = c(NA, 2876, 10732, 2605, 1345, 2223, 2585, 2146, 4038, 
    1917), C79 = c(NA, 1727, NA, NA, NA, NA, NA, NA, NA, NA), 
    C89 = c(NA, 423, 6388, 497, NA, NA, 422, 705, 9871, NA), 
    E81 = c(NA, 128, 11868, 826, NA, NA, 267, NA, 2640, 202), 
    U11 = c(NA, 568, 3682, NA, NA, NA, 6, 843, 120, NA), T29 = c(NA, 
    15, 443, 185, 16, 126, 10, 156, 743, NA), C30 = c(NA, NA, 
    1068, NA, NA, NA, 431, 185, 1001, NA), C94 = c(NA, NA, 2779, 
    NA, NA, NA, NA, NA, NA, NA), E31 = c(NA, NA, 4084, NA, NA, 
    239, 202, 271, NA, 158), E94 = c(NA, NA, 2761, NA, NA, NA, 
    NA, NA, NA, NA), F81 = c(NA, NA, 131, NA, NA, NA, NA, NA, 
    6, NA), F94 = c(NA, NA, 4, NA, NA, NA, NA, NA, NA, NA), G31 = c(NA, 
    NA, 65, NA, NA, NA, NA, NA, NA, NA), G81 = c(NA, NA, 2180, 
    NA, NA, NA, NA, NA, NA, NA), G94 = c(NA, NA, 801, NA, NA, 
    NA, NA, NA, NA, NA), F25 = c(NA, NA, 4, NA, NA, NA, NA, NA, 
    NA, NA), D80 = c(NA, NA, NA, 1, NA, NA, NA, NA, NA, NA), 
    E24 = c(NA, NA, NA, 285, NA, 184, NA, 42, NA, 135), E12 = c(NA, 
    NA, NA, NA, 13, NA, 10, 144, 94, NA), T99 = c(NA, NA, NA, 
    NA, 121, NA, 242, NA, NA, NA), D30 = c(NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_), A80 = c(NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_), E80 = c(NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_
    ), A91 = c(NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_), E45 = c(NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_), C50 = c(NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_), I91 = c(NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_), M12 = c(NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_
    ), D79 = c(NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_), F91 = c(NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_), C42 = c(NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_), D94 = c(NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_), L59 = c(NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_
    ), M24 = c(NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_), M91 = c(NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_), U01 = c(NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_), M18 = c(NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_), M81 = c(NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_
    ), C80 = c(NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_), G24 = c(NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_), F80 = c(NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_), D91 = c(NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_), L89 = c(NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_
    ), L44 = c(NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_), T21 = c(NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_), C91 = c(NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_), C21 = c(NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_), D50 = c(NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_
    ), G80 = c(NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_), A03 = c(NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_), A93 = c(NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_), A92 = c(NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_), E92 = c(NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_
    ), E93 = c(NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_), F92 = c(NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_), G93 = c(NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_), I92 = c(NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_), I93 = c(NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_
    ), E03 = c(NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_), F24 = c(NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_), F66 = c(NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_), F23 = c(NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_), F93 = c(NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_
    ), G92 = c(NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_), G91 = c(NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_), D92 = c(NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_), F52 = c(NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_), G03 = c(NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_
    ), T11 = c(NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_), G45 = c(NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_), F79 = c(NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_), A94 = c(NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_), M93 = c(NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_
    ), M92 = c(NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_), F60 = c(NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_), C92 = c(NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_), C93 = c(NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_), L25 = c(NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_
    ), D93 = c(NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_), M80 = c(NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_), F50 = c(NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_), L32 = c(NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_), L81 = c(NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_
    ), B94 = c(NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_), F03 = c(NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_), M94 = c(NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_), A60 = c(NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_), E60 = c(NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_
    ), L52 = c(NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_), L62 = c(NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_), A12 = c(NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_), I94 = c(NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_), L92 = c(NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_
    ), L91 = c(NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_), L93 = c(NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_), L36 = c(NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_), L01 = c(NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_), L50 = c(NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_
    ), A09 = c(NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_), F12 = c(NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_), G12 = c(NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_), A10 = c(NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_), D11 = c(NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_
    ), B92 = c(NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_), A45 = c(NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_), F22 = c(NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_), F45 = c(NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_), M67 = c(NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_
    ), M87 = c(NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_), Y11 = c(NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_), Y14 = c(NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_), Y21 = c(NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_), G90 = c(NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_
    ), L12 = c(NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_), F77 = c(NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_), F85 = c(NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_), Q18 = c(NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_), U95 = c(NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_
    ), L04 = c(NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_), L67 = c(NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_), G60 = c(NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_), L79 = c(NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_), L05 = c(NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_
    ), S89 = c(NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_), T50 = c(NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_), X42 = c(NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_), Y12 = c(NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_), L18 = c(NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_
    ), L29 = c(NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_), L80 = c(NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_), L23 = c(NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_), L87 = c(NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_), S67 = c(NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_
    ), Y51 = c(NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_), Y53 = c(NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_), Y61 = c(NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_), F05 = c(NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_), M60 = c(NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_
    ), L61 = c(NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_), L66 = c(NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_), Y52 = c(NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_), J85 = c(NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_), L94 = c(NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_
    ), B93 = c(NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_), L60 = c(NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_), M68 = c(NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_), F90 = c(NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_), `Total Revenue` = c(4524495818, 4524495818, 4524495818, 
    4524495818, 4524495818, 4524495818, 4524495818, 4524495818, 
    4524495818, 4524495818), na.rm = c(TRUE, TRUE, TRUE, TRUE, 
    TRUE, TRUE, TRUE, TRUE, TRUE, TRUE)), row.names = c(NA, 10L
), class = "data.frame")

and the code I used here

agg2017 %>% 
  rowwise() %>% 
  mutate(`Total revenue` = sum(c(B01, B21, B22, B30, B42, B46, B50, B59, B79, B80, B89, B91, B92, B93, B94, C21,
                                 C30, C42, C46, C50, C79, C80, C89, C91, C92, C93, C94, D21, D30, D42, D46, D50, D79,
                                 D80, D89, D91, D92, D93, D94, T01, T09, T10, T11, T12, T13, T14, T15, T16, T19, T20, 
                                 T21, T22, T23, T24, T25, T27, T28, T29, T40, T41, T50, T51, T53, T99, A01, A03, A09, A10,
                                 A12, A16, A18, A21, A36, A44, A45, A50, A59, A60, A61, A80, A81, A87, A89, U01, 
                                 U11, U20, U30, U40, U41, U50, U95, U99, A90, A91,
                                 A92, A93, A94, X01, X02, X05, X08, Y01, Y02, Y11, Y12, Y51, Y52))) %>% 
  ungroup() -> tst

head(tst$`Total revenue`, 30)
 [1] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA

I know the data has a lot of NA's but it does have values for some of the columns listed
What's going on?

Thank you so much!!!

Hlynur · August 3, 2020, 8:10am

Hi, really glad that my answer has been of some help.

Using na.rm = TRUE works for me on the subset of the data you posted. Is it possible that you accidentally put na.rm = TRUE within the call to c()? (I ask because that's what I did when first trying to fix the problem ). If I make sure that na.rm = TRUE comes as an argument to sum(), right after the parenthesis closes on c(), it works

library(tidyverse)

agg2017 %>% 
  as_tibble() %>% 
  rowwise() %>% 
  mutate(`Total revenue` = sum(c(B01, B21, B22, B30, B42, B46, B50, B59, B79, B80, B89, B91, B92, B93, B94, C21,
                                 C30, C42, C46, C50, C79, C80, C89, C91, C92, C93, C94, D21, D30, D42, D46, D50, D79,
                                 D80, D89, D91, D92, D93, D94, T01, T09, T10, T11, T12, T13, T14, T15, T16, T19, T20, 
                                 T21, T22, T23, T24, T25, T27, T28, T29, T40, T41, T50, T51, T53, T99, A01, A03, A09, A10,
                                 A12, A16, A18, A21, A36, A44, A45, A50, A59, A60, A61, A80, A81, A87, A89, U01, 
                                 U11, U20, U30, U40, U41, U50, U95, U99, A90, A91,
                                 A92, A93, A94, X01, X02, X05, X08, Y01, Y02, Y11, Y12, Y51, Y52), na.rm = TRUE)) %>% 
  ungroup() %>% 
  select(`Total revenue`)
## A tibble: 10 x 1
#`Total revenue`
#<dbl>
#  1        33248411
#  2           15516
#  3          120783
#  4           10684
#  5            6352
#  6           15782
#  7            7088
#  8            9052
#  9           46127
# 10            8608

Erica · August 3, 2020, 12:53pm

You;re right, I had put it inside the parentheses. Thank you so much!!!

system · August 10, 2020, 12:53pm

This topic was automatically closed 7 days after the last reply. New replies are no longer allowed.