【求助】按条件过滤数据

相似文件 换一批

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

现有数据如下，一共两列： data.xlsx (8.46 KB)

2018-2-28 14:56:03 上传

       当数据处于同一个chr时，计算每个数据之间的差值，当相邻数据差值之和小于等于150且数据不小于三个时，保留，不符合条件的删除。例如，第2行与第3行、第3行与第4行的position项差距较大，所以1、2行的数据应该删除，第4-10行差值之和小于150，且包含7个数据，所以保留。按照这种规律标有颜色的会被保留，剩下的被删除：

    请问这个用R怎么写，试了半天还是没有思路，不知道我描述清楚没有。。。。。。。

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏1 回帖

关键词：position 不知道 TIO CHR POS

1. library(tidyverse)
   
2. mydata %>%
   
3.   group_by(chr) %>%
   
4.   arrange(position, .by_group = TRUE) %>%
   
5.   mutate(diff = c(0, diff(position)),
   
6.          diff_flag = cumsum(diff >= 150)) %>%
   
7.   group_by(chr, diff_flag) %>%
   
8.   mutate(num = n()) %>%
   
9.   ungroup() %>%
   
10.   filter(num >= 3) %>%
    
11.   select(1:2)

复制代码

结果和你的有一个差异，chr10 28033134，你没有保留，但我的算法保留了，我觉得应该是要保留的，请检查一下。

cheetahfly 发表于 2018-2-28 16:57
结果和你的有一个差异，chr10 28033134，你没有保留，但我的算法保留了，我觉得应该是要保留的，请检查一 ...

非常感谢，检查了一下，是应该保留。这个包我还没有用过，去研究研究

cheetahfly 发表于 2018-2-28 16:57
结果和你的有一个差异，chr10 28033134，你没有保留，但我的算法保留了，我觉得应该是要保留的，请检查一 ...

我安装了这个包，导致我dplyr包是用不了，怎么办

xing035400 发表于 2018-3-1 11:12
我安装了这个包，导致我dplyr包是用不了，怎么办

是不是因为tidyverse 覆盖了 dplyr 的一部分函数， 你可以在调用 dplyr 的函数方法时在前面加上 dplyr::
例如 dplyr::select

菊花冰糖水 发表于 2018-3-1 14:57
是不是因为tidyverse 覆盖了 dplyr 的一部分函数， 你可以在调用 dplyr 的函数方法时在前面加上 dplyr::
...

后来我卸载了tidyverse包，还是有问题

xing035400 发表于 2018-3-1 17:21
后来我卸载了tidyverse包，还是有问题

按说调用 dplyr 函数的时候加上 dplyr::就不会出现错误调用函数的问题。如果结果和之前不一致我觉得要从其他方向考虑一下。

菊花冰糖水 发表于 2018-3-2 14:08
按说调用 dplyr 函数的时候加上 dplyr::就不会出现错误调用函数的问题。如果结果和之前不一致我觉得要从其 ...

已经解决了，还是中文这一块问题