r - 是否有一种优雅的方式来划分两个摘要的结果？

Question

我有两个数据框 ww 和 pp，它们看起来像这样:

> dput(ww)
structure(list(W_GEW = c(10.2353040126845, 0.831884443723177, 
0.0575388010226499, 1.32932249946475, 2.13545328158608, 1.50612800756816, 
2.41411532240417, 0.272920801750955, 3.79195336455684, 0.450693179969236, 
0.883654166531408, 0.603534664687679, 0.198110485099713, 0.342652280675876, 
0.289928028243725, 6.48435880991905, 0.95407351540347, 8.05653638220527, 
2.28316474679988, 0.487760497903453, 0.224996354463327, 0.553435733400469, 
0.330192354360254, 0.0876946359196857, 1.70203606555473), HP_SEX = c(2, 
2, 1, 2, 1, 1, 1, 2, 1, 2, 2, 2, 2, 1, 1, 1, 2, 2, 1, 1, 2, 1, 
2, 1, 1)), row.names = c(NA, -25L), class = c("tbl_df", "tbl", 
"data.frame"))

> dput(pp)
structure(list(P_GEW = c(0.576128482078767, 1.25180927581909, 
0.102502993629535, 1.35359813678368, 1.56322227704677, 0.692845851279203, 
0.387225241732726, 6.42061013929398, 1.06395274428815, 1.17730952115625, 
1.57374838709105, 6.32338670691465, 1.12358902975024, 0.559852743148847, 
0.244938631158, 0.342884568546086, 0.166101271235082, 1.32857499401104, 
0.898931543979382, 0.830611859630472, 0.822265711424057, 0.183075471978817, 
3.67290388261628, 0.149252785037704, 0.523294953428612), HP_SEX = c(1, 
1, 1, 2, 2, 2, 2, 2, 2, 1, 2, 1, 1, 1, 1, 2, 1, 2, 2, 2, 2, 1, 
2, 1, 1)), row.names = c(NA, -25L), class = c("tbl_df", "tbl", 
"data.frame"))

我的第一步是按 HP_SEX 对数据帧进行分组，并像这样总结它们:

> ww %>% 
+   group_by(HP_SEX) %>%
+   summarise(W_GEW = sum(W_GEW))
# A tibble: 2 x 2
  HP_SEX W_GEW
   <dbl> <dbl>
1      1  22.1
2      2  24.4

> pp %>% 
+   group_by(HP_SEX) %>%
+   summarise(P_GEW = sum(P_GEW)) 
# A tibble: 2 x 2
  HP_SEX P_GEW
   <dbl> <dbl>
1      1  12.4
2      2  21.0

我的下一步是将 W_GEW 与由 HP_SEX 分隔的 P_GEW 分开，如下所示:

22.1/12.4 = 1.78

24.4/21.0 = 1.16

我很好奇，如果有一种优雅的方法可以一步或使用一个功能。

谢谢你。

Answer 1

一种方法可能是:

ww %>% 
 group_by(HP_SEX) %>%
 summarise(W_GEW = sum(W_GEW)) %>%
 left_join(pp %>% 
            group_by(HP_SEX) %>%
            summarise(P_GEW = sum(P_GEW)), by = c("HP_SEX" = "HP_SEX")) %>%
 mutate(GEW = W_GEW/P_GEW)

  HP_SEX W_GEW P_GEW   GEW
   <dbl> <dbl> <dbl> <dbl>
1      1  22.1  12.4  1.79
2      2  24.4  21.0  1.16