作者celestialgod (天)
看板R_Language
标题Re: [问题] 指定group
时间Sun Aug 16 08:51:27 2015
我不知道这样有没有快多少... 参考看看
library(data.table) # setnames
library(plyr) # dlply, laply, mapvalues
library(dplyr) # mutate, summarise, group_by
library(magrittr) # %>%, %$%
dt <- data.table(
series = c(1,2,3,4,5,6,1,1,2,2,2,2,2),
cate = c(2,3,4,2,2,5,1,2,2,3,4,4,4)
)
## method 1 (应该会较慢，而且程式比较丑)
maptable = dlply(dt, .(series), table) %>% laply(function(x){
as.integer(colnames(x)[which.max(x)])
}) %>% cbind(as.integer(names(.)), .)
dt2 = dt %>% mutate(cate = mapvalues(series, maptable[,1], maptable[,2]))
## method 2
-maptable = dt %>% group_by(series, cate) %>% summarise(n()) %>%
group_by(series) %>% summarise(maxCate = max(cate))
dt2 = dt %>% mutate(cate = mapvalues(series, maptable$series,
maptable$maxCate))
## 把上面连在一起的版本 (没有maptable这个暂存变数)
dt2 = dt %>% group_by(series, cate) %>% summarise(n()) %>%
group_by(series) %>% summarise(maxCate = max(cate)) %>%
setnames("series", "maxSeries") %$%
mutate(dt, cate = mapvalues(series, maxSeries, maxCate))
: 例如
: series cate
: [1,] 1 2
: [2,] 2 3
: [3,] 3 4
: [4,] 4 2
: [5,] 5 2
: [6,] 6 5
: [7,] 1 1
: [8,] 1 2
: [9,] 2 2
: [10,] 2 3
: [11,] 2 4
: [12,] 2 4
: [13,] 2 4
: 在series中1出现三次,其cate分别为2,2,1 频率最高的为2
: 想将所有series为1的族群 其cate栏位接指定为2
: 又例如 series中为2的族群 其cate 分别为 3,2,3,4,4,4 频率最高的4
: 想将所有series为2的族群 其cate栏位皆指定为4
: 请问除了用for loop外有其他方法吗?
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 123.205.27.107
※ 文章网址: https://webptt.com/cn.aspx?n=bbs/R_Language/M.1439686290.A.B89.html
1F：推 lambking: 谢谢你! 有快很多 08/16 18:17
2F：→ lambking: 另外想请问如果是series中有missing data 08/16 18:18
3F：→ lambking: 能针对missing data 保留其原本cate 值吗 08/16 18:19
在maptable产生之後加上na.omit，然後加上控制项去控制NA不要覆盖到原本的cate就好
PS: na.omit是把有包含NA的那一个ROW移除，是base的函数
dt <- data.table(
series = c(1,2,3,4,5,6,1,1,2,2,2,2,2,NA,NA,NA),
cate = c(2,3,4,2,2,5,1,2,2,3,4,4,4,3,3,1)
)
dt2 = dt %>% group_by(series, cate) %>% summarise(n()) %>%
group_by(series) %>% summarise(maxCate = max(cate)) %>%
setnames("series", "maxSeries") %>% na.omit %$%
mutate(dt, cate = ifelse(is.na(series), cate,
mapvalues(series, maxSeries, maxCate)))
# series cate
# 1: 1 2
# 2: 2 4
# 3: 3 4
# 4: 4 2
# 5: 5 2
# 6: 6 5
# 7: 1 2
# 8: 1 2
# 9: 2 4
#10: 2 4
#11: 2 4
#12: 2 4
#13: 2 4
#14: NA 3
#15: NA 3
#16: NA 1
4F：推 lambking: 谢谢! 了解了! 08/16 20:59
5F：→ lambking: 但在实际执行时 会出现以下错误讯息 08/16 21:01
6F：→ lambking: Error in n() : This function should not be called 08/16 21:01
7F：→ lambking: directly 08/16 21:03
8F：→ lambking: 查到是因为plyr 和dplyr皆有summarise功能 08/16 21:04
我的版本没问题
你可以指定dplyr:::summarise
9F：推 lambking: 谢谢, 但改完之後反而出现setnames的问题 08/16 22:42
10F：→ lambking: Error in setnames(., "series", "maxSeries") : 08/16 22:42
11F：→ lambking: Items of 'old' not found in column names: series 08/16 22:42
12F：→ lambking: 请问你用的是哪一个版本呢? 谢谢 08/16 22:43
我都没有问题...
我用的R: Revolution R Open 3.2.0
我的session information:
> sessionInfo()
R version 3.2.0 (2015-04-16)
Platform: x86_64-w64-mingw32/x64 (64-bit)
Running under: Windows 7 x64 (build 7601) Service Pack 1
locale:
[1] LC_COLLATE=Chinese (Traditional)_Taiwan.950
[2] LC_CTYPE=Chinese (Traditional)_Taiwan.950
[3] LC_MONETARY=Chinese (Traditional)_Taiwan.950
[4] LC_NUMERIC=C
[5] LC_TIME=Chinese (Traditional)_Taiwan.950
attached base packages:
[1] stats graphics grDevices utils datasets methods base
other attached packages:
[1] magrittr_1.5 dplyr_0.4.2 plyr_1.8.3 data.table_1.9.4
loaded via a namespace (and not attached):
[1] lazyeval_0.1.10 R6_2.0.1 assertthat_0.1 parallel_3.2.0
[5] DBI_0.3.1 tools_3.2.0 reshape2_1.4.1 Rcpp_0.12.0
[9] stringi_0.4-1 stringr_1.0.0 chron_2.3-45
13F：推 lambking: 谢谢 我再试试看 08/17 00:14
14F：→ lambking: 另外想请问如果series 和 cate 都是 character 08/17 00:15
这个应该不用改太多，你自己改一下就好了
※ 编辑: celestialgod (123.205.27.107), 08/17/2015 09:23:53