Tenho um DF
com 14 colunas: item
, grupo
e anos de 2010 a 2021. Eu necessito saber o percentual de cada item
dentro do grupo
ao qual ele faz parte. Tentei utilizar um group_by
com o summarise
, mas como preciso do cálculo para cada um dos anos, não deu certo.
Tentei usar um mutate_at
, mas também não deu certo.
Meu dput
:
structure(list(Item = c("Carvão mineral", "Minerais não-metálicos",
"Petróleo, gás natural e serviços de apoio", "Minério de ferro",
"Minerais metálicos não-ferrosos", "Carne de bovinos e outros prod. de carne",
"Carne de suíno", "Carne de aves", "Pescado industrializado",
"Leite resfriado, esterilizado e pasteurizado", "Outros produtos do laticínio",
"Açúcar", "Conservas de frutas, legumes, outros vegetais e sucos de frutas",
"Óleos e gorduras vegetais e animais", "Café beneficiado",
"Arroz beneficiado e produtos derivados do arroz", "Produtos derivados do trigo, mandioca ou milho",
"Rações balanceadas para animais", "Outros produtos alimentares",
"Bebidas", "Produtos do fumo", "Fios e fibras têxteis beneficiadas",
"Tecidos", "Art. têxteis de uso doméstico e outros têxteis",
"Artigos do vestuário e acessórios", "Calçados e artefatos de couro",
"Produtos de madeira, exclusive móveis", "Celulose", "Papel, papelão, embalagens e artefatos de papel",
"Serviços de impressão e reprodução"), Atividade = c("Extrativas ex-petróleo e gás",
"Extrativas ex-petróleo e gás", "Combustíveis e lubrificantes básicos",
"Extrativas ex-petróleo e gás", "Extrativas ex-petróleo e gás",
"3.10 Fabricação de produtos alimentícios", "3.10 Fabricação de produtos alimentícios",
"3.10 Fabricação de produtos alimentícios", "3.10 Fabricação de produtos alimentícios",
"3.10 Fabricação de produtos alimentícios", "3.10 Fabricação de produtos alimentícios",
"3.10 Fabricação de produtos alimentícios", "3.10 Fabricação de produtos alimentícios",
"3.10 Fabricação de produtos alimentícios", "3.10 Fabricação de produtos alimentícios",
"3.10 Fabricação de produtos alimentícios", "3.10 Fabricação de produtos alimentícios",
"3.10 Fabricação de produtos alimentícios", "3.10 Fabricação de produtos alimentícios",
"3.11 Fabricação de bebidas", "3.12 Fabricação de produtos do fumo",
"3.13 Fabricação de produtos têxteis", "3.13 Fabricação de produtos têxteis",
"3.13 Fabricação de produtos têxteis", "3.14 Confecção de artigos do vestuário e acessórios",
"3.15 Preparação de couros e fabricação de artefatos de couro, artigos para viagem e calçados",
"3.16 Fabricação de produtos de madeira", "3.17 Fabricação de celulose, papel e produtos de papel",
"3.17 Fabricação de celulose, papel e produtos de papel", "3.18 Impressão e reprodução de gravações"
), `2010` = c(0, 3800, 0, 5705, 1540, 15471, 1203, 4387, 909,
2045, 5336, 4713, 5285, 3850, 2291, 3668, 5273, 3974, 23901,
14037, 1459, 1110, 2513, 6913, 27945, 10840, 4671, 1059, 9586,
5218), `2011` = c(0, 4514, 0, 7229, 2054, 14430, 1048, 3753,
861, 2644, 7965, 5350, 5701, 3700, 2356, 3645, 5988, 4355, 26759,
12830, 1754, 2102, 4234, 8279, 33938, 12963, 5309, 1145, 11274,
6862), `2012` = c(0, 3846, 0, 6400, 2316, 21339, 1195, 5872,
1098, 2727, 8653, 5985, 6920, 5022, 2870, 4973, 7719, 5868, 29907,
15592, 1965, 2095, 4347, 10057, 39308, 14406, 6012, 967, 14133,
8790), `2013` = c(0, 4414, 0, 7347, 2580, 27309, 1502, 7518,
1681, 3905, 10728, 6090, 7918, 5326, 3961, 6340, 10684, 6448,
36141, 17583, 2430, 2473, 4684, 12889, 41084, 16638, 6535, 969,
15818, 9260), `2014` = c(0, 4367, 0, 6346, 2209, 31058, 1720,
7881, 1850, 5661, 14518, 6624, 10580, 5920, 4441, 6648, 11421,
7070, 47139, 21731, 2607, 2295, 4921, 13610, 50464, 17915, 6599,
972, 14611, 10414), `2015` = c(0, 4776, 0, 4921, 2606, 35003,
1890, 8601, 1780, 5157, 15002, 6797, 10865, 6057, 5158, 7244,
12818, 7217, 50349, 22886, 2924, 2354, 4893, 14240, 50223, 19992,
6800, 913, 15571, 11804), `2016` = c(0, 3822, 0, 4858, 2102,
37319, 2022, 9592, 1871, 5742, 15526, 9216, 10579, 5864, 5419,
8239, 13446, 8808, 53770, 23388, 3052, 1965, 4533, 14555, 51787,
20728, 6116, 1124, 16783, 12900), `2017` = c(0, 4231, 0, 7002,
2506, 41952, 2295, 10938, 2117, 6115, 16304, 8685, 10980, 6003,
6808, 9077, 14351, 8411, 58195, 26806, 3469, 2080, 4966, 17092,
53900, 20718, 7553, 1324, 18283, 12907), `2018` = c(0, 4231,
0, 7002, 2506, 41952, 2295, 10938, 2117, 6115, 16304, 8685, 10980,
6003, 6808, 9077, 14351, 8411, 58195, 26806, 3469, 2080, 4966,
17092, 53900, 20718, 7553, 1324, 18283, 12907), `2019` = c(0,
4231, 0, 7002, 2506, 41952, 2295, 10938, 2117, 6115, 16304, 8685,
10980, 6003, 6808, 9077, 14351, 8411, 58195, 26806, 3469, 2080,
4966, 17092, 53900, 20718, 7553, 1324, 18283, 12907), `2020` = c(0,
4231, 0, 7002, 2506, 41952, 2295, 10938, 2117, 6115, 16304, 8685,
10980, 6003, 6808, 9077, 14351, 8411, 58195, 26806, 3469, 2080,
4966, 17092, 53900, 20718, 7553, 1324, 18283, 12907), `2021` = c(0,
4231, 0, 7002, 2506, 41952, 2295, 10938, 2117, 6115, 16304, 8685,
10980, 6003, 6808, 9077, 14351, 8411, 58195, 26806, 3469, 2080,
4966, 17092, 53900, 20718, 7553, 1324, 18283, 12907)), row.names = c(NA,
-30L), class = c("tbl_df", "tbl", "data.frame"), na.action = structure(c(`1` = 1L), class = "omit"))
Para os itens do primeiro grupo (Extrativas ex-petróleo e gás), respectivamente, carvão mineral, minerais não-metálicos, minério de ferro e minerais metálicos não-ferrosos, eu precisaria, para 2010, encontrar os seguintes valores: 0.0, 0.344, 0.517 e 0.139.
A ideia geral é encontrar esses valores percentuais para todos os anos.
Códito que utilizei:
library(tidyverse)
variaveis <- c("Carvão mineral", "Minerais não-metálicos", "Minério de ferro", "Minerais metálicos não-ferrosos")
dados %>%
group_by(item, atividade) %>%
mutate_at(vars(variaveis), funs(. / sum(.)))