使用库(tidyverse)。
我有一个非常大的数据帧结构像这样(数以千计的观察按站分组为31个级别):
A tibble: 31 x 5
STATION NAME LONGITUDE LATITUDE ELEVATION
<fct> <list> <list> <list> <list>
1 US1ORLA0001 <chr [31]> <NULL> <NULL> <NULL>
2 US1ORLA0003 <chr [31]> <NULL> <NULL> <NULL>
3 US1ORLA0004 <chr [31]> <NULL> <NULL> <NULL>
4 US1ORLA0005 <chr [31]> <NULL> <NULL> <NULL>
5 US1ORLA0031 <chr [31]> <NULL> <NULL> <NULL>
6 US1ORLA0040 <chr [31]> <NULL> <NULL> <NULL>
7 US1ORLA0050 <chr [31]> <NULL> <NULL> <NULL>
8 US1ORLA0076 <chr [31]> <NULL> <NULL> <NULL>
9 US1ORLA0088 <chr [31]> <NULL> <NULL> <NULL>
10 US1ORLA0089 <chr [31]> <NULL> <NULL> <NULL>
... with 21 more rows
我还有一个使用以下代码创建的此数据集的摘要版本:
summary <- raw %>%
group_by(STATION) %>%
summarise(count = n()) %>%
arrange(desc(count))
导致这样:
# A tibble: 31 x 2
STATION count
<fct> <int>
1 USC00353995 17136
2 USC00351324 14110
3 USR0000OHIG 11610
4 USR0000OVIL 11500
5 USR0000OGOO 10855
6 USC00352972 10068
7 USR0000ODUN 8286
8 USC00355204 5880
9 USR0000OCLY 4514
10 USC00352973 4295
# ... with 21 more rows
我的目标是以left_join的方式将原始数据tibble连接到摘要数据tibble上,这将导致每个级别有一行(STATION因子)的数据框、一个COUNT列以及原始数据中的其余列(name、long、lat、ele)。像这样:
# A tibble: 31 x 6
# Groups: STATION [31]
STATION count NAME LONGITUDE LATITUDE ELEVATION
<fct> <int> <fct> <dbl> <dbl> <dbl>
1 USC00353995 17136 HONEYMAN STATE PARK, OR US -124. 43.9 35.1
# ... with 30 more rows
我尝试使用以下命令执行此操作:
desired_output <- left_join(summary, raw, by = c("STATION" = "STATION"))
我希望这会保留x中的行(每个级别为1)和y中的列。问题是,它也保留了y中的所有行,所以我最终得到了这样的结果:
# A tibble: 135,905 x 6
# Groups: STATION [31]
STATION count NAME LONGITUDE LATITUDE ELEVATION
<fct> <int> <fct> <dbl> <dbl> <dbl>
1 USC00353995 17136 HONEYMAN STATE PARK, OR US -124. 43.9 35.1
2 USC00353995 17136 HONEYMAN STATE PARK, OR US -124. 43.9 35.1
3 USC00353995 17136 HONEYMAN STATE PARK, OR US -124. 43.9 35.1
4 USC00353995 17136 HONEYMAN STATE PARK, OR US -124. 43.9 35.1
5 USC00353995 17136 HONEYMAN STATE PARK, OR US -124. 43.9 35.1
6 USC00353995 17136 HONEYMAN STATE PARK, OR US -124. 43.9 35.1
7 USC00353995 17136 HONEYMAN STATE PARK, OR US -124. 43.9 35.1
8 USC00353995 17136 HONEYMAN STATE PARK, OR US -124. 43.9 35.1
9 USC00353995 17136 HONEYMAN STATE PARK, OR US -124. 43.9 35.1
10 USC00353995 17136 HONEYMAN STATE PARK, OR US -124. 43.9 35.1
# ... with 135,895 more rows
这对我的目的没有用。
实现这个和类似结果的最佳方法是什么?我愿意完全不同地做这件事(不使用总结等)。
我们在原始数据集上创建一个列“count”,然后在列上应用不同的
,这样我们每个组只能得到一个观察值
library(dplyr)
raw %>%
group_by(STATION) %>%
mutate(count = n()) %>%
ungroup %>%
distinct(STATION, .keep_all = TRUE)
注意:在输入数据集中,某些列显示为list
。根据显示的预期输出不清楚
data. table方式:
library(data.table)
setDT(raw)
raw[ , .(count=.N, lat=unique(LATITUDE), lng=unique(LONGITUDE)), by=STATION]