R_Language 板


LINE

※ 引述《angela79979 (mini)》之铭言: : [问题叙述]: : 有两个list : a.list , b.list : >head(a.list) : $'1' : [1] 3 4 5 8 15 : $'3' : [1] 2 3 6 9 12 14 : ... : >head(b.list) : $'2' : [1] 2 3 5 13 24 : $'1' : [1] 2 3 5 6 7 8 9 12 : ... : 想比较a.list和b.list中 相同编号的list的重复element个数 : 例如: : similarity<-sum(table(a.list$'1'[a.list$'1' %in% b.list$'1'])) : 欲对每一个编号的list都作去最比对 : 但使用loop却没办法对a.list$'i' 或 b.list$'i'作回圈 : 想请问有没有其他的方法 ## data generating - assume no duplicated names and element library(magrittr) list_1 = replicate(50, rbinom(1, 15, .5) %>% sample(1:20, .) %>% sort, simplify = FALSE) %>% set_names(sample(1:50, 50) %>% as.character) list_2 = replicate(50, rbinom(1, 15, .5) %>% sample(1:20, .) %>% sort, simplify = FALSE) %>% set_names(sample(1:50, 50) %>% as.character) ## method 1 - using for intersectNames = intersect(names(list_1), names(list_2)) %>% sort similarity = vector('numeric', length(intersectNames)) %>% set_names(as.character(intersectNames)) for (i in seq_along(intersectNames)) similarity[i] = sum(list_1[[intersectNames[i]]] %in% list_2[[intersectNames[i]]]) print(head(similarity)) # 1 10 11 12 13 14 # 1 2 4 1 3 1 ## method 2 - using map2 (or mapply) library(purrr) intersectNames = intersect(names(list_1), names(list_2)) %>% sort similarity = map2(list_1[intersectNames], list_2[intersectNames], ~ sum(.x %in% .y)) %>% do.call(c, .) %>% set_names(intersectNames) print(head(similarity)) # 1 10 11 12 13 14 # 1 2 4 1 3 1 PS: 其实map2就是mapply... 只是 function(x,y) sum(x %in% y) 改成 ~ sum(.x %in% .y) --



※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 180.218.154.163
※ 文章网址: https://webptt.com/cn.aspx?n=bbs/R_Language/M.1446739230.A.4F1.html ※ 编辑: celestialgod (140.109.73.190), 11/06/2015 14:45:46
1F:推 cywhale: method 2~~笔记.. 11/10 12:38







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

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

TOP