1

我想像这样拆分数据框

chr.pos nt.pos  CNV
1   74355   0
1   431565  0
1   675207  0
1   783605  1
1   888149  1
1   991311  1
1   1089305 1
1   1177669 1
1   1279886 0
1   1406311 0
1   1491385 0
1   1579761 0
2   1670488 1
2   1758800 1
2   1834256 0
2   1902924 1
2   1978088 1
2   2063124 0

关键是要获取 chr 相同且 CNV=1 列的区间列表，但要考虑到它们之间的 0 inervals

[[1]]

   1    783605  1
   1    888149  1
   1    991311  1
   1    1089305 1
   1    1177669 1
[[2]]
   2    1670488 1
   2    1758800 1
[[3]]
   2    1902924 1
   2    1978088 1

有任何想法吗？

4

1 回答 1

5

您可以使用rle创建一个变量来使用split

  # create a group identifier
  DF$GRP <- with(rle(DF$CNV), rep(seq_along(lengths),lengths))


  # split a subset of DF which contains only CNV==1
  split(DF[DF$CNV==1,],DF[DF$CNV==1,'GRP'] )


$`2`
  chr.pos  nt.pos CNV GRP
4       1  783605   1   2
5       1  888149   1   2
6       1  991311   1   2
7       1 1089305   1   2
8       1 1177669   1   2

$`4`
   chr.pos  nt.pos CNV GRP
13       2 1670488   1   4
14       2 1758800   1   4

$`6`
   chr.pos  nt.pos CNV GRP
16       2 1902924   1   6
17       2 1978088   1   6
于 2012-10-30T11:04:04.633 回答