R_Language 板


LINE

※ 引述《locka (locka)》之铭言: : [问题叙述]: : 版上前辈大家晚安~ : 假设我的原始资料栏位有year,month,weekday,y等栏位 : 我想要对他重复做一样的事情 : (根据不同的栏位grouping,计算每组的数量,组内y的平均然後画图) : 因为差别只在於grouping的栏位不同,所以在想说可不可以用函数包起来 : ex: : df_group_fn(df,"year","month") >>> 回传以year,month栏位grouping後计算的结果 : df_group_fn(df,"month","weekday") >>> 回传以month,weekday分组後计算的结果 : 也就是只要输入该data frame跟要grouping的栏位 : 就可以直接回传整理好的结果 : 原本想要用dplyr做,大概像下面这样: : df_group_fn <- function(df,col_1,col_2){ : df %>% group_by(col_1,col_2) %>% summarise(count=n(),avg=mean(y)) %>% : ggplot(aes(mean,n)+geom_point() : } : 不过会卡在指定栏位参数因为是字串的关系, 在group_by那边会有问题 : 所以试着改用data.table的写法: : df_group_fn <- function(df,col_1,col_2){ : df <- as.data.table(df) : df[,`:=`(count=.N, avg=mean(y)),by=c(col_1,col_2)] : ... : } : 可是data.table不会像dplyr一样 : 产生只留下grouping跟summarise栏位的dataframe : 他是在原始的data里面新增栏位,这样我就不知道怎麽画图了... : 总结我的问题: : 1. 希望有高手可以指点用dplyr跟data.table把function写得更有弹性的方法 : 2. 如果我今天不想把grouping的栏位数量写死, : (例如我输入"year"它就只根据year栏位分组, : 输入"year","month","weekday"就根据那三个栏位分组,该怎麽做呢? : 3. 最後想问大家实务上会这麽做吗? 很希望可以听到版上大家分享!! : 先谢谢各位版上先进了 m(_ _)m : [关键字]: : : function, data.table, grouping : 好读版:http://pastebin.com/Yxres7jy 我会建议用wrapr去做这件事情 下面先把一般写法列出给原PO参考 library(dplyr) library(pipeR) library(ggplot2) library(data.table) data("diamonds", package = "ggplot2") # 一般写法 (dplyr) df_group_fn <- function(df, meanCol, col_1, col_2){ df %>>% group_by_(.dots = c(col_1, col_2)) %>>% summarise_(.dots = c(n = "n()", mean = paste0("mean(", meanCol, ")"))) %>>% {ggplot(., aes(mean,n)) + geom_point()} } df_group_fn(diamonds, "price", "cut", "color") # 一般写法 (data.table) dt_group_fn <- function(dt, meanCol, col_1, col_2){ dt[ , .(n = .N, mean = eval(parse(text = paste0("mean(", meanCol, ")")))), by = c(col_1, col_2)] %>>% {ggplot(., aes(mean,n)) + geom_point()} } dt_group_fn(data.table(diamonds), "price", "cut", "color") # wrapr + dplyr library(wrapr) df_group_fn2 <- function(df, meanCol, col_1, col_2){ let(list(y = meanCol, c1 = col_1, c2 = col_2), { df %>>% group_by(c1, c2) %>>% summarise(n = n(), mean = mean(y)) }) %>>% {ggplot(., aes(mean,n)) + geom_point()} } df_group_fn2(diamonds, "price", "cut", "color") # wrapr + data.table dt_group_fn2 <- function(dt, meanCol, col_1, col_2){ let(list(y = meanCol, c1 = col_1, c2 = col_2), { dt[ , .(n = .N, mean = mean(y)), by = .(c1, c2)] }) %>>% {ggplot(., aes(mean,n)) + geom_point()} } dt_group_fn2(data.table(diamonds), "price", "cut", "color") # 进阶,不把栏位给死的方法: # dplyr df_group_fn3 <- function(df, meanCol, groupByCols){ let(list(y = meanCol), { df %>>% group_by_(.dots = groupByCols) %>>% summarise(n = n(), mean = mean(y)) }) %>>% {ggplot(., aes(mean,n)) + geom_point()} } df_group_fn3(diamonds, "price", c("cut", "color")) # data.table dt_group_fn3 <- function(dt, meanCol, groupByCols){ let(list(y = meanCol), { dt[ , .(n = .N, mean = mean(y)), by = groupByCols] }) %>>% {ggplot(., aes(mean,n)) + geom_point()} } dt_group_fn3(data.table(diamonds), "price", c("cut", "color")) 实务上,我自己是做比较接近data engineer的工作 基本上user就会有类似需要,这时候弹性的函数就显得很重要 所以怎麽去利用eval, parse以及...就变得很重要 (没看错就是三个.) 除非全部都是处理data.frame,就可用dplyr透过lazyeval去做 不用wrapr,写起来最漂亮的应该是下面这样: (更正,应该是user用起来最爽XD) # data.table + ... + substitute dt_group_fn3 <- function(dt, meanCol, ...){ groupByCols <- as.character(as.list(substitute(list(...)))[-1L]) y <- substitute(meanCol) dt[ , .(n = .N, mean = mean(y)), by = groupByCols] %>>% {ggplot(., aes(mean,n)) + geom_point()} } dt_group_fn3(data.table(diamonds), price, cut, color) -- R资料整理套件系列文: magrittr #1LhSWhpH (R_Language) https://goo.gl/72l1m9 data.table #1LhW7Tvj (R_Language) https://goo.gl/PZa6Ue dplyr(上.下) #1LhpJCfB,#1Lhw8b-s (R_Language) https://goo.gl/I5xX9b tidyr #1Liqls1R (R_Language) https://goo.gl/i7yzAz pipeR #1NXESRm5 (R_Language) https://goo.gl/zRUISx --



※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 36.235.41.96
※ 文章网址: https://webptt.com/cn.aspx?n=bbs/R_Language/M.1489058786.A.745.html
1F:推 locka: 先感谢C大详细的说明!第一次看到wrapr,我再研究看看怎麽用 03/09 20:19
2F:→ locka: 说到'...'我常常看到套件里面使用这个参数,藉机请问它的处 03/09 20:20
3F:→ celestialgod: andrew板主有写过一篇 请参考#1LV4sfXT 03/09 20:22
4F:→ locka: 谢谢我有看到了,所以用list处理。那麽substitute跟parse呢? 03/09 20:29
5F:→ locka: 不清楚甚麽时候用eval(parse(text=xx))甚麽时候substitute 03/09 20:31
这个我也不会教XD,去看advanced R吧,看你能学到多少了Orz
6F:推 locka: 不过最後一个做法里面,为什麽price等参数不需要引号啊? 03/09 21:30
7F:→ locka: 哈哈好~我会去找来看 谢谢版主大大~ 03/09 21:32
透过substitute转成symbol / name,所以不需要quote~~
8F:推 ginseng21: 这篇收获良多 03/09 21:37
※ 编辑: celestialgod (36.235.41.96), 03/09/2017 21:49:54
9F:推 locka: 真的!! 03/09 22:28
10F:推 cywhale: 推~~ 03/09 22:43







like.gif 您可能会有兴趣的文章
icon.png[问题/行为] 猫晚上进房间会不会有憋尿问题
icon.pngRe: [闲聊] 选了错误的女孩成为魔法少女 XDDDDDDDDDD
icon.png[正妹] 瑞典 一张
icon.png[心得] EMS高领长版毛衣.墨小楼MC1002
icon.png[分享] 丹龙隔热纸GE55+33+22
icon.png[问题] 清洗洗衣机
icon.png[寻物] 窗台下的空间
icon.png[闲聊] 双极の女神1 木魔爵
icon.png[售车] 新竹 1997 march 1297cc 白色 四门
icon.png[讨论] 能从照片感受到摄影者心情吗
icon.png[狂贺] 贺贺贺贺 贺!岛村卯月!总选举NO.1
icon.png[难过] 羡慕白皮肤的女生
icon.png阅读文章
icon.png[黑特]
icon.png[问题] SBK S1安装於安全帽位置
icon.png[分享] 旧woo100绝版开箱!!
icon.pngRe: [无言] 关於小包卫生纸
icon.png[开箱] E5-2683V3 RX480Strix 快睿C1 简单测试
icon.png[心得] 苍の海贼龙 地狱 执行者16PT
icon.png[售车] 1999年Virage iO 1.8EXi
icon.png[心得] 挑战33 LV10 狮子座pt solo
icon.png[闲聊] 手把手教你不被桶之新手主购教学
icon.png[分享] Civic Type R 量产版官方照无预警流出
icon.png[售车] Golf 4 2.0 银色 自排
icon.png[出售] Graco提篮汽座(有底座)2000元诚可议
icon.png[问题] 请问补牙材质掉了还能再补吗?(台中半年内
icon.png[问题] 44th 单曲 生写竟然都给重复的啊啊!
icon.png[心得] 华南红卡/icash 核卡
icon.png[问题] 拔牙矫正这样正常吗
icon.png[赠送] 老莫高业 初业 102年版
icon.png[情报] 三大行动支付 本季掀战火
icon.png[宝宝] 博客来Amos水蜡笔5/1特价五折
icon.pngRe: [心得] 新鲜人一些面试分享
icon.png[心得] 苍の海贼龙 地狱 麒麟25PT
icon.pngRe: [闲聊] (君の名は。雷慎入) 君名二创漫画翻译
icon.pngRe: [闲聊] OGN中场影片:失踪人口局 (英文字幕)
icon.png[问题] 台湾大哥大4G讯号差
icon.png[出售] [全国]全新千寻侘草LED灯, 水草

请输入看板名称,例如:Tech_Job站内搜寻

TOP