我有2个向量;
a<-c(1,2,4,3,5,6)
b<-c(1,4,3,2,5,6)
A和B之间的断点数为3 ([1,2] [2,4] [3,5])
B和A之间的断点数为3 ([1,4] [2,3] [2,5])
我基本上想检查数字是否彼此相邻并将其打印出来,因此我不必每次手动执行
这是两个解决方案。首先让
a <- c(1, 2, 4, 3, 5, 6)
b <- c(1, 4, 3, 2, 6, 5)
注意我切换了b
的最后两个元素。然后
# Creating a vector of pairs of subsequent elements of x
fun <- function(x) paste(head(x, -1), tail(x, -1))
BP <- function(x, y) fun(a)[!fun(a) %in% fun(b)]
BP(a, b)
# [1] "1 2" "2 4" "3 5" "5 6"
考虑了此订单。也就是说,它包括5 6
作为断点。另外,如果订单不重要,并且仅是彼此相邻的:
BP <- function(x, y) fun(a)[!(fun(a) %in% fun(b) | fun(a) %in% fun(rev(b)))]
BP(a, b)
# [1] "1 2" "2 4" "3 5"
然后length(BP(a, b))
给出距离。