- html - 出于某种原因,IE8 对我的 Sass 文件中继承的 html5 CSS 不友好?
- JMeter 在响应断言中使用 span 标签的问题
- html - 在 :hover and :active? 上具有不同效果的 CSS 动画
- html - 相对于居中的 html 内容固定的 CSS 重复背景?
我有一个非常慢且无法正常工作的 for 循环,它在 1 个 data.frame 中查找条形码,然后在另一个 data.frame 中搜索该条形码。第二个 data.frame 的 bar_code 可以出现多次。每次找到条形码时,计数器都应计算条形码出现的次数,并将条形码的数量写入第一个数据框。
我的尝试:
for(i in 1:length(tcgadataUniek$Tumor_Sample_Barcode)){
for(j in 1:length(hprdDataSorted$Samples.Int1)){
count<-0
if(i==j){
count<-count+1
} else {
count<-count+0
}
hprdDataSorted$Samples.Int2<-count[j]
}
}
第一个 Data.Frame 如下所示 (csv):
HUGO.Int1,HUGO.Int2,barcode.Int1
A1CF,APOBEC1,TCGA-B6-A0RS-01A-11D-A099-09
A1CF,TNPO2,TCGA-B6-A0RS-01A-11D-A099-09
A1CF,SYNCRIP,TCGA-B6-A0RS-01A-11D-A099-09
A1CF,KHSRP,TCGA-B6-A0RS-01A-11D-A099-09
A2M,SHBG,TCGA-D8-A1JK-01A-11D-A13L-09
A2M,C11orf58,TCGA-D8-A1JK-01A-11D-A13L-09
A2M,ATF7IP,TCGA-D8-A1JK-01A-11D-A13L-09
AAMP,TH1L,TCGA-A8-A08S-01A-11W-A050-09
AARS,EEF1B2,TCGA-AO-A0JC-01A-11W-A071-09
保存重复条形码 (csv) 的第二个 Data.frame
Sample_Barcode
TCGA-A8-A08G-01A-11W-A019-09
TCGA-AO-A03O-01A-11W-A019-09
TCGA-AO-A03O-01A-11W-A019-09
TCGA-B6-A0RS-01A-11D-A099-09
TCGA-BH-A0HP-01A-12D-A099-09
TCGA-BH-A0HP-01A-12D-A099-09
TCGA-BH-A18H-01A-11D-A12B-09
TCGA-BH-A18H-01A-11D-A12B-09
TCGA-BH-A18J-01A-11D-A12B-09
TCGA-D8-A1JK-01A-11D-A13L-09
TCGA-E2-A1BC-01A-11D-A14G-09
TCGA-E2-A1BC-01A-11D-A14G-09
TCGA-E9-A1NH-01A-11D-A14G-09
TCGA-E9-A22B-01A-11D-A159-09
如果来自 barcode.Int1(dataframe 1)的条形码在 Sample_barcode 中出现了 3 次,脚本应该在脚本正在寻找的 barcode.Int1 旁边添加一个 3。例如:
HUGO.Int1,HUGO.Int2,barcode.Int1, number_of_times
A1CF,APOBEC1,TCGA-B6-A0RS-01A-11D-A099-09,5
最佳答案
Paul 的评论非常恰当,它会大大加快合并步骤。我会使用 table
获取第二个 data.frame 中唯一条形码的计数,并将其merge
到第一个 data.frame 中,见下文:
dat <- structure(list(HUGO.Int1 = c("A1CF", "A1CF", "A1CF", "A1CF",
"A2M", "A2M", "A2M", "AAMP", "AARS"), HUGO.Int2 = c("APOBEC1",
"TNPO2", "SYNCRIP", "KHSRP", "SHBG", "C11orf58", "ATF7IP", "TH1L",
"EEF1B2"), barcode.Int1 = c("TCGA-B6-A0RS-01A-11D-A099-09", "TCGA-B6-A0RS-01A-11D-A099-09",
"TCGA-B6-A0RS-01A-11D-A099-09", "TCGA-B6-A0RS-01A-11D-A099-09",
"TCGA-D8-A1JK-01A-11D-A13L-09", "TCGA-D8-A1JK-01A-11D-A13L-09",
"TCGA-D8-A1JK-01A-11D-A13L-09", "TCGA-A8-A08S-01A-11W-A050-09",
"TCGA-AO-A0JC-01A-11W-A071-09")), .Names = c("HUGO.Int1", "HUGO.Int2",
"barcode.Int1"), class = "data.frame", row.names = c(NA, -9L))
dat2 <- structure(list(Sample_Barcode = c("TCGA-A8-A08G-01A-11W-A019-09",
"TCGA-AO-A03O-01A-11W-A019-09", "TCGA-AO-A03O-01A-11W-A019-09",
"TCGA-B6-A0RS-01A-11D-A099-09", "TCGA-BH-A0HP-01A-12D-A099-09",
"TCGA-BH-A0HP-01A-12D-A099-09", "TCGA-BH-A18H-01A-11D-A12B-09",
"TCGA-BH-A18H-01A-11D-A12B-09", "TCGA-BH-A18J-01A-11D-A12B-09",
"TCGA-D8-A1JK-01A-11D-A13L-09", "TCGA-E2-A1BC-01A-11D-A14G-09",
"TCGA-E2-A1BC-01A-11D-A14G-09", "TCGA-E9-A1NH-01A-11D-A14G-09",
"TCGA-E9-A22B-01A-11D-A159-09")), .Names = "Sample_Barcode", class = "data.frame", row.names = c(NA,
-14L))
foo <- as.data.frame(table(dat2))
merge(dat, foo, by.x='barcode.Int1', by.y='dat2', all.x=TRUE)
# barcode.Int1 HUGO.Int1 HUGO.Int2 Freq
# 1 TCGA-A8-A08S-01A-11W-A050-09 AAMP TH1L NA
# 2 TCGA-AO-A0JC-01A-11W-A071-09 AARS EEF1B2 NA
# 3 TCGA-B6-A0RS-01A-11D-A099-09 A1CF TNPO2 1
# 4 TCGA-B6-A0RS-01A-11D-A099-09 A1CF SYNCRIP 1
# 5 TCGA-B6-A0RS-01A-11D-A099-09 A1CF APOBEC1 1
# 6 TCGA-B6-A0RS-01A-11D-A099-09 A1CF KHSRP 1
# 7 TCGA-D8-A1JK-01A-11D-A13L-09 A2M C11orf58 1
# 8 TCGA-D8-A1JK-01A-11D-A13L-09 A2M ATF7IP 1
# 9 TCGA-D8-A1JK-01A-11D-A13L-09 A2M SHBG 1
data.table版本:
library(data.table)
foo <- data.table(as.data.frame(table(dat2)))
setnames(foo, c('barcode.Int1', 'Freq'))
setkey(foo, barcode.Int1)
dat <- data.table(dat, key='barcode.Int1')
foo[dat]
# barcode.Int1 Freq HUGO.Int1 HUGO.Int2
# 1: NA NA AAMP TH1L
# 2: NA NA AARS EEF1B2
# 3: TCGA-B6-A0RS-01A-11D-A099-09 1 A1CF APOBEC1
# 4: TCGA-B6-A0RS-01A-11D-A099-09 1 A1CF TNPO2
# 5: TCGA-B6-A0RS-01A-11D-A099-09 1 A1CF SYNCRIP
# 6: TCGA-B6-A0RS-01A-11D-A099-09 1 A1CF KHSRP
# 7: TCGA-D8-A1JK-01A-11D-A13L-09 1 A2M SHBG
# 8: TCGA-D8-A1JK-01A-11D-A13L-09 1 A2M C11orf58
# 9: TCGA-D8-A1JK-01A-11D-A13L-09 1 A2M ATF7IP
在“纯”数据表中:
dat <- data.table(dat, key='barcode.Int1')
dat2 <- data.table(dat2)
setnames(dat2, 'barcode.Int1')
setkey(dat2, barcode.Int1)
counts <- dat2[, list(count= .N), by=barcode.Int1]
counts[dat]
关于r - 在数据框中查找条形码并计算这种情况发生的次数,for 循环非常慢?,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/14342661/
我在 GlassFish (J2EE_1.4) 上的 NetBeans 中开发企业项目。我的项目中有一些实体 bean、一些 session bean 和消息驱动 bean。我以如下方式使用 serv
什么在速度方面更好...... 我正在尝试确定用户是否已将某个 URL 添加到他们的快捷方式列表中。如果他们添加了 URL,页面上就会有一个链接,用于从快捷方式中删除该页面,否则他们可以将其添加到快捷
我的问题如下: 我打开一个Excel-File,但我不知道我的客户在模板文件中使用了哪些可能的标头变量。它们可以是:#DATE,#TIME,#NAME等。因此,我需要查找这些变量,以及是否已使用过:替
我有一堆以“-e”结尾的文件要删除。 $ find . -name "*-e" exec rm {} \; find: exec: unknown primary or operator 正则表达式是
我有一个简单的问题:是否可以在 TypeScript 中获取联合的一部分的类型? 例如,您可以经常使用如下查找类型: interface Person { name: string; } type
我正在尝试设置 Point Cloud Library启用 CUDA 选项的主干构建。 我相信我已经按照 these instructions 正确安装了 CUDA . 在 PCL 构建的 cmake
我将首先说我所知道的唯一 VBA 是操作录制的宏的反复试验。我是一名注册会计师,试图以艰难的方式学习 VBA(并希望我去学校学习计算机编程!)。 我有带有多个工作表的大型工作簿。 G 列中以黄色突出显
当文件数达到阈值时,我试图删除目录中最旧的文件。 list_of_files = os.listdir('log') if len([name for name in list_of_files
我有一个数组,它有一些重复的值。 我必须计算每个重复项的数量及其索引。 打印如: Index of b: 1 Index of b: 4 Index of c: 2 Index of c: 3 Ind
我已经搜索了我的问题的解决方案,但没有成功。热键 ctrl+F 找到的 eclipse 查找/替换功能不起作用。注意:通过 Eclipse 菜单 Edit>Find Replace(不工作我的意思是
我想检查 div 是否包含类为“error”的子级,但条件是错误类显示不等于无。 (意味着错误类必须可见。 如何更改我的以下代码: $(".related_field").each(function
这个问题已经有答案了: 已关闭13 年前。 Possible Duplicate: Can jQuery provide the tag name? 嗨! 这个问题太基础了,我不好意思问,但我尝试了
我一直听说这是 cygwin 的路径问题。它阻止了 emacs 在我的 cygwin 中工作。当我在 cli(不是 bash/cygwin)上执行 find 时,无论我输入什么,我都会得到同样的错误。
我正在使用此变量来获取一个或多个与我需要的值相匹配的值。 var mail = $("#dat").contents().find("td:contains('" + name + "')" ).si
请原谅这个长问题。我只是不确定解决这个问题的最佳方法是什么。 我有一个电子表格(Google 表格),其中包含用户和地址列表,我需要从中创建邮寄标签。该电子表格是从我们的学生信息系统导出的。这些地址应
我正在 Excel VBA 中创建一个公式,以解析单元格中以逗号分隔的“部分”列表。在另一个工作表中查找具有该零件名称的单元格,然后使用找到的该单元格的地址来获取同一行不同列的零件成本。我为此工作了数
我被要求在网络应用程序上实现一些电子邮件地址验证 - 我确信我们都已经经历过一千次了...但是,这一次我被要求在域上进行 MX 查找查看它是否接受电子邮件。 有人知道这样做有任何潜在的问题吗? mx
我有一个切换按钮,可读取.wave文件,并且字符串更改为暂停,然后..... 我的问题是,当用户播放声音时,按钮字符串更改为暂停,结束声音后,该字符串仍为暂停状态,我想将其更改为播放。但是我不知道如何
对于令人困惑的标题提前表示歉意。我的问题如下,我在大约 600 个文件中有以下文本: $_REQUEST['FOO'] 我想将其替换为以下内容: $this->input->post('FOO') 为
我正在使用 Ansible 的查找功能查找 INI 文件中的值。这是文档中的示例: - debug: msg="User in integration is {{ lookup('ini', 'use
我是一名优秀的程序员,十分优秀!