我想知道什么是定义所有范围的最佳方法,这些范围不在给定的一组范围内。举例来说,如果我有已知坐标一组基因:查找定义的一组范围之外的所有范围
dtGenes <- fread(
"id,start,end
1,1000,1300
2,1200,1500
3,1600,2600
4,3000,4000
")
比方说,我知道,染色体的总长度(为简单起见,假设它们都是同一个染色体上)是10000。所以,最后我期望有间隔区下面的列表:
"startR,endR
0,1000
1500,1600
2600,3000
4000,10000
"
可以Bioconductor的的IRange
有用吗?或者还有其他一些好方法来解决这个问题?
完美的解决方案,非常感谢! –