我编写了一个代码,用于将某些位置(数字)映射到一系列位置(数字范围)
library(sqldf)
library(tcltk)
setwd('G:\\Ann\\op\\')
a <- sqldf("select * from file1test f1 inner join file2 f2 on (f1.POS > f2.Start and f1.POS<= f2.Stop) ")
write.table(a, file = "op",sep="\t")
file1test 包含由值
组成的列位置POSPOS
5757
7456
1212
4567
9870
2350
和 file2 包含两列Start和Stop
Start Stop
5751 5758
7440 7460
1201 1224
4560 4569
并且生成的代码的输出类似于
POS Start Stop
5757 5751 5758
1212 1201 1224
4567 4560 4569
7456 7440 7460
但现在我需要打印范围之外的位置,如下所示
POS Start Stop
5757 5751 5758
7456 7440 7460
1212 1201 1224
4567 4560 4569
9870
2350