R_Language 板


LINE

library(plyr) library(dplyr) library(data.table) library(magrittr) set.seed(10) tmp = llply(LETTERS[1:4], function(x) paste0(x, 1:40)) out = expand.grid(tmp[[1]], tmp[[2]], tmp[[3]], tmp[[4]]) %>% tbl_dt() out %<>% .[rep(1:nrow(.), sample(1:5, nrow(.), TRUE)), ] %>% mutate(Value = runif(nrow(.))) dim(out) # 7677529 5 s = proc.time() res = out %>% group_by(Var1, Var2, Var3, Var4) %>% summarise(num_rep = length(Value), value = mean(Value)) %>% filter(num_rep > 1) proc.time() - s # user system elapsed # 3.23 0.11 3.66 结论:七百六十七万列 三秒 收工了 ※ 引述《Godkin (山里的人)》之铭言: : [问题类型]: : : 效能谘询(我想让R 跑更快) : : [软体熟悉度]: : 请把以下不需要的部份删除 : 使用者(已经有用R 做过不少作品) : : [问题叙述]: : 将data frame当中重复出现二次以上的资料抓出,并求取某个栏位的平均值 : 我有一个含有接近四百万比资料的data frame, 格式大概长这个样子 : ID1 ID2 Name1 Name2 Value : ------------------------------------- : m001 g001 mm1 gg1 0.5 : m001 g001 mm1 gg1 0.6 : m002 g001 mm2 gg1 0.3 : m002 g002 mm2 gg2 0.5 : m002 g002 mm2 gg2 0.7 : m002 g002 mm2 gg2 0.2 : m003 g002 mm3 gg2 0.1 : .... .... .... ... ... : ------------------------------------ : 我想把当中重复出现二次以上的资料抓出来,并计算value栏位的平均 : 使得资料变成下面这个样子 : ID1 ID2 Name1 Name2 Value : ------------------------------------- : m001 g001 mm1 gg1 0.55 : m002 g002 mm2 gg2 0.47 : .... .... .... .... .... : ------------------------------------- : 用过for loop + aggregate, foreach + doParallel以及用Rcpp写for loop : 的方式去处理过,但是实在是相当慢,foreach+doParallel跟Rcpp for的版本 : 几乎都跑了接近快12小时, 不知道版上的各位大大有没有什麽好的建议? : : [环境叙述]: : : 请提供 sessionInfo() 的输出结果, : 里面含有所有你使用的作业系统、R 的版本和套件版本资讯, : 让版友更容易找出错误 : : [关键字]: : : data frame, subsetting : --



※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 180.218.154.163
※ 文章网址: https://webptt.com/cn.aspx?n=bbs/R_Language/M.1444297906.A.0B6.html
1F:推 Godkin: 强大!! 感恩!!! 10/08 19:17







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灯, 水草

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

TOP