个人中心
文章发布
退出
作者热门文章
- html - 出于某种原因，IE8 对我的 Sass 文件中继承的 html5 CSS 不友好？
- JMeter 在响应断言中使用 span 标签的问题
- html - 在 :hover and :active? 上具有不同效果的 CSS 动画
- html - 相对于居中的 html 内容固定的 CSS 重复背景？
25
4
我有两个数据集 Data 和 Data1。我想合并这些保持所有差异，同时在新表中为所有公共(public)行添加数值。有什么简单的工具吗？
head(Data)
contig position variantID refAllele altAllele refCount altCount totalCount lowMAPQDepth lowBaseQDepth rawDepth otherBases improperPairs
1 chr1 905373 . T C 2 4 6 0 0 6 0 0
2 chr1 911428 . C T 1 2 3 0 0 3 0 0
3 chr1 953279 . T C 146 126 272 0 0 273 1 0
4 chr1 962184 . T C 14 15 29 0 0 29 0 0
5 chr1 1024129 . T G 1 0 1 0 0 1 0 0
6 chr1 1039514 . C T 1 1 2 0 0 2 0 0
head(Data1)
contig position variantID refAllele altAllele refCount altCount totalCount lowMAPQDepth lowBaseQDepth rawDepth otherBases improperPairs
1 chr1 905373 . T C 2 3 5 0 0 5 0 0
2 chr1 933024 . C T 1 0 1 0 0 1 0 0
3 chr1 953279 . T C 122 124 246 0 0 248 2 0
4 chr1 962184 . T C 17 21 38 0 0 38 0 0
5 chr1 1022518 . G T 0 1 1 0 0 1 0 0
6 chr1 1024129 . T G 1 2 3 0 0 3 0 0
想要的输出示例
contig position variantID refAllele altAllele refCount altCount totalCount lowMAPQDepth lowBaseQDepth rawDepth otherBases improperPairs
1 chr1 905373 . T C 4 7 11 0 0 11 0 0
2 chr1 911428 . C T 1 2 3 0 0 3 0 0
2 chr1 933024 . C T 1 0 1 0 0 1 0 0
4 chr1 953279 . T C 268 150 518 0 0 521 3 0
正如我们在 column position site 905373 中看到的，common 是从 refCount 列添加在一起的。而站点 911428 和 933024 对于它们的数据集都是唯一的，但已插入到新数据集中。他们创建输出表的方式不痛苦吗？
Data <- structure(list(contig = c("chr1", "chr1", "chr1", "chr1", "chr1",
"chr1"), position = c(905373L, 911428L, 953279L, 962184L, 1024129L,
1039514L), variantID = c(".", ".", ".", ".", ".", "."), refAllele = c("T",
"C", "T", "T", "T", "C"), altAllele = c("C", "T", "C", "C", "G",
"T"), refCount = c(2L, 1L, 146L, 14L, 1L, 1L), altCount = c(4L,
2L, 126L, 15L, 0L, 1L), totalCount = c(6L, 3L, 272L, 29L, 1L,
2L), lowMAPQDepth = c(0L, 0L, 0L, 0L, 0L, 0L), lowBaseQDepth = c(0L,
0L, 0L, 0L, 0L, 0L), rawDepth = c(6L, 3L, 273L, 29L, 1L, 2L),
otherBases = c(0L, 0L, 1L, 0L, 0L, 0L), improperPairs = c(0L,
0L, 0L, 0L, 0L, 0L)), row.names = c(NA, 6L), class = "data.frame")
Data1 <- structure(list(contig = c("chr1", "chr1", "chr1", "chr1", "chr1",
"chr1"), position = c(905373L, 933024L, 953279L, 962184L, 1022518L,
1024129L), variantID = c(".", ".", ".", ".", ".", "."), refAllele = c("T",
"C", "T", "T", "G", "T"), altAllele = c("C", "T", "C", "C", "T",
"G"), refCount = c(2L, 1L, 122L, 17L, 0L, 1L), altCount = c(3L,
0L, 124L, 21L, 1L, 2L), totalCount = c(5L, 1L, 246L, 38L, 1L,
3L), lowMAPQDepth = c(0L, 0L, 0L, 0L, 0L, 0L), lowBaseQDepth = c(0L,
0L, 0L, 0L, 0L, 0L), rawDepth = c(5L, 1L, 248L, 38L, 1L, 3L),
otherBases = c(0L, 0L, 2L, 0L, 0L, 0L), improperPairs = c(0L,
0L, 0L, 0L, 0L, 0L)), row.names = c(NA, 6L), class = "data.frame")
最佳答案
这里有一个可能性:
Data
#> contig position variantID refAllele altAllele refCount altCount totalCount
#> 1 chr1 905373 . T C 2 4 6
#> 2 chr1 911428 . C T 1 2 3
#> 3 chr1 953279 . T C 146 126 272
#> 4 chr1 962184 . T C 14 15 29
#> 5 chr1 1024129 . T G 1 0 1
#> 6 chr1 1039514 . C T 1 1 2
#> lowMAPQDepth lowBaseQDepth rawDepth otherBases improperPairs
#> 1 0 0 6 0 0
#> 2 0 0 3 0 0
#> 3 0 0 273 1 0
#> 4 0 0 29 0 0
#> 5 0 0 1 0 0
#> 6 0 0 2 0 0
Data1
#> contig position variantID refAllele altAllele refCount altCount totalCount
#> 1 chr1 905373 . T C 2 3 5
#> 2 chr1 933024 . C T 1 0 1
#> 3 chr1 953279 . T C 122 124 246
#> 4 chr1 962184 . T C 17 21 38
#> 5 chr1 1022518 . G T 0 1 1
#> 6 chr1 1024129 . T G 1 2 3
#> lowMAPQDepth lowBaseQDepth rawDepth otherBases improperPairs
#> 1 0 0 5 0 0
#> 2 0 0 1 0 0
#> 3 0 0 248 2 0
#> 4 0 0 38 0 0
#> 5 0 0 1 0 0
#> 6 0 0 3 0 0
aggregate(. ~ contig + position + variantID + refAllele + altAllele, rbind(Data, Data1), sum)
#> contig position variantID refAllele altAllele refCount altCount totalCount
#> 1 chr1 905373 . T C 4 7 11
#> 2 chr1 953279 . T C 268 250 518
#> 3 chr1 962184 . T C 31 36 67
#> 4 chr1 1024129 . T G 2 2 4
#> 5 chr1 911428 . C T 1 2 3
#> 6 chr1 933024 . C T 1 0 1
#> 7 chr1 1039514 . C T 1 1 2
#> 8 chr1 1022518 . G T 0 1 1
#> lowMAPQDepth lowBaseQDepth rawDepth otherBases improperPairs
#> 1 0 0 11 0 0
#> 2 0 0 521 3 0
#> 3 0 0 67 0 0
#> 4 0 0 4 0 0
#> 5 0 0 3 0 0
#> 6 0 0 1 0 0
#> 7 0 0 2 0 0
#> 8 0 0 1 0 0
关于r - 合并两个数据集，保留所有行差异并添加相似行，我们在Stack Overflow上找到一个类似的问题： https://stackoverflow.com/questions/75176056/
25
4
0
仅使用POSIX shell 的特性，是否有一个“简单命令”，它什么也不做，也不改变$? 的值。人们通常描述:作为 shell 的无操作命令，但这总是设置 $?为零，所以这不是我想要的。 这是生成 s
我需要保留当前的 GlBlendFunc 以便在我做一些工作后可以恢复它。似乎这不是可以用 GLPushAttrib 保存的属性之一，还有其他类似的方法可以用来保存状态吗？ 最佳答案 glGet与
我目前正致力于创建 Fuseki 三重存储浏览器的项目。我需要可视化 TripleStore 中的所有数据并使应用程序可浏览。唯一的问题是 QuerySolution 遗漏了三元组中的“”。 如果我使
我知道没有必要保留委托(delegate)，以避免保留循环。我在一次采访中碰巧遇到了一个问题，“如果保留了 appDelegate 会怎样？”。我对此没有答案，并根据我的知识在这里寻求答案。谢谢 最佳
我正在用 C 编写第一个更严肃的程序，但我陷入了困境。我需要将此列表排序为不同的单独文件，因此它看起来像这样: BE30B Berlin 2014-04-02 Gale 02 BE30B Berlin
我有 4 个页面，我使用 ajax 和 historyPopState 在它们之间进行更改。但是有一个问题，因为页面非常不同，它们使用不同的样式表和脚本。我可以一遍又一遍地下载它们，但我想问一下是否有
我有一个表单，其中包含从一个表创建的多行(与其他表没有关系)。当我保存表单时，我所做的每项更改都会保存，但数据库中确实有一个额外的空行。请参阅下文(希望)了解所有必要的信息。 PropertyAddi
我正在编写一个可以在 Canvas 上绘制气泡的应用程序。我有 MainActivity，它的布局是一个简单的 LinearLayout，我将其用作 fragment 的容器。当我在 Canvas 上
我想弄清楚为什么我通过这段代码得到 EXC_BAD_ACESS。我没有线索。谁能帮帮我。 - (void)loadJsonFile:(NSString*)fileName { NSError
我正在编写一个代码，它将遍历单词中的每个单词，在字典中查找它们，然后将字典值附加到计数器。但是，如果我打印计数器，我只会从我的 if 语句中获取最后一个数字(如果有的话)。如果我将 print cou
抱歉，标题有点困惑。 我想隐藏这些选择选项，只在用户选择按钮时显示。一旦用户选择了一个按钮，我希望选择字段将继续显示......但是如果用户选择了不同的按钮，我想从以前的选择中重置所选的选项并炫耀。
>>> t = "first%s\n" >>> t = t %("second") >>> print t firstsecond 无论如何我可以保留最后的“\n”并得到“firstsecond\n”
我试图弄清楚如何更改 if 语句中的变量并使其保持全局不变。 用户输入:!change Hi var A = "Hello" if (msg.content.includes ('!change'))
您好，我有一个关于在重新创建 Activity 时保留 fragment 的问题。我听说一种方法是在 onCreate 方法中使用 setRetainFragment(true) 。问题是——这与跟踪
我想知道是否有人可以帮助我解决下面代码中的内存管理问题。我对 rootController 特别感兴趣，它是在我执行 initWithRootViewController 时被保留还是在窗口 addS
我想确保我在这里正确理解了内存管理。是否有任何特殊原因在这里使用其中一种 assignCurrentDate 方法而不是另一种方法？此外，所有这些都不会导致内存泄漏，对吗？ 在 .h 中我们有: NS
我对 Angular2 删除尾部斜杠有疑问。我已经设置了我的 dotnet 核心应用程序来添加它们，但是一旦加载了 js，它们就会被删除。 在 Angular2 中甚至可能吗？ (我的客户需要它，所以
local_settings.py 反模式的原因之一是把 SECRET_KEY, AWS设置文件中的键等值有问题: secret 通常应该是这样的: secret !将它们保存在版本控制中意味着拥有存
import('./A'); import('./B'); import('./C'); export class Person {}; A、B 和 C 是纯 JS (es5) 库，它们使用全局 wi
df 是一个测试数据框，其中我只想保留 Hits 列中包含 | 字符的行和行其中包含 hits(即 Hit1、Hit2 等)以及 Hits 列中的空白单元格。 df 从开始 (^) 到结束 ($)
我是一名优秀的程序员,十分优秀！