R_Language 板


LINE

※ 引述《po5113 (π)》之铭言: : [问题类型]: : : 程式谘询(我想用R 做某件事情,但是我不知道要怎麽用R 写出来) : [软体熟悉度]: : 新手(没写过程式,R 是我的第一次) : [问题叙述]: : 想找出row里面的最大值并挑出该值的colnames : 资料与前几篇matrix找最大值的column name类似 : 部分资料如下: : [,1] [,2] [,3] : [1,] 3.490245e-11 8.579553e-05 9.999142e-01 : [2,] 3.020358e-11 7.424585e-05 9.999258e-01 : [3,] 7.690743e-12 1.890628e-05 9.999811e-01 : [4,] 8.779051e-10 2.153561e-03 9.978464e-01 : [5,] 6.732860e-20 1.655181e-13 1.000000e+00 : [6,] 9.544444e-15 2.346370e-08 1.000000e+00 : [7,] 4.519803e-13 1.111130e-06 9.999989e-01 : ... : ... : 不过目的在找出各row中最大值所在的column name : 因此output中会有跟row number一样的个数 : 有想过用apply加上colnames : 但似乎是遇到了点瓶颈 : 想请问在处理各列单行内的资料时 : 是否有适合使用的函数或其他写法呢? : 感谢各位提供宝贵的意见! : : [关键字]: : row, max, column names 根据我的测试,十万列,apply都还很快,0.57秒而已XD numCols <- 1e2 numRows <- 1e5 mat <- matrix(rnorm(numCols*numRows), numRows, dimnames = list(NULL, paste0(sample(LETTERS, numCols, TRUE), 1:numCols))) st <- proc.time() loc_max <- apply(mat, 1, which.max) out <- colnames(mat)[loc_max] proc.time() - st # user system elapsed # 0.53 0.04 0.57 如果觉得太慢,需要加速什麽的可以试试看matrixStats: library(matrixStats) st <- proc.time() loc_max <- rowOrderStats(mat, which = 1) out2 <- colnames(mat)[loc_max] proc.time() - st # user system elapsed # 0.12 0.00 0.12 all.equal(out, out2) # TRUE 只要原本的20%时间 -- R资料整理套件系列文: magrittr #1LhSWhpH (R_Language) http://tinyurl.com/j3ql84c data.table #1LhW7Tvj (R_Language) http://tinyurl.com/hr77hrn dplyr(上) #1LhpJCfB (R_Language) http://tinyurl.com/jtg4hau dplyr(下) #1Lhw8b-s (R_Language) tidyr #1Liqls1R (R_Language) http://tinyurl.com/jq3o2g3 --



※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 211.76.63.212
※ 文章网址: https://webptt.com/cn.aspx?n=bbs/R_Language/M.1466080833.A.F5D.html ※ 编辑: celestialgod (211.76.63.212), 06/16/2016 21:02:53
1F:推 po5113: 自己没想到这样写就可以了...QQ 06/16 22:28
2F:→ po5113: 还有很多不熟悉的地方,感谢!! 06/16 22:28







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